Molecular packing and intermolecular interactions in two structural polymorphs of N-palmitoylethanolamine, a type 2 cannabinoid receptor agonist

The molecular structure, packing properties, and intermolecular interactions of two structural polymorphs of N-palmitoylethanolamine (NPEA) have been determined by single-crystal X-ray diffraction. Polymorphs alpha and beta crystallized in monoclinic space group P2(1)/c and orthorhombic space group Pbca, respectively. In both polymorphs, NPEA molecules are organized in a tail-to-tail manner, resembling a bilayer membrane. Although the molecular packing in polymorph alpha is similar to that in N-myristoylethanolamine and N-stearoylethanolamine, polymorph beta is a new form. The acyl chains in both polymorphs are tilted by approximately 35 degrees with respect to the bilayer normal, with their hydrocarbon moieties packed in an orthorhombic subcell. In both structures, the hydroxy group of NPEA forms two hydrogen bonds with the hydroxy groups of molecules in the opposite leaflet, resulting in extended, zig-zag type H-bonded networks along the b-axis in polymorph alpha and along the a-axis in polymorph beta. Additionally, the amide N-H and carbonyl groups of adjacent molecules are involved in N-H...O hydrogen bonds that connect adjacent molecules along the b-axis and a-axis, respectively, in alpha and beta. Whereas in polymorph alpha the L-shaped NPEA molecules in opposite layers are arranged to yield a Z-like organization, in polymorph beta one of the two NPEA molecules is rotated 180 degrees , leading to a W-like arrangement. Lattice energy calculations indicate that polymorph alpha is more stable than polymorph beta by approximately 2.65 kcal/mol.     

Conformational flexibility in a triazole nucleoside derivative: 4-cyano-5-cyanomethyl-1-(2,3,5-tri-O-acetyl-beta-D-ribofuranosyl)-1,2,3-triazole

