Buried Hydrogen Bonds Research Articles

Effect of the Lys62Ala Mutation on the Thermal Stability of BstHPr Protein by Molecular Dynamics.

We analyzed the thermal stability of the BstHPr protein through the site-directed point mutation Lys62 replaced by Ala residue using molecular dynamics simulations at five different temperatures: 298, 333, 362, 400, and 450 K, for periods of 1 μs and in triplicate. The results from the mutant thermophilic BstHPrm protein were compared with those of the wild-type thermophilic BstHPr protein and the mesophilic BsHPr protein. Structural and molecular interaction analyses show that proteins lose stability as temperature increases. Mutant and wild-type proteins behave similarly up to 362 K. However, at 400 K the mutant protein shows greater structural instability, losing more buried hydrogen bonds and exposing more of its non-polar residues to the solvent. Therefore, in this study, we confirmed that the salt bridge network of the Glu3-Lys62-Glu36 triad, made up of the Glu3-Lys62 and Glu36-Lys62 ion pairs, provides thermal stability to the thermophilic BstHPr protein.

De novo design of obligate ABC-type heterotrimeric proteins

The de novo design of three protein chains that associate to form a heterotrimer (but not any of the possible two-chain heterodimers) and that can drive the assembly of higher-order branching structures is an important challenge for protein design. We designed helical heterotrimers with specificity conferred by buried hydrogen bond networks and large aromatic residues to enhance shape complementary packing. We obtained ten designs for which all three chains cooperatively assembled into heterotrimers with few or no other species present. Crystal structures of a helical bundle heterotrimer and extended versions, with helical repeat proteins fused to individual subunits, showed all three chains assembling in the designed orientation. We used these heterotrimers as building blocks to construct larger cyclic oligomers, which were structurally validated by electron microscopy. Our three-way junction designs provide new routes to complex protein nanostructures and enable the scaffolding of three distinct ligands for modulation of cell signaling.

Potent Th2 Cytokine Bias of Natural Killer T Cell by CD1d Glycolipid Ligands: Anchoring Effect of Polar Groups in the Lipid Component.

Th2-biasing CD1d ligands are attractive potential candidates for adjuvants and therapeutic drugs. However, the number of potent ligands is limited, and their biasing mechanism remain unclear. Herein, a series of novel Th2-biasing CD1d glycolipid ligands, based on modification of their lipid part, have been identified. These have shown high binding affinities and efficient Th2 cytokine production. Importantly, the truncated acyl chain containing variants still retain their binding affinities and agonistic activities, which can be associated with an "anchoring effect," that is, formation of a buried hydrogen bond between a polar group on the acyl chain and the CD1d lipid-binding pocket. The analysis indicated that the appearance rates of ligand-CD1d complexes on the cell surface were involved in Th2-biasing responses. The designed ligands, having the anchor in the shorter lipid part, are one of the most potent Th2-biasing ligands among the known ligands.

Potent Th2 Cytokine Bias of Natural Killer T Cell by CD1d Glycolipid Ligands: Anchoring Effect of Polar Groups in the Lipid Component

AbstractTh2‐biasing CD1d ligands are attractive potential candidates for adjuvants and therapeutic drugs. However, the number of potent ligands is limited, and their biasing mechanism remain unclear. Herein, a series of novel Th2‐biasing CD1d glycolipid ligands, based on modification of their lipid part, have been identified. These have shown high binding affinities and efficient Th2 cytokine production. Importantly, the truncated acyl chain containing variants still retain their binding affinities and agonistic activities, which can be associated with an “anchoring effect,” that is, formation of a buried hydrogen bond between a polar group on the acyl chain and the CD1d lipid‐binding pocket. The analysis indicated that the appearance rates of ligand–CD1d complexes on the cell surface were involved in Th2‐biasing responses. The designed ligands, having the anchor in the shorter lipid part, are one of the most potent Th2‐biasing ligands among the known ligands.

Buried Hydrogen Bond Interactions Contribute to the High Potency of Complement Factor D Inhibitors.

