Strength of Calpha-H...O=C hydrogen bonds in transmembrane proteins

A large number of Calpha-H...O contacts are present in transmembrane protein structures, but contribution of such interactions to protein stability is still not well understood. According to previous ab initio quantum calculations, the stabilization energy of a Calpha-H...O contact is about 2-3 kcal/mol. However, experimental studies on two different Calpha-H...O hydrogen bonds present in transmembrane proteins lead to conclusions that one contact is only weakly stabilizing and the other is not even stabilizing. We note that most previous computational studies were on optimized geometries of isolated molecules, but the experimental measurements were on those in the structural context of transmembrane proteins. In the present study, 263 Calpha-H...O=C contacts in alpha-helical transmembrane proteins were extracted from X-ray crystal structures, and interaction energies were calculated with quantum mechanical methods. The average stabilization energy of a Calpha-H...O=C interaction was computed to be 1.4 kcal/mol. About 13% of contacts were stabilizing by more than 3 kcal/mol, and about 11% were destabilizing. Analysis of the relationships between energy and structure revealed four interaction patterns: three types of attractive cases in which additional Calpha-H...O or N-H...O contact is present and a type of repulsive case in which repulsion between two carbonyl oxygen atoms occur. Contribution of Calpha-H...O=C contacts to protein stability is roughly estimated to be greater than 5 kcal/mol per helix pair for about 16% of transmembrane helices but for only 3% of soluble protein helices. The contribution would be larger if Calpha-H...O contacts involving side chain oxygen were also considered.     

超嗜热酯酶APE1547中特殊位置氢键对酶活力和热稳定性的影响

探讨了APE1547蛋白的β-推进器结构中第3和第4“叶片”间的侧链氢键(Thr127-Gly154, Leu182-Arg145-Glu122)对蛋白质的作用.     

Experimental measurement of the strength of a C alpha-H...O bond in a lipid bilayer

Due to the apolar nature of the lipid bilayer, the weak Calpha-H...O H-bond is thought to contribute significantly toward the stability of transmembrane helical bundles such as glycophorin A (GPA). Here for the first time we measured the strength of such a bond, using vibrational frequency shifts of a dimeric and nondimeric variants of GPA containing a Gly CD2 label. Although the resulting estimated bond strength of 0.88 kcal/mol is relatively weak, several such bonds could contribute significantly toward bundle stabilization.     

Solution and solid-state models of peptide CH...O hydrogen bonds

Fumaramide derivatives were analyzed in solution by (1)H NMR spectroscopy and in the solid state by X-ray crystallography in order to characterize the formation of CH...O interactions under each condition and to thereby serve as models for these interactions in peptide and protein structure. Solutions of fumaramides at 10 mM in CDCl(3) were titrated with DMSO-d(6), resulting in chemical shifts that moved downfield for the CH groups thought to participate in CH...O=S(CD(3))(2) hydrogen bonds concurrent with NH...O=S(CD(3))(2) hydrogen bonding. In this model, nonparticipating CH groups under the same conditions showed no significant change in chemical shifts between 0.0 and 1.0 M DMSO-d(6) and then moved upfield at higher DMSO-d(6) concentrations. At concentrations above 1.0 M DMSO-d(6), the directed CH...O=S(CD(3))(2) hydrogen bonds provide protection from random DMSO-d(6) contact and prevent the chemical shifts for participating CH groups from moving upfield beyond the original value observed in CDCl(3). X-ray crystal structures identified CH...O=C hydrogen bonds alongside intermolecular NH...O=C hydrogen bonding, a result that supports the solution (1)H NMR spectroscopy results. The solution and solid-state data therefore both provide evidence for the presence of CH...O hydrogen bonds formed concurrent with NH...O hydrogen bonding in these structures. The CH...O=C hydrogen bonds in the X-ray crystal structures are similar to those described for antiparallel beta-sheet structure observed in protein X-ray crystal structures.     

