Synthesis of peptides containing oxo amino acids and their crystallographic analysis

An isolated uncharged hydrogen bond acceptor such as the carbonyl functionality of an aldehyde or a keto group is absent in natural amino acids. Although glutamine and asparagine are known to hydrogen bond through the amide carbonyl group in their side chains, they also possess the amide ? NH₂ group, which can act as a hydrogen bond donor. This makes the structural study of peptides containing an oxo residue, with an isolated carbonyl group in the side chain, interesting. Here, we report the synthesis of δ‐ and ε‐oxo amino acids and their incorporation into oligopeptides as the N‐terminal residue. The resultant oxo peptides were extensively studied using X‐ray crystallography to understand the interactions offered by the oxo group in peptide crystals. We find that the oxo groups are capable of providing additional hydrogen bonding opportunities to the peptides, resulting in increased intermolecular interactions in crystals. The study thus offers avenues for the utilization of oxo residues to introduce intermolecular interactions in synthetic peptides.     

Interaction of the peptide bond with solvent water: a vapor phase analysis

A dynamic technique, using radioactivity as a means of detection, makes it possible to measure the partial pressures of highly polar compounds in dilute aqueous solution. The results can be expressed in terms of the dimensionless distribution coefficient for transfer of a compound from dilute aqueous solution to the vapor phase. For acetic acid this coefficient is 1.1 X 10(-5), for acetamide 7.6 X 10(-8), for N-methylacetamide 4.1 X 10(-8), and for N,N-dimethylacetamide 5.4 X 10(-7). Thus acetamide is much more strongly solvated than the uncharged acetic acid molecule. The results suggest: (1) that the peptide bond represents an extreme among uncharged functional groups in the degree to which it is stabilized by solvent water; (2) that the very great hydrophilic character of the peptide bond may be associated mainly with hydrogen bonding of the solvent to the carbonyl oxygen atom (rather than the N-H group); and (3) that the observed equilibria of biosynthesis and hydrolysis of peptide bonds in aqueous solution are largely determined by differences between reactants and products in their free energies of solvation. It is anticipated that where "bound" water is found in proteins, it will often be found to be associated with peptide bonds, and will tend to be associated with the C-O group rather than with the N-H group.     

Structural basis of DNA flexibility

It has been recently shown by us, on the basis of crystal structure database that the flexibility of B-DNA double helices depends significantly on their base sequence. Our model building studies further indicated that the existence of bifurcated cross-strand hydrogen bonds between successive base pairs is possibly the main factor behind the sequence directed DNA flexibility. These cross-strand hydrogen bonds are, of course, weaker than the usual Watson-Crick hydrogen bonds and their bond geometry is characterized by relatively larger bond lengths and smaller bond angles. We have tried to improve our model structures by incorporating non-planarity of the amino groups in DNA bases due to the presence of lone pair electrons at the nitrogen atoms. Energy minimization studies have been carried out by using different quantum chemical methods, whereby it is found that in all cases of N-H....O type cross-strand hydrogen bonds, the bond geometry improves significantly. In the cases of N-H....N type hydrogen bonds, however, no such consistent improvements can be noticed. Perhaps the true picture would emerge only if all the other interactions present in the DNA macromolecule could be appropriately taken into account.     

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.     

Xanthate sulfur as a hydrogen bond acceptor: the free xanthate anion and ligand sulfur in nickel tris ethylxanthate

Xanthates, like thiolates, form a variety of complexes with metals in which coordinating sulfur can serve as a hydrogen bond acceptor. Nickel tris xanthate complexes [Ni(xan)₃]⁻, (xan = o-ethylxanthate, N-(carbamoylmethyl)ethylxanthate) have been synthesized and compared by a combination of X-ray crystallographic and spectroscopic measurements. Recent results from our studies of N-H?S hydrogen bonding interactions in metal-xanthate complexes shows N-S distances to be longer than those in related thiolate complexes, indicative of weaker hydrogen bonds for the xanthates. The complex (Et₄N)[N-(carbamoylmethyl)ethylxanthate)] adopts an extended conformation in both the solid state and solution and lacks either intraligand or intermolecular N-H?S hydrogen bonds. The complex (CTA)[Ni(exa)₃] exhibits N-H?S hydrogen bonds between the amide group of the counterion and the ligand sulfur. The amide-sulfur N-H?S distance is 3.567 Å.     

