Contribution of intra- and intermolecular hydrogen bonds to the conformational stability of human lysozyme(,).

In globular proteins, there are intermolecular hydrogen bonds between protein and water molecules, and between water molecules, which are bound with the proteins, in addition to intramolecular hydrogen bonds. To estimate the contribution of these hydrogen bonds to the conformational stability of a protein, the thermodynamic parameters for denaturation and the crystal structures of five Thr to Val and five Thr to Ala mutant human lysozymes were determined. The denaturation Gibbs energy (DeltaG) of Thr to Val and Thr to Ala mutant proteins was changed from 4.0 to -5.6 kJ/mol and from 1.6 to -6.3 kJ/mol, respectively, compared with that of the wild-type protein. The contribution of hydrogen bonds to the stability (DeltaDeltaG(HB)) of the Thr and other mutant human lysozymes previously reported was extracted from the observed stability changes (DeltaDeltaG) with correction for changes in hydrophobicity and side chain conformational entropy between the wild-type and mutant structures. The estimation of the DeltaDeltaG(HB) values of all mutant proteins after removal of hydrogen bonds, including protein-water hydrogen bonds, indicates a favorable contribution of the intra- and intermolecular hydrogen bonds to the protein stability. The net contribution of an intramolecular hydrogen bond (DeltaG(HB[pp])), an intermolecular one between protein and ordered water molecules (DeltaG(HB[pw])), and an intermolecular one between ordered water molecules (DeltaG(HB[ww])) could be estimated to be 8. 5, 5.2, and 5.0 kJ/mol, respectively, for a 3 A long hydrogen bond. This result shows the different contributions to protein stability of intra- and intermolecular hydrogen bonds. The entropic cost due to the introduction of a water molecule (DeltaG(H)()2(O)) could be also estimated to be about 8 kJ/mol.     

Role of surface hydrophobic residues in the conformational stability of human lysozyme at three different positions

To evaluate the contribution of the amino acid residues on the surface of a protein to its stability, a series of hydrophobic mutant human lysozymes (Val to Gly, Ala, Leu, Ile, Met, and Phe) modified at three different positions on the surface, which are located in the alpha-helix (Val 110), the beta-sheet (Val 2), and the loop (Val 74), were constructed. Their thermodynamic parameters of denaturation and crystal structures were examined by calorimetry and by X-ray crystallography at 100 K, respectively. Differences in the denaturation Gibbs energy change between the wild-type and the hydrophobic mutant proteins ranged from 4.6 to -9.6 kJ/mol, 2.7 to -1.5 kJ/mol, and 3.6 to -0.2 kJ/mol at positions 2, 74, and 110, respectively. The identical substitution at different positions and different substitutions at the same position resulted in different degrees of stabilization. Changes in the stability of the mutant proteins could be evaluated by a unique equation considering the conformational changes due to the substitutions [Funahashi et al. (1999) Protein Eng. 12, 841-850]. For this calculation, secondary structural propensities were newly considered. However, some mutant proteins were not adapted to the equation. The hydration structures around the mutation sites of the exceptional mutant proteins were affected due to the substitutions. The stability changes in the exceptional mutant proteins could be explained by the formation or destruction of the hydration structures. These results suggest that the hydration structure mediated via hydrogen bonds covering the protein surface plays an important role in the conformational stability of the protein.     

Contribution of polar groups in the interior of a protein to the conformational stability

It has been generally believed that polar residues are usually located on the surface of protein structures. However, there are many polar groups in the interior of the structures in reality. To evaluate the contribution of such buried polar groups to the conformational stability of a protein, nonpolar to polar mutations (L8T, A9S, A32S, I56T, I59T, I59S, A92S, V93T, A96S, V99T, and V100T) in the interior of a human lysozyme were examined. The thermodynamic parameters for denaturation were determined using a differential scanning calorimeter, and the crystal structures were analyzed by X-ray crystallography. If a polar group had a heavy energy cost to be buried, a mutant protein would be remarkably destabilized. However, the stability (Delta G) of the Ala to Ser and Val to Thr mutant human lysozymes was comparable to that of the wild-type protein, suggesting a low-energy penalty of buried polar groups. The structural analysis showed that all polar side chains introduced in the mutant proteins were able to find their hydrogen bond partners, which are ubiquitous in protein structures. The empirical structure-based calculation of stability change (Delta Delta G) [Takano et al. (1999) Biochemistry 38, 12698--12708] revealed that the mutant proteins decreased the hydrophobic effect contributing to the stability (Delta G(HP)), but this destabilization was recovered by the hydrogen bonds newly introduced. The present study shows the favorable contribution of polar groups with hydrogen bonds in the interior of protein molecules to the conformational stability.     