Modeling competitive interactions in proteins: vibrational spectroscopy of M+(n-methylacetamide)1(H2O)n=0-3, M=Na and K, in the 3 microm region

To properly understand the preferred structures and biological properties of proteins, it is important to understand how they are influenced by their immediate environment. Competitive intrapeptide, peptide...water, ion...water, and ion...peptide interactions, such as hydrogen bonding, play a key role in determining the structures, properties, and functionality of proteins. The primary types of hydrogen bonding involving proteins are intramolecular amide...amide (N-H...O=C) and intermolecular amide...water (O-H...O=C and H-O...H-N). n-Methylacetamide (NMA) is a convenient model for investigating these competitive interactions. An analysis of the IR photodissociation (IRPD) spectra of M+(n-methylacetamide)1(H2O)n=0-3 (M=Na and K) in the O-H and N-H spectral regions is presented. Ab initio calculations (MP2/cc-pVDZ) are used as a guide in identifying both the type and location of hydrogen bonds present. In larger clusters, where several structural isomers may be present in the molecular beam, ab initio calculations are also used to suggest assignments for the observed spectral features. The results presented offer insight to the nature of ion...NMA interactions in an aqueous environment and reveal how different ion...ligand pairwise interactions direct the extent of water...water and water...NMA hydrogen bonding observed.     

Calculation of relative binding affinities of fructose 1,6-bisphosphatase mutants with adenosine monophosphate using free energy perturbation method

The free energy perturbation (FEP) methodology is the most accurate means of estimating relative binding affinities between inhibitors and protein variants. In this article, the importance of hydrophobic and hydrophilic residues to the binding of adenosine monophosphate (AMP) to the fructose 1,6-bisphosphatase (FBPase), a target enzyme for type-II diabetes, was examined by FEP method. Five mutations were made to the FBPase enzyme with AMP inhibitor bound: 113Tyr --> 113Phe, 31Thr --> 31Ala, 31Thr --> 31Ser, 177Met --> 177Ala, and 30Leu --> 30Phe. These mutations test the strength of hydrogen bonds and van der Waals interactions between the ligand and enzyme. The calculated relative free energies indicated that: 113Tyr and 31Thr play an important role, each via two hydrogen bonds affecting the binding affinity of inhibitor AMP to FBPase, and any changes in these hydrogen bonds due to mutations on the protein will have significant effect on the binding affinity of AMP to FBPase, consistent to experimental results. Also, the free energy calculations clearly show that the hydrophilic interactions are more important than the hydrophobic interactions of the binding pocket of FBPase.     

Correlation of protein structure and dynamics to scalar couplings across hydrogen bonds

NMR-observable scalar couplings across hydrogen bonds in nucleic acids and proteins present a quantitative measure for the geometry and--by the implicit experimental time averaging--dynamics of hydrogen bonds. We have carried out in-depth molecular dynamics (MD) simulations with various force fields on three proteins: ubiquitin, the GB1 domain of protein G, and the SMN Tudor domain, for which experimental h3JNC' scalar couplings of backbone hydrogen bonds and various high-resolution X-ray structures are available. Theoretical average values for h3JNC' were calculated from the snapshots of these MD simulations either by density functional theory or by a geometric parametrization (Barfield, M. J. Am. Chem. Soc. 2002, 124, 4158-4168). No significant difference was found between the two methods. The results indicate that time-averaging using explicit water solvation in the MD simulations improves significantly the agreement between experimental and theoretical values for the lower resolution structures ubiquitin (1.8 A), Tudor domain (1.8 A), and protein G (2.1 A). Only marginal improvement is found for the high-resolution structure (1.1 A) of protein G. Hence, experimental h3JNC' values are compatible with a static, high-resolution structural model. The MD averaging of the low-resolution structures moves the averages of the rHO distance and the H...O=C angle theta closer to their respective values in the high-resolution structures, thereby improving the agreement using experimental h3JNC' data. In contrast, MD averaging with implicit water models deteriorates the agreement with experiment for all proteins. The differing behavior can be explained by an artifactual lengthening of H-bonds caused by the implicit water models.     