Molecular packing and intermolecular interactions in N-acylethanolamines: crystal structure of N-myristoylethanolamine.

N-Acylethanolamines elicited much interest in recent years owing to their occurrence in biological membranes under conditions of stress as well as under normal conditions. The molecular conformation, packing properties and intermolecular interactions of N-myristoylethanolamine (NMEA) have been determined by single crystal X-ray diffraction analysis. The lipid crystallized in the space group P21/a with unit cell dimensions: a=9.001, b=4.8761, c=39. 080. There are four symmetry-related molecules in the monoclinic unit cell. The molecules are organized in a tail-to-tail fashion, similar to the arrangement in a bilayer membrane. The hydrophobic acyl chain of the NMEA molecule is tilted with respect to the bilayer normal by an angle of 37 degrees. Each hydroxy group forms two hydrogen bonds, one as a donor and the other as an acceptor, with the hydroxy groups of molecules in the opposing leaflet. These O-H...O hydrogen bonds form an extended, zig-zag type network along the b-axis. In addition, the N-H and C=O groups of adjacent molecules are involved in N-H...O hydrogen bonds, which also connect adjacent molecules along the b-axis.     

Lectin-carbohydrate interactions. Studies of the nature of hydrogen bonding between D-galactose and certain D-galactose-specific lectins, and between D-mannose and concanavalin A

The binding of galactose-specific lectins from Erythrina indica (EIL), Erythrina arborescens (EAL), Ricinus communis (agglutinin; RCA-I), Abrus precatorius (agglutinin; APA), and Bandeiraea simplicifolia (lectin I; BSL-I) to fluoro-, deoxy-, and thiogalactoses were studied in order to determine the strength of hydrogen bonds between the hydroxyl groups of galactose and the binding sites of the proteins. The results have allowed insight into the nature of the donor/acceptor groups in the lectins that are involved in hydrogen bonding with the sugar. The data indicate that the C-2 hydroxyl group of galactose is involved in weak interactions as a hydrogen-bond acceptor with uncharged groups of EIL and EAL. With RCA-I, the C-2 hydroxyl group forms two weak hydrogen bonds in the capacity of a hydrogen-bond acceptor and a donor. On the other hand, there is a strong hydrogen bond between the C-2 hydroxyl group of galactose, which acts as a donor, and a charged group on BSL-I. The C-2 hydroxyl group of the sugar is also a hydrogen-bond donor to APA. The lectins are involved in strong hydrogen bonds through charged groups with the C-3 and C-4 hydroxyl groups of galactose, with the latter serving as hydrogen-bond donors. The C-6 hydroxyl group of the sugar is weakly hydrogen bonded with neutral groups of EIL, EAL, and APA. With BSL-I, however, a strong hydrogen bond is formed at this position with a charged group of the lectin. The C-6 hydroxyl groups is a hydrogen-bond acceptor for EIL and EAL, a hydrogen-bond donor for APA and BSL-I, and appears not to be involved in binding to RCA-I. The data with the thiosugars indicate the involvement of the C-1 hydroxyl group of galactose in binding to EIL, EAL, and BSL-I, but not to RCA-I and APA. We have also performed a similar analysis of the binding data of fluoro- and deoxysugars to concanavalin A [Poretz, R. D. and Goldstein, I. J. (1970) Biochemistry 9, 2890-2896]. This has allowed comparison of the donor/acceptor properties and free energies of hydrogen bonding of the hydroxyl groups of methyl alpha-D-mannopyranoside to concanavalin A with the results in the present study. On the basis of this analysis, new assignments are suggested for amino acid residues of concanavalin A [corrected] that may be involved in hydrogen bonding to the sugar.     

