Main chain hydrogen bond interactions in the binding of proline-rich gluten peptides to the celiac disease-associated HLA-DQ2 molecule

Binding of peptide epitopes to major histocompatibility complex proteins involves multiple hydrogen bond interactions between the peptide main chain and major histocompatibility complex residues. The crystal structure of HLA-DQ2 complexed with the alphaI-gliadin epitope (LQPFPQPELPY) revealed four hydrogen bonds between DQ2 and peptide main chain amides. This is remarkable, given that four of the nine core residues in this peptide are proline residues that cannot engage in amide hydrogen bonding. Preserving main chain hydrogen bond interactions despite the presence of multiple proline residues in gluten peptides is a key element for the HLA-DQ2 association of celiac disease. We have investigated the relative contribution of each main chain hydrogen bond interaction by preparing a series of N-methylated alphaI epitope analogues and measuring their binding affinity and off-rate constants to DQ2. Additionally, we measured the binding of alphaI-gliadin peptide analogues in which norvaline, which contains a backbone amide hydrogen bond donor, was substituted for each proline. Our results demonstrate that hydrogen bonds at P4 and P2 positions are most important for binding, whereas the hydrogen bonds at P9 and P6 make smaller contributions to the overall binding affinity. There is no evidence for a hydrogen bond between DQ2 and the P1 amide nitrogen in peptides without proline at this position. This is a unique feature of DQ2 and is likely a key parameter for preferential binding of proline-rich gluten peptides and development of celiac disease.     

Hydrogen bonding in helical polypeptides from molecular dynamics simulations and amide hydrogen exchange analysis: alamethicin and melittin in methanol.

Molecular dynamics simulations of ion channel peptides alamethicin and melittin, solvated in methanol at 27 degrees C, were run with either regular alpha-helical starting structures (alamethicin, 1 ns; melittin 500 ps either with or without chloride counterions), or with the x-ray crystal coordinates of alamethicin as a starting structure (1 ns). The hydrogen bond patterns and stabilities were characterized by analysis of the dynamics trajectories with specified hydrogen bond angle and distance criteria, and were compared with hydrogen bond patterns and stabilities previously determined from high-resolution NMR structural analysis and amide hydrogen exchange measurements in methanol. The two alamethicin simulations rapidly converged to a persistent hydrogen bond pattern with a high level of 3(10) hydrogen bonding involving the amide NH's of residues 3, 4, 9, 15, and 18. The 3(10) hydrogen bonds stabilizing amide NH's of residues C-terminal to P2 and P14 were previously proposed to explain their high amide exchange stabilities. The absence, or low levels of 3(10) hydrogen bonds at the N-terminus or for A15 NH, respectively, in the melittin simulations, is also consistent with interpretations from amide exchange analysis. Perturbation of helical hydrogen bonding in the residues before P14 (Aib10-P14, alamethicin; T11-P14, melittin) was characterized in both peptides by variable hydrogen bond patterns that included pi and gamma hydrogen bonds. The general agreement in hydrogen bond patterns determined in the simulations and from spectroscopic analysis indicates that with suitable conditions (including solvent composition and counterions where required), local hydrogen-bonded secondary structure in helical peptides may be predicted from dynamics simulations from alpha-helical starting structures. Each peptide, particularly alamethicin, underwent some large amplitude structural fluctuations in which several hydrogen bonds were cooperatively broken. The recovery of the persistent hydrogen bonding patterns after these fluctuations demonstrates the stability of intramolecular hydrogen-bonded secondary structure in methanol (consistent with spectroscopic observations), and is promising for simulations on extended timescales to characterize the nature of the backbone fluctuations that underlie amide exchange from isolated helical polypeptides.     

Correlation of the sweetness of variants of the protein brazzein with patterns of hydrogen bonds detected by NMR spectroscopy