Aberrant activation of the complement system is associated with diseases, including paroxysmal nocturnal hemoglobinuria and age-related macular degeneration. Complement factor D is the rate-limiting enzyme for activating the alternative pathway in the complement system. Recent development led to a class of potent amide containing pyrrolidine derived factor D inhibitors. Here, we used biochemical enzymatic and biolayer interferometry assays to demonstrate that the amide group improves the inhibitor potency by more than 80-fold. Our crystal structures revealed buried hydrogen bond interactions are important. Molecular orbital analysis from quantum chemistry calculations dissects the chemical groups participating in these interactions. Free energy calculation supports the differential contributions of the amide group to the binding affinities of these inhibitors. Cell-based hemolysis assay confirmed these compounds inhibit factor D mediated complement activation via the alternative pathway. Our study highlights the important interactions contributing to the high potency of factor D inhibitors reported recently.

Revealing buried hydrogen bonds

Surface Science One way to stabilize self-assembly monolayers of alkane thiols on gold surfaces is to introduce amide groups into the chain near the thiol groups to form hydrogen-bonding networks. This subsurface location has made direct structural studies challenging, but Thomas et al. have used scanning tunneling microscopy and local tunneling barrier height studies of 3-mercapto- N -nonylpropionamide on Au{111} (gold films on mica) at cryogenic temperatures to map the spatial orientation of the hydrogen bond. These networks can cross the boundaries between the domains that form with different rotational orientations. ACS Nano 10.1021/acsnano.6b01717 (2016).

A comparison of successful and failed protein interface designs highlights the challenges of designing buried hydrogen bonds

A comparison of successful and failed protein interface designs highlights the challenges of designing buried hydrogen bonds

A comparison of successful and failed protein interface designs highlights the challenges of designing buried hydrogen bonds

The accurate design of new protein-protein interactions is a longstanding goal of computational protein design. However, most computationally designed interfaces fail to form experimentally. This investigation compares five previously described successful de novo interface designs with 158 failures. Both sets of proteins were designed with the molecular modeling program Rosetta. Designs were considered a success if a high-resolution crystal structure of the complex closely matched the design model and the equilibrium dissociation constant for binding was less than 10 μM. The successes and failures represent a wide variety of interface types and design goals including heterodimers, homodimers, peptide-protein interactions, one-sided designs (i.e., where only one of the proteins was mutated) and two-sided designs. The most striking feature of the successful designs is that they have fewer polar atoms at their interfaces than many of the failed designs. Designs that attempted to create extensive sets of interface-spanning hydrogen bonds resulted in no detectable binding. In contrast, polar atoms make up more than 40% of the interface area of many natural dimers, and native interfaces often contain extensive hydrogen bonding networks. These results suggest that Rosetta may not be accurately balancing hydrogen bonding and electrostatic energies against desolvation penalties and that design processes may not include sufficient sampling to identify side chains in preordered conformations that can fully satisfy the hydrogen bonding potential of the interface.

Molecular Determinants of Self-Association and Rearrangement of a Trimeric Intermediate during the Assembly of a Parvovirus Capsid

The minute virus of mice (MVM) provides a simple model for the dissection of the molecular determinants of the self-assembly, stability, and dynamics of a biological supramolecular complex. MVM assembly involves the trimerization of capsid subunits in the cytoplasm; trimers are transported to the nucleus, where they suffer a conformational change and are made competent for capsid formation. Our previous study revealed that capsid assembly from trimers is dependent on stronger intertrimer interactions that are equally spaced in an equatorial belt surrounding each trimer. We have now targeted the interfaces between monomers within each trimer to identify the molecular determinants of trimerization and the rearrangement needed for capsid assembly. Twenty-eight amino acid residues per monomer were individually mutated to alanine to remove most of the stronger intersubunit interactions. The effects on trimer and capsid assembly and virus infectivity in cells were analyzed. No side chain was individually required for trimer assembly in the cytoplasm; in contrast, half of them were required to make the trimers competent for nuclear capsid assembly, even though none was close to intertrimer interfaces. These critical side chains are conserved and participate in extensive hydrophobic contacts, buried hydrogen bonds, or salt bridges between subunits. This study on MVM capsid assembly reveals that: (i) trimerization is a robust process, insensitive to removal of individual intersubunit interactions; and (ii) the rearrangement of the trimer intermediate required for capsid assembly is a global process that depends on the establishment of many interactions along the protein–protein interfaces within each trimer.

