Effects of hydrated water on protein unfolding

The conformational stability of a protein in aqueous solution is described in terms of the thermodynamic properties such as unfolding Gibbs free energy, which is the difference in the free energy (Gibbs function) between the native and random conformations in solution. The properties are composed of two contributions, one from enthalpy due to intramolecular interactions among constituent atoms and chain entropy of the backbone and side chains, and the other from the hydrated water around a protein molecule. The hydration free energy and enthalpy at a given temperature for a protein of known three-dimensional structure can be calculated from the accessible surface areas of constituent atoms according to a method developed recently. Since the hydration free energy and enthalpy for random conformations are computed from those for an extended conformation, the thermodynamic properties of unfolding are evaluated quantitatively. The evaluated hydration properties for proteins of known transition temperature (Tm) and unfolding enthalpy (delta Hm) show an approximately linear dependence on the number of constituent heavy atoms. Since the unfolding free energy is zero at Tm, the enthalpy originating from interatomic interactions of a polypeptide chain and the chain entropy are evaluated from an experimental value of delta Hm and computed properties due to the hydrated water around the molecule at Tm. The chain enthalpy and entropy thus estimated are largely compensated by the hydration enthalpy and entropy, respectively, making the unfolding free energy and enthalpy relatively small. The computed temperature dependences of the unfolding free energy and enthalpy for RNase A, T4 lysozyme, and myoglobin showed a good agreement with the experimental ones.(ABSTRACT TRUNCATED AT 250 WORDS)     

Distributions of water around amino acid residues in proteins

The atomic co-ordinates from 16 high-resolution (less than or equal to 1.7 A = 0.1 nm), non-homologous proteins have been used to study the distributions of water molecule sites around the 20 different amino acid residues. The proportion of residues whose main-chain atoms are in contact with water molecules was fairly constant (between 40% and 60%), irrespective of the nature of the side-chain. However, the proportion of residues whose side-chain atoms were in contact with water molecules showed a clear (inverse) correlation with the hydrophobicity of the residue, being as low as 14% for leucine and isoleucine but greater than 80% for asparagine and arginine. Despite the problems in determining accurate water molecule sites from X-ray diffraction data and the complexity of the protein surface, distinct non-random distributions of water molecules were found. These hydration patterns are consistent with the expected stereochemistry of the potential hydrogen-bonding sites on the polar side-chains. The water molecules around apolar side-chains lie predominantly at van der Waals' contact distances, but most of these have a primary, shorter contact with a neighbouring polar atom. Further analysis of these distributions, combined with energy minimization techniques, should lead to improved modelling of protein structures, including their primary shells of hydration.     

H-bonding in protein hydration revisited

H-bonding between protein surface polar/charged groups and water is one of the key factors of protein hydration. Here, we introduce an Accessible Surface Area (ASA) model for computationally efficient estimation of a free energy of water-protein H-bonding at any given protein conformation. The free energy of water-protein H-bonds is estimated using empirical formulas describing probabilities of hydrogen bond formation that were derived from molecular dynamics simulations of water molecules at the surface of a small protein, Crambin, from the Abyssinian cabbage (Crambe abyssinica) seed. The results suggest that atomic solvation parameters (ASP) widely used in continuum hydration models might be dependent on ASA for polar/charged atoms under consideration. The predictions of the model are found to be in qualitative agreement with the available experimental data on model compounds. This model combines the computational speed of ASA potential, with the high resolution of more sophisticated solvation methods.     

Significant role of electrostatic interactions for stabilization of protein assemblies