Structural scaffold of 18-crown-6 tetracarboxylic acid for optical resolution of chiral amino acid: X-ray crystal analyses and energy calculations of complexes of D- and L-isomers of tyrosine, isoleucine, methionine and phenylglycine

To clarify the structural scaffold of (+)-18-crown-6 tetracarboxylic acid ((+)-18C6H4) for the optical resolution of a chiral amino acid, the crystal structures of its equimolar complexes with L- and D-isomers of tyrosine (Tyr), isoleucine (Ile), methionine (Met) and phenylglycine (PheG) were analysed by X-ray diffraction methods. (+)-18C6H4 took very similar conformations for all complexes. Although the chemical structure of (+)-18C6H4 is C2-symmetric, it took a similar asymmetric ring conformation of radius ca. 6.0 A. In all complexes, the amino group of chiral amino acids was located near the center of the ring and formed three hydrogen bonds and five electrostatic interactions with eight oxygen atoms of the ether ring and carboxyl groups. Also, the Calpha atom of chiral amino acids participated in Calpha-H...O interaction with the oxygen atom of (+)-18C6H4. In contrast, the carboxyl group of chiral amino acids did not directly interact with (+)-18C6H4. These results indicate that the structural scaffold of (+)-18C6H4 for the optical resolution of chiral amino acids is mainly based on the mode of interaction of (+)-18C6H4 with the amino and Calpha-H groups of chiral amino acids. The differences in interaction pattern and binding energy between the L- and D-isomers of each amino acid are discussed in relation to the chiral recognition of (+)-18C6H4.     

High-field 1H MAS and 15N CP-MAS NMR studies of alanine tripeptides and oligomers: distinction of antiparallel and parallel beta-sheet structures and two crystallographically independent molecules

beta-Strand peptides are known to assemble into either antiparallel (AP) or parallel (P) beta-sheet forms which are very important motifs for protein folding and fibril formations occurring in silk fibroin or amyloid proteins. Well-resolved 1H NMR signals including NH protons were observed for alanine tripeptides (Ala)3 with the AP and P structures as well as (Ala)n (n = 4-6) by high-field/fast magic-angle spinning NMR. Amide NH and amino NH3+ 1H signals of (Ala)3 with the P structure were well resonated at 7.5 and 8.9 ppm, respectively, whereas they were not resolved for the AP structure. Notably, NH 1H signals of (Ala)3 and (Ala)4 taking the P structure are resonated at higher field than those of the AP structure by 1.0 and 1.1 ppm, respectively. Further, NH 15N signals of (Ala)3 with the AP structure were resonated at lower field by 2 to 5 ppm than those of (Ala)3 with the P structure. These relative 1H and 15N hydrogen bond shifts of the P structure with respect to those of the AP structure are consistent with the relative hydrogen bond lengths of the interstrand N-H...O=C bonds. Distinction between the two crystallographically independent chains present in the AP and P structures was feasible by 15N chemical shifts but not by 1H chemical shifts because of insufficient spectral resolution in the latter. Calculated 1H and 15N shielding constants by density functional theory are generally consistent with the experimental data, although some discrepancies remain depending upon the models used.     

Hydrogen Bond Fluctuations Control Photochromism in a Reversibly Photo‐Switchable Fluorescent Protein

Reversibly switchable fluorescent proteins (RSFPs) are essential for high‐resolution microscopy of biological samples, but the reason why these proteins are photochromic is still poorly understood. To address this problem, we performed molecular dynamics simulations of the fast switching Met159Thr mutant of the RSFP Dronpa. Our simulations revealed a ground state structural heterogeneity in the chromophore pocket that consists of three populations with one, two, or three hydrogen bonds to the phenolate moiety of the chromophore. By means of non‐adiabatic quantum mechanics/molecular dynamics simulations, we demonstrated that the subpopulation with a single hydrogen bond is responsible for off‐switching through photo‐isomerization of the chromophore, whereas two or more hydrogen bonds inhibit the isomerization and promote fluorescence instead. While rational design of new RSFPs has so far focused on structure alone, our results suggest that structural heterogeneity must be considered as well.     