The crystal-structure determination of the molecular structure of the hydrophobic compound, 4-cyano-5-cyanomethyl-1-(2,3,5-tri-O-acetyl-beta-D-ribofuranosyl)-1,2,3-triazole, C16H17N5O7, provides us with two different conformations of ribofuranosyl moieties [(C2'-exo, C3'-endo) and C2'-exo] with two markedly different N-glycosidic angles. There are two molecules in the asymmetric unit of the crystal and the overall stereochemistry of the molecules are influenced predominantly by weak intramolecular bifurcated and trifurcated hydrogen bonds of the type C-H...O and C-H...N, where endo-H atoms attached to C2' and C3' are involved. The molecules in the crystal are interconnected with longer intermolecular bonds of the same type. There are empty channels (occupying 14.0% of the whole volume of the unit cell), which are extended along b-axis of the entire crystal.     

Strong and weak hydrogen bonds in protein-ligand complexes of kinases: a comparative study

Strong and weak hydrogen bonds between protein and ligand are analyzed in a group of 233 X-ray crystal structures of the kinase family. These kinases are from both eukaryotic and prokaryotic organisms. The dataset comprises of 44 sub-families, out of which 35 are of human origin and the rest belong to other organisms. Interaction analysis was carried out in the active sites, defined here as a sphere of 10 A radius around the ligand. A majority of the interactions are observed between the main chain of the protein and the ligand atoms. As a donor, the ligand frequently interacts with amino acid residues like Leu, Glu and His. As an acceptor, the ligand interacts often with Gly, and Leu. Strong hydrogen bonds N-H...O, O-H...O, N-H...N and weak bonds C-H...O, C-H...N are common between the protein and ligand. The hydrogen bond donor capacity of Gly in N-H...O and C-H...O interactions is noteworthy. Similarly, the acceptor capacity of main chain Glu is ubiquitous in several kinase sub-families. Hydrogen bonds between protein and ligand form characteristic hydrogen bond patterns (supramolecular synthons). These synthon patterns are unique to each sub-family. The synthon locations are conserved across sub-families due to a higher percentage of conserved sequences in the active sites. The nature of active site water molecules was studied through a novel classification scheme, based on the extent of exposure of water molecules. Water which is least exposed usually participates in hydrogen bond formation with the ligand. These findings will help structural biologists, crystallographers and medicinal chemists to design better kinase inhibitors.     

Occurrence of bifurcated three-center hydrogen bonds in proteins

Analysis of 13 high-resolution protein X-ray crystal structures shows that 1204 (24%) of all the 4974 hydrogen bonds are of the bifurcated three-center type with the donor X-H opposing two acceptors A1, A2. They occur systematically in alpha-helices where 90% of the hydrogen bonds are of this type; the major component is (n + 4)N-H ... O = C(n) as expected for a 3.6(13) alpha-helix, and the minor component is (n + 4)N-H ... O = C(n + 1), as observed in 3(10) helices; distortions at the C-termini of alpha-helices are stabilized by three-center bonds. In beta-sheets 40% of the hydrogen bonds are three-centered. The frequent occurrence of three-center hydrogen bonds suggests that they should not be neglected in protein structural studies.     

Estimation of strength in different extra Watson-Crick hydrogen bonds in DNA double helices through quantum chemical studies

It was shown earlier, from database analysis, model building studies, and molecular dynamics simulations that formation of cross-strand bifurcated or Extra Watson-Crick hydrogen (EWC) bonds between successive base pairs may lead to extra rigidity to DNA double helices of certain sequences. The strengths of these hydrogen bonds are debatable, however, as they do not have standard linear geometry criterion. We have therefore carried out detailed ab initio quantum chemical studies using RHF/6-31G(2d,2p) and B3LYP/6-31G(2p,2d) basis sets to determine strengths of several bent hydrogen bonds with different donor and acceptors. Interaction energy calculations, corrected for the basis set superposition errors, suggest that N-H...O type bent EWC hydrogen bonds are possible along same strands or across the strands between successive base pairs, leading to significant stability (ca. 4-9 kcal/mol). The N-H...N and C-H...O type interactions, however, are not so stabilizing. Hence, consideration of EWC N-H...O H-bonds can lead to a better understanding of DNA sequence directed structural features.     

Crystal and molecular structure of cyclo(L-prolyl-glycyl)3. A cyclic hexapeptide with a cis peptide bond

The crystal structure of cyclo(L-Pro-Gly)3 was solved using X-ray crystallographic techniques. The backbone of the peptide is asymmetric and is made up of five trans peptide units and one cis peptide. There is a hydrogen bonded water bridge that links the carbonyl oxygens, O1 and O4. The molecules exist as dimers in the crystal lattice. The two molecules of the dimer are related by crystallographic twofold symmetry and are linked by two N-H ... O hydrogen bonds. The crystals are trigonal, space group P3(2)12 with a = 11.379(3), c = 32.93(1) and z = 6. The structure was solved by multisolution methods and refined by least squares technique to an R of 0.083.     

Molecular and crystal structures of chitosan/HI type I salt determined by X-ray fiber diffraction

The three-dimensional structure of chitosan/HI type I salt was determined by the X-ray fiber diffraction technique and linked-atom least-squares refinement method. Two polymer chains and four iodide ions (I(-)) crystallized in a monoclinic unit cell with dimensions a = 9.46(2), b = 9.79(2)], c (fiber axis)=10.33(2)A, beta = 105.2(2) degrees and a space group P2(1). Chitosan chains adopted an extended twofold helical conformation that was stabilized by O-3...O-5 hydrogen bonds, and the O-6 atom adopted nearly gt orientation. Polymer chains zigzag along the b-axis and directly connect to each other by N-2...O-6 hydrogen bonds. Two columns of iodide ions were shown to pack at the bending points of the zigzag sheets, and their locations are closely related to those of water columns in the hydrated chitosan. The iodide ions stabilized the salt structure by forming hydrogen bonds with the N-2 and O-6 atoms of the polymer chains together with an electrostatic interaction between N-2 and the iodide ions.     

