The pseudomolecule method and the structure of globular proteins. II. The example of ribonuclease F1 and T1

The pseudomolecule approach to the structure of globular proteins in which a small number of water molecules are incorporated into the "molecule" is tested again by comparing the ribbon of hydrogen bonds in two proteins, ribonuclease F1 and T1. These two molecules are 59% homologous and have the same backbone conformation both globally and locally. The two ribbons of hydrogen bonds that cover the whole of the backbone are conserved with an accuracy of some 95% providing that allowance is made for the intrusion into one of the pair of such extra factors as the presence of adducts or metal ions, the insertions and the absence of a few water molecules from one of the x-ray data sets. Without these corrections, the conservation of the ribbon is some 85%. There are 35 conserved hydrogen-bonding residues, nearly all of which show many unions to the backbone or interactions with the active site. There are 36 point mutations that involve one or two hydrogen-bonding side chains and nearly all of these have either none or one hydrogen bond to the backbone. These are minor contributors to the ribbon of hydrogen bonds. Of the 71 residues involved in these two categories, all but six fit into the pseudomolecular picture of the structure of globular proteins. The remaining 30 residues almost all contain conserved hydrocarbon side chains that may have a second order effect on the structure through their space filling effects.     

The ribbon of hydrogen bonds and the pseudomolecule in the three-dimensional structure of globular proteins. III. Bovine pancreatic ribonuclease A and bovine seminal ribonuclease

The model of the three-dimensional structure of globular proteins, which is based on a ribbon of hydrogen bonds along the whole of the backbone, is now applied to the comparison between monomeric bovine pancreatic ribonuclease A and dimeric bovine seminal ribonuclease. Some waters are involved in the hydrogen bonding of the ribbon, and the protein molecule plus these waters forms a pseudomolecule. The conformations of the three backbones are essentially identical and the three ribbons of hydrogen bonds are conserved with greater than 90% accuracy. We suggest that the conservation of the backbone conformations of the two molecules is a consequence of the conservation of the ribbons of hydrogen bonds. There are 16 simple mutations between the two molecules, of which 15 involve only side-chain groups with no more than one hydrogen bond to the backbone. Such mutations are not sufficient to change the ribbon of hydrogen bonds and hence there is no change in the backbone conformation. Generalizing this result, we suggest that the conservation of the ribbon is the reason why single point mutations rarely change the conformation of the backbone of the globular proteins.     

The Ribbon of Hydrogen Bonds in the Three-Dimensional Structure of Globular Proteins