A review of adsorption of amino acids on minerals: Was it important for origin of life?

Minerals more readily adsorb amino acids with charged R groups than uncharged R groups, so that the incorporation of amino acids with charged R groups into peptides would be more frequent than for amino acids with uncharged R groups. However, 74% of the amino acids in the proteins of modern organisms contain uncharged R groups. Thus, what could have been the mechanisms that produced peptides/proteins with more amino acids with uncharged R groups than precursors with charged R groups? Should we expect the composition of amino acids adsorbed on minerals to be similar to those of present proteins? Was the adsorption of amino acids on minerals important for the origin of life? The lipid world offers an alternative view of origin of life. Liposomes contributed to elongation of peptides as well as select hydrophobic amino acids and peptides. These experiments could be showing the mechanism, which hydrophobic amino acids have been selected. However, liposomes have no influence on the stereoselectivity in the oligomerization of amino acids. In the present paper, several other mechanisms are also discussed that could produce peptides with a greater proportion of amino acids with uncharged R groups.     

Crystal structure of oxidized flavodoxin from a red alga Chondrus crispus refined at 1.8 A resolution. Description of the flavin mononucleotide binding site.

In order to describe the detailed conformation of the oxidized flavodoxin from a eukaryotic red alga, Chondrus crispus, the crystal structure has been refined by a restrained least-squares method. The crystallographic R factor is 0.168 for 13,899 reflections with F greater than 2 sigma F between 6.0 and 1.8 A resolution. The refined model includes 173 amino acid residues, flavin mononucleotide (FMN) and 110 water molecules. The root-mean-square deviation in bond lengths from ideal values is 0.015 A, and the mean co-ordinate error is estimated to be 0.2 A. The FMN is located at the periphery of the molecule. The orientation of the isoalloxazine ring is such that the C-7 and C-8 methyl groups are exposed to solvent and the pyrimidine moiety is buried in the protein. Three peptide segments, T8-T13, T55-T58 and D94-C103, are involved in FMN binding. The first segment of T8-T13 enfolds the phosphate group of the FMN. The three oxygen atoms in the phosphate group form extensive hydrogen bonds with amide groups of the main chain and the O gamma atoms of the side-chains in this segment. T55 O and W56 N epsilon 1 in the second segment form hydrogen bonds with O-2 in the ribityl moiety and one of the oxygen atoms in the phosphate group, respectively. The O gamma H of T58 forms a hydrogen bond with the N-5 atom in the isoalloxazine ring, which is expected to be protonated in the semiquinone form. The third segment is in contact with the isoalloxazine ring. It appears that the hydrogen bond acceptor of the NH of Asp94 in the third segment is O-2 rather than N-1 in the isoalloxazine ring. The isoalloxazine ring is flanked by the side-chains of Trp56 and Tyr98; it forms an angle of 38 degrees with the indole ring of Trp56 and is almost parallel to the benzene ring of Tyr98. The environment of the phosphate group is conserved as in other flavodoxins whereas that of the isoalloxazine ring differs. The relationship between the hydrogen bond to the N-5 in the ring and the redox potential for the oxidized/semiquinone couple is discussed.     

Helix stability of oligoglycine,oligoalanine, and oligo‐β‐alanine dodecamers reflected by hydrogen‐bond persistence