Computational design of second‐site suppressor mutations at protein–protein interfaces

The importance of a protein-protein interaction to a signaling pathway can be established by showing that amino acid mutations that weaken the interaction disrupt signaling, and that additional mutations that rescue the interaction recover signaling. Identifying rescue mutations, often referred to as second-site suppressor mutations, controls against scenarios in which the initial deleterious mutation inactivates the protein or disrupts alternative protein-protein interactions. Here, we test a structure-based protocol for identifying second-site suppressor mutations that is based on a strategy previously described by Kortemme and Baker. The molecular modeling software Rosetta is used to scan an interface for point mutations that are predicted to weaken binding but can be rescued by mutations on the partner protein. The protocol typically identifies three types of specificity switches: knob-in-to-hole redesigns, switching hydrophobic interactions to hydrogen bond interactions, and replacing polar interactions with nonpolar interactions. Computational predictions were tested with two separate protein complexes; the G-protein Galpha(i1) bound to the RGS14 GoLoco motif, and UbcH7 bound to the ubiquitin ligase E6AP. Eight designs were experimentally tested. Swapping a buried hydrophobic residue with a polar residue dramatically weakened binding affinities. In none of these cases were we able to identify compensating mutations that returned binding to wild-type affinity, highlighting the challenges inherent in designing buried hydrogen bond networks. The strongest specificity switches were a knob-in-to-hole design (20-fold) and the replacement of a charge-charge interaction with nonpolar interactions (55-fold). In two cases, specificity was further tuned by including mutations distant from the initial design. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Development and validation of an empirical free energy function for calculating protein–protein binding free energy surfaces

In a previous paper, we described a novel empirical free energy function that was used to accurately predict experimental binding free energies for a diverse test set of 31 protein–protein complexes to within ≈ 1.0 kcal. Here, we extend that work and show that an updated version of the function can be used to (1) accurately predict native binding free energies and (2) rank crystallographic, native-like and non-native binding modes in a physically realistic manner. The modified function includes terms designed to capture some of the unfavorable interactions that characterize non-native interfaces. The function was used to calculate one-dimensional binding free energy surfaces for 21 protein complexes. In roughly 90% of the cases tested, the function was used to place native-like and crystallographic binding modes in global free energy minima. Our analysis further suggests that buried hydrogen bonds might provide the key to distinguishing native from non-native interactions. To the best of our knowledge our function is the only one of its kind, a single expression that can be used to accurately calculate native and non-native binding free energies for a large number of proteins. Given the encouraging results presented in this paper, future work will focus on improving the function and applying it to the protein–protein docking problem.

Modest stabilization by most hydrogen-bonded side-chain interactions in membrane proteins

Understanding the energetics of molecular interactions is fundamental to all of the central quests of structural biology including structure prediction and design, mapping evolutionary pathways, learning how mutations cause disease, drug design, and relating structure to function. Hydrogen-bonding is widely regarded as an important force in a membrane environment because of the low dielectric constant of membranes and a lack of competition from water. Indeed, polar residue substitutions are the most common disease-causing mutations in membrane proteins. Because of limited structural information and technical challenges, however, there have been few quantitative tests of hydrogen-bond strength in the context of large membrane proteins. Here we show, by using a double-mutant cycle analysis, that the average contribution of eight interhelical side-chain hydrogen-bonding interactions throughout bacteriorhodopsin is only 0.6 kcal mol(-1). In agreement with these experiments, we find that 4% of polar atoms in the non-polar core regions of membrane proteins have no hydrogen-bond partner and the lengths of buried hydrogen bonds in soluble proteins and membrane protein transmembrane regions are statistically identical. Our results indicate that most hydrogen-bond interactions in membrane proteins are only modestly stabilizing. Weak hydrogen-bonding should be reflected in considerations of membrane protein folding, dynamics, design, evolution and function.