Contribution of electrostatic interactions to stability of BPTI orthorhombic, pig-insulin cubic crystals, and horse L ferritin crystals was evaluated with numerical calculation of Poisson-Boltzmann equation based on a dielectric model. The stability of a ferritin molecule (24-mer) composed of 24 subunits was also evaluated. It was found that the surface charge-charge interactions at separation distances (< 5 Å) were insensitive to variations in the ionic strength, and thus stabilized assembled states of the proteins (i.e., crystalline state and oligomeric state). It was also revealed that the charge density and the packing of the protein crystals were largely responsible for the ionic strength dependence of the crystal stability. The stability of the 5PTI crystalline state with a high charge density drastically increased as the concentration of the solvent ions increased. In contrast, that of the insulin crystal with a low charge density and large solvent region was insensitive to changes in the ionic concentration. The electrostatic interaction between ferritin 24-mers was attributed to two salt bridges mediated by Cd ion. For the stability of the ferritin 24-mer, which is evolutionally designed, the electrostatic stabilization between the subunits was attributed to polar bonds such as buried salt bridges or hydrogen bonds, which occasionally yielded more than 5 kcal/mol and were numerous and very strong compared with the bonds between molecules in the 5PTI and 9INS crystals.By analyzing the atomic charge-charge interactions in detail, it was found that charge pairs separated by less than 3 Å, such as hydrogen bonds, dominantly stabilize the assembled states, and that pairs 3 to 5 Å apart were also important. The stability of the assembled states evaluated by the total EET was determined by the fine balance between the two competing contributions arising from the stabilizing atoms and the destabilizing atoms.Changes of the ASA and hydration free energy were also evaluated in accordance with the process of the subunit assembly. The change of hydration free energy, which was very large (i.e., ~+ 100 kcal/mol/subunit) and unfavorable for the assembly, was proportional to the electrostatic hydration energy (i.e., Born energy change in the hydration process). Hydrophobic groups were likely to appear more frequently than hydrophilic groups at the interfaces.This study offers a method which can improve the stability of protein crystals by introducing polar or charged residues that are properly designed to form specific hydrogen bonds or salt bridges between neighboring protein molecules. This method is also applicable to crystallography, because it improves refinement of protein structures in crystals by taking the inter-protein interactions into account.     

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.     

Excluded volume approximation to protein-solvent interaction. The solvent contact model.

Important properties of globular proteins, such as the stability of its folded state, depend sensitively on interactions with solvent molecules. Existing methods for estimating these interactions, such as the geometrical surface model, are either physically misleading or too time consuming to be applied routinely in energy calculations. As an alternative, we derive here a simple model for the interactions between protein atoms and solvent atoms in the first hydration layer, the solvent contact model, based on the conservation of the total number of atomic contacts, a consequence of the excluded-volume effect. The model has the conceptual advantage that protein-protein contacts and protein-solvent contacts are treated in the same language and the technical advantage that the solvent term becomes a particularly simple function of interatomic distances. The model allows rapid calculation of any physical property that depends only on the number and type of protein-solvent nearest-neighbor contacts. We propose use of the method in the calculation of protein solvation energies, conformational energy calculations, and molecular dynamics simulations.     

Protein-water interactions in ribonuclease A and angiogenin: a molecular dynamics study

It is known that water molecules play an important role in the biological functioning of proteins. The members of the ribonuclease A (RNase A) family of proteins, which are sequentially and structurally similar, are known to carry out the obligatory function of cleaving RNA and individually perform other diverse biological functions. Our focus is on elucidating whether the sequence and structural similarity lead to common hydration patterns, what the common hydration sites are and what the differences are. Extensive molecular dynamics simulations followed by a detailed analysis of protein-water interactions have been carried out on two members of the ribonuclease A superfamily-RNase A and angiogenin. The water residence times are analyzed and their relationship with the characteristic properties of the protein polar atoms, such as their accessible surface area and mean hydration, is studied. The capacity of the polar atoms to form hydrogen bonds with water molecules and participate in protein-water networks are investigated. The locations of such networks are identified for both proteins.     

Electrostatic interactions during electron transfer reactions between c-type cytochromes and flavodoxin

The interaction of three different c-type cytochromes with flavodoxin has been studied by computer graphics modelling and computational methods. Flavodoxin and each cytochrome can make similar hypothetical electron transfer complexes that are characterized by nearly coplanar arrangement of the prosthetic groups, close intermolecular contacts at the protein-protein interface, and complementary intermolecular salt linkages. Computation of the electrostatic free energy of each complex showed that all were electrostatically stable. However, both the magnitude and behavior of the electrostatic stabilization as a function of solution ionic strength differed for the three cytochrome c-flavodoxin complexes. Variation in the computed electrostatic stabilization appears to reflect differences in the surface distribution of all charged groups in the complex, rather than differences localized at the site of intermolecular contact. The computed electrostatic association constants for the complexes and the measured kinetic rates of electron transfer in solution show a remarkable similarity in their ionic strength dependence. This correlation suggests electrostatic interactions influence electron transfer rates between protein molecules at the intermolecular association step. Comparative calculations for the three cytochrome c-flavodoxin complexes show that these ionic strength effects also involve all charged groups in both redox partners.     

Hydrogen bonding of adenine derivatives to tyrosine side chain.