In sequence-function investigations, approaches are needed for rapidly screening protein variants for possible changes in conformation. Recent NMR methods permit direct detection of hydrogen bonds through measurements of scalar couplings that traverse hydrogen bonds (trans-hydrogen bond couplings). We have applied this approach to screen a series of five single site mutants of the sweet protein brazzein with altered sweetness for possible changes in backbone hydrogen bonding with respect to wild-type. Long range, three-dimensional data correlating connectivities among backbone 1HN, 15N, and 13C' atoms were collected from the six brazzein proteins labeled uniformly with carbon-13 and nitrogen-15. In wild-type brazzein, this approach identified 17 backbone hydrogen bonds. In the mutants, altered magnitudes of the couplings identified hydrogen bonds that were strengthened or weakened; missing couplings identified hydrogen bonds that were broken, and new couplings indicated the presence of new hydrogen bonds. Within the series of brazzein mutants investigated, a pattern was observed between sweetness and the integrity of particular hydrogen bonds. All three "sweet" variants exhibited the same pattern of hydrogen bonds, whereas all three "non-sweet" variants lacked one hydrogen bond at the middle of the alpha-helix, where it is kinked, and one hydrogen bond in the middle of beta-strands II and III, where they are twisted. Two of the non-sweet variants lack the hydrogen bond connecting the N and C termini. These variants showed greater mobility in the N- and C-terminal regions than wild-type brazzein.     

BACE1 inhibitors: optimization by replacing the P1' residue with non-acidic moiety

Recently, we reported potent BACE1 inhibitors KMI-429, -684, and -574 possessing a hydroxymethylcarbonyl isostere as a substrate transition-state mimic. These inhibitors showed potent inhibitory activities in enzymatic and cell assays, especially, KMI-429 was confirmed to significantly inhibit Abeta production in vivo. However, acidic moieties at the P(4) and P(1)' positions of KMI-compounds were thought to be unfavorable for membrane permeability across the blood-brain barrier. Herein, we replaced acidic moieties at the P(4) position with other hydrogen bond acceptor groups, and these inhibitors exhibited improved BACE1 inhibitory activities in cultured cells. In this study, we replaced the acidic moieties at the P(1)' position with non-acidic and low molecular sized moieties.     

Water molecules coupled to the redox-active tyrosine Y(D) in photosystem II as detected by FTIR spectroscopy

The redox-active tyrosine YD (D2-Tyr160) in photosystem II (PSII) serves as a side-path electron donor to P680. When YD is oxidized, a proton is released from phenolic OH, and a neutral radical YD* is formed. A hydrogen bond network around YD must be deeply involved in the mechanism of the YD reaction. In this study, we have detected water molecules structurally coupled to YD by means of Fourier transform infrared (FTIR) spectroscopy. Light-induced YD*/YD FTIR difference spectrum of a hydrated film of the PSII core complexes from Thermosynechococcus elongatus showed major signals at 3636(-)/3617(+) and 3594(+)/3585(-) cm-1 in the weakly hydrogen bonded OH stretching region. These peaks downshifted by 11-12 cm-1 upon H218O substitution and almost disappeared upon H/D exchange, and hence, they were definitely assigned to the water OH vibrations. Small intramolecular couplings of 3-6 cm-1 estimated from the OH frequencies of residual HOD species in a deuterated film indicate that these OH signals arise from two different water molecules that have significantly asymmetric hydrogen bond structures. Similar OH signals were observed in PSII-enriched membranes from spinach, suggesting that two water molecules commonly exist near YD irrespective of biological species. These water molecules are coupled to YD most probably through a hydrogen bond network or one of them possibly interacts directly with YD, and thus, they may play crucial roles in the YD reaction by forming a proton-transfer pathway and tuning the redox potential of YD.     

Penicillopepsin-JT2, a recombinant enzyme from Penicillium janthinellum and the contribution of a hydrogen bond in subsite S3 to k(cat)