Experimental verification of the 'stability profile of mutant protein' (SPMP) data using mutant human lysozymes.

The stability profile of mutant protein (SPMP) (Ota,M., Kanaya,S. and Nishikawa,K., 1995, J. Mol. Biol., 248, 733-738) estimates the changes in conformational stability due to single amino acid substitutions using a pseudo-energy potential developed for evaluating structure-sequence compatibility in the structure prediction method, the 3D-1D compatibility evaluation. Nine mutant human lysozymes expected to significantly increase in stability from SPMP were constructed, in order to experimentally verify the reliability of SPMP. The thermodynamic parameters for denaturation and crystal structures of these mutant proteins were determined. One mutant protein was stabilized as expected, compared with the wild-type protein. However, the others were not stabilized even though the structural changes were subtle, indicating that SPMP overestimates the increase in stability or underestimates negative effects due to substitution. The stability changes in the other mutant human lysozymes previously reported were also analyzed by SPMP. The correlation of the stability changes between the experiment and prediction depended on the types of substitution: there were some correlations for proline mutants and cavity-creating mutants, but no correlation for mutants related to side-chain hydrogen bonds. The present results may indicate some additional factors that should be considered in the calculation of SPMP, suggesting that SPMP can be refined further.     

Large-scale networks of hydration water molecules around proteins investigated by cryogenic X-ray crystallography.

The structures at protein-water interface, i.e. the hydration structure of proteins, have been investigated by cryogenic X-ray crystal structure analyses. Hydration structures appeared far clearer at cryogenic temperature than at ambient temperature, presumably because the motions of hydration water molecules were quenched by cooling. Based on the structural models obtained, the hydration structures were systematically analyzed with respect to the amount of water molecules, the interaction modes between water molecules and proteins, the local and the global distribution of them on the surface of proteins. The standard tetrahedral interaction geometry of water in bulk retained at the interface and enabled the three-dimensional chain connection of hydrogen bonds between hydration water molecules and polar protein atoms. Large-scale networks of hydrogen bonds covering the entire surface of proteins were quite flexible to accommodate to the large-scale conformational changes of proteins and seemed to have great influences on the dynamics and function of proteins. The present observation may provide a new concept for discussing the dynamics of proteins in aqueous solution.     

Are the parameters of various stabilization factors estimated from mutant human lysozymes compatible with other proteins?

The various factors which contribute to protein stability have been extensively examined using mutant proteins, but the same kinds of substitutions have given different results depending on the substitution sites. Recently, the contributions of some stabilization factors have been quantitatively derived as parameters by a unique equation, considering the conformational changes due to the mutations using mutant human lysozymes [Funahashi et al. (1999) Protein ENG: 12, 841-850]. To evaluate these parameters estimated from the mutant human lysozymes, stability-structure datasets for the mutant T4 lysozymes were selected. The stabilities for the mutant T4 lysozymes could be roughly estimated using these parameters. Notable differences between the estimated and experimental stabilities were caused by the uncertainty in part of the structures due to some Arg and Lys residues fluctuating on the surface of the T4 lysozyme. Excluding these atoms from the estimation gave a good correlation between the estimated and experimental stabilities. These results suggest that the parameters of the various stabilization factors derived from the mutant human lysozymes are compatible with the mutant T4 lysozymes, although they should be improved with respect to some points using more information.     

Role of amino acid residues at turns in the conformational stability and folding of human lysozyme