Contribution of multicentered short hydrogen bond arrays to potency of active site-directed serine protease inhibitors

We describe and compare the pH dependencies of the potencies and of the bound structures of two inhibitor isosteres that form multicentered short hydrogen bond arrays at the active sites of trypsin, thrombin, and urokinase type plasminogen activator (urokinase or uPA) over certain ranges of pH. Depending on the pH, short hydrogen bond arrays at the active site are mediated by two waters, one in the oxyanion hole (H(2)O(oxy)) and one on the other (S2) side of the inhibitor (H(2)O(S2)), by one water (H(2)O(oxy)), or by no water. The dramatic variation in the length of the active site hydrogen bonds as a function of pH, of inhibitor, and of enzyme, along with the involvement or absence of ordered water, produces a large structural manifold of active site hydrogen bond motifs. Diverse examples of multicentered and two-centered short hydrogen bond arrays, both at and away from the active site, recently discovered in several protein crystal systems, suggest that short hydrogen bonds in proteins may be more common than has been recognized. The short hydrogen bond arrays resemble one another with respect to ionic nature, highly polar environment, multitude of associated ordinary hydrogen bonds, and disparate pK(a) values of participating groups. Comparison of structures and K(i) values of trypsin complexes at pH values where the multicentered short hydrogen bond arrays mediating inhibitor binding are present or absent indicate that these arrays have a minor effect on inhibitor potency. These features suggest little covalent nature within the short hydrogen bonds, despite their extraordinary shortness (as short as 2.0 A).     

Dynamics of noncovalent interactions in all-α and all-β class proteins: implications for the stability of amyloid aggregates

A fully folded functional protein is stabilized by several noncovalent interactions. When a protein undergoes conformational motions, the existing noncovalent interactions may be maintained. They may also break or new interactions may be formed. Knowledge of the dynamical nature of the different types of noncovalent interactions is extremely important to understand the structural stability, function, and folding of a protein. There are experimental limitations to investigate the dynamics of different noncovalent interactions simultaneously in a biomolecule. We have carried out molecular dynamics simulations on four different proteins, two belonging to all-α class proteins and the other two are representatives of all-β class proteins. The dynamical nature of eight different noncovalent interactions was studied by monitoring the maximum residence time (MRT) and lifetime (LT). The conventional hydrogen bonds are the dominant interactions in all four proteins, and the majority of those formed between the main-chain atoms were maintained during most of the simulation time with MRT greater than 10 ns. Such interactions with more than 1 ns lifetime provide stability to the secondary structures, and hence they are responsible for the overall stability of the protein. The weak C-H···O hydrogen bond is the next major type of interactions. However, a large number of such interactions are observed between the main-chain atoms only in all-β proteins as interstrand interactions, and, surprisingly, they are observed during most part of the simulation although their average lifetime is only about 20 to 30 ps. The strong cation···π and salt-bridge interactions are present few in number. However, in many cases they are almost uninterrupted indicating the higher strength of these interactions. Four other interactions involving the π-electron cloud of aromatic rings are very small in number, and, in many cases, their presence is not maintained throughout the simulation. Our results clearly indicate that the weak C-H···O interactions between the main-chain atoms are the distinguishing factor between the all-α and all-β class of proteins, and these interstrand interactions can provide additional stability to all-β protein structures. Based on these results, we hypothesize that such weak C-H···O interstrand interactions could play a major role in providing stability to amyloid type of aggregates that are responsible for the pathological state of many proteins.     

Microsolvation of formamide: a rotational study