Crystal structure of N-(1-deoxy-beta-D-fructopyranos-1-yl)-L-proline-an Amadori compound

Here we report the crystal structure data on N-(1-deoxy-beta-D-fructopyranos-1-yl)-L-proline (Fru-Pro)-an Amadori compound. X-ray crystal and molecular structures of its two isomorphous crystalline forms, (Fru-Pro)xMeOH, C(11)H(19)NO(7)xCH(4)O (1a) and (Fru-Pro)x2H(2)O, C(11)H(19)NO(7)x2H(2)O (1b) were determined. In 1a and 1b the compound crystallizes as the beta-anomer with the overall geometry of Fru-Pro zwitterions being very similar. Fructose ring adopts the chair (2)C(5) conformation with the proline moiety bonded to equatorial C-1 atom and remaining in a trans-gauche (tg) orientation with respect to the sugar ring. The five-membered pyrrolidine ring adopts an envelope conformation, with the Cbeta atom puckered. Fructosyl and carboxylate groups are in bisectional and axial positions of pyrrolidine ring, respectively. The overall molecular geometry of Fru-Pro zwitterions, especially the relative orientation of sugar and amino acid moieties, is stabilized by intramolecular, three-centred N-H...O(Fru)/O(Pro) hydrogen bonds (with bifurcated acceptor) formed between aminium and hydroxyl/carboxylate groups. The packing diagrams are very similar in both 1a and 1b with the adjacent zwitterions linked to each other by the extensive network of O-H...O and C-H...O hydrogen bonds to form channels along the a-axis, filled up with solvent molecules.     

Channel‐like crystal structure of cinchoninium L‐O‐phosphoserine salt dihydrate

Studies on the interactions between L ‐O‐ phosphoserine, as one of the simplest fragments of membrane components, and the Cinchona alkaloid cinchonine, in the crystalline state were performed. Cinchoninium L ‐O‐phosposerine salt dihydrate (PhSerCin) crystallizes in a monoclinic crystal system, space group P2₁, with unit cell parameters: a = 8.45400(10) Å, b = 7.17100(10) Å, c = 20.7760(4) Å, α = 90°, β = 98.7830(10)°, γ = 90°, Z = 2. The asymmetric unit consists of the cinchoninium cation linked by hydrogen bonds to a phosphoserine anion and two water molecules. Intermolecular hydrogen bonds connecting phosphoserine anions via water molecules form chains extended along the b axis. Two such chains symmetrically related by twofold screw axis create a “channel.” On both sides of this channel cinchonine cations are attached by hydrogen bonds in which the atoms N1, O12, and water molecules participate. This arrangement mimics the system of bilayer biological membrane. Chirality 2010. © 2009 Wiley‐Liss, Inc.     

Structure of myohemerythrin in the azidomet state at 1.7/1.3 A resolution

The molecular model of myohemerythrin, an oxygen-carrying protein from sipunculan worms, has been refined by stereochemically restrained least-squares minimization at 1.7/1.3 A resolution to a conventional R-value of 0.158. The estimated positional standard deviation is better than 0.15 A for most of the 979 protein atoms. The average isotropic displacement parameter, B, for the protein atoms is 23.1 A2. This high average B parameter appears to be due to the overall motion of the molecule, which correlates with the observed anisotropic diffraction. The side-chains of seven residues were modeled in two conformations, i.e. the side-chains were discretely disordered, and B parameters for several lysine and glutamate side-chains indicate that they are poorly localized. Of the residues in myohemerythrin, 66% are helical, with 62% occurring in four long alpha-helices with mean values for the backbone torsion angles of phi = -65 degrees, psi = -42 degrees, and for the hydrogen bonds distances of N ... O, 3.0 A and H ... O, 2.1 A, and angles of N ... O = C, 153 degrees, N-H ... O, 157 degrees, and H ... O = C, 147 degrees. For two-thirds of the alpha-helical residues, the torsional rotation of the C alpha-C beta bond, chi 1, is approximately -60 degrees, and for one-third chi 1 is approximately 180 degrees. Although most turns in myohemerythrin are well-categorized by previous classification, two do not fit in established patterns. Also included in the refined model are three sulfate ions, all partially occupied, and 157 water molecules, 40% of which are modeled fully occupied. Only one water molecule is internal to the protein, the remainder occur on the surface and are observed principally between symmetry-related molecules contributing, along with van der Waals' contacts, most of the interactions between molecules. There are eight intermolecular protein-protein hydrogen bonds, of which only four are between well-located atoms.     