The nucleotide sequence of the gene (pepA) of a zymogen of an aspartic proteinase from Penicillium janthinellum with a 71% identity in the deduced amino acid sequence to penicillopepsin (which we propose to call penicillopepsin-JT1) has been determined. The gene consists of 60 codons for a putative leader sequence of 20 amino acid residues, a sequence of about 150 nucleotides that probably codes for an activation peptide and a sequence with two introns that codes for the active aspartic proteinase. This gene, inserted into the expression vector pGPT-pyrG1, was expressed in an aspartic proteinase-free strain of Aspergillus niger var. awamori in high yield as a glycosylated form of the active enzyme that we call penicillopepsin-JT2. After removal of the carbohydrate component with endoglycosidase H, its relative molecular mass is between 33,700 and 34,000. Its kinetic properties, especially the rate-enhancing effects of the presence of alanine residues in positions P3 and P2' of substrates, are similar to those of penicillopepsin-JT1, endothiapepsin, rhizopuspepsin, and pig pepsin. Earlier findings suggested that this rate-enhancing effect was due to a hydrogen bond between the -NH- of P3 and the hydrogen bond accepting oxygen of the side chain of the fourth amino acid residue C-terminal to Asp215. Thr219 of penicillopepsin-JT2 was mutated to Ser, Val, Gly, and Ala. Thr219Ser showed an increase in k(cat) when a P3 residue was present in the substrate, which was similar to that of the wild-type, whereas the mutants Thr219Val, Thr219Gly, and Thr219Ala showed no significant increase when a P3 residue was added. The results show that the putative hydrogen bond alone is responsible for the increase. We propose that by locking the -NH- of P3 to the enzyme, the scissile peptide bond between P1 and P1' becomes distorted toward a tetrahedral conformation and becomes more susceptible to nucleophilic attack by the catalytic apparatus without the need of a conformational change in the enzyme.     

The structural basis of a conserved P2 threonine in canonical serine proteinase inhibitors

Bowman-Birk inhibitors (BBIs) are a well-studied family of canonical inhibitor proteins of serine proteinases. In nature, the active region of BBIs possesses a highly conserved Thr at the P2 position. The importance of this residue has been reemphasized by synthetic BBI reactive site loop proteinomimetics. In particular, this residue was exclusively identified for active chymotrypsin inhibitors selected from a BBI template-assisted combinatorial peptide library. A further kinetic analysis of 26 P2 variant peptides revealed that Thr provides both optimal binding affinity and optimal resistance against enzymatic turnover by chymotrypsin. Herein, we report the (1)H-NMR spectroscopic study of a 5-membered sub-set of these reactive site loop peptides representing a stepwise elimination of the Thr side-chain functionalities and inversion of its side-chain chirality. The P2 Thr variant adopts a three-dimensional structure that closely mimics the one of the corresponding region of the complete protein. This validates the use of this template for the investigation of structure-function relationships. While the overall backbone geometry is similar in all studied variants, conformational changes induced by the modification of the P2 side chain have now been identified and provide a rational explanation of the kinetically observed functional differences. Eliminating the gamma-methyl group has little structural effect, whereas the elimination of the gamma-oxygen atom or the inversion of the side-chain chirality results in characteristic changes to the intramolecular hydrogen bond network. We conclude that the transannular hydrogen bond between the P2 Thr side-chain hydroxyl and the P5' backbone amide is an important conformational constraint and directs the hydrophobic contact of the P2 Thr side chain with the enzyme surface in a functionally optimal geometry, both in the proteinomimetic and the native protein. In at least four canonical inhibitor protein families similar structural arrangements for a conserved P2 Thr have been observed, which suggests an analogous functional role. Substitutions at P2 of the proteinomimetic also affect the conformational balance between cis and trans isomers at a distant Pro-Pro motif (P3'-P4'). Presented with a mixture of cis/trans isomers chymotrypsin appears to interact preferably with the conformer that retains the cis-P3' Pro-trans-P4' Pro geometry found in the parent BBI protein.     

Direct analysis of backbone-backbone hydrogen bond formation in protein folding transition states