High resolution proton magnetic resonance measurements provide evidence for the formation of hydrogen-bonded complexes between 9-ethyladenine and p-cresol used as a model of tyrosine side chain in CDCl3. We have calculated the sum of the association constants corresponding to the three existing 1:1 complexes: K=6.3+/-0.15. By methylation of the amino group of adenine, we were able to calculate the ratio of the two strongest equilibrium constants K7/K1=1.6+/-0.3. Theoretical computations by the complete neglect of differential overlap (CNDO/2) method indicate that several hydrogen-bonded planar complexes can form between 9-methyladenine and phenol. The computed energy of the complexes with 6-dimethylamino adenine removes some ambiguity concerning the computed ratio of the association constants. Comparison of the calculated energies with free energies experimentally determined in organic solvent shows that despite the competition with CDCl3, which associates with both solute molecules, the preferential order of association is conserved. The small variations of charge density of adenine carbon atoms when complexed with phenol are in agreement with very small chemical shifts observed by 13C-nuclear magnetic resonance.     

Stochastic modeling of noninteracting probes in the protein structure space for construction of knowledge-based potentials for atom-atom interactions

Knowledge-based potentials are extensively used to represent atomic interactions in modeling the protein structure. We consider a number of problems in constructing efficient knowledge-based potentials for biopolymer modeling. We show that some limitations can be overcome by normalizing estimated interactions through the distribution of distances between noninteracting random probes in protein structure space. We demonstrate that knowledge-based potentials thus constructed can be efficiently applied for analysis of the hydration state of proteins atoms. With this approach, one can predict the locations of structural water molecules in a protein globule. We have also succeeded in recognizing the correctly folded protein structure among many misfolded decoys in cases when the interaction with water solvent is dominant for structure formation.     

Frictional models for stochastic simulations of proteins

As a first step toward a systematic parametrization of friction constants of atoms in proteins, a model in which frictional resistance is placed explicitly on each atom accessible to solvent is used to calculate overall translational and rotational diffusion constants. It is found that these quantities are relatively insensitive to the precise value of the atomic friction constant, as long as the effective hydrodynamic radius of the surface atoms is approximately 1 Å. However, if only protein atoms are included in the calculation, no reasonable range atomic of radii can reproduce the experimental translational diffusion constant to better than 20% for lysozyme and 5% for ribonuclease. When a hydration shell of approximately 70% coverage for lysozyme and 20% for ribonuclease is included, there is quantitative agreement with experimental results. The sensitivity of peptide diffusion to levels of hydration is also investigated; it is found that for glycine, two bound waters are required to provide agreement with experiment. These findings imply that the effects of solvent damping will be underestimated in stochastic simulations of proteins and peptides unless bound waters are taken into account.     

A Connected-cluster of hydration around myoglobin: Correlation between molecular dynamics simulations and experiment

An analysis of a molecular dynamics simulation of metmyoglobin in an explicit solvent environment of 3,128 water molecules has been performed. Both statics and dynamics of the protein-solvent interface are addressed in a comparison with experiment. Three-dimensional density distributions, temperature factors, and occupancy weights are computed for the solvent by using the trajectory coordinates. Analysis of the hydration leads to the localization of more than 500 hydration sites distributed into multiple layers of solvation located between 2.6 and 6.8 Å from the atomic protein surface. After locating the local solvent density maxima or hydration sites we conclude that water molecules of hydration positions and hydration sites are distinct concepts. Both global and detailed properties of the hydration cluster around myoglobin are compared with recent neutron and X-ray data on myoglobin. Questions arising from differences between X-ray and neutron data concerning the locations of the protein-bound water are investigated. Analysis of water site differences found from X-ray and neutron experiments compared with our simulation shows that the simulation gives a way to unify the hydration picture given by the two experiments. © 1994 John Wiley & Sons, Inc.     

Empirical solvation models can be used to differentiate native from near-native conformations of bovine pancreatic trypsin inhibitor.

Several hydration models for peptides and proteins based on solvent accessible surface area have been proposed previously. We have evaluated some of these models as well as four new ones in the context of near-native conformations of a protein. In addition, we propose an empirical site-site distance-dependent correction that can be used in conjunction with any of these models. The set of near-native structures consisted of 39 conformations of bovine pancreatic trypsin inhibitor (BPTI) each of which was a local minimum of an empirical energy function (ECEPP) in the absence of solvent. Root-mean-square (rms) deviations from the crystallographically determined structure were in the following ranges: 1.06-1.94 A for all heavy atoms, 0.77-1.36 A for all backbone heavy atoms, 0.68-1.33 A for all alpha-carbon atoms, and 1.41-2.72 A for all side-chain heavy atoms. We have found that there is considerable variation among the solvent models when evaluated in terms of concordance between the solvation free energy and the rms deviations from the crystallographically determined conformation. The solvation model for which the best concordance (0.939) with the rms deviations of the C alpha atoms was found was derived from NMR coupling constants of peptides in water combined with an exponential site-site distance dependence of the potential of mean force. Our results indicate that solvation free energy parameters derived from nonpeptide free energies of hydration may not be transferrable to peptides. Parameters derived from peptide and protein data may be more applicable to conformational analysis of proteins. A general approach to derive parameters for free energy of hydration from ensemble-averaged properties of peptides in solution is described.     

The kinetics of interactions of bilirubin with lipid bilayers and with serum albumin.

Rate constants for the hydration of bilirubin bound to unilamellar bilayers of dioleoylphosphatidylcholine and albumin were measured by stopped-flow methods. Rate constants for association of bilirubin with these vesicles and albumin were calculated from measured rate constants for dissociation and the equilibrium binding constants of bilirubin and lipids or albumin. Rate constants for hydration (dissociation) for bilirubin bound to dioleoylphosphatidylcholine and albumin were 71 s-1 and 1.8 s-1 respectively. Rate constants for association were 4.0 10(7) s-1 and 1.1 10(9) M-1 s-1, respectively. Both rates for interactions of bilirubin with bilayers were essentially independent of temperature in the range 0-40 degrees C, indicating that barriers to entry and exit of bilirubin from bilayers were entropic. Rates of transbilayer movement of bilirubin in dioleoylphosphatidylcholine were too fast to resolve by measuring rates of hydration of bilirubin. Rate constants for hydration of bilirubin bound to bilayers with less avidity for bilirubin as compared with dioleoylphosphatidylcholine also were too fast to measure with stopped-flow methods. In addition to providing details of the energetic basis for interactions between bilirubin and membranes, the data allow for calculating the maximal rates at which bilirubin could transfer spontaneously from sites on albumin in blood to the interior of cells. The data show, in this regard, that this rate is 10-50 fold faster than measured rates of uptake of bilirubin by intact liver.     

Analysis of protein main-chain solvation as a function of secondary structure

We have analysed the hydration of main-chain carbonyl and amide groups in 24 high-resolution well-refined protein structures as a function of the secondary structure in which these polar groups occur. We find that main-chain atoms in beta-sheets are as hydrated as those in alpha-helices, with most interactions involving "free" amide and carbonyl groups that do not participate in secondary structure hydrogen bonds. The distributions of water molecules around these non-bonded carbonyl groups reflect specific steric interactions due to the local secondary structure. Approximately 20% and 4%, respectively of bonded carbonyl and amide groups interact with solvent. These include interactions with carbonyl groups on the exposed faces of alpha-helices that have been correlated previously with bending of the helix. Water molecules interacting with alpha-helices occur mainly at the amino and carbonyl termini of the helices, in which case the solvent sites maintain the hydrogen bonding by bridging between residues i and i-3 or i-4 at the amino terminus and between i and i+3 or i+4 at the carbonyl terminus. We also see a number of solvent-mediated Ncap and Ccap interactions. The water molecules interacting with beta-sheets occur mainly at the edges, in which case they extend the sheet structure, or at the ends of strands, in which case they extend the beta-ladder. In summary, the solvent networks appear to extend the hydrogen-bonding structure of the secondary structures. In beta-turns, which usually occur at the surface of a protein, exposed amide and carbonyl groups are often hydrated, especially close to glycine residues. Occasionally water molecules form a bridge between residues i and i+3 in the turn and this may provide extra stabilization.     

Fatty acid binding to plasma albumin.

A review of the available information about fatty acid binding to plasma albumin is presented. Albumin is composed of a single polypeptide chain, folded so as to form three or four spherical units. The strong fatty acid binding sites probably are located in crevices between these spherical regions. The anionic form of the fatty acid binds to albumin. Most of the binding energy comes from nonpolar interactions between the fatty acid hydrocarbon chain and uncharged amino acid side chains that line the binding sites. The binding sites are somewhat pliable, and their configuration can adapt to fit the incoming fatty acid. Stepwise association constants for binding to human albumin of fatty acids containing 6-18 carbon atoms are presented. These data indicate that each mole of fatty acid binds with a different affinity and that the association constants for multiple binding diminish sequentially, i.e., kappa 1 greater than kappa 2 greater than kappa 3 greater ... greater kappan. Because of uncertainties concerning fatty acid association in aqueous solutions, the constants for the 14-18 carbon acids probably are not definitive. In the usual physiological concentration range, free fatty acids do not displace appreciable amounts of a second organic compound from albumin. Sensitive spectrophotometric analyses revealed, however, that even small increases in free fatty acid concentration alter the molecular interaction between human albumin and another organic compound.     

Subunit interactions and the allosteric response in phosphorylase.

The contribution of intersubunit interactions to allosterically induced conformational changes in phosphorylase are considered. Phosphorylase a, Pa (phosphorylated at Ser-14), is significantly in the active (R) conformation, while phosphorylase b, Pb (nonphosphorylated), is predominantly in the inactive (T) conformation. The structure of glucose-inhibited (T) Pa has been determined at 2.5-A resolution and atomic coordinates have been measured. These data have been used to calculate the solvent accessible surface area at the subunit interface and map noncovalent interactions between protomers. The subunit contact involves only 6% of the Pa monomer surface, but withdraws an area of 4,600 A2 from solvent. The contact region is confined to the N-terminal (regulatory) domain of the subunit. Half of the residues involved are among the 70 N-terminal peptides. A total of approximately 100 atoms take part in polar or nonpolar contacts of less than 4.0 A with atoms of the symmetry-related monomer. The contact surface surrounds a central cavity at the core of the interface of sufficient volume to accommodate 150-180 solvent molecules. There are four intersubunit salt bridges. Two of these (Arg 10/Asp 32, Ser-14-P/Arg 43) are interactions between the N-terminus of one protomer with an alpha-helix loop segment near the N-terminus of the symmetry-related molecule. These two are relatively solvent accessible. The remainder (Arg 49/Glu 195, Arg 184/Asp 251) are nearer the interface core and are less accessible. The salt bridges at the N-terminus are surrounded by the polar and nonpolar contacts which may contribute to their stability. Analysis of the difference electron density between the isomorphous Pa and Pb crystal structures reveals that the N-terminal 17 residues of Pb are disordered. Pb thus lacks two intermolecular and one intersubunit (Ser-14-P/Arg 69) salt linkage present in Pa. The absence of these interactions in Pb is manifested in the difference in the free energy of T leads to R activation, which is 4 kcal more than that for Pa. Difference Fourier analysis of the T leads to R transition in substrate-activated crystals of Pa suggests that the 70 N-terminal residues undergo a concerted shift towards the molecular core; salt bridges are probably conserved in the transition. It is proposed that the N-terminus, when "activated" by phosphorylation (via a specific kinase) behaves as an intramolecular "effector" of the R state in phosphorylase and serves as the vehicle of homotropic cooperativity between subunits of the dimer.     

Dipole–dipole interactions between tryptophan side chains and hydration water molecules dominate the observed dynamic stokes shift of lysozyme

The fluorescence intensity of tryptophan residues in hen egg-white lysozyme was measured up to 500 ps after the excitation by irradiation pulses at 290 nm. From the time-dependent variation of fluorescence intensity in a wavelength range of 320–370 nm, the energy relaxation in the dynamic Stokes shift was reconstructed as the temporal variation in wavenumber of the estimated fluorescence maximum. The relaxation was approximated by two exponential curves with decay constants of 1.2 and 26.7 ps. To interpret the relaxation, a molecular dynamics simulation of 75 ns was conducted for lysozyme immersed in a water box. From the simulation, the energy relaxation in the electrostatic interactions of each tryptophan residue was evaluated by using a scheme derived from the linear response theory. Dipole–dipole interactions between each of the Trp62 and Trp123 residues and hydration water molecules displayed an energy relaxation similar to that experimentally observed regarding time constants and magnitudes. The side chains of these residues were partly or fully exposed to the solvent. In addition, by inspecting the variation in dipole moments of the hydration water molecules around lysozyme, it was suggested that the observed relaxation could be attributed to the orientational relaxation of hydration water molecules participating in the hydrogen-bond network formed around each of the two tryptophan residues.     

Quantification of helix-helix binding affinities in micelles and lipid bilayers

A theoretical approach for estimating association free energies of alpha-helices in nonpolar media has been developed. The parameters of energy functions have been derived from DeltaDeltaG values of mutants in water-soluble proteins and partitioning of organic solutes between water and nonpolar solvents. The proposed approach was verified successfully against three sets of published data: (1) dissociation constants of alpha-helical oligomers formed by 27 hydrophobic peptides; (2) stabilities of 22 bacteriorhodopsin mutants, and (3) protein-ligand binding affinities in aqueous solution. It has been found that coalescence of helices is driven exclusively by van der Waals interactions and H-bonds, whereas the principal destabilizing contributions are represented by side-chain conformational entropy and transfer energy of atoms from a detergent or lipid to the protein interior. Electrostatic interactions of alpha-helices were relatively weak but important for reproducing the experimental data. Immobilization free energy, which originates from restricting rotational and translational rigid-body movements of molecules during their association, was found to be less than 1 kcal/mole. The energetics of amino acid substitutions in bacteriorhodopsin was complicated by specific binding of lipid and water molecules to cavities created in certain mutants.