N-H...O, O-H...O, and C-H...O hydrogen bonds in protein-ligand complexes: strong and weak interactions in molecular recognition

The characteristics of N-H...O, O-H...O, and C-H...O hydrogen bonds are examined in a group of 28 high-resolution crystal structures of protein-ligand complexes from the Protein Data Bank and compared with interactions found in small-molecule crystal structures from the Cambridge Structural Database. It is found that both strong and weak hydrogen bonds are involved in ligand binding. Because of the prevalence of multifurcation, the restrictive geometrical criteria set up for hydrogen bonds in small-molecule crystal structures may need to be relaxed in macromolecular structures. For example, there are definite deviations from linearity for the hydrogen bonds in protein-ligand complexes. The formation of C-H...O hydrogen bonds is influenced by the activation of the C(alpha)-H atoms and by the flexibility of the side-chain atoms. In contrast to small-molecule structures, anticooperative geometries are common in the macromolecular structures studied here, and there is a gradual lengthening as the extent of furcation increases. C-H...O bonds formed by Gly, Phe, and Tyr residues are noteworthy. The numbers of hydrogen bond donors and acceptors agree with Lipinski's "rule of five" that predicts drug-like properties. Hydrogen bonds formed by water are also seen to be relevant in ligand binding. Ligand C-H...O(w) interactions are abundant when compared to N-H...O(w) and O-H...O(w). This suggests that ligands prefer to use their stronger hydrogen bond capabilities for use with the protein residues, leaving the weaker interactions to bind with water. In summary, the interplay between strong and weak interactions in ligand binding possibly leads to a satisfactory enthalpy-entropy balance. The implications of these results to crystallographic refinement and molecular dynamics software are discussed.     

Variability in the backbone conformation of cyclic pentapeptides. Crystal structure of cyclic(Gly-L-Pro-D-Phe-Gly-L-Ala)

All the peptide bonds in cyclic(Gly-L-Pro-D-Phe-Gly-L-Ala) are in the trans conformation; however, the peptide bond C'5-N1 is twisted by 19 degrees from planarity (omega 5 = -161 degrees). A Type II beta-turn encompasses the L-Pro-D-Phe residues. Carbonyl oxygens O2, O4 and O5 are directed to the same side of the average plane through the backbone ring and they form hydrogen bonds with N3, N5 and N1, respectively, in adjacent molecules in a stacked column where the adjacent molecules are related by one translational unit. The conformation of the backbone is different from that established in other molecules with the DLDDL chirality sequence. The P21 cell contains two molecules of C21H26N5O5 with a = 4.836(2) A, b = 18.346(8) A, c = 12.464(5) A and beta = 100.05(4) degrees. The R factor for 1382 data with [F0[ greater than 1 sigma is 7.0%.     

Observation of a unique pattern of bifurcated hydrogen bonds in the crystal structures of the N-glycoprotein linkage region models