Continuous Molecular Evolution of Protein-Domain Structures by Single Amino Acid Changes

Protein structures cluster into families of folds that can result from extremely different amino acid sequences [1]. Because the enormous amount of genetic information generates a limited number of protein folds [2], a particular domain structure often assumes numerous functions. How new protein structures and new functions evolve under these limitations remains elusive. Molecular evolution may be driven by the ability of biomacromolecules to adopt multiple conformations as a bridge between different folds [3-6]. This could allow proteins to explore new structures and new tasks while part of the structural ensemble retains the initial conformation and function as a safeguard [7]. Here we show that a global structural switch can arise from single amino acid changes in cysteine-rich domains (CRD) of cnidarian nematocyst proteins. The ability of these CRDs to form two structures with different disulfide patterns from an identical cysteine pattern is distinctive [8]. By applying a structure-based mutagenesis approach, we demonstrate that a cysteine-rich domain can interconvert between two natively occurring domain structures via a bridge state containing both structures. Comparing cnidarian CRD sequences leads us to believe that the mutations we introduced to stabilize each structure reflect the birth of new protein folds in evolution.

Mutation of Asparagine 52 to Glycine Promotes the Alkaline Form of Iso-1-cytochrome c and Causes Loss of Cooperativity in Acid Unfolding

The kinetics and thermodynamics of the alkaline and acid conformational transitions of a Lys 79 --> Ala/Asn 52 --> Gly (A79G52) variant of iso-1-cytochrome c are studied. The Lys 79 --> Ala mutation is designed to limit heme ligation in the alkaline conformer to Lys 73. The Asn 52 --> Gly mutation is intended to shift the population of the alkaline conformer to physiological pH based on the hierarchical nature of the cooperative substructures of this protein. The midpoint pH for formation of the alkaline conformer is approximately 7.45. The kinetics for the alkaline conformational transition of the A79G52 variant are consistent with the ionization constant, pK(H), for the trigger group controlling formation of the alkaline conformer being approximately 9.5. This pK(H) is low for alkaline conformers involving lysine-heme ligation but is consistent with the pK(a) of the highest of three ionizable groups which modulate formation of the histidine-heme alkaline conformer of a His 73 variant of iso-1-cytochrome c [Martinez, R. E., and Bowler, B. E. (2004) J. Am. Chem. Soc. 126, 6751-6758]. The acid transition of the A79G52 variant is split into two phases. Both the Lys 79 --> Ala and Asn 52 --> Gly mutations are expected to affect the buried hydrogen bond network of cytochrome c, suggesting that this network is an important modulator of the acid unfolding of cytochrome c.

Structural Basis for Preferential Recognition of Diaminopimelic Acid-type Peptidoglycan by a Subset of Peptidoglycan Recognition Proteins

Drosophila peptidoglycan recognition protein (PGRP)-LCx and -LCa are receptors that preferentially recognize meso-diaminopimelic acid (DAP)-type peptidoglycan (PGN) present in Gram-negative bacteria over lysine-type PGN of gram-positive bacteria and initiate the IMD signaling pathway, whereas PGRP-LE plays a synergistic role in this process of innate immune defense. How these receptors can distinguish the two types of PGN remains unclear. Here the structure of the PGRP domain of Drosophila PGRP-LE in complex with tracheal cytotoxin (TCT), the monomeric DAP-type PGN, reveals a buried ionic interaction between the unique carboxyl group of DAP and a previously unrecognized arginine residue. This arginine is conserved in the known DAP-type PGN-interacting PGRPs and contributes significantly to the affinity of the protein for the ligand. Unexpectedly, TCT induces infinite head-to-tail dimerization of PGRP-LE, in which the disaccharide moiety, but not the peptide stem, of TCT is positioned at the dimer interface. A sequence comparison suggests that TCT induces heterodimerization of the ectodomains of PGRP-LCx and -LCa in a closely analogous manner to prime the IMD signaling pathway, except that the heterodimer formation is nonperpetuating.