To clarify the role of amino acid residues at turns in the conformational stability and folding of a globular protein, six mutant human lysozymes deleted or substituted at turn structures were investigated by calorimetry, GuHCl denaturation experiments, and X-ray crystal analysis. The thermodynamic properties of the mutant and wild-type human lysozymes were compared and discussed on the basis of their three-dimensional structures. For the deletion mutants, Delta47-48 and Delta101, the deleted residues are in turns on the surface and are absent in human alpha-lactalbumin, which is homologous to human lysozyme in amino acid sequence and tertiary structure. The stability of both mutants would be expected to increase due to a decrease in conformational entropy in the denatured state; however, both proteins were destabilized. The destabilizations were mainly caused by the disappearance of intramolecular hydrogen bonds. Each part deleted was recovered by the turn region like the alpha-lactalbumin structure, but there were differences in the main-chain conformation of the turn between each deletion mutant and alpha-lactalbumin even if the loop length was the same. For the point mutants, R50G, Q58G, H78G, and G37Q, the main-chain conformations of these substitution residues located in turns adopt a left-handed helical region in the wild-type structure. It is thought that the left-handed non-Gly residue has unfavorable conformational energy compared to the left-handed Gly residue. Q58G was stabilized, but the others had little effect on the stability. The structural analysis revealed that the turns could rearrange the main-chain conformation to accommodate the left-handed non-Gly residues. The present results indicate that turn structures are able to change their main-chain conformations, depending upon the side-chain features of amino acid residues on the turns. Furthermore, stopped-flow GuHCl denaturation experiments on the six mutants were performed. The effects of mutations on unfolding-refolding kinetics were significantly different among the mutant proteins. The deletion/substitutions in turns located in the alpha-domain of human lysozyme affected the refolding rate, indicating the contribution of turn structures to the folding of a globular protein.     

Contribution of hydrogen bonds to the conformational stability of human lysozyme: calorimetry and X-ray analysis of six Ser --> Ala mutants.

To further examine the contribution of hydrogen bonds to the conformational stability of the human lysozyme, six Ser to Ala mutants were constructed. The thermodynamic parameters for denaturation of these six Ser mutant proteins were investigated by differential scanning calorimetry (DSC), and the crystal structures were determined by X-ray analysis. The denaturation Gibbs energy (DeltaG) of the Ser mutant proteins was changed from 2.0 to -5.7 kJ/mol, compared to that of the wild-type protein. With an analysis in which some factors that affected the stability due to mutation were considered, the contribution of hydrogen bonds to the stability (Delta DeltaGHB) was extracted on the basis of the structures of the mutant proteins. The results showed that hydrogen bonds between protein atoms and between a protein atom and a water bound with the protein molecule favorably contribute to the protein stability. The net contribution of one intramolecular hydrogen bond to protein stability (DeltaGHB) was 8.9 +/- 2.6 kJ/mol on average. However, the contribution to the protein stability of hydrogen bonds between a protein atom and a bound water molecule was smaller than that for a bond between protein atoms.     

Contribution of amino acid substitutions at two different interior positions to the conformational stability of human lysozyme.

To elucidate correlative relationships between structural change and thermodynamic stability in proteins, a series of mutant human lysozymes modified at two buried positions (Ile56 and Ile59) were examined. Their thermodynamic parameters of denaturation and crystal structures were studied by calorimetry and X-ray crystallography. The mutants at positions 56 and 59 exhibited different responses to a series of amino acid substitutions. The changes in stability due to substitutions showed a linear correlation with changes in hydrophobicity of substituted residues, having different slopes at each mutation site. However, the stability of each mutant was found to be represented by a unique equation involving physical properties calculated from mutant structures. By fitting present and previous stability data for mutant human lysozymes substituted at various positions to the equation, the magnitudes of the hydrophobicity of a carbon atom and the hydrophobicity of nitrogen and neutral oxygen atoms were found to be 0.178 and -0.013 kJ/mol.A(2), respectively. It was also found that the contribution of a hydrogen bond with a length of 3.0 A to protein stability was 5.1 kJ/mol and the entropy loss of newly introduction of a water molecules was 7.8 kJ/mol.     

Hydration in proteins observed by high-resolution neutron crystallography

It is well known that water molecules surrounding a protein play important roles in maintaining its structural stability. Water molecules are known to participate in several physiological processes through the formation of hydrogen bonds. However, the hydration structures of most proteins are not known well at an atomic level at present because X-ray protein crystallography has difficulties to localize hydrogen atoms. In contrast, neutron crystallography has no problem in determining the position of hydrogens with high accuracy.1 In this article, the hydration structures of three proteins are described- myoglobin, wild-type rubredoxin, and a mutant rubredoxin-the structures of which were solved at 1.5- or 1.6-A resolution by neutron structure determination. These hydration patterns show fascinating features and the water molecules adopt a variety of shapes in the neutron Fourier maps, revealing details of intermolecular hydrogen bond formation and dynamics of hydration. Our results further show that there are strong relationships between these shapes and the water environments.     

A new scale for side-chain contribution to protein stability based on the empirical stability analysis of mutant proteins