Microsolvated formamide clusters have been generated in a supersonic jet expansion and characterized using Fourier transform microwave spectroscopy. Three conformers of the monohydrated cluster and one of the dihydrated complex have been observed. Seven monosubstituted isotopic species have been measured for the most stable conformer of formamide...H(2)O, which adopts a closed planar ring structure stabilized by two intermolecular hydrogen bonds (N-H...O(H)-H...O=C). The two higher energy forms of formamide...H(2)O have been observed for the first time. The second most stable conformer is stabilized by a O-H...O=C and a weak C-H...O hydrogen bond, while, in the less stable form, water accepts a hydrogen bond from the anti hydrogen of the amino group. For formamide...(H(2)O)(2), the parent and nine monosubstituted isotopic species have been observed. In this cluster the two water molecules close a cycle with the amide group through three intermolecular hydrogen bonds (N-H...O(H)-H...O(H)-H...O=C), the nonbonded hydrogen atoms of water adopting an up-down configuration. Substitution (r(s)) and effective (r(0)) structures have been determined for formamide, the most stable form of formamide...H(2)O and formamide...(H(2)O)(2). The results on monohydrated formamide clusters can help to explain the observed preferences of bound water in proteins. Clear evidence of sigma-bond cooperativity effects emerges when comparing the structures of the mono- and dihydrated formamide clusters. No detectable structural changes due to pi-bond cooperativity are observed on formamide upon hydration.     

Crystallographic characterization of helical secondary structures in alpha/beta-peptides with 1:1 residue alternation

Oligomers that contain both alpha- and beta-amino acid residues in a 1:1 alternating pattern have recently been shown by several groups to adopt helical secondary structures in solution. The beta-residue substitution pattern has a profound effect on the type of helix formed and the stability of the helical conformation. On the basis of two-dimensional NMR data, we have previously proposed that beta-residues with a five-membered ring constraint promote two different types of alpha/beta-peptide helix. The "11-helix" contains i, i+3 CO...H-N hydrogen bonds between backbone amide groups; these hydrogen bonds occur in 11-atom rings. The alpha/beta-peptide "14/15-helix" contains i, i+4 CO...H-N hydrogen bonds, which occur in alternating 14- and 15-atom rings. Here we provide crystallographic data for 14 alpha/beta-peptides that form the 11-helix and/or the 14/15-helix. These results were obtained for a series of oligomers containing beta-residues derived from ( S,S)- trans-2-aminocyclopentanecarboxylic acid (ACPC) and alpha-residues derived from alpha-aminoisobutyric acid (Aib) or l-alanine (Ala). The crystallized alpha/beta-peptides range in length from 4 to 10 residues. Nine of the alpha/beta-peptides display the 11-helix in the solid state, three display the 14/15-helix, and two display conformations that contain both i, i+3 and i, i+4 CO...H-N hydrogen bonds, but not bifurcated hydrogen bonds. Only 3 of the 14 crystal structures presented here have been previously described. These results suggest that longer alpha/beta-peptides prefer the 14/15-helix over the 11-helix, a conclusion that is consistent with previously reported NMR data obtained in solution.     

Recurrence of carboxylic acid-pyridine supramolecular synthon in the crystal structures of some pyrazinecarboxylic acids.

X-ray crystal structures of pyrazinic acid 1 and isomeric methylpyrazine carboxylic acids 2-4 are analyzed to examine the occurrence of carboxylic acid-pyridine supramolecular synthon V in these heterocyclic acids. Synthon V, assembled by (carboxyl)O-H...N(pyridine) and (pyridine)C-H...O(carbonyl) hydrogen bonds, controls self-assembly in the crystal structures of pyridine and pyrazine monocarboxylic acids. The recurrence of acid-pyridine heterodimer V compared to the more common acid-acid homodimer I in the crystal structures of pyridine and pyrazine monocarboxylic acids is explained by energy computations in the RHF 6-31G* basis set. Both the O-H.N and the C-H...O hydrogen bonds in synthon V result from activated acidic donor and basic acceptor atoms in 1-4. Pyrazine 2,3- and 2,5-dicarboxylic acids 10 and 11 crystallize as dihydrates with a (carboxyl)O-H...O(water) hydrogen bond in synthon VII, a recurring pattern in the diacid structures. In summary, the carboxylic acid group forms an O-H...N hydrogen bond in pyrazine monocarboxylic acids and an O-H...O hydrogen bond in pyrazine dicarboxylic acids. This structural analysis correlates molecular features with supramolecular synthons in pyridine and pyrazine carboxylic acids for future crystal engineering strategies.     