Thermodynamics of Hydrogen Bonding in Hydrophilic and Hydrophobic Media

The thermodynamics of hydrogen bond breaking and formation was studied in solutions of alcohol (methanol, ethanol, 1-propanol) molecules. An extensive series of over 400 molecular dynamics simulations with an aggregate length of over 900 ns was analyzed using an analysis technique in which hydrogen bond (HB) breaking is interpreted as an Eyring process, for which the Gibbs energy of activation DeltaG can be determined from the HB lifetime. By performing simulations at different temperatures, we were able to determine the enthalpy of activation DeltaH and the entropy of activation TDeltaS for this process from the Van't Hoff relation. The equilibrium thermodynamics was determined separately, based on the number of donor hydrogens that are involved in hydrogen bonds. Results (DeltaH) are compared to experimental data from Raman spectroscopy and found to be in good agreement for pure water and methanol. The DeltaG as well as the DeltaG are smooth functions of the composition of the mixtures. The main result of the calculations is that DeltaG is essentially independent of the environment (around 5 kJ/mol), suggesting that buried hydrogen bonds (e.g., in proteins) do not contribute significantly to protein stability. Enthalpically HB formation is a downhill process in all substances; however, for the alcohols there is an entropic barrier of 6-7 kJ/mol, at 298.15 K, which cannot be detected in pure water.

Thermodynamics and Kinetics of Formation of the Alkaline State of a Lys 79→Ala/Lys 73→His Variant of Iso-1-cytochrome c

The alkaline transition kinetics of a Lys 73-->His (H73) variant of iso-1-cytochrome c are triggered by three ionizable groups [Martinez, R. E., and Bowler, B. E. (2004) J. Am. Chem. Soc. 126, 6751-6758]. To eliminate ambiguities caused by overlapping phases due to formation of the Lys 79 alkaline conformer and proline isomerization associated with the His 73 alkaline conformer, we mutated Lys 79 to Ala in the H73 variant (A79H73). The stability and guanidineHCl m-values of the A79H73 and H73 variants at pH 7.5 are the same. The Ala 79 mutation causes formation of the alkaline conformer to depend on [NaCl]. The salt dependence saturates at 500 mM NaCl, and the thermodynamics of alkaline state formation for the A79H73 and H73 variants become identical. The salt dependence is consistent with loss of an electrostatic contact between Lys 79 and heme propionate D in the A79H73 variant. The kinetics of alkaline state formation for the A79H73 variant support the three trigger group model developed for the H73 variant, with the primary trigger, pK(HL), being ionization of His 73. The low pH ionization, pK(H1), is perturbed by the Ala 79 mutation indicating that this ionization is modulated by the buried hydrogen bond network involving heme propionate D. The A79H73 variant has a high spin heme above pH 9 suggesting that the high pH ionization, pK(H2), involves a high spin heme conformer. The proline isomerization phase is modulated by both pK(HL) and pK(H2) indicating that it is sensitive to protein conformation.

Role of Hydrogen Bond Networks and Dynamics in Positive and Negative Cooperative Stabilization of a Protein