Elucidation of the intra- and intermolecular carbohydrate-protein interactions would greatly contribute toward obtaining a better understanding of the structure-function correlations of the protein-linked glycans. The weak interactions involving C-H...O have recently been attracting immense attention in the domain of biomolecular recognition. However, there has been no report so far on the occurrence of C-H...O hydrogen bonds in the crystal structures of models and analogs of N-glycoproteins. We present herein an analysis of C-H...O interactions in the crystal structures of all N-glycoprotein linkage region models and analogs. The study reveals a cooperative network of bifurcated hydrogen bonds consisting of N-H...O and C-H...O interactions seen uniquely for the models. The cooperative network consists of two antiparallel chains of bifurcated hydrogen bonds, one involving N1-H, C2'-H and O1' of the aglycon moiety and the other involving N2-H, C1-H and O1' of the sugar. Such bifurcated hydrogen bonds between the core glycan and protein are likely to play an important role in the folding and stabilization of proteins.     

A two-stranded template-based approach to G.(C-A) triad formation: designing novel structural elements into an existing DNA framework

We have designed a DNA sequence, d(G-G-G-T-T-C-A-G-G), which dimerizes to form a 2-fold symmetric G-quadruplex in which G(syn). G(anti).G(syn).G(anti) tetrads are sandwiched between all trans G. (C-A) triads. The NMR-based solution structural analysis was greatly aided by monitoring hydrogen bond alignments across N-H...N and N-H...O==C hydrogen bonds within the triad and tetrad, in a uniformly ((13)C,(15)N)-labeled sample of the d(G-G-G-T-T-C-A-G-G) sequence. The solution structure establishes that the guanine base-pairs with the cytosine through Watson-Crick G.C pair formation and with adenine through sheared G.A mismatch formation within the G.(C-A) triad. A model of triad DNA was constructed that contains the experimentally determined G.(C-A) triad alignment as the repeating stacked unit.     

Behavior of cholesterol and its effect on head group and chain conformations in lipid bilayers: a molecular dynamics study.

Cholesterol molecules were put into a computer-modeled hydrated bilayer of dimyristoyl phosphatidyl choline molecules, and molecular dynamics simulations were run to characterize the effect of this important molecule on membrane structure and dynamics. The effect was judged by observing differences in order parameters, tilt angles, and the fraction of gauche bonds along the hydrocarbon chains between lipids adjacent to cholesterol molecules and comparing them with those further away. It was observed that cholesterol causes an increase in the fraction of trans dihedrals and motional ordering of chains close to the rigid steroid ring system with a decrease in the kink population. The hydrogen-bonding interactions between cholesterol and lipid molecules were determined from radial distribution calculations and showed the cholesterol hydroxyl groups either solvated by water, or forming hydrogen bond contacts with the oxygens of lipid carbonyl and phosphate groups. The dynamics and conformation of the cholesterol molecules were investigated and it was seen that they had a smaller tilt with respect to the bilayer normal than the lipid chains and furthermore that the hydrocarbon tail of the cholesterol was conformationally flexible.     

Helix packing of leucine-rich peptides: a parallel leucine ladder in the structure of Boc-Aib-Leu-Aib-Aib-Leu-Leu-Leu-Aib-Leu-Aib-OMe.

The packing of peptide helices in crystals of the leucine-rich decapeptide Boc-Aib-Leu-Aib-Aib-Leu-Leu-Leu-Aib-Leu-Aib-OMe provides an example of ladder-like leucylleucyl interactions between neighboring molecules. The peptide molecule forms a helix with five 5----1 hydrogen bonds and two 4----1 hydrogen bonds near the C terminus. Three head-to-tail NH ... O = C hydrogen bonds between helices form continuous columns of helices in the crystal. The helicial columns associate in an antiparallel fashion, except for the association of Leu ... Leu side chains, which occurs along the diagonal of the cell where the peptide helices are parallel. The peptide, with formula C56H102N10O13, crystallizes in space group P2(1)2(1)2(1) with Z = 4 and cell parameters a = 16.774(3) A, b = 20.032(3) A and c = 20.117(3) A; overall agreement factor R = 10.7% for 2014 data with magnitude of F(obs) greater than 3 sigma (F); resolution 1.0 A.     