Helices are important structural/recognition elements in proteins and peptides. Stability and conformational differences between helices composed of α‐ and β‐amino acids as scaffolds for mimicry of helix recognition has become a theme in medicinal chemistry. Furthermore, helices formed by β‐amino acids are experimentally more stable than those formed by α‐amino acids. This is paradoxical because the larger sizes of the hydrogen‐bonding rings required by the extra methylene groups should lead to entropic destabilization. In this study, molecular dynamics simulations using the second‐generation force field, AMOEBA (Ponder, J.W., et al., Current status of the AMOEBA polarizable force field. J Phys Chem B, 2010. 114 (8): p. 2549–64.) explored the stability and hydrogen‐bonding patterns of capped oligo‐β‐alanine, oligoalanine, and oligoglycine dodecamers in water. The MD simulations showed that oligo‐β‐alanine has strong acceptor+2 hydrogen bonds, but surprisingly did not contain a large content of 3₁₂‐helical structures, possibly due to the sparse distribution of the 3₁₂‐helical structure and other structures with acceptor+2 hydrogen bonds. On the other hand, despite its backbone flexibility, the β‐alanine dodecamer had more stable and persistent <3.0 Å hydrogen bonds. Its structure was dominated more by multicentered hydrogen bonds than either oligoglycine or oligoalanine helices. The 3₁ (PII) helical structure, prevalent in oligoglycine and oligoalanine, does not appear to be stable in oligo‐β‐alanine indicating its competition with other structures (stacking structure as indicated by MD analyses). These differences are among the factors that shape helical structural preferences and the relative stabilities of these three oligopeptides. Proteins 2014; 82:3043–3061. © 2014 Wiley Periodicals, Inc.     

Hydrogen bonding in nucleosides and nucleotides

An analysis of the hydrogen bonding in 76 nucleoside and 11 nucleotide crystal structures shows that the hydrogen bond lengths fall into well-defined categories according to the nature of the donor or acceptor groups. The shortest bonds are those involving P---OH or O=P groups. For donor groups, the sequence in bond lengths is

P—OH<C—OH< N−H<O_w(H)—H<N(H)—H<C−H

There are ten examples of two centre

H—H…O

bonds, which are comparable in length with P---OH …O bonds. The acceptor seqeunce is

O=P<OH₂<OH₂<O=CO(H)C<N N(H₂)C<Cl^…<O<S=C

The number of three-centre bonds, about 24%, is comparable to that observed in the carbohydrates and the amino acids. Most hydrogen bonds are involved in short finite chains. Only in the nucleotides are cyclic hydrogen bonding schemes observed.     

Synthesis, X-ray crystallographic structures of thio substituted N-acetyl N'-methylamide alanine and evaluation of sp sulfur parameters of the CFF91 force field.

Acetyl thioalanine N-methyl (Ac-Alat-NHMe) and thioacetyl alanine N-methyl (Act-Ala-NHMe) were synthesized, crystallized and their X-ray diffraction structures determined for the first time. Both molecules adopted beta-sheet conformations and showed similar hydrogen bonding patterns with one molecular surface forming two oxo hydrogen bonds and the other forming two thio hydrogen bonds. The crystal structure data for the two thioamides provided a validation of the thioamide parameters for the newly derived CFF91 force field because the observed crystal (phi, psi) angles were situated in the global minimum regions of the theoretical (phi, psi) map predicted using the parameters. In addition, the parameters were further validated because conformational energy minimization of the crystal structure produced low deviations in unit cell dimensions, bond lengths, bond angles and torsional angles, and a 120-ps molecular dynamics simulation also gave a low deviation for the most probable N-H...S=C bond distance.     

How alcohol chain-length and concentration modulate hydrogen bond formation in a lipid bilayer