The hydrophobicity scales for amino acid side chains based on the transfer Gibbs energy (DeltaG(trans)) of amino acids from non-aqueous phases to water have been widely used to estimate the contribution of buried side chains to the conformational stability of proteins. In this paper, we propose a new scale for the side-chain contribution to protein stability, which is derived from data on protein denaturation experiments using systematic and comprehensive mutant proteins. In the experiments, the contribution of some physical properties were quantitatively determined as parameters in a unique equation representing the stability change (DeltaDeltaG) of mutant proteins as a function of the structural changes due to the mutations. These parameters are able conveniently to provide a scale for the side-chain contribution to protein stability. This new scale also has the advantage over the previously reported hydrophobicity scales of residues with the contributions of hydrogen bonds or secondary structural propensity. It may find practical application in algorithms for the prediction of protein structures.     

Structural studies of mutants of T4 lysozyme that alter hydrophobic stabilization

Multiple replacements at amino acid position 3 of bacteriophage T4 lysozyme have shown that the conformational stability of the protein is directly governed by the hydrophobicity of the residue substituted (Matsumura, M., Becktel, W. J., and Matthews, B. W. (1988) Nature 334, 406-410). Of the 13 mutant lysozymes made by site-directed mutagenesis, two variants, one with valine (I3V) and the other with tyrosine (I3Y), were crystallized and their structures solved. In this report we describe the crystal structures of these variants at 1.7 A resolution. While the structure of the I3V mutant is essentially the same as that of wild-type lysozyme, the I3Y mutant has substantial changes in its structure. The most significant of these are that the side chain of the tyrosine is not accommodated within the interior of the protein and the amino-terminal polypeptide (residues 1-9) moves 0.6-1.1 A relative to the wild-type structure. Using coordinates based on the wild-type and available mutant structures, solvent accessible surface area of residue 3 as well as the adjacent 9 residues in the folded form were calculated. The free energy of stabilization based on the transfer of these residues from a fully extended form to the interior to the folded protein was found to correlate well with the protein stability determined by thermodynamic analysis. The enhanced thermostability of the variant Ile-3----Leu, relative to wild-type lysozyme, can also be rationalized by surface-area calculations based on a model-built structure. Noncrystallization of most lysozyme variants at position 3 appears to be due to disruption of intermolecular contacts in the crystal. The Ile-3----Val variant is closely isomorphous with wild-type and maintains the same crystal contacts. In the Ile-3----Tyr variant, however, a new set of contacts is made in which direct protein-protein hydrogen bonds are replaced by protein-water-protein hydrogen bonds as well as a novel hydrogen bond involving the phenolic hydroxyl of the substituted tyrosine.     

Conservative mutation Met8 --> Leu affects the folding process and structural stability of squash trypsin inhibitor CMTI-I

Protein molecules can accommodate a large number of mutations without noticeable effects on their stability and folding kinetics. On the other hand, some mutations can have quite strong effects on protein conformational properties. Such mutations either destabilize secondary structures, e.g., alpha-helices, are incompatible with close packing of protein hydrophobic cores, or lead to disruption of some specific interactions such as disulfide cross links, salt bridges, hydrogen bonds, or aromatic-aromatic contacts. The Met8 --> Leu mutation in CMTI-I results in significant destabilization of the protein structure. This effect could hardly be expected since the mutation is highly conservative, and the side chain of residue 8 is situated on the protein surface. We show that the protein destabilization is caused by rearrangement of a hydrophobic cluster formed by side chains of residues 8, Ile6, and Leu17 that leads to partial breaking of a hydrogen bond formed by the amide group of Leu17 with water and to a reduction of a hydrophobic surface buried within the cluster. The mutation perturbs also the protein folding. In aerobic conditions the reduced wild-type protein folds effectively into its native structure, whereas more then 75% of the mutant molecules are trapped in various misfolded species. The main conclusion of this work is that conservative mutations of hydrophobic residues can destabilize a protein structure even if these residues are situated on the protein surface and partially accessible to water. Structural rearrangement of small hydrophobic clusters formed by such residues can lead to local changes in protein hydration, and consequently, can affect considerably protein stability and folding process.     

Conformational transitions in eosinophil cationic protein: a molecular dynamics study in aqueous environment

Extensive molecular dynamics simulations have been performed on eosinophil cationic protein (ECP). The two structures found in the crystallographic dimer (ECPA and ECPB) have been independently simulated. A small difference in the pattern of the sidechain hydrogen bonds in the starting structure has resulted in interesting differences in the conformations accessed during the simulations. In one simulation (ECPB), a stable equilibrium conformation was obtained and in the other (ECPA), conformational transitions at the level of sidechain interactions were observed. The conformational transitions exhibit the involvement of the solvent (water) molecules with a pore-like construct in the equilibrium conformation and an opening for a large number of water molecules during the transition phase. The details of these transitions are examined in terms of intra-protein hydrogen bonds, protein-water networks and the residence times of water molecules on the polar atoms of the protein. These properties show some significant differences in the region between the N-terminal helix and the loop before the C-terminal strand as a function of different conformations accessed during the simulations. However, the stable hydrogen bonds, the protein-water networks, and the hydration patterns in most part of the protein including the active site are very much similar in both the simulations, indicating the fact that these are intrinsic properties of proteins.     

Water-protein interactions from high-resolution protein crystallography

To understand the role of water in life at molecular and atomic levels, structures and interactions at the protein-water interface have been investigated by cryogenic X-ray crystallography. The method enabled a much clearer visualization of definite hydration sites on the protein surface than at ambient temperature. Using the structural models of proteins, including several hydration water molecules, the characteristics in hydration structures were systematically analysed for the amount, the interaction geometries between water molecules and proteins, and the local and global distribution of water molecules on the surface of proteins. The tetrahedral hydrogen-bond geometry of water molecules in bulk solvent was retained at the interface and enabled the extension of a three-dimensional chain connection of a hydrogen-bond network among hydration water molecules and polar protein atoms over the entire surface of proteins. Networks of hydrogen bonds were quite flexible to accommodate and/or to regulate the conformational changes of proteins such as domain motions. The present experimental results may have profound implications in the understanding of the physico-chemical principles governing the dynamics of proteins in an aqueous environment and a discussion of why water is essential to life at a molecular level.     

Crystal structures of the apo- and holomutant human lysozymes with an introduced Ca2+ binding site

The three-dimensional structures of apo- and holomutant human lysozymes (D86/92 lysozyme), in which a calcium binding site was designed and created for enhancing molecular stability by replacing both Gln86 and Ala92 with aspartic acids, were refined at 1.8-A resolution by x-ray crystallography. The overall structures and crystallographic thermal factors of all three proteins, the apo-, holo-D86/92, and the wild-type human lysozymes, were essentially identical; these results showed that the introduction of the calcium binding site did not affect either the overall structure or molecular rigidity of the proteins. However, structure analyses of the apo-D86/92 lysozyme revealed that the mutations affected the side chain conformation of residue 86 and hydrogen networks between the protein and the internal solvent molecules. In the structure of the holo-D86/92 lysozyme, seven oxygen ligands formed a slightly distorted pentagonal bipyramid around the calcium ion, indicating that the coordination around the calcium ion was quite similar to that in baboon alpha-lactalbumin. The pentagonal bipyramid coordination could be one of the most widely found and appropriate calcium binding schemes in proteins.     

Contribution of salt bridges near the surface of a protein to the conformational stability

Salt bridges play important roles in the conformational stability of proteins. However, the effect of a surface salt bridge on the stability remains controversial even today; some reports have shown little contribution of a surface salt bridge to stability, whereas others have shown a favorable contribution. In this study, to elucidate the net contribution of a surface salt bridge to the conformational stability of a protein, systematic mutant human lysozymes, containing one Glu to Gln (E7Q) and five Asp to Asn mutations (D18N, D49N, D67N, D102N, and D120N) at residues where a salt bridge is formed near the surface in the wild-type structure, were examined. The thermodynamic parameters for denaturation between pH 2.0 and 4.8 were determined by use of a differential scanning calorimeter, and the crystal structures were analyzed by X-ray crystallography. The denaturation Gibbs energy (DeltaG) of all mutant proteins was lower than that of the wild-type protein at pH 4, whereas there was little difference between them near pH 2. This is caused by the fact that the Glu and Asp residues are ionized at pH 4 but protonated at pH 2, indicating a favorable contribution of salt bridges to the wild-type structure at pH 4. Each contribution was not equivalent, but we found that the contributions correlate with the solvent inaccessibility of the salt bridges; the salt bridge contribution was small when 100% accessible, while it was about 9 kJ/mol if 100% inaccessible. This conclusion indicates how to reconcile a number of conflicting reports about role of surface salt bridges in protein stability. Furthermore, the effect of salts on surface salt bridges was also examined. In the presence of 0.2 M KCl, the stability at pH 4 decreased, and the differences in stability between the wild-type and mutant proteins were smaller than those in the absence of salts, indicating the compensation to the contribution of salt bridges with salts. Salt bridges with more than 50% accessibility did not contribute to the stability in the presence of 0.2 M KCl.     

Effects of disulfide bonds on compactness of protein molecules revealed by volume, compressibility, and expansibility changes during reduction

To elucidate the effects of disulfide bonds on the compactness of protein molecules, the partial specific volume (v(o)) and coefficients of adiabatic compressibility (beta(s)(o)) and thermal expansibility (alpha) of five globular proteins (ovalbumin, beta-lactoglobulin, lysozyme, ribonuclease A, and bovine serum albumin) were measured in aqueous solutions with pH values of 7 and 2 at 25 degrees C when their disulfide bonds were totally reduced by carboxamidomethylation. Circular dichroism and fluorescence spectra show that the secondary and tertiary structures are partly disrupted by reduction, depending on the number of disulfide bonds in the proteins and the pH of the medium. The conformational changes are accompanied by decreases in v(o) and beta(s)(o) and by an increase in alpha, indicating that reduction decreases the internal cavity and increases surface hydration. The beta(s)(o) values of native or oxidized proteins decrease, and the effects of reduction on the volumetric parameters become more significant as the number of disulfide bonds increases and as they are formed over a larger distance in the primary structure. These results demonstrate that disulfide bonds play an important role, mainly via entropic forces, in the three-dimensional structure and compactness of protein molecules.     

The mechanism of the type III antifreeze protein action: a computational study

The random network model of water quantitatively describes the different hydration heat capacities of polar and apolar solutes in terms of differential distortions of the water-water hydrogen bonding angle in the first hydration shell. This method of hydration analysis is applied here to study the hydration of the wild type III thermal hysteresis protein from eel pout and three mutations at residue 16. Wild type and one mutant have full activity, the other two mutants have little or no anti-freeze (thermal hysteresis) activity. The analysis reveals significant differences in the hydration structure of the ice-binding site (centered on residue 16) among four proteins. For the A16T and A16Y mutants with reduced activity, polar groups have a typical polar-like hydration. For the wild type and mutant A16C with 100% of the wild type activity, polar groups have unusual, very apolar-like hydration. In the latter case, hydrating water molecules form a more ice-like pattern of hydrogen bonding on the ice-binding face, while in the former case water-water H-bonds are more distorted and more heterogenous. Overall, the binding surface of active protein strongly enhances the water tetrahedral structure, i.e. promotes ice-like hydration. It is concluded that the specific shape, residue size and clustering of both polar/apolar groups are essential for the binding surface to recognize, and preferentially interact with nascent ice crystals forming in liquid water.     

Increasing Protein Stability by Polar Surface Residues: Domain-Wide Consequences of Interactions Within a Loop

We have examined the influence of surface hydrogen bonds on the stability of proteins by studying the effects of mutations of human immunoglobulin light chain variable domain (V_L). In addition to the variants Y27dD, N28F, and T94H of protein κIV Len that were previously described, we characterized mutants M4L, L27cN, L27cQ, and K39T, double mutant M4L/Y27dD, and triple mutant M4L/Y27dD/T94H. The triple mutant had an enhanced thermodynamic stability of 4.2 kcal/mol. We determined the structure of the triple mutant by x-ray diffraction and correlated the changes in stability due to the mutations with changes in the three-dimensional structure. Y27dD mutant had increased stability of Len by 2.7 kcal/mol, a large value for a single mutation. Asp27d present in CDR1 formed hydrogen bonds with the side-chain and main-chain atoms within the loop. In the case of the K39T mutant, which reduces stability by 2 kcal/mol, Lys39 in addition to forming a hydrogen bond with a carbonyl oxygen of a neighboring loop may also favorably influence the surface electrostatics of the molecule. We showed that hydrogen bonds between residues in surface loops can add to the overall stability of the V_L domains. The contribution to stability is further increased if the surface residue makes more than one hydrogen bond or if it forms a hydrogen bond between neighboring turns or loops separated from each other in the amino acid sequence. Based on our experiments we suggest that stabilization of proteins might be systematically accomplished by introducing additional hydrogen bonds on the surface. These substitutions are more straightforward to predict than core-packing interactions and can be selected to avoid affecting the protein’s function.