Structure of macrocyclic K+, Rb+-complexon of meso-valinomycin monohydrate, cyclo[-(D-Val-Hyi-Val-D-Hyi)3-].H2O, in a crystalline complex with dioxane by x-ray structural data]

The crystal structure of a valinomycin analogue, cyclo[-(D-Val-Hyi-Val-D-Hyi)3-]x(C60H102N6O18) crystallized with dioxane and water molecules, has been solved by X-ray direct methods. The conformation found is analogous to one established for free meso-valinomycin crystallized from other organic solvents. It is characterized by a centrosymmetric bracelet form, stabilized by six intramolecular 4----1 type hydrogen bonds between amide N-H and C = O groups. One water molecule is fixed asymmetrically by hydrogen bonds in the internal negatively charged cavity of the complexon. The meso-valinomycin molecule "bracelets" in the crystal form stacks alternatively with dioxane molecules.     

Simulation study of a gramicidin/lipid bilayer system in excess water and lipid. I. Structure of the molecular complex

This paper reports on a simulation of a gramicidin channel inserted into a fluid phase DMPC bilayer with 100 lipid molecules. Two lipid molecules per leaflet were removed to insert the gramicidin, so the resulting preparation had 96 lipid molecules and 3209 water molecules. Constant surface tension boundary conditions were employed. Like previous simulations with a lower lipid/gramicidin ratio (Woolf, T. B., and B. Roux. 1996. Proteins: Struct., Funct., Genet. 24:92-114), it is found that tryptophan-water hydrogen bonds are more common than tryptophan-phospholipid hydrogen bonds. However, one of the tryptophan NH groups entered into an unusually long-lived hydrogen bonding pattern with two glycerol oxygens of one of the phospholipid molecules. Comparisons were made between the behavior of the lipids adjacent to the channel with those farther away. It was found that hydrocarbon chains of lipids adjacent to the channel had higher-order parameters than those farther away. The thickness of the lipid bilayer immediately adjacent to the channel was greater than it was farther away. In general, the lipids adjacent to the membrane had similar orientations to those seen by Woolf and Roux, while those farther away had similar orientations to those pertaining before the insertion of the gramicidin. A corollary to this observation is that the thickness of the hydrocarbon region adjacent to the gramicidin was much thicker than what other studies have identified as the "hydrophobic length" of the gramicidin channel.     

Hydrogen bond supramolecular self-assembly in nickel(II) dithiophosphates, Ni[S2P(OR)2]2, R = sec-Bu, iso-Bu, and their bis(pyrazole) adducts

The reaction between nickel(II) nitrate and potassium phosphorus-1,1-dithiolates (di-sec-butyl and di-iso-butyl) in methanol yields 2:1 complexes which were characterized by FT-IR and NMR spectroscopy. 2:1 pyrazole adducts of both compounds were also obtained.The X-ray diffraction analysis of the compounds reveals square planar, four-coordination geometry for the homoleptic compounds and a six-coordinated distorted octahedral geometry for the adducts. In Ni[S₂P(OBu^s)₂]₂ the molecules are associated through C-H?O hydrogen bonds (2.652 Å), and in Ni[S₂P(OBuⁱ)₂]₂ the molecules are associated through C-H?S hydrogen bonds (2.948 Å). The pyrazole adducts are associated through N-H?O bonds and N-H?S bonds from the pyrazole nitrogen atoms, to form supramolecular assemblies. Thus, Ni[S₂P(OBu^s)₂(Pz)₂]₂ (Pz = pyrazole) forms bi-dimensional layers through N-H?O and N-H?S bonds (2.502 and 2.965 Å, respectively), whereas Ni[S₂P(OBuⁱ)₂(Pz)₂]₂ forms linear chains with N-H?S bonds 2.728 Å. The dithiophosphato groups behave as isobidentate chelating ligands.