Hydrogen bond strength modulates the mechanical strength of ferric-thiolate bonds in rubredoxin

It has long been recognized that hydrogen bonds formed by protein backbone amides with cysteinyl S(γ) atoms play important roles in modulating the functional and structural properties of the iron-sulfur centers in proteins. Here we use single molecule atomic force microscopy, cyclic voltammetry, and protein engineering techniques to investigate directly how the strength of N-H···S(γ) hydrogen bonds in the secondary coordination sphere affects the mechanical stability of Fe(III)-thiolate bonds of rubredoxin. Our results show that the mechanical stability of Fe(III)-thiolate bonds in rubredoxin correlates with the strength of N-H···S(γ) hydrogen bonds as reflected by the midpoint reduction potential, providing direct evidence that N-H···S(γ) hydrogen bonds play important roles in modulating the mechanical and kinetic properties of the Fe(III)-thiolate bonds of iron-sulfur proteins and corroborating the important roles of the protein environment in tuning the properties of metal-thiolate bonds.     

A versatile template for the formation of [2]pseudorotaxanes. 1,2-Bis(pyridinium)ethane axles and 24-crown-8 ether wheels

Linear 1,2-bis(pyridinium)ethane 'axles' and macrocyclic 24-membered crown ether 'wheels' (, and ) combine to form [2]pseudorotaxanes. These interpenetrated adducts are held together by N+...O ion-dipole interactions, a series of C-H...O hydrogen bonds and pi-stacking between electron-poor pyridinium rings of the axle and electron-rich catechol rings of the wheel. 1H NMR spectroscopy was used to identify the structural details of the interaction and to determine the thermodynamics of the binding process in solution. Analysis of nine of these adducts by single crystal X-ray crystallography allowed a detailed study of the non-covalent interactions in the solid state. A wide variety of structural changes could be made to the system. The versatility and potential of the template for the construction of permanently interlocked structures such as rotaxanes and catenanes is discussed.     

Peptide bond isosteres: ester or (E)-alkene in the backbone of the collagen triple helix

[structure: see text] Collagen is the most abundant protein in animals. Interstrand N-H...O=C hydrogen bonds between backbone amide groups form a ladder in the middle of the collagen triple helix. Isosteric replacement of the hydrogen-bond-donating amide with an ester or (E)-alkene markedly decreases the conformational stability of the triple helix. Thus, this recurring hydrogen bond is critical to the structural integrity of collagen. In this context, an ester isostere confers more stability than does an (E)-alkene.     

Intramolecular Hydrogen Bonds in Low‐Molecular‐Weight Polyethylene Glycol

We used static DFT calculations to analyze, in detail, the intramolecular hydrogen bonds formed in low‐molecular‐weight polyethylene glycol (PEG) with two to five repeat subunits. Both red‐shifted O?H???O and blue‐shifting C?H???O hydrogen bonds, which control the structural flexibility of PEG, were detected. To estimate the strength of these hydrogen bonds, the quantum theory of atoms in molecules was used. Car–Parrinello molecular dynamics simulations were used to mimic the structural rearrangements and hydrogen‐bond breaking/formation in the PEG molecule at 300 K. The time evolution of the H???O bond length and valence angles of the formed hydrogen bonds were fully analyzed. The characteristic hydrogen‐bonding patterns of low‐molecular‐weight PEG were described with an estimation of their lifetime. The theoretical results obtained, in particular the presence of weak C?H???O hydrogen bonds, could serve as an explanation of the PEG structural stability in the experimental investigation.