Protein-solvent Interface Research Articles

Hydrogen bond analysis by MD simulation of copper plastocyanin at different hydration levels

Abstract Hydrogen bond interactions in copper plastocyanin aqueous solutions have been investigated in details by computer simulation. Molecular dynamics simulations of the protein in the presence of a different number of water molecules (80, 228, 682, 2516) have been performed for 110 picoseconds. Intraprotein hydrogen bonds occurring on the MD simulations are investigated and compared to the hydrogen bonds detected in the crystal structure as obtained from X-ray data. Protein-water hydrogen bonds are analyzed in terms of the average number of hydrogen bonds formed by each aminoacid residue and in terms of the number of different water molecules forming hydrogen bonds with each amino-acid residue during the simulation time. A strong influence on the protein-solvent interface of the electrical character of the aminoacid residues involved in the protein-water hydrogen bonds is evidenced. The results point out the crucial role played by the hydration level on the organization of the hydrogen bond network surrounding a protein molecule and on the protein-solvent coupling. The behaviour of the protein-solvent hydrogen bonds is discussed in connection with the protein dynamics.

New methods for the analysis of the protein-solvent interface

The protein-solvent interface is complex and may include solvent channels and cavities as well as the normal surface water molecules. We describe several algorithms for investigating the intra- and inter-molecular interactions of proteins in general but with the aim of developing methods to accurately and definitively characterise the interactions of water and other small ligands with proteins. Specifically, we present the methods which underlie three programs (AQUARIUS2, HOLE and PRO_ACT) which can be used to to look at different aspects of these interactions.

A global model of the protein-solvent interface

The solvent structure and dynamics around myoglobin is investigated at the microscopic level of detail by computer simulation. We analyze a molecular dynamics trajectory in terms of solvent mobility and probability distribution. Local events, occurring in the protein-solvent interfacial region, which are often masked by other approaches are thus revealed. Specifically, the local solvent mobility is greatly enhanced for certain locations at the protein surface and in its interior. In addition, a strong correlation between the solvent mobility and density emerges on both global and local scales. We propose a simple model where the solvent distribution measured perpendicularly to the protein surface is utilized to reconstruct the simulated network of hydration within 6 A from the protein surface with a relative error of only 17%. The global precision of this solvation model matches results obtained with more complicated models usually used in refinement procedures in x-ray and neutron experiments but with far fewer parameters. The dramatically improved correspondence between observed and calculated x-ray intensities at low resolution relative to other methods both confirms the validity of the approach used in the MD (molecular dynamics) simulations and allows the results of this study to be implemented in solvent studies on real systems.

Distribution function implied dynamics versus residence times and correlations: Solvation shells of myoglobin

The dynamics of water at the protein-solvent interface is investigated through the analysis of a molecular dynamics simulation of metmyoglobin in explicit aqueous environment. Distribution implied dynamics, harmonic and quasi-harmonic, are compared with the simulated macroscopic dynamics. The distinction between distinguishable solvent molecules and hydration sites developed in the previous paper is used. The simulated hydration region within 7 A from the protein surface is analyzed using a set of 551 hydration sites characterized by occupancy weights and temperature B-factors determined from the simulation trajectory. The precision of the isotropic harmonic and anisotropic harmonic models for the description of proximal solvent fluctuations is examined. Residence times and dipole reorientation times of water around the protein surface are compared with NMR and ESR results. A correlation between diffraction experiment quantities such as the occupancy weights and temperature factors and the residence and correlation times resulting from magnetic resonance experiments is found via comparison with simulation.

Improved Continuum Electrostatic Modelling in Proteins, with Comparison to Experiment

Electrostatic interactions in macromolecules can be calculated with the method of finite differences applied to a continuum model. The accuracy of dielectric and counterion continuum modelling has been tested for long-range interactions by comparison with available experimental data over a range of ionic strengths. Various model parameters have been adjusted. Some have little effect, such as protein dielectric and the selection of Van der Waals radii. It is shown that the reduction in interaction due to dielectric effects is overestimated when a dielectric constant of 80 is assigned to all solvent accessible regions. Improved agreement is seen when the effects of the Kirkwood correlation sphere and dielectric saturation are included. Further support for the use of dielectric saturation arises from a correlation of solvent polarization saturation with crystallographic ordered water structure. Calculations over the medium ionic strength range indicate that requiring counterions to maintain a solvent layer places too great a restriction on their approach to the protein-solvent interface. However, counterion accessibility that coincides with the solvent accessible region gives too much interaction damping. Modelling of observed ion binding sites suggests that a counterion response which includes an ion desolvation term, obtained by difference calculation, will improve the computation of ionic strength effects. This study demonstrates that there is scope for improvement in continuum electrostatics calculations, and shows that progress is possible with the inclusion of physically realistic solvent and counterion properties at the protein surface.

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 A 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.

Fatty acid interactions with rat intestinal and liver fatty acid-binding proteins expressed in Escherichia coli: A comparative 13C NMR study

Abstract Enterocytes in the small intestinal mucosa contain abundant quantities of two homologous cytosolic proteins known as intestinal and liver fatty acid-binding proteins (I- and L-FABP, respectively). To elucidate structure-function relationships for these proteins, the interactions between 13C-enriched palmitate and oleate and Escherichia coli-expressed rat I- and L-FABP were systematically compared using 13C NMR spectroscopy. NMR spectra of samples containing fatty acids (FA) and I-FABP at different molar ratios (all at pH 7.2 and 37 degrees C) exhibited a single carboxyl resonance corresponding to FA bound to I-FABP (181.4 ppm, peak I) and an additional carboxyl resonance corresponding to unbound FA in a bilayer phase (179.6 ppm). Peak I reached a maximum intensity corresponding to 1 mol of bound FA/mol of I-FABP under all sample conditions examined. NMR spectra for samples containing FA and L-FABP also exhibited a single carboxyl resonance corresponding to FA bound to L-FABP but at a different chemical shift value (182.2 ppm, peak L). Its maximum intensity varied depending on the physical state of the unbound FA (liquid crystalline or crystalline), the FA used (palmitate or oleate), and the sample pH. In the presence of a liquid crystalline (bilayer) phase, up to 1 (oleate) or 2 (palmitate) mol of FA were bound/mol of L-FABP, but in the presence of a crystalline phase (1:1 acid-soap), up to 3 mol of palmitate were bound/mol of L-FABP (all at pH 7.2). Peak I exhibited little or no ionization shift over a wide pH range (pH 3.0-11.0), and its chemical shift was unaffected by the ionization of Lys and His residues. Hence, the carboxylate group of FA bound to I-FABP was solvent inaccessible and most likely involved in an ion-pair electrostatic interaction with the delta-guanidinium moiety of an Arg residue. In contrast, peak L exhibited an ionization shift and an estimated apparent pKa value similar to that obtained for monomeric FA in water, suggesting that the carboxylate groups of FA bound to L-FABP were solvent accessible and located at or near the protein solvent interface. With decreasing pH, FA dissociated from L-FABP but not I-FABP, as monitored by NMR peak intensities. Concurrently, a large decrease in circular dichroism molar ellipticity was observed with L-FABP but not I-FABP. In conclusion, I-FABP and L-FABP are distinct with regards to their FA-binding stoichiometries, binding mechanisms, and sensitivity to pH.(ABSTRACT TRUNCATED AT 400 WORDS)

Electrostatic interactions in sperm whale myoglobin. Site specificity, roles in structural elements, and external electrostatic potential distributions.

The electrostatic free energy contribution to the stability of sperm whale ferrimyoglobin was evaluated according to the static accessibility modified Tanford-Kirkwood model. The electrostatic free energy contribution of each distinct structural element was divided into one term arising from interactions between it and other elements (interelemental) and another from interactions within the particular element itself (intraelemental). At pH 7 the majority of the terms were found to be stabilizing. The interelemental terms are the dominant ones for most structural elements. The small interelemental terms of the C and D helices are compensated by large intraelemental interactions which stabilize these short helices. Perturbations in pH can be accommodated by the structural elements through a redistribution of stabilizing and destabilizing interactions. The electrostatic potentials calculated at the surface of the protein indicate that the internal compensation of local potentials achieved during folding results in a generally neutral protein-solvent interface save for two distinct areas of nonzero potential. The accessibility of each charged atom to solvent was analyzed in terms of the surface area lost to charged, polar and nonpolar atoms separately. The net solvent accessibility lost parallels closely that lost to nonpolar atoms alone, indicating a specific role for nonpolar atoms in defining dielectric shielding of charged atoms, aside from their participation in the well-known hydrophobic interactions.

Mechanism of surface peptide proton exchange in bovine pancreatic trypsin inhibitor salt effects and O-protonation

The acid-catalyzed hydrogen exchange rate constants k H, and the base-catalyzed rate constants k OH, have been determined (in the preceding paper) for the 25 most rapidly exchanging NH groups of bovine pancreatic trypsin inhibitor. Most of these NH groups are at the protein-solvent interface. The correlation of k H, but not k OH, with the static accessibility and hydrogen bonding of the peptide carbonyl O atom indicates that the mechanism of acid catalysis in proteins involves O-protonation. Agreement between the ionic strength dependence observed for k H and k OH and the ionic strength dependence calculated for an O-protonation mechanism supports this conclusion. N-protonation for acid catalysis, as well as N-deprotonation for base catalysis, have traditionally been assumed in the mechanism of the chemical step in peptide amide proton exchange. A preference for the alternative O-protonation mechanism has far-reaching implications in the interpretation of protein hydrogen exchange kinetics. With an O-protonation mechanism, acid-catalyzed rates of surface NH groups are primarily a function of the average solvent accessibility of the carbonyl O atoms in the dynamic solution structure, while base-catalyzed rates of surface NH groups measure solvent accessibility of the peptide N. The relative dynamic accessibilities of peptide O atoms, as measured by relative values of k H (corrected for electrostatic effects), correlate with O static accessibilities in the crystal structure. A lower correlation of static accessibility of N atoms with k OH is observed for surface NH groups in peptide groups in which the carbonyl O is not hydrogen bonded. For some surface NH groups, the observed pH of minimum rate, pH min, deviates widely from the pH min of model compounds. This is explained as the combined result of electrostatic effects and of the differences in accessibility of the carbonyl O and N atoms that result in a change in the relative values of k H and k OH as compared to those of model peptides. A mechanism whereby exchange of interior sites is catalyzed by interactions of catalysis ions with protein surface atoms via charge transfer is suggested.

Hydrogen kinetics of peptide amide protons at the bovine pancreatic trypsin inhibitor protein-solvent interface

Hydrogen exchange rate constants of the 25 most rapidly exchanging peptide amide protons in bovine pancreatic trypsin inhibitor have been determined over a range of pH that spans pH min, the pH of minimum rate. Most of these are on the protein surface, exposed to solvent and not hydrogen bonded in the crystal structure. Contrary to commonly held assumptions, the exchange kinetics of surface NH groups are not equivalent to the kinetics of NH groups in peptides in the extended configuration. All surface NH groups exchange more slowly than NH groups in model peptides, with rate constants distributed over a range of more than two orders of magnitude. In addition, their pH min values vary widely. For most of the surface NH groups, pH min is lower than in model compounds and, for several, pH min is < 1. These results indicate that the local environment of the surface peptide groups when the exchange event occurs is very different from that of extended peptides. Analysis based on consideration of an O-protonation mechanism for acid catalysis and of electrostatic effects on exchange kinetics further indicates (see the accompanying paper) that, in general, exchange of surface NH groups occurs from a conformation of the protein approximated by the crystal structure. The 1H- 2H exchange rate constants were measured from 300 MHz nuclear magnetic resonance spectra in which assigned surface N 1H resonances are resolved by the use of partially deuterated protein samples. A marked pH dependence of the chemical shifts observed in the pH range 1 to 4.5 for several surface NH groups reflects the titration of nearby carboxyl groups.

Anion binding and pH-dependent electrostatic effects in ribonuclease.

The solvent-accessibility-modified, Tanford-Kirkwood, discrete charge model for electrostatic effects is applied to both ribonuclease A and ribonuclease S. The behavior of individual titratable sites and the pH-dependent free energy of denaturation are correctly predicted. The use of the solvent-accessibility factor in reducing charge-site interactions introduces a higher Coulombic shielding for solvent-exposed sites. This shielding is interpreted as a higher local strength or alternatively a higher effective dielectric constant. Specific anion binding sites are determined by locating areas of high positive electrostatic potential at the protein solvent interface. The potential and thus the anion affinity of a given site are calculated and shown to vary with the pH-dependent charge array. pH-dependent anion binding constants are calculated for the ribonuclease S active site. These binding constants and the predicted response of the active-site histidine pK1/2 values to anion binding are shown to agree with experimental determinations.

Magnetic cross-relaxation among protons in protein solutions

The magnetic spin-lattice relaxation rates of solvent water nuclei are known to increase upon addition of diamagnetic solute protein. This enhancement of the relaxation rate is a function of magnetic field, and the orientational relaxation time of the protein molecules can be deduced from analysis of the field-dependent relaxation rates. Although the nature of the interactions that convey information about the dynamics of protein motion to the solvent molecules is not established, it is known that there is a contribution to the relaxation rates of solvent protons that plays no role in the relaxation of solvent deuterons and 17O nuclei. We show here that the additional interaction arises from a cross-relaxation process between solvent and solute protons. We introduce a heuristic three-parameter model in which protein protons and solvent protons are considered as two separate thermodynamic systems that interact across the protein-solvent interface. The three parameters are the intrinsic relaxation rates of each system and a cross-relaxation term. The sign of the latter term must always be positive, for all values of magnetic field, in order for magnetization energy to flow from the hotter to the cooler system. We find that the magnetic field-dependence of the cross-relaxation contribution is much like that of the remaining solvent proton relaxation, i.e., about the same as the deuteron relaxation field dependence. This finding is not compatible with the predictions of expressions for the cross-relaxation that have been used by other authors, but not applied to data over a wide range of magnetic field strength. The model predicts that the relaxation behavior of both the protein protons and the solvent protons is the sum of two exponentials, the relative contributions of which would vary with protein concentration and solvent isotopic composition in a fashion suggestive of the presence of two classes of protein protons, when there is in reality only one. This finding has immediate implications for the interpretation of published proton relaxation rates in complex systems such as tissues; these data should be reexamined with cross-relaxation taken into account.

Hydrocarbon binding by particles of bovine casein

Casein, isolated from milk by high-speed centrifugation, was washed with water and dehydrated by serial transfer through liquids of decreasing polarity to provide a system consisting of protein dispersed in anhydrous, liquid hydrocarbon. Heats of fusion of solvents associated with casein were measured to determine how much hydrocarbon is free and can be frozen and what fraction is bound to the protein and immobilized. Ring substitution decreased binding in the aromatic compounds studied, i.e. : 8.4 mmoles benzene, 4.7 mmoles p-xylene, and 3.3 mmoles mesitylene bound per gram casein. Binding of even n-paraffins was significantly less than that of the aromatic compounds and decreased with increasing chain length, i.e. : 2.0 mmoles C 8H 18, 1.9 mmoles C 10H 22, 1.6 mmoles C 12H 26, 1.1 mmoles C 14H 30, and 1.0 mmoles C 16H 34 bound per gram casein. The observation of unfreezable hydrocarbon indicates that some solvents, commonly considered as hydrophobic, attain preferred orientation at the protein-solvent interface.