Molecular dynamics simulations are used to measure the change in properties of a hydrated dipalmitoylphosphatidylcholine bilayer when solvated with ethanol, propanol, and butanol solutions. There are eight oxygen atoms in dipalmitoylphosphatidylcholine that serve as hydrogen bond acceptors, and two of the oxygen atoms participate in hydrogen bonds that exist for significantly longer time spans than the hydrogen bonds at the other six oxygen atoms for the ethanol and propanol simulations. We conclude that this is caused by the lipid head group conformation, where the two favored hydrogen-bonding sites are partially protected between the head group choline and the sn-2 carbonyl oxygen. We find that the concentration of the alcohol in the ethanol and propanol simulations does not have a significant influence on the locations of the alcohol/lipid hydrogen bonds, whereas the concentration does impact the locations of the butanol/lipid hydrogen bonds. The concentration is important for all three alcohol types when the lipid chain order is examined, where, with the exception of the high-concentration butanol simulation, the alcohol molecules having the longest hydrogen-bonding relaxation times at the favored carbonyl oxygen acceptor sites also have the largest order in the upper chain region. The lipid behavior in the high-concentration butanol simulation differs significantly from that of the other alcohol concentrations in the order parameter, head group rotational relaxation time, and alcohol/lipid hydrogen-bonding location and relaxation time. This appears to be the result of the system being very near to a phase transition, and one occurrence of lipid flip-flop is seen at this concentration.     

Atoms-in-molecules study of the genetically encoded amino acids. II. Computational study of molecular geometries

The geometries of the 20 genetically encoded amino acids were optimized at the restricted Hartree-Fock level of theory using the 6-31+G* basis set. A detailed comparison showed the calculated geometries to be in excellent agreement with those determined by X-ray crystallography. The study demonstrated that the geometric parameters for the main-chain group and for the bonds and common functional groups of the side-chains exhibit a high degree of transferability among the members of this set of molecules. This geometric transferability is a necessary prerequisite for the corresponding transferability of their electron density distributions and hence of their bond and atomic properties. The transferability of the electron distributions will be demonstrated and exploited in the following paper of this series, which uses the topology of the electron density to define an atom within the quantum theory of atoms in molecules. Particular features of the geometries of the amino acids are discussed. It has been shown, for example, how the apparent anomaly of the Calpha-N bond length in a peptide being shorter than in the charged species Calpha-NH3+ is resolved when the charge separation is gauged by the differences in the charges of the Calpha and N atoms as opposed to the use of formal charges. A compilation of literature sources on experimental geometries covering each member of the 20 amino acids is presented. A set of rules for labeling the atoms and bonds, complementing the generally accepted IUPAC-IUB rules, is proposed to uniquely identify every atom and bond in the amino acids.     

A survey of hydrogen bond geometries in the crystal structures of amino acids

An analysis of the geometries of the hydrogen bonds observed by neutron diffraction in thirt-two crystal structures of amino acids shows the following results. Of the 168 hydrogen bonds in the data set, 64 involve the zwitterion groups 

and $\text{C}$O₂. Another 18 are from

to sulphate or carbonyl oxygens. The majority, 46, of these

H … O bonds are three-centered (bifurcated). Nine are four-centered (trifurcated). The geometry in which the three-centered hydrogen bond involves both oxygens of the same carboxylate group is not especially favoured. When it does occur, one hydrogen bond is generally shorter and the other longer, than when the bonding involves oxygens on different carboxylate groups. The shortest hydrogen bonds are the OH … O $\text{C}$, from a carboxylic acid hydroxyl to a carboxylate oxygen, and NH … O$\text{C}$ when the nitrogen is the ring atom in histidine or proline. Carboxylate groups, on average, accept six hydrogen bonds, with no examples of less than four bonds. The reason for the large number of three-centered

H … O$\text{C}$ bonds is therefore a proton deficiency arising from the disparity between the tripled donor property of the

groups and the sextuple, on average, acceptor property of the carboxylate groups. There is good geometrical evidence for the existence of H … O and H … Cl^? hydrogen bonds, especially involving the hydrogen atoms on α-atoms.     

Phylogenetic analysis of condensation domains in the nonribosomal peptide synthetases