Here we investigate the role of backbone-backbone hydrogen bonding interactions in stabilizing the protein folding transition states of two model protein systems, the B1 domain of protein L (ProtL) and the P22 Arc repressor. A backbone modified analogue of ProtL containing an amide-to-ester bond substitution between residues 105 and 106 was prepared by total chemical synthesis, and the thermodynamic and kinetic parameters associated with its folding reaction were evaluated. Ultimately, these parameters were used in a Phi-value analysis to determine if the native backbone-backbone hydrogen bonding interaction perturbed in this analogue (i.e. a hydrogen bond in the first beta-turn of ProtL's beta-beta-alpha-beta-beta fold) was formed in the transition state of ProtL's folding reaction. Also determined were the kinetic parameters associated with the folding reactions of two Arc repressor analogues, each containing an amide-to-ester bond substitution in the backbone of their polypeptide chains. These parameters were used together with previously established thermodynamic parameters for the folding of these analogues in Phi-value analyses to determine if the native backbone-backbone hydrogen bonding interactions perturbed in these analogues (i.e. a hydrogen bond at the end of the intersubunit beta-sheet interface and hydrogen bonds at the beginning of the second alpha-helix in Arc repressor's beta-alpha-alpha structure) were formed in the transition state of Arc repressor's folding reaction. Our results reveal that backbone-backbone hydrogen bonding interactions are formed in the beta-turn and alpha-helical transition state structures of ProtL and Arc repressor, respectively; and they were not formed in the intersubunit beta-sheet interface of Arc repressor, a region of Arc repressor's polypeptide chain previously shown to have other non-native-like conformations in Arc's protein folding transition state.     

Learning about protein hydrogen bonding by minimizing contrastive divergence

Defining the strength and geometry of hydrogen bonds in protein structures has been a challenging task since early days of structural biology. In this article, we apply a novel statistical machine learning technique, known as contrastive divergence, to efficiently estimate both the hydrogen bond strength and the geometric characteristics of strong interpeptide backbone hydrogen bonds, from a dataset of structures representing a variety of different protein folds. Despite the simplifying assumptions of the interatomic energy terms used, we determine the strength of these hydrogen bonds to be between 1.1 and 1.5 kcal/mol, in good agreement with earlier experimental estimates. The geometry of these strong backbone hydrogen bonds features an almost linear arrangement of all four atoms involved in hydrogen bond formation. We estimate that about a quarter of all hydrogen bond donors and acceptors participate in these strong interpeptide hydrogen bonds.     

Oxygen activation by cytochrome P450 monooxygenase

Unlike photosystem II (PSII) that catalyzes formation of the O-O bond, the cytochromes P450 (P450), members of a superfamily of hemoproteins, catalyze the scission of the O-O bond of dioxygen molecules and insert a single oxygen atom into unactivated hydrocarbons through a hydrogen abstraction-oxygen rebound mechanism. Hydroxylation of the unactivated hydrocarbons at physiological temperatures is vital for many cellar processes such as the biosynthesis of many endogenous compounds and the detoxification of xenobiotics in humans and plants. Even though it carries out the opposite of the water splitting reaction, P450 may share similarities to PSII in proton delivery networks, oxygen and water access channels, and consecutive electron transfer processes. In this article, we review recent advances in understanding the molecular mechanisms by which P450 activates dioxygen.     

NMR study of the inhibition of pepsin by glyoxal inhibitors: mechanism of tetrahedral intermediate stabilization by the aspartyl proteases

Z-Ala-Ala-Phe-glyoxal (where Z is benzyloxycarbonyl) has been shown to be a competitive inhibitor of pepsin with a Ki = 89 +/- 24 nM at pH 2.0 and 25 degrees C. Both the ketone carbon (R13COCHO) and the aldehyde carbon (RCO13CHO) of the glyoxal group of Z-Ala-Ala-Phe-glyoxal have been 13C-enriched. Using 13C NMR, it has been shown that when the inhibitor is bound to pepsin, the glyoxal keto and aldehyde carbons give signals at 98.8 and 90.9 ppm, respectively. This demonstrates that pepsin binds and preferentially stabilizes the fully hydrated form of the glyoxal inhibitor Z-Ala-Ala-Phe-glyoxal. From 13C NMR pH studies with glyoxal inhibitor, we obtain no evidence for its hemiketal or hemiacetal hydroxyl groups ionizing to give oxyanions. We conclude that if an oxyanion is formed its pKa must be >8.0. Using 1H NMR, we observe four hydrogen bonds in free pepsin and in pepsin/Z-Ala-Ala-Phe-glyoxal complexes. In the pepsin/pepstatin complex an additional hydrogen bond is formed. We examine the effect of pH on hydrogen bond formation, but we do not find any evidence for low-barrier hydrogen bond formation in the inhibitor complexes. We conclude that the primary role of hydrogen bonding to catalytic tetrahedral intermediates in the aspartyl proteases is to correctly orientate the tetrahedral intermediate for catalysis.     

Temperature-dependence of protein hydrogen bond properties as studied by high-resolution NMR

The temperature-dependence of a large number of NMR parameters describing hydrogen bond properties in the protein ubiquitin was followed over a range from 5 to 65 degrees C. The parameters comprise hydrogen bond (H-bond) scalar couplings, h3JNC', chemical shifts, amide proton exchange rates, 15N relaxation parameters as well as covalent 1JNC' and 1JNH couplings. A global weakening of the h3JNC' coupling with increasing temperature is accompanied by a global upfield shift of the amide protons and a decrease of the sequential 1JNC' couplings. If interpreted as a linear increase of the N...O distance, the change in h3JNC' corresponds to an average linear thermal expansion coefficient for the NH-->O hydrogen bonds of 1.7 x 10(-4)/K, which is in good agreement with overall volume expansion coefficients observed for proteins. A residue-specific analysis reveals that not all hydrogen bonds are affected to the same extent by the thermal expansion. The end of beta-sheet beta1/beta5 at hydrogen bond E64-->Q2 appears as the most thermolabile, whereas the adjacent hydrogen bond I3-->L15 connecting beta-strands beta1 and beta2 is even stabilized slightly at higher temperatures. Additional evidence for the stabilization of the beta1/beta2 beta-hairpin at higher temperatures is found in reduced hydrogen exchange rates for strand end residue V17. This reduction corresponds to a stabilizing change in free energy of 9.7 kJ/mol for the beta1/beta2 hairpin. The result can be linked to the finding that the beta1/beta2 hairpin behaves as an autonomously folding unit in the A-state of ubiquitin under changed solvent conditions. For several amide groups the temperature-dependencies of the amide exchange rates and H-bond scalar couplings are uncorrelated. Therefore, amide exchange rates are not a sole function of the hydrogen bond "strength" as given by the electronic overlap of donors and acceptors, but are clearly dependent on other blocking mechanisms.     

Differential stability of the triple helix of (Pro-Pro-Gly)10 in H2O and D2O: thermodynamic and structural explanations

(Pro-Pro-Gly)10 [(PPG10)], a collagen-like polypeptide, forms a triple-helical, polyproline-II structure in aqueous solution at temperatures somewhat lower than physiological, with a melting temperature of 24.5 degrees C. In this article, we present circular dichroism spectra that demonstrate an increase of the melting temperature with the addition of increasing amounts of D2O to an H2O solution of (PPG)10, with the melting temperature reaching 40 degrees C in pure D2O. A thermodynamic analysis of the data demonstrates that this result is due to an increasing enthalpy of unfolding in D2O vs. H2O. To provide a theoretical explanation for this result, we have used a model for hydration of (PPG)10 that we developed previously, in which inter-chain water bridges are formed between sterically crowded waters and peptide bond carbonyls. Energy minimizations were performed upon this model using hydrogen bond parameters for water, and altered hydrogen bond parameters that reproduced the differences in carbonyl oxygen-water oxygen distances found in small-molecule crystal structures containing oxygen-oxygen hydrogen bonds between organic molecules and H2O or D2O. It was found that using hydrogen bond parameters that reproduced the distance typical of hydrogen bonds to D2O resulted in a significant lowering of the potential energy of hydrated (PPG)10. This lowering of the energy involved energetic terms that were only indirectly related to the altered hydrogen bond parameters, and were therefore not artifactual; the intra-(PPG10) energy, plus the water-(PPG10) van der Waals energy (not including hydrogen bond interactions), were lowered enough to qualitatively account for the lower enthalpy of the triple-helical conformation, relative to the unfolded state, in D2O vs. H2O. This result indicates that the geometry of the carbonyl-D2O hydrogen bonds allows formation of good hydrogen bonds without making as much of an energetic sacrifice from other factors as in the case of hydration by H2O.     

Structures of thiolate- and carboxylate-ligated ferric H93G myoglobin: models for cytochrome P450 and for oxyanion-bound heme proteins

Crystal structures of the ferric H93G myoglobin (Mb) cavity mutant containing either an anionic proximal thiolate sulfur donor or a carboxylate oxygen donor ligand are reported at 1.7 and 1.4 A resolution, respectively. The crystal structure and magnetic circular dichroism spectra of the H93G Mb beta-mercaptoethanol (BME) thiolate adduct reveal a high-spin, five-coordinate complex. Furthermore, the bound BME appears to have an intramolecular hydrogen bond involving the alcohol proton and the ligated thiolate sulfur, mimicking one of the three proximal N-H...S hydrogen bonds in cytochrome P450. The Fe is displaced from the porphyrin plane by 0.5 A and forms a 2.41 A Fe-S bond. The Fe(3+)-S-C angle is 111 degrees , indicative of a covalent Fe-S bond with sp(3)-hybridized sulfur. Therefore, the H93G Mb.BME complex provides an excellent protein-derived structural model for high-spin ferric P450. In particular, the Fe-S bond in high-spin ferric P450-CAM has essentially the same geometry despite the constraints imposed by covalent linkage of the cysteine to the protein backbone. This suggests that evolution led to the geometric optimization of the proximal Fe-S(cysteinate) bond in P450. The crystal structure and spectral properties of the H93G Mb acetate adduct reveal a high-spin, six-coordinate complex with proximal acetate and distal water axial ligands. The distal His-64 forms a hydrogen bond with the bound water. The Fe-acetate bonding geometry is inconsistent with an electron pair along the Fe-O bond as the Fe-O-C angle is 152 degrees and the Fe is far from the plane of the acetate. Thus, the Fe-O bonding is ionic. The H93G Mb cavity mutant has already been shown to be a versatile model system for the study of ligand binding to heme proteins; this investigation affords the first structural evidence that nonimidazole exogenous ligands bind in the proximal ligation site.     

Successful prediction of the hydrogen bond network of the 3-oxo-12alpha-hydroxy-5beta-cholan-24-oic acid crystal from resolution of the crystal structure of deoxycholic acid and its three 3,12-dihydroxy epimers

Crystal structures of p-xylene-crystallized deoxycholic acid (3alpha,12alpha-dihydroxy-5beta-cholan-24-oic acid) and its three epimers (3beta,12alpha-; 3alpha,12beta-; and 3beta,12beta-) have been solved. Deoxycholic acid forms a crystalline (P21) complex with the solvent with a 2:1 stoichiometry whereas crystals of the three epimers do not form inclusion compounds. Crystals of the 3beta,12beta-epimer are hexagonal, whereas the 3alpha,12beta-and 3beta,12alpha-epimers crystallize in the P2(1)2(1)2(1) orthorhombic space group. The three hydrogen bond sites (two hydroxy groups, i. e. O3-H, and O12-H, and the carboxylic acid group of the side chain, O24bO24a-H) simultaneously act as hydrogen bond donors and acceptors. The hydrogen bond network in the crystals was analyzed and the following sequences have been observed: two chains (abcabc... or acbacb... ) and two rings (abc or acb), which constitute a complete set of all the possible sequences which can be drawn for an intermolecular hydrogen bond network formed by three hydrogen bond donor/acceptor sites forming crossing hydrogen bonds. The orientation of O3-H (alpha or beta) determines the sequence of the acceptor and the donor groups involved in the pattern: O24a --> O12 --> O3 --> O24b when it is alpha and O24a --> O3 --> O12--> O24B when it is beta. These observations were used to predict the hydrogen bond network of p-xylene-crystallized 3-oxo,12alpha-hydroxy-5beta-cholan-24-oic acid. This compound has two hydrogen bond donor and three potential hydrogen bond acceptor sites. According to the previous sequence set, this compound should crystallize in the monoclinic P21 system, should form a complex with the solvent, O24b should not participate in the hydrogen bond network, and the chain sequence O24a --> O12 --> O3 would be followed. All predictions were confirmed experimentally.     

Influence of the environment in the conformation of alpha-helices studied by protein database search and molecular dynamics simulations

The influence of the solvent on the main-chain conformation (phi and Psi dihedral angles) of alpha-helices has been studied by complementary approaches. A first approach consisted in surveying crystal structures of both soluble and membrane proteins. The residues of analysis were further classified as exposed to either the water (polar solvent) or the lipid (apolar solvent) environment or buried to the core of the protein (intermediate polarity). The statistical results show that the more polar the environment, the lower the value of phi(i) and the higher the value of Psi(i) are. The intrahelical hydrogen bond distance increases in water-exposed residues due to the additional hydrogen bond between the peptide carbonyl oxygen and the aqueous environment. A second approach involved nanosecond molecular dynamics simulations of poly-Ala alpha-helices in environments of different polarity: water to mimic hydrophilic environments that can form hydrogen bonds with the peptide carbonyl oxygen and methane to mimic hydrophobic environments without this hydrogen bond capabilities. These simulations reproduce similar effects in phi and Psi angles and intrahelical hydrogen bond distance and angle as observed in the protein survey analysis. The magnitude of the intrahelical hydrogen bond in the methane environment is stronger than in the water environment, suggesting that alpha-helices in membrane-embedded proteins are less flexible than in soluble proteins. There is a remarkable coincidence between the phi and Psi angles obtained in the analysis of residues exposed to the lipid in membrane proteins and the results from computer simulations in methane, which suggests that this simulation protocol properly mimic the lipidic cell membrane and reproduce several structural characteristics of membrane-embedded proteins. Finally, we have compared the phi and Psi torsional angles of Pro kinks in membrane protein crystal structures and in computer simulations.     

Structure-based design and synthesis of novel P2/P3 modified,non-peptidic β-secretase (BACE-1) inhibitors

Starting from peptidomimetic BACE-1 inhibitors, the P2 amino acid including the P2/P3 peptide bond was replaced by a rigid 3-aminomethyl cyclohexane carboxylic acid. Co-crystallization revealed an unexpected binding mode with the P3/P4 amide bond placed into the S3 pocket resulting in a new hydrogen bond interaction pattern. Further optimization based on this structure resulted in highly potent BACE-1 inhibitors with selectivity over BACE-2 and cathepsin D.     

A double-deletion method to quantifying incremental binding energies in proteins from experiment: example of a destabilizing hydrogen bonding pair

The contribution of a specific hydrogen bond in apoflavodoxin to protein stability is investigated by combining theory, experiment and simulation. Although hydrogen bonds are major determinants of protein structure and function, their contribution to protein stability is still unclear and widely debated. The best method so far devised to estimate the contribution of side-chain interactions to protein stability is double mutant cycle analysis, but the interaction energies so derived are not identical to incremental binding energies (the energies quantifying net contributions of two interacting groups to protein stability). Here we introduce double-deletion analysis of 'isolated' residue pairs as a means to precisely quantify incremental binding. The method is exemplified by studying a surface-exposed hydrogen bond in a model protein (Asp96/Asn128 in apoflavodoxin). Combined substitution of these residues by alanines slightly destabilizes the protein due to a decrease in hydrophobic surface burial. Subtraction of this effect, however, clearly indicates that the hydrogen-bonded groups in fact destabilize the native conformation. In addition, molecular dynamics simulations and classic double mutant cycle analysis explain quantitatively that, due to frustration, the hydrogen bond must form in the native structure because when the two groups get approximated upon folding their binding becomes favorable. We would like to remark that 1), this is the first time the contribution of a specific hydrogen bond to protein stability has been measured by experiment; and 2), more hydrogen bonds need to be analyzed to draw general conclusions on protein hydrogen bond energetics. To that end, the double-deletion method should be of help.