We report quantum mechanical computations and experimental evidence which suggest that the backbone conformation of globular proteins depends generally on the conservation of that part of the hydrogen bond network or ribbon which is joined, in general, directly to the backbone and is largely independent of the remainder of this whole network of hydrogen bonds. The familiar hydrogen bonds of the helix and the sheet form about one-half of this ribbon of hydrogen bonds. Both water molecules and hydrogen bonding side chain groups are involved in the formation of the ribbon.This view of the three-dimensional structure of globular proteins in terms of the `molecule' allows us to deal with the non-secondary structure as well as with the familiar secondary structure. It also suggests that the ribbon contains approximately the same number of hydrogen bonds within all three structures – the helix, the sheet and the coil – and that this is the reason for the ease of interconversion of these three structures.The quantum mechanical computations on hydrogen bonding suggest that delocalised water molecules which have substantial mobility are an essential part of the ribbon. This situation arises because the hydrogen bonding groups of the protein molecule are not free to move to optimise the hydrogen bonding geometries as are the oxygen atoms in the waters and ices. Such delocalised water molecules either have high B values or are invisible in the X-ray data and yet are able to form a structure which is as strong as a normal hydrogen bond.The experimental data on the point mutations of the THRI57 residue of the T4 phage lysome provides an initial test of this model. Both the local backbone conformation and the ribbon of hydrogen bonds are conserved throughout all the mutations of residue 157,providing that the delocalised water molecules are accepted as a genuine part of the structure. These mutations include the introduction of hydrocarbon side chains at position 157 when water molecules or other side chain groups take over the formation of the hydrogen bonds.We suggest that, provided steric effects are not important, many point mutations succeed because they leave the ribbon of hydrogen bonds (and so the backbone conformation) largely unchanged.     

The role of local hydrogen bonds in formation of irregular regions of globular protein polypeptide chains]

An assumption is made on the substantial role of local hydrogen bonds in formation of irregular regions of globular protein polypeptide chains. The statistics of the amino acid composition of irregular regions is examined from this point of view. A statistical analysis of side group-backbone hydrogen bonds is carried out for three proteins: alpha-chy-motrypsin, lysozyme and myoglobin. It is shown that short side groups participate in formation of local hydrogen bonds more often than long ones. Conformations of amino acid residues in the first and the last positions are studied in beta-bends of 9 proteins. It is shown that over 70% of these residues are in conformations corresponding to the formation of local hydrogen bonds of three types: backbone-backbone, side groupbackbone, backbone-water molecule-backbone. Thus, the participation of the cooperative hydrogen-bonding network in stabilization of beta-bends is demonstrated.     

Protein folding: could hydrophobic collapse be coupled with hydrogen-bond formation?

A judicious examination of an exhaustive PDB sample of soluble globular proteins of moderate size (N<102) reveals a commensurable relationship between hydrophobic surface burial and number of backbone hydrogen bonds. An analysis of 50,000 conformations along the longest all-atom MD trajectory allows us to infer that not only the hydrophobic collapse is concurrent with the formation of backbone amide-carbonyl hydrogen bonds, they are also dynamically coupled processes. In statistical terms, hydrophobic clustering of the side chains is inevitably conducive to backbone burial and the latter process becomes thermodynamically too costly and kinetically unfeasible without amide-carbonyl hydrogen-bond formation. Furthermore, the desolvation of most hydrogen bonds is exhaustive along the pathway, implying that such bonds guide the collapse process.     

Importance of hydrogen-bonding interactions involving the side chain of Asp158 in the catalytic mechanism of papain

In a previous study, it was shown that replacing Asp158 in papain by Asn had little effect on activity and that the negatively charged carboxylate of Asp158 does not significantly stabilize the active site thiolate-imidazolium ion pair of papain (Ménard et al., 1990). In this paper, we report the kinetic characterization of three more mutants at this position: Asp158Gly, Asp158Ala, and Asp158Glu. From the pH-activity profiles of these and other mutants of papain, it has been possible to develop a model that enables us to dissect out the contribution of the various mutations toward (i) intrinsic activity, (ii) ion pair stability, and (iii) the electrostatic potential at the active site. Results obtained with mutants that place either Gly or Ala at position 158 indicate that the hydrogen bonds involving the side chain of Asp158 in wild-type papain are indirectly important for enzyme activity. When CBZ-Phe-Arg-MCA is used as a substrate, the (kcat/KM)obs values at pH 6.5 are 3650 and 494 M-1 s-1 for Asp158Gly and Asp158Ala, respectively, as compared to 119,000 M-1 s-1 for papain. Results with the Asp158Glu mutant suggest that the side chain of Glu moves closer to the active site and cannot form hydrogen bonds similar to those involving Asp158 in papain. From the four mutations introduced at position 158 in papain, we can conclude that it is not the charge but the hydrogen-bonding interactions involving the side chain of Asp158 that contribute the most to the stabilization of the thiolate-imidazolium ion pair in papain. However, the charge and the hydrogen bonds of Asp158 both contribute to the intrinsic activity of the enzyme.     

Hydrogen bonding stabilizes globular proteins.

It is clear that intramolecular hydrogen bonds are essential to the structure and stability of globular proteins. It is not clear, however, whether they make a net favorable contribution to this stability. Experimental and theoretical studies are at odds over this important question. Measurements of the change in conformational stability, delta (delta G), for the mutation of a hydrogen bonded residue to one incapable of hydrogen bonding suggest a stabilization of 1.0 kcal/mol per hydrogen bond. If the delta (delta G) values are corrected for differences in side-chain hydrophobicity and conformational entropy, then the estimated stabilization becomes 2.2 kcal/mol per hydrogen bond. These and other experimental studies discussed here are consistent and compelling: hydrogen bonding stabilizes globular proteins.     

Relationship between intramolecular hydrogen bonding and solvent accessibility of side-chain donors and acceptors in proteins

This study shows that intramolecular hydrogen bonding in proteins depends on the accessibility of donors and acceptors to water molecules. The frequency of occurrence of H-bonded side chains in proteins is inversely proportional to the solvent accessibility of their donors and acceptors. Estimates of the notional free energy of hydrogen bonding suggest that intramolecular hydrogen-bonding interactions of buried and half-buried donors and acceptors can contribute favorably to the stability of a protein, whereas those of solvent-exposed polar atoms become less favorable or unfavorable.     

A view of the three dimensional structure of globular proteins

A view of the three dimensional structure of globular proteins based on continuous networks of hydrogen bonds is proposed. Active sites of enzymes and ion sites are prominent and, within the networks, there are islands of hydrophobic regions giving an overall piebald effect to the appearance of the molecule. This point of view was originally suggested by the results of quantum mechanical computations on the coupling between hydrogen bonds. A formalism for the total energy of a globular protein in water is also suggested.The study of five lines of experimental evidence supports this suggestion. The analysis of the experimental X-ray data for ten globular proteins, using the NETWORK program, revealed the existence of these hydrogen bond networks; X-ray data showed that water molecules tend to occupy fixed positions relative to the protein molecule; a survey has shown that water molecules tend to occupy specific positions relative to the hydrogen bonding side chains; experimental evidence on the bulk properties of lysozyme showed that there exist tightly bound water molecules; graphics studies of the ribonucleaseA molecule demonstrated the networks and the piebald effect. This point of view is pictorially simple and, to illustrate the use of such networks, we discuss the simple ion pairs which occur as substructures within the networks.     

Amino-acid solvation structure in transmembrane helices from molecular dynamics simulations

Understanding the solvation of amino acids in biomembranes is an important step to better explain membrane protein folding. Several experimental studies have shown that polar residues are both common and important in transmembrane segments, which means they have to be solvated in the hydrophobic membrane, at least until helices have aggregated to form integral proteins. In this work, we have used computer simulations to unravel these interactions on the atomic level, and classify intramembrane solvation properties of amino acids. Simulations have been performed for systematic mutations in poly-Leu helices, including not only each amino acid type, but also every z-position in a model helix. Interestingly, many polar or charged residues do not desolvate completely, but rather retain hydration by snorkeling or pulling in water/headgroups--even to the extent where many of them exist in a microscopic polar environment, with hydration levels corresponding well to experimental hydrophobicity scales. This suggests that even for polar/charged residues a large part of solvation cost is due to entropy, not enthalpy loss. Both hydration level and hydrogen bonding exhibit clear position-dependence. Basic side chains cause much less membrane distortion than acidic, since they are able to form hydrogen bonds with carbonyl groups instead of water or headgroups. This preference is supported by sequence statistics, where basic residues have increased relative occurrence at carbonyl z-coordinates. Snorkeling effects and N-/C-terminal orientation bias are directly observed, which significantly reduces the effective thickness of the hydrophobic core. Aromatic side chains intercalate efficiently with lipid chains (improving Trp/Tyr anchoring to the interface) and Ser/Thr residues are stabilized by hydroxyl groups sharing hydrogen bonds to backbone oxygens.     

The transmembrane 7-alpha-bundle of rhodopsin: distance geometry calculations with hydrogen bonding constraints.

A 3D model of the transmembrane 7-alpha-bundle of rhodopsin-like G-protein-coupled receptors (GPCRs) was calculated using an iterative distance geometry refinement with an evolving system of hydrogen bonds, formed by intramembrane polar side chains in various proteins of the family and collectively applied as distance constraints. The alpha-bundle structure thus obtained provides H bonding of nearly all buried polar side chains simultaneously in the 410 GPCRs considered. Forty evolutionarily conserved GPCR residues form a single continuous domain, with an aliphatic "core" surrounded by six clusters of polar and aromatic side chains. The 7-alpha-bundle of a specific GPCR can be calculated using its own set of H bonds as distance constraints and the common "average" model to restrain positions of the helices. The bovine rhodopsin model thus determined is closely packed, but has a few small polar cavities, presumably filled by water, and has a binding pocket that is complementary to 11-cis (6-s-cis, 12-s-trans, C = N anti)-retinal or to all-trans-retinal, depending on conformations of the Lys296 and Trp265 side chains. A suggested mechanism of rhodopsin photoactivation, triggered by the cis-trans isomerization of retinal, involves rotations of Glu134, Tyr223, Trp265, Lys296, and Tyr306 side chains and rearrangement of their H bonds. The model is in agreement with published electron cryomicroscopy, mutagenesis, chemical modification, cross-linking, Fourier transform infrared spectroscopy, Raman spectroscopy, electron paramagnetic resonance spectroscopy, NMR, and optical spectroscopy data. The rhodopsin model and the published structure of bacteriorhodopsin have very similar retinal-binding pockets.     

Altered hydrogen bonding of Arg82 during the proton pump cycle of bacteriorhodopsin: a low-temperature polarized FTIR spectroscopic study

Light-driven proton transport in bacteriorhodopsin (BR) is achieved by dynamic rearrangement of the hydrogen-bonding network inside the membrane protein. Arg82 is located between the Schiff base region and proton release group, and has a major influence on the pK(a) values of these groups. It is believed that Arg82 changes its hydrogen-bonding acceptors during the pump cycle of BR, stages of which are correlated with proton movement along the transport pathway. In this study, we compare low-temperature polarized FTIR spectra of [eta(1,2)-(15)N]arginine-labeled BR in the 2750-2000 cm(-1) region with those of unlabeled BR for the K, L, M, and N intermediates. In the K-minus-BR difference spectra, (15)N-shifted modes were found at 2292 (-)/2266 (+) cm(-1) and at 2579 (-)/2567 (+) cm(-1). The former corresponds to strong hydrogen bonding, while the latter corresponds to very weak hydrogen bonding. Both N-D stretches probably originate from Arg82, the former oriented toward water 406 and the latter toward the extracellular side, and both hydrogen bonds are somewhat strengthened upon retinal photoisomerization. This perturbation of arginine hydrogen bonding is entirely relaxed in the L intermediate where no (15)N-isotope shifts are observed in the difference spectrum. In the M intermediate, the frequency is not significantly altered from that in BR. However, the polarized FTIR spectra strongly suggest that the dipolar orientation of the strongly hydrogen bonded N-D group of Arg82 is changed from perpendicular to parallel to the membrane plane. Such a change is presumably related to the motion of the Arg82 side chain from the Schiff base region to the extracellular proton release group. Additional bands corresponding to weak hydrogen bonding were observed in both the M-minus-BR and N-minus-BR spectra. Changes in hydrogen-bonding structures involving Arg82 are discussed on the basis of these FTIR observations.     

The mechanism of glycerol conduction in aquaglyceroporins.

BACKGROUND: The E. coli glycerol facilitator, GlpF, selectively conducts glycerol and water, excluding ions and charged solutes. The detailed mechanism of the glycerol conduction and its relationship to the characteristic secondary structure of aquaporins and to the NPA motifs in the center of the channel are unknown. RESULTS: Molecular dynamics simulations of GlpF reveal spontaneous glycerol and water conduction driven, on a nanosecond timescale, by thermal fluctuations. The bidirectional conduction, guided and facilitated by the secondary structure, is characterized by breakage and formation of hydrogen bonds for which water and glycerol compete. The conduction involves only very minor changes in the protein structure, and cooperativity between the GlpF monomers is not evident. The two conserved NPA motifs are strictly linked together by several stable hydrogen bonds and their asparagine side chains form hydrogen bonds with the substrates passing the channel in single file. CONCLUSIONS: A complete conduction of glycerol through the GlpF was deduced from molecular dynamics simulations, and key residues facilitating the conduction were identified. The nonhelical parts of the two half-membrane-spanning segments expose carbonyl groups towards the channel interior, establishing a curve-linear pathway. The conformational stability of the NPA motifs is important in the conduction and critical for selectivity. Water and glycerol compete in a random manner for hydrogen bonding sites in the protein, and their translocations in single file are correlated. The suggested conduction mechanism should apply to the whole family.     

Hydrogen bonds in protein-DNA complexes: where geometry meets plasticity

Recognition of a DNA sequence by a protein is achieved by interface-coupled chemical and shape complementation. This complementation between the two molecules is clearly directional and is determined by the specific chemical contacts including mainly hydrogen bonds. Directionality is an instrumental property of hydrogen bonding as it influences molecular conformations, which also affects DNA-protein recognition. The prominent elements in the recognition of a particular DNA sequence by a protein are the hydrogen-bond donors and acceptors of the base pairs into the grooves of the DNA that must interact with complementary moieties of the protein partner. Protein side chains make most of the crucial contacts through bidentate and complex hydrogen-bonding interactions with DNA base edges hence conferring remarkable specificity.     

Forces stabilizing proteins

The goal of this article is to summarize what has been learned about the major forces stabilizing proteins since the late 1980s when site-directed mutagenesis became possible. The following conclusions are derived from experimental studies of hydrophobic and hydrogen bonding variants. (1) Based on studies of 138 hydrophobic interaction variants in 11 proteins, burying a –CH₂− group on folding contributes 1.1 ± 0.5 kcal/mol to protein stability. (2) The burial of non-polar side chains contributes to protein stability in two ways: first, a term that depends on the removal of the side chains from water and, more importantly, the enhanced London dispersion forces that result from the tight packing in the protein interior. (3) Based on studies of 151 hydrogen bonding variants in 15 proteins, forming a hydrogen bond on folding contributes 1.1 ± 0.8 kcal/mol to protein stability. (4) The contribution of hydrogen bonds to protein stability is strongly context dependent. (5) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (6) Polar group burial can make a favorable contribution to protein stability even if the polar group is not hydrogen bonded. (7) Hydrophobic interactions and hydrogen bonds both make large contributions to protein stability.     

Side-chain conformations in 4-alpha-helical bundles.

The distribution of the chi(1), chi(2) dihedral angles in a dataset consisting of 12 unrelated 4-alpha-helical bundle proteins was determined and qualitatively compared with that observed in globular proteins. The analysis suggests that the 4-alpha-helical bundle motif could occasionally impose steric constraints on side chains: (i) the side-chain conformations are limited to only a subset of the conformations observed in globular proteins and for some amino acids they are sterically more constrained than those in helical regions of globular proteins; (ii) aspartic acid and asparagine occasionally adopt rotamers that have not been previously reported for globular or helical proteins; (iii) some rotamers of tyrosine and isoleucine are predominantly or exclusively associated with hydrophobic core positions (a, d); (iv) mutations in the hydrophobic core occur preferentially between residue types which among other physicochemical properties also share a predominant rotamer.     

A simple and novel interpretation of the three-dimensional structure of globular proteins based on quantum-mechanical computations on small model molecules. I

It is suggested that the three-dimensional structure of globular proteins is partly determined by a framework of strengthened hydrogen bonds that involves both ionic side chains and water molecules in addition to the polypeptide backbone. This conclusion follows from a combination of the results of ab initio molecular-orbital computations on small model molecules and high-accuracy x-ray data on the rubredoxin molecule. The computations yield the idea of hydrogen-bonded bridges that are built from tens of atoms, and the experimental information yields the idea that the bridges are assembled into clusters, each of which is built from hundreds of atoms. Some 10 such clusters then form a globular protein.     

Interhelical hydrogen bonds in the CFTR membrane domain.

Critical mutations in the membrane-spanning domains of proteins cause many human diseases. We report the expression in Escherichia coli of helix-loop-helix segments of the cystic fibrosis transmembrane conductance regulator (CFTR) chloride channel domain in milligram quantities. Analysis of gel migration patterns of these constructs, in conjunction with circular dichroism spectroscopy, demonstrate that a neutral-to-charged, CF-phenotypic point mutation of a hydrophobic residue (V232D) in the CFTR transmembrane (TM) helix 4 induces a hydrogen bond with neighboring wild type Gln 207 in TM helix 3. As an electrostatic crosslink within a hydrocarbon phase, such a hydrogen bond could alter the normal assembly and alignment of CFTR TM helices and/or impede their movement in response to substrate transport. Our results imply that membrane proteins may be vulnerable to loss of function through formation of membrane-buried interhelical hydrogen bonds by partnering of proximal polar side chains.     

Structural and mutational studies of the carboxylate cluster in iron-free ribonucleotide reductase R2

The R2 protein of ribonucleotide reductase features a di-iron site deeply buried in the protein interior. The apo form of the R2 protein has an unusual clustering of carboxylate side chains at the empty metal-binding site. In a previous study, it was found that the loss of the four positive charge equivalents of the diferrous site in the apo protein appeared to be compensated for by the protonation of two histidine and two carboxylate side chains. We have studied the consequences of removing and introducing charged residues on the local hydrogen-bonding pattern in the region of the carboxylate cluster of Corynebacterium ammoniagenes and Escherichia coli protein R2 using site-directed mutagenesis and X-ray crystallography. The structures of the metal-free forms of wild-type C. ammoniagenes R2 and the mutant E. coli proteins D84N, S114D, E115A, H118A, and E238A have been determined and their hydrogen bonding and protonation states have been structurally assigned as far as possible. Significant alterations to the hydrogen-bonding patterns, protonation states, and hydration is observed for all mutant E. coli apo proteins as compared to wild-type apo R2. Further structural variations are revealed by the wild-type apo C. ammoniagenes R2 structure. The protonation and hydration effects seen in the carboxylate cluster appear to be due to two major factors: conservation of the overall charge of the site and the requirement of electrostatic shielding of clustered carboxylate residues. Very short hydrogen-bonding distances between some protonated carboxylate pairs are indicative of low-barrier hydrogen bonding.     

Crystal structure of Bordetella pertussis BugD solute receptor unveils the basis of ligand binding in a new family of periplasmic binding proteins

Periplasmic binding proteins of a new family particularly well represented in Bordetella pertussis have been called Bug receptors. One B.pertussis Bug protein is part of a tripartite tricarboxylate transporter while the functions of the other 77 are unknown. We report the first structure of a Bug receptor, BugD. It adopts the characteristic Venus flytrap motif observed in other periplasmic binding proteins, with two globular domains bisected by a deep cleft. BugD displays a closed conformation resulting from the fortuitous capture of a ligand, identified from the electron density as an aspartate. The structure reveals a distinctive alpha carboxylate-binding motif, involving two water molecules that bridge the carboxylate oxygen atoms to the protein. Both water molecules are hydrogen bonded to a common carbonyl group from Ala14, and each forms a hydrogen bond with one carboxylate oxygen atom of the ligand. Additional hydrogen bonds are found between the ligand alpha carboxylate oxygen atoms and protein backbone amide groups and with a threonine hydroxyl group. This specific ligand-binding motif is highly conserved in Bug proteins, indicating that they may all be receptors of amino acids or other carboxylated solutes, with a similar binding mode. The present structure thus unveils the bases of ligand binding in this large family of periplasmic binding proteins, several hundred members of which have been identified in various bacterial species.