Condensation (C) domains in the nonribosomal peptide synthetases are capable of catalyzing peptide bond formation between two consecutively bound various amino acids. C-domains coincide in frequency with the number of peptide bonds in the product peptide. In this study, a phylogenetic approach was used to investigate structural diversity of bacterial C-domains. Phylogenetic trees show that the C-domains are clustered into three functional groups according to the types of substrate donor molecules. They are l-peptidyl donors, d-peptidyl donors, and N-acyl donors. The fact that C-domain structure is not subject to optical configuration of amino acid acceptor molecules supports an idea that the conversion from l to d-form of incorporating amino acid acceptor occurs during or after peptide bond formation. l-peptidyl donors and d-peptidyl donors are suggested to separate before separating the lineage of Gram-positive and Gram-negative bacteria in the evolution process.     

Simulation of Interactions between Nucleic Acid Bases by Refined Atom-Atom Potential Functions

Abstract

Energy of interaction between nitrogen bases of nucleic acids has been calculated as a function of parameters determining the mutual position of two bases. Refined atom-atom potential functions are suggested. These functions contain terms proportional to the first (electrostatics), sixth (or tenth for the atoms forming a hydrogen bond) and twelfth (repulsion of all atoms) powers of interatomic distance. Calculations have shown that there are two groups of minima of the base interaction energy. The minima of the first group correspond to coplanar arrangement of the base pairs and hydrogen bond formation. The minima of the second group correspond to the position of bases one above the other in almost parallel planes. There are 28 energy minima corresponding to the formation of coplanar pairs with two (three for the G:C pair) almost linear N-H … O and (or) N-H … N hydrogen bonds. The position of nitrogen bases paired by two such H-bonds in any crystal of nucleic acid component, in polynucleotide complexes and in tRNA is close to the position in one of these minima. Besides, for each pair there are energy minima corresponding to the formation of a single N-H … O or N-H … N and one C-H … O or C-H … N hydrogen bond. The form of potential surface in the vicinity of minima has been characterized. The results of calculations agree with the experimental data and with more rigorous calculations based on quantum- mechanical approach.     

The specific guanine binding site of the ribonuclease T1 family enzymes and of G-proteins is modeled in the cocrystal formed by 7-methylguanosine-5'-phosphate and phenylalanine

In the cocrystal formed by 7-methylguanosine-5'-phosphate.phenylalanine.6H2O, the interactions between guanine and phenylalanine are similar to those observed in the complex of ribonuclease T1 with 2'-guanylic acids, and those of the two G-proteins, Elongation Factor-Tu and ras oncogene p21, with GDP. They are similar in the following three points: (a) guanine N(1)H and N(2)H donate cyclic N-H...O hydrogen bonds to the carboxylate group of phenylalanine in the former cocrystal and to the side chain carboxylate group of Asp or Glu in the latter proteins, (b) O(6) of guanine accepts hydrogen bond(s) from main-chain NH group(s), and (c) the purine moiety is sandwiched between aromatic (or hydrophobic) amino acid side chains.     

Exhaustive Metropolis Monte Carlo sampling and analysis of polyalanine conformations adopted under the influence of hydrogen bonds

We propose a novel Metropolis Monte Carlo procedure for protein modeling and analyze the influence of hydrogen bonding on the distribution of polyalanine conformations. We use an atomistic model of the polyalanine chain with rigid and planar polypeptide bonds, and elastic alpha carbon valence geometry. We adopt a simplified energy function in which only hard-sphere repulsion and hydrogen bonding interactions between the atoms are considered. Our Metropolis Monte Carlo procedure utilizes local crankshaft moves and is combined with parallel tempering to exhaustively sample the conformations of 16-mer polyalanine. We confirm that Flory's isolated-pair hypothesis (the steric independence between the dihedral angles of individual amino acids) does not hold true in long polypeptide chains. In addition to 3(10)- and alpha-helices, we identify a kink stabilized by 2 hydrogen bonds with a shared acceptor as a common structural motif. Varying the strength of hydrogen bonds, we induce the helix-coil transition in the model polypeptide chain. We compare the propensities for various hydrogen bonding patterns and determine the degree of cooperativity of hydrogen bond formation in terms of the Hill coefficient. The observed helix-coil transition is also quantified according to Zimm-Bragg theory.