Cooperativity mediated through hydrogen bond networks in yeast iso-1-cytochrome c was studied using a thermodynamic triple mutant cycle. Three known stabilizing mutations, Asn 26 to His, Asn 52 to Ile, and Tyr 67 to Phe, were used to construct the triple mutant cycle. The side chain of His 26, a wild-type residue, forms two hydrogen bonds that bridge two substructures of the wild-type protein, and Tyr 67 and Asn 52 are part of an extensive buried hydrogen bond network. The stabilities of all variants in the triple mutant cycle were determined by guanidine hydrochloride denaturation methods and used to determine the pairwise, Delta(2)G(int), and triple interaction energies. His 26 and Ile 52 interact cooperatively (Delta(2)G(int) is 1-2 kcal/mol), whereas the two other pairs of mutations interact anticooperatively (Delta(2)G(int) is -0.5 to -1.5 kcal/mol). Previously reported structural data for iso-1-cytochrome c variants containing these mutations show that changes in the strength of the His 26 to Glu 44 hydrogen bond, apparently caused by changes in main chain dynamics, provide a mechanism for the long distance (His 26 to Phe 67 and His 26 to Ile 52) propagation of pairwise interaction energies. Opposing changes in the strength of the His 26 to Glu 44 hydrogen bond caused by the N52I and Y67F mutations generate a negative triple interaction energy (-0.9 +/-0.7 kcal/mol) that combined with cancellation of cooperative and anticooperative pairwise interactions produce apparent additivity for the stabilizing effects of the single mutations in the triple mutant variant.

Modulation of the midpoint potential of the [2Fe-2S] Rieske iron sulfur center by Qo occupants in the bc1 complex.

Following addition of myxothiazol to antimycin-treated chromatophores from Rhodobacter sphaeroides poised at an ambient redox potential (E(h)) of approximately 300 mV, the amplitude of the flash-induced cytochrome c(1) oxidation in the ms range increased, indicating a decrease in the availability of electrons from the immediate donor to c(1), the Rieske iron-sulfur protein (ISP). Because the effect was seen only over the limited E(h) range, we conclude that it is due to a decrease in the apparent midpoint redox potential (E(m)) of the ISP by about 40 mV on addition of myxothiazol. This is in line with the change in E(m) previously seen in direct redox titrations. Our results show that the reduced ISP binds with quinone at the Q(o) site with a higher affinity than does the oxidized ISP. The displacement of ubiquinone by myxothiazol leads to elimination of this preferential binding of the ISP reduced form and results in a shift in the midpoint potential of ISP to a more negative value. A simple hypothesis to explain this effect is that myxothiazol prevents formation of hydrogen bond of ubiquinone with the reduced ISP. We conclude that all Q(o) site occupants (ubiquinone, UHDBT, stigmatellin) that form hydrogen bonds with the reduced ISP shift the apparent E(m) of the ISP in the same direction to more positive values. Inhibitors that bind in the domain of the Q(o) site proximal to heme b(L) (myxothiazol, MOA-stilbene) and displace ubiquinone from the site cause a decrease in E(m) of ISP. We present a new formalism for treatment of the relation between E(m) change and the binding constants involved, which simplifies analysis. Using this formalism, we estimated that binding free energies for hydrogen bond formation with the Q(o) site occupant, range from the largest value of approximately 23 kJ mol(-1) in the presence of stigmatellin (appropriate for the buried hydrogen bond shown by structures), to a value of approximately 3.5 kJ mol(-1) in the native complex. We discuss this range of values in the context of a model in which the native structure constrains the interaction of ISP with the Q(o) site occupant so as to favor dissociation and the faster kinetics of unbinding necessary for rapid turnover.

The Role of Buried Hydrogen Bonds in Self-Assembled Mixed Composition Thiols on Au{111}

We have investigated the role of internal functionality in self-assembled monolayers of a family of amide-containing alkanethiol molecules on Au{111} using scanning tunneling microscopy. In addition to van der Waals interactions that are present within n-alkanethiol self-assembled monolayers, hydrogen bonding between adjacent buried amide groups contributes to the stability of the amide-containing molecules on the surface and causes spontaneous phase separation upon coadsorption with an n-alkanethiol. A deposition solution concentration dependence study reveals that this is an observed trend across a range of examined solution compositions. Additionally, hydrogen bonding affects the packing structure of the amide-containing alkanethiol self-assembled monolayers. Although they adopt the same (√3×√3)R30° base lattice as n-alkanethiolate self-assembled monolayers, the amide-containing molecules form superlattice structures that are more linear than n-alkanethiol monolayers due to the hydrogen bonds they form. The internal functionality of monolayers can be used to control their formation and stability.