The role of solvent polarity in the free energy of transfer of amino acid side chains from water to organic solvents

The transfer free energies of amino acid side chains from water to N-methylacetamide have been determined and compared with those obtained from other model systems. Although the process of transfer from water to N-methylacetamide represents transfer from a lower dielectric phase to a higher dielectric phase, the transfer free energies of most of the amino acid side chains are nearly the same as those obtained from the water to ethanol system. Among the apolar side chains studied, only the transfer free energies of methionine and the aromatic side chains are apparently influenced to some extent by the polarity of the organic solvent phase. The transfer free energies of the neutral polar side chains also exhibit significant dependence on solvent polarity. The van't Hoff plots for most of the apolar side chains exhibit nonlinear curves, indicating that the enthalpy of transfer from water to N-methylacetamide is temperature-dependent. It is suggested that to assess the contribution of the hydrophobic free energy to the stability of globular proteins, it is probably not necessary to account for variation in the internal environment of the protein.     

On the dependence of molecular conformation on the type of solvent environment: a molecular dynamics study of cyclosporin A

The dependence of the conformation of cyclosporin A (CPA), a cyclic undecapeptide with potent immunosuppressive activity, on the type of solvent environment is examined using the computer simulation method of molecular dynamics (MD). Conformational and dynamic properties of CPA in aqueous solution are obtained from MD simulations of a CPA molecule dissolved in a box with water molecules. Corresponding properties of CPA in apolar solution are obtained from MD simulations of CPA in a box with carbontetrachloride. The results of these simulations in H2O and in CCl4 are compared to each other and to those of previous simulations of crystalline CPA and of an isolated CPA molecule. The conformation of the backbone of the cyclic polypeptide is basically independent of the type of solvent. In aqueous solution the beta-pleated sheet is slightly weaker and the gamma-turn is a bit less pronounced than in apolar solution. Side chains may adopt different conformations in different solvents. In apolar solution the hydrophobic side chain of the MeBmt residue is in an extended conformation with its hydroxyl group hydrogen bonded to the backbone carbonyl group. In aqueous solution this hydrophobic side chain folds over the core of the molecule and the mentioned hydrogen bond is broken in favor of hydrogen bonding to water molecules. The conformation obtained from the MD simulation in CCl4 nicely agrees with experimental atom-atom distance data as obtained from nmr experiments in chloroform. In aqueous solution the relaxation of atomic motion tends to be slower than in apolar solution.     

A structural model of the acetylcholine receptor channel based on partition energy and helix packing calculations.

A structural model of the transmembrane portion of the acetylcholine receptor was developed from sequences of all its subunits by using transfer energy calculations to locate transmembrane alpha-helices and to calculate which helical side chains should be in contact with water inside the channel, with portions of other transmembrane helices, or with lipid hydrocarbon chains. "Knobs-into-holes" side chain packing calculations were used with other factors to stack the transmembrane alpha-helices together. In the model each subunit has the following structures in order along the sequence from the NH2 terminus: a large extracellular domain of undetermined structure, a short apolar alpha-helix that lies on the extracellular lipid surface of the membrane; three apolar transmembrane alpha-helices (I, II, and III), a cytoplasmic domain of undetermined structure, an amphipathic transmembrane alpha-helix (L) that forms the channel lining, a short extracellular alpha-helix, another apolar transmembrane alpha-helix (IV), and a small cytoplasmic domain formed by the COOH-terminal end of the chain. Three concentric layers form the pore. A bundle of five amphipathic L helices forms the channel lining. This bundle is surrounded by a bundle of 10 alternating II and III helices. Helices I and IV cover portions of the outer surface of the bundle formed by helices II and III. Positions of disulfide bridges are predicted and a mechanism for opening and closing conformational changes is proposed that requires tilting transmembrane helices and possibly a thiol-disulfide interchange reaction.     

Packed protein bilayers in the 0.90 A resolution structure of a designed alpha helical bundle.

A 12-residue peptide designed to form an alpha-helix and self-associate into an antiparallel 4-alpha-helical bundle yields a 0.9 A crystal structure revealing unanticipated features. The structure was determined by direct phasing with the "Shake-and-Bake" program, and contains four crystallographically distinct 12-mer peptide molecules plus solvent for a total of 479 atoms. The crystal is formed from nearly ideal alpha-helices hydrogen bonded head-to-tail into columns, which in turn pack side-by-side into sheets spanning the width of the crystal. Within each sheet, the alpha-helices run antiparallel and are closely spaced (9-10 A center-to-center). The sheets are more loosely packed against each other (13-14 A between helix centers). Each sheet is amphiphilic: apolar leucine side chains project from one face, charged lysine and glutamate side chains from the other face. The sheets are stacked with two polar faces opposing and two apolar faces opposing. The result is a periodic biomaterial composed of packed protein bilayers, with alternating polar and apolar interfaces. All of the 30 water molecules in the unit cell lie in the polar interface or between the stacked termini of helices. A section through the sheet reveals that the helices packed at the apolar interface resemble the four-alpha-helical bundle of the design, but the helices overhang parts of the adjacent bundles, and the helix crossing angles are less steep than intended (7-11 degrees rather than 18 degrees).     

The refinement and the structure of the dimer of alpha-chymotrypsin at 1.67-A resolution

The two molecules of the asymmetric unit of the pH 3.5 conformer of alpha-chymotrypsin have been refined at 1.67-A resolution using restrained least squares methods with Hendrickson's program (PROLSQ). The final R factor is 0.179 (including 247 water molecules). The folding of the main chain of the independent molecules is the same within experimental error but the same does not generally apply to the side chain stereochemistry. From this we conclude that the folding of a protein structure is basically independent of most of the detailed stereochemistry of its side chains. The side chains of the interface region between the independent molecules display pronounced asymmetry. This asymmetry suggests that dynamic and asymmetrical structural changes take place at the time of oligomerization leading to more energetically favorable interactions for the dimer. Comparison of the structures of the independent molecules of alpha-chymotrypsin with the structure of monomeric gamma-chymotrypsin revealed that although the folding of the three molecules is essentially the same, numerous and significant differences pervade the side chain stereochemistry attributable to general flexibility. The specificity site of alpha-chymotrypsin is occupied by ordered water molecules in a similar way to gamma-chymotrypsin and other proteins. Some of these water molecules are displaced when substrate binds to the enzyme, while the others appear to help identify and position the aromatic side chain in catalysis.     

Computer simulation of hydration networks around amino acids

The networks of solvent hydrogen bonds around polar and apolar amino acids have been studied by computer simulation techniques using a non-pair additive model for the water molecules interactions. Analysis of the simulated aqueous solutions has shown the presence of water molecules which (a) form a bridge around individual polar solute atoms (self-bridging loops) and (b) form chains between different polar solute atoms (polar bridging chains). Some of these networks associated with polar solute atoms from pentagons but 4, 6 and 7 sided polygons are also seen. The water molecule close to apolar solute atoms (<4.0 Å) also form irregular networks with polygons of 4, 5, 6 and 7 sides. These networks are compared with those found experimentally in ice, clathrates and crystal hydrates of macromolecules.     

Stabilization of internal charges in a protein: water penetration or conformational change?

The ionizable amino acid side chains of proteins are usually located at the surface. However, in some proteins an ionizable group is embedded in an apolar internal region. Such buried ionizable groups destabilize the protein and may trigger conformational changes in response to pH variations. Because of the prohibitive energetic cost of transferring a charged group from water to an apolar medium, other stabilizing factors must be invoked, such as ionization-induced water penetration or structural changes. To examine the role of water penetration, we have measured the 17O and 2H magnetic relaxation dispersions (MRD) for the V66E and V66K mutants of staphylococcal nuclease, where glutamic acid and lysine residues are buried in predominantly apolar environments. At neutral pH, where these residues are uncharged, we find no evidence of buried water molecules near the mutation site. This contrasts with a previous cryogenic crystal structure of the V66E mutant, but is consistent with the room-temperature crystal structure reported here. MRD measurements at different pH values show that ionization of Glu-66 or Lys-66 is not accompanied by penetration of long-lived water molecules. On the other hand, the MRD data are consistent with a local conformational change in response to ionization of the internal residues.     

Statistical and conformational analysis of the electron density of protein side chains

Protein side chains make most of the specific contacts between proteins and other molecules, and their conformational properties have been studied for many years. These properties have been analyzed primarily in the form of rotamer libraries, which cluster the observed conformations into groups and provide frequencies and average dihedral angles for these groups. In recent years, these libraries have improved with higher resolution structures and using various criteria such as high thermal factors to eliminate side chains that may be misplaced within the crystallographic model coordinates. Many of these side chains have highly non-rotameric dihedral angles. The origin of side chains with high B-factors and/or with non-rotameric dihedral angles is of interest in the determination of protein structures and in assessing the prediction of side chain conformations. In this paper, using a statistical analysis of the electron density of a large set of proteins, it is shown that: (1) most non-rotameric side chains have low electron density compared to rotameric side chains; (2) up to 15% of chi1 non-rotameric side chains in PDB models can clearly be fit to density at a single rotameric conformation and in some cases multiple rotameric conformations; (3) a further 47% of non-rotameric side chains have highly dispersed electron density, indicating potentially interconverting rotameric conformations; (4) the entropy of these side chains is close to that of side chains annotated as having more than one chi(1) rotamer in the crystallographic model; (5) many rotameric side chains with high entropy clearly show multiple conformations that are not annotated in the crystallographic model. These results indicate that modeling of side chains alternating between rotamers in the electron density is important and needs further improvement, both in structure determination and in structure prediction.     

Crystallographic study of hydration of an internal cavity in engineered proteins with buried polar or ionizable groups

Although internal water molecules are essential for the structure and function of many proteins, the structural and physical factors that govern internal hydration are poorly understood. We have examined the molecular determinants of internal hydration systematically, by solving the crystal structures of variants of staphylococcal nuclease with Gln-66, Asn-66, and Tyr-66 at cryo (100 K) and room (298 K) temperatures, and comparing them with existing cryo and room temperature structures of variants with Glu-66, Asp-66, Lys-66, Glu-92 or Lys-92 obtained under conditions of pH where the internal ionizable groups are in the neutral state. At cryogenic temperatures the polar moieties of all these internal side chains are hydrated except in the cases of Lys-66 and Lys-92. At room temperature the internal water molecules were observed only in variants with Glu-66 and Tyr-66; water molecules in the other variants are probably present but they are disordered and therefore undetectable crystallographically. Each internal water molecule establishes between 3 and 5 hydrogen bonds with the protein or with other internal water molecules. The strength of interactions between internal polar side chains and water molecules seems to decrease from carboxylic acids to amides to amines. Low temperature, low cavity volume, and the presence of oxygen atoms in the cavity increase the positional stability of internal water molecules. This set of structures and the physical insight they contribute into internal hydration will be useful for the development and benchmarking of computational methods for artificial hydration of pockets, cavities, and active sites in proteins.     

Conformational characteristics of peptides and unanticipated results from crystal structure analyses

Preferred conformation and types of molecular folding are some of the topics that can be addressed by structure analysis using x-ray diffraction of single crystals. The conformations of small linear peptide molecules with 2-6 residues are affected by polarity of solvent, presence of water molecules, hydrogen bonding with neighboring molecules, and other packing forces. Larger peptides, both cyclic and linear, have many intramolecular hydrogen bonds, the effect of which outweighs any intermolecular attractions. Numerous polymorphs of decapeptides grown from a variety of solvents, with different cocrystallized solvents, show a constant conformation for each peptide. Large conformational changes occur, however, upon complexation with metal ions. A new form of free valinomycin grown from DMSO exhibits near three-fold symmetry with only three intramolecular hydrogen bonds. The peptide is in the form of a shallow bowl with a hydrophobic exterior. Near the bottom of the interior of the bowl are three carbonyl oxygens, spaced and directed so that they are in position to form three ligands to a K+, e.g., complexation can be completed by the three lobes containing the beta-bends closing over and encapsulating the K+ ion. In another example, free antamanide and the biologically inactive perhydro analogue, in which four phenyl groups become cyclic hexyl groups, have essentially the same folding of backbone and side chains. The conformation changes drastically upon complexation with Li+ or Na+. However, the metal ion complex of natural antamanide has a hydrophobic globlar form whereas the metal ion complex of the inactive perhydro analogue has a polar band around the middle. The structure results indicate that the antamanide molecule is in a complexed form during its biological activity. Single crystal x-ray diffraction structure analyses have identified the manner in which water molecules are essential to creating minipolar areas on apolar helices. Completely apolar peptides, such as membrane-active peptides, can acquire amphiphilic character by insertion of a water molecule into the helical backbone of Boc-Aib-Ala-Leu-Aib-Ala-Leu-Aib-Ala-Leu-Aib-OMe, for example. The C-terminal half assumes an alpha-helix conformation, whereas the N-terminal half is distorted by an insertion of a water molecule W(1) between N(Ala5) and O(Ala2), forming hydrogen bonds N(5)H...W(1) and W(1)...O(2). The distortion of the helix exposes C = O(Aib1) and C = O(Aib4) to the outside environment with the consequence of attracting additional water molecules. The leucyl side chains are on the other side of the molecule. Thus a helix with an apolar sequence can mimic an amphiphilic helix.     

Pseudo-dynamic contact surface areas: estimation of apolar bonding

A new Monte Carlo based algorithm has been written for the computation of pseudo-dynamic contact surface areas. The linear correlation of this contact area with solute transfer free energies (water leads to organic liquid) is established for apolar amino acid side chains. The slope of these linear plots, deltaGosp, is a unitary free energy which has potential use in the estimation of apolar bond free energies in proteins. The magnitude of deltaGosp is dependent upon the nature of the organic solvent involved in the transfer process, varying from 86 to 130 cal/A2. Analogues linear correlations with the same range of deltaGosp values are observed for inhibitors of protein-catalyzed reactions.     

Structural Characteristics in Protein Hydration Investigated by Cryogenic X-ray Crystal Structure Analyses

Cryogenic X-ray crystallography has heen applied to investigate thehydration structures of proteins. The amount of hydration water moleculesidentified at cryogenic temperature is more than twice those at ambienttemperature, and the structural models of proteins with a lot of hydrationwater molecules have provided much information to elucidate the static anddynamical characteristics of hydration structures of proteins. On proteinsurface, hydration water molecules distribute non-randomly and stillretain the tetrahedral hydrogen-bond geometry as well as in bulk solvent.In addition, water molecules form clathrate-like arrangements to cover thehydrophobic residues exposed to solvent. The standard interaction geometryenables the three-dimensional extension of hydrogen-bond networks aroundprotein molecules and, simultaneously, ensures the concerted reorganizationof hydration structures during the dynamical motion of proteins at work.The hydration structure analyses at cryogenic temperatures may contributeto understanding physical principles governing the dynamics of `molecularmachines' in aqueous environment.     

Position-dependence of stabilizing polar interactions of asparagine in transmembrane helical bundles

Recent studies with model peptides and statistical analyses of the crystal structures of membrane proteins have shown that buried polar interactions contribute significantly to the stabilization of the three-dimensional structures of membrane proteins. Here, we probe how the location of these polar groups along the transmembrane helices affect their free energies of interaction. Asn residues were placed singly and in pairs at three positions within a model transmembrane helix, which had previously been shown to support the formation of trimers in micelles. The model helix was designed to form a transmembrane coiled coil, with Val side chains at the "a" positions of the heptad repeat. Variants of this peptide were prepared in which an Asn residue was introduced at one or more of the "a" positions, and their free energies of association were determined by analytical ultracentrifugation. When placed near the middle of the transmembrane helix, the formation of trimers was stabilized by at least -2.0 kcal/mol per Asn side chain. When the Asn was placed at the interface between the hydrophobic and polar regions of the peptide, the substitution was neither stabilizing nor destabilizing (0.0 +/- 0.5 kcal/mol of monomer). Finally, it has previously been shown that a Val-for-Asn mutation in a water-soluble coiled coil destabilizes the structure by approximately 1.5 kcal/mol of monomer [Acharya, A., et al. (2002) Biochemistry 41, 14122-14131]. Thus, the headgroup region of a micelle appears to have a conformational impact intermediate between that of bulk water and the apolar region of micelle. A similarly large dependence on the location of the polar residues was found in a statistical survey of helical transmembrane proteins. The tendency of different types of residues to be buried in the interiors versus being exposed to lipids was analyzed. Asn and Gln show a very strong tendency to be buried when they are located near the middle of a transmembrane helix. However, when placed near the ends of transmembrane helices, they show little preference for the surface versus the interior of the protein. These data show that Asn side chains within the apolar region of the transmembrane helix provide a significantly larger driving force for association than Asn residues near the apolar/polar interface. Thus, although polar interactions are able to strongly stabilize the folding of membrane proteins, the energetics of association depend on their location within the hydrophobic region of a transmembrane helix.     

Artificial antifreeze polypeptides: alpha-helical peptides with KAAK motifs have antifreeze and ice crystal morphology modifying properties.

Antifreeze polypeptides from fish are generally thought to inhibit ice crystal growth by specific adsorption onto ice surfaces and preventing addition of water molecules to the ice lattice. Recent studies have suggested that this adsorption results from hydrogen bonding through the side chains of polar amino acids as well as hydrophobic interactions between the non-polar domains on the ice-binding side of antifreeze polypeptides and the clathrate-like surfaces of ice. In order to better understand the activity of one of the antifreeze polypeptide families, namely the alpha-helical type I antifreeze polypeptides, four alpha-helical peptides having sequences not directly analogous to those of known antifreeze polypeptides and containing only positively charged and non-polar side chains were synthesized. Two peptides with regularly spaced lysine residues, GAAKAAKAAAAAAAKAAKAAAAAAAKAAKAAGGY-NH2 and GAALKAAKAAAAAALKAAKAAAAAALKAAKAAGGY-NH2, showed antifreeze activity, albeit weaker than seen in natural antifreeze polypeptides, by the criteria of freezing point depression (thermal hysteresis) and ice crystal modification to a hexagonal trapezohedron. Peptides with irregular spacing of Lys residues were completely inactive. Up to now, lysine residues have not been generally associated with antifreeze activity, though they have been implicated in some antifreeze polypeptides. This work also shows that lysine residues in themselves, when properly positioned on an alpha-helical polyalanine scaffold, have all the requisite properties needed for such an activity.     

Waterproofing the heme pocket. Role of proximal amino acid side chains in preventing hemin loss from myoglobin

The ability of myoglobin to bind oxygen reversibly depends critically on retention of the heme prosthetic group. Globin side chains at the Leu(89)(F4), His(97)(FG3), Ile(99)(FG5), and Leu(104)(G5) positions on the proximal side of the heme pocket strongly influence heme affinity. The roles of these amino acids in preventing heme loss have been examined by determining high resolution structures of 14 different mutants at these positions using x-ray crystallography. Leu(89) and His(97) are important surface amino acids that interact either sterically or electrostatically with the edges of the porphyrin ring. Ile(99) and Leu(104) are located in the interior region of the proximal pocket beneath ring C of the heme prosthetic group. The apolar amino acids Leu(89), Ile(99), and Leu(104) "waterproof" the heme pocket by forming a barrier to solvent penetration, minimizing the size of the proximal cavity, and maintaining a hydrophobic environment. Substitutions with smaller or polar side chains at these positions result in exposure of the heme to solvent, the appearance of crystallographically defined water molecules in or near the proximal pocket, and large increases in the rate of hemin loss. Thus, the naturally occurring amino acid side chains at these positions serve to prevent hydration of the His(93)-Fe(III) bond and are highly conserved in all known myoglobins and hemoglobins.     

Hydration of oligosaccharides: anomalous hydration ability of trehalose.

The disaccharide trehalose extensively exists in anhydrobiotic organism and is considered to play an important role in preserving the integrity of biomembrane. However, the preserving mechanism remains unclear. In this report, we examine the hydration abilities of trehalose and several oligosaccharides composed of alpha-D-glucopyranosyl residues. The unfrozen water fraction per molecule was determined from differential scanning calorimetry measurements of their aqueous solutions. Further, the NMR relaxation time of the natural abundance 17O of water is measured for several saccharide solutions. These results indicate that trehalose has the highest hydration ability among the saccharides studied. In other words, trehalose can effectively lower the mobility of water molecules hydrogen-bonded with saccharides. It is thus reasonable that, among the disaccharides studied, trehalose exhibits the maximum stabilizing effect on the bilayer structure of lipid whose acyl chains are bonded with each other by the apolar interaction, because the apolar interaction is strengthened with the stabilization of the surrounding water structure.     

Refined atomic model of glutamine synthetase at 3.5 A resolution

An atomic model of 43,692 non-hydrogen atoms has been determined for the 12-subunit enzyme glutamine synthetase from Salmonella typhimurium, by methods of x-ray diffraction including restrained least-squares atomic refinement against 65,223 unique reflections. At 3.5 A resolution the crystallographic R-factor (on 2 sigma data) is 25.8%. As reported earlier for the unrefined structure, the 12 subunits are arranged in two layers of six; at the interface of pairs of subunits within each layer, cylindrical active sites are formed by six anti-parallel beta strands contributed by one subunit and two strands by the neighboring subunit. This interpretation of the electron density map has now been supported by comparison with glutamine synthetase from Escherichia coli by the Fourier difference method. Each active site cylinder holds two Mn2+ ions, with each ion having as ligands three protein side chains and two water molecules (one water shared by both metals), as well as a histidyl side chain just beyond liganding distance. The protein ligands to Mn2+ 469 are Glu-131, Glu-212, and Glu-220; those to Mn2+ 470 are Glu-129, His-269, and Glu-357. The two layers of subunits are held together largely by the apolar COOH terminus, a helical thong, which inserts into a hydrophobic pocket formed by two neighboring subunits on the opposite ring. Also between layers, there is a hydrogen-bonded beta sheet interaction, as there is between subunits within a ring, but hydrophobic interactions account for most of the intersubunit stability. The central loop, which extends into the central aqueous channel, is subject to attack by at least five enzymes and is discussed as an enzyme "passive site."