1
|
Kruchinin SE, Kislinskaya EE, Chuev GN, Fedotova MV. Protein 3D Hydration: A Case of Bovine Pancreatic Trypsin Inhibitor. Int J Mol Sci 2022; 23:ijms232314785. [PMID: 36499117 PMCID: PMC9737982 DOI: 10.3390/ijms232314785] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2022] [Revised: 11/20/2022] [Accepted: 11/23/2022] [Indexed: 11/29/2022] Open
Abstract
Characterization of the hydrated state of a protein is crucial for understanding its structural stability and function. In the present study, we have investigated the 3D hydration structure of the protein BPTI (bovine pancreatic trypsin inhibitor) by molecular dynamics (MD) and the integral equation method in the three-dimensional reference interaction site model (3D-RISM) approach. Both methods have found a well-defined hydration layer around the protein and revealed the localization of BPTI buried water molecules corresponding to the X-ray crystallography data. Moreover, under 3D-RISM calculations, the obtained positions of waters bound firmly to the BPTI sites are in reasonable agreement with the experimental results mentioned above for the BPTI crystal form. The analysis of the 3D hydration structure (thickness of hydration shell and hydration numbers) was performed for the entire protein and its polar and non-polar parts using various cut-off distances taken from the literature as well as by a straightforward procedure proposed here for determining the thickness of the hydration layer. Using the thickness of the hydration shell from this procedure allows for calculating the total hydration number of biomolecules properly under both methods. Following this approach, we have obtained the thickness of the BPTI hydration layer of 3.6 Å with 369 water molecules in the case of MD simulation and 3.9 Å with 333 water molecules in the case of the 3D-RISM approach. The above procedure was also applied for a more detailed description of the BPTI hydration structure near the polar charged and uncharged radicals as well as non-polar radicals. The results presented for the BPTI as an example bring new knowledge to the understanding of protein hydration.
Collapse
Affiliation(s)
- Sergey E. Kruchinin
- G.A. Krestov Institute of Solution Chemistry, The Russian Academy of Sciences, Akademicheskaya St., 1, 153045 Ivanovo, Russia
| | - Ekaterina E. Kislinskaya
- Department of Fundamental and Applied Chemistry, Institute of Mathematics, Information Technology and Science, Ivanovo State University, Ermak St., 39, 153025 Ivanovo, Russia
| | - Gennady N. Chuev
- Institute of Theoretical and Experimental Biophysics, The Russian Academy of Sciences, Institutskaya St., Pushchino, 142290 Moscow, Russia
- Correspondence: (G.N.C.); (M.V.F.)
| | - Marina V. Fedotova
- G.A. Krestov Institute of Solution Chemistry, The Russian Academy of Sciences, Akademicheskaya St., 1, 153045 Ivanovo, Russia
- Correspondence: (G.N.C.); (M.V.F.)
| |
Collapse
|
2
|
Abstract
Much of biology happens at the protein-water interface, so all dynamical processes in this region are of fundamental importance. Local structural fluctuations in the hydration layer can be probed by 17O magnetic relaxation dispersion (MRD), which, at high frequencies, measures the integral of a biaxial rotational time correlation function (TCF)-the integral rotational correlation time. Numerous 17O MRD studies have demonstrated that this correlation time, when averaged over the first hydration shell, is longer than in bulk water by a factor 3-5. This rotational perturbation factor (RPF) has been corroborated by molecular dynamics simulations, which can also reveal the underlying molecular mechanisms. Here, we address several outstanding problems in this area by analyzing an extensive set of molecular dynamics data, including four globular proteins and three water models. The vexed issue of polarity versus topography as the primary determinant of hydration water dynamics is resolved by establishing a protein-invariant exponential dependence of the RPF on a simple confinement index. We conclude that the previously observed correlation of the RPF with surface polarity is a secondary effect of the correlation between polarity and confinement. Water rotation interpolates between a perturbed but bulk-like collective mechanism at low confinement and an exchange-mediated orientational randomization (EMOR) mechanism at high confinement. The EMOR process, which accounts for about half of the RPF, was not recognized in previous simulation studies, where only the early part of the TCF was examined. Based on the analysis of the experimentally relevant TCF over its full time course, we compare simulated and measured RPFs, finding a 30% discrepancy attributable to force field imperfections. We also compute the full 17O MRD profile, including the low-frequency dispersion produced by buried water molecules. Computing a local RPF for each hydration shell, we find that the perturbation decays exponentially with a decay "length" of 0.3 shells and that the second and higher shells account for a mere 3% of the total perturbation measured by 17O MRD. The only long-range effect is a weak water alignment in the electric field produced by an electroneutral protein (not screened by counterions), but this effect is negligibly small for 17O MRD. By contrast, we find that the 17O TCF is significantly more sensitive to the important short-range perturbations than the other two TCFs examined here.
Collapse
Affiliation(s)
- Filip Persson
- Division of Biophysical Chemistry, Department of Chemistry, Lund University, P.O. Box 124, SE-22100 Lund, Sweden
| | - Pär Söderhjelm
- Division of Biophysical Chemistry, Department of Chemistry, Lund University, P.O. Box 124, SE-22100 Lund, Sweden
| | - Bertil Halle
- Division of Biophysical Chemistry, Department of Chemistry, Lund University, P.O. Box 124, SE-22100 Lund, Sweden
| |
Collapse
|
3
|
Statistical survey of the buried waters in the Protein Data Bank. Amino Acids 2015; 48:193-202. [PMID: 26315961 DOI: 10.1007/s00726-015-2064-4] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2015] [Accepted: 07/29/2015] [Indexed: 12/18/2022]
Abstract
The structures of buried water molecules were studied in an ensemble of high-quality and non-redundant protein crystal structures. Buried water molecules were clustered and classified in lake-like clusters, which are completely isolated from the bulk solvent, and bay-like clusters, which are in contact with the bulk solvent through a surface water molecule. Buried water molecules are extremely common: lake-like clusters are found in 89 % of the protein crystal structures and bay-like clusters in 93 %. Clusters with only one water molecule are much more common than larger clusters. Both cluster types incline to be surrounded by loop residues, and to a minor extent by residues in extended secondary structure. Helical residues on the contrary do not tend to surround clusters of buried water molecules. One buried water molecule is found every 30-50 amino acid residues, depending on the secondary structures that are more abundant in the protein. Both main- and side-chain atoms are in contact with buried waters; they form four hydrogen bonds with the first water and 1-1.5 additional hydrogen bond for each additional water in the cluster. Consequently, buried water molecules appear to be firmly packed and rigid like the protein atoms. In this regard, it is remarkable to observe that prolines often surround water molecules buried in the protein interior. Interestingly, clusters of buried water molecules tend to be just beneath the protein surface. Moreover, water molecules tend to form a one-dimensional wire rather than more compact arrangements. This agrees with recent evidence of the mechanisms of solvent exchange between internal cavities and bulk solvent.
Collapse
|
4
|
Persson F, Halle B. Transient access to the protein interior: simulation versus NMR. J Am Chem Soc 2013; 135:8735-48. [PMID: 23675835 DOI: 10.1021/ja403405d] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Many proteins rely on rare structural fluctuations for their function, whereby solvent and other small molecules gain transient access to internal cavities. In magnetic relaxation dispersion (MRD) experiments, water molecules buried in such cavities are used as intrinsic probes of the intermittent protein motions that govern their exchange with external solvent. While this has allowed a detailed characterization of exchange kinetics for several proteins, little is known about the exchange mechanism. Here, we use a millisecond all-atom MD trajectory produced by Shaw et al. (Science2010, 330, 341) to characterize water exchange from the four internal hydration sites in the protein bovine pancreatic trypsin inhibitor. Using a recently developed stochastic point process approach, we compute the survival correlation function probed by MRD experiments as well as other quantities designed to validate the exchange-mediated orientational randomization (EMOR) model used to interpret the MRD data. The EMOR model is found to be quantitatively accurate, and the simulation reproduces the experimental mean survival times for all four sites with activation energy discrepancies in the range 0-3 kBT. On the other hand, the simulated hydration sites are somewhat too flexible, and the water flip barrier is underestimated by up to 6 kBT. The simulation reveals that water molecules gain access to the internal sites by a transient aqueduct mechanism, migrating as single-file water chains through transient (<5 ns) tunnels or pores. The present study illustrates the power of state-of-the-art molecular dynamics simulations in validating and extending experimental results.
Collapse
Affiliation(s)
- Filip Persson
- Biophysical Chemistry, Lund University, P.O. Box 124, SE-22100 Lund, Sweden
| | | |
Collapse
|
5
|
Xue Y, Ward JM, Yuwen T, Podkorytov IS, Skrynnikov NR. Microsecond Time-Scale Conformational Exchange in Proteins: Using Long Molecular Dynamics Trajectory To Simulate NMR Relaxation Dispersion Data. J Am Chem Soc 2012; 134:2555-62. [PMID: 22206299 DOI: 10.1021/ja206442c] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Yi Xue
- Department of Chemistry,
Purdue University, 560 Oval
Drive, West Lafayette, Indiana 47907-2084, United States
| | - Joshua M. Ward
- Department of Chemistry,
Purdue University, 560 Oval
Drive, West Lafayette, Indiana 47907-2084, United States
| | - Tairan Yuwen
- Department of Chemistry,
Purdue University, 560 Oval
Drive, West Lafayette, Indiana 47907-2084, United States
| | - Ivan S. Podkorytov
- Department of Chemistry,
Purdue University, 560 Oval
Drive, West Lafayette, Indiana 47907-2084, United States
| | - Nikolai R. Skrynnikov
- Department of Chemistry,
Purdue University, 560 Oval
Drive, West Lafayette, Indiana 47907-2084, United States
| |
Collapse
|
6
|
Korb JP, Goddard Y, Pajski J, Diakova G, Bryant RG. Extreme-Values Statistics and Dynamics of Water at Protein Interfaces. J Phys Chem B 2011; 115:12845-58. [DOI: 10.1021/jp2053426] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Jean-Pierre Korb
- Physique de la Matière Condensée, Ecole Polytechnique, CNRS, 91128 Palaiseau, France
| | - Yanina Goddard
- Chemistry Department, University of Virginia, Charlottesville, Virginia 22904, United States
| | - Jason Pajski
- Chemistry Department, University of Virginia, Charlottesville, Virginia 22904, United States
| | - Galina Diakova
- Chemistry Department, University of Virginia, Charlottesville, Virginia 22904, United States
| | - Robert G. Bryant
- Chemistry Department, University of Virginia, Charlottesville, Virginia 22904, United States
| |
Collapse
|
7
|
Oikawa M, Yonetani Y. Molecular dynamics free energy calculations to assess the possibility of water existence in protein nonpolar cavities. Biophys J 2010; 98:2974-83. [PMID: 20550910 DOI: 10.1016/j.bpj.2010.01.029] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2009] [Revised: 01/13/2010] [Accepted: 01/14/2010] [Indexed: 11/19/2022] Open
Abstract
Are protein nonpolar cavities filled with water molecules? Although many experimental and theoretical investigations have been done, particularly for the nonpolar cavity of IL-1 beta, the results are still conflicting. To study this problem from the thermodynamic point of view, we calculated hydration free energies of four protein nonpolar cavities by means of the molecular dynamics thermodynamic integration method. In addition to the IL-1 beta cavity (69 A(3)), we selected the three largest nonpolar cavities of AvrPphB (81 A(3)), Trp repressor (87 A(3)), and hemoglobin (108 A(3)) from the structural database, in view of the simulation result from another study that showed larger nonpolar cavities are more likely to be hydrated. The calculations were performed with flexible and rigid protein models. The calculated free energy changes were all positive; hydration of the nonpolar cavities was energetically unfavorable for all four cases. Because hydration of smaller cavities should happen more rarely, we conclude that existing protein nonpolar cavities are not likely to be hydrated. Although a possibility remains for much larger nonpolar cavities, such cases are not found experimentally. We present a hypothesis to explain this: hydrated nonpolar cavities are quite unstable and the conformation could not be maintained.
Collapse
Affiliation(s)
- Masataka Oikawa
- Computational Biology Group, Quantum Beam Science Directorate, Japan Atomic Energy Agency, Kyoto, Japan.
| | | |
Collapse
|
8
|
Brownian dynamics of water confined in AOT reverse micelles: A field-cycling deuteron NMR relaxometry study. CR CHIM 2010. [DOI: 10.1016/j.crci.2009.09.004] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
9
|
Persson E, Halle B. Nanosecond to Microsecond Protein Dynamics Probed by Magnetic Relaxation Dispersion of Buried Water Molecules. J Am Chem Soc 2008; 130:1774-87. [DOI: 10.1021/ja0775873] [Citation(s) in RCA: 83] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Erik Persson
- Department of Biophysical Chemistry, Center for Molecular Protein Science, Lund University, SE-22100 Lund, Sweden
| | - Bertil Halle
- Department of Biophysical Chemistry, Center for Molecular Protein Science, Lund University, SE-22100 Lund, Sweden
| |
Collapse
|
10
|
Bueno M, Camacho CJ, Sancho J. SIMPLE estimate of the free energy change due to aliphatic mutations: Superior predictions based on first principles. Proteins 2007; 68:850-62. [PMID: 17523191 DOI: 10.1002/prot.21453] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
The bioinformatics revolution of the last decade has been instrumental in the development of empirical potentials to quantitatively estimate protein interactions for modeling and design. Although computationally efficient, these potentials hide most of the relevant thermodynamics in 5-to-40 parameters that are fitted against a large experimental database. Here, we revisit this longstanding problem and show that a careful consideration of the change in hydrophobicity, electrostatics, and configurational entropy between the folded and unfolded state of aliphatic point mutations predicts 20-30% less false positives and yields more accurate predictions than any published empirical energy function. This significant improvement is achieved with essentially no free parameters, validating past theoretical and experimental efforts to understand the thermodynamics of protein folding. Our first principle analysis strongly suggests that both the solute-solute van der Waals interactions in the folded state and the electrostatics free energy change of exposed aliphatic mutations are almost completely compensated by similar interactions operating in the unfolded ensemble. Not surprisingly, the problem of properly accounting for the solvent contribution to the free energy of polar and charged group mutations, as well as of mutations that disrupt the protein backbone remains open.
Collapse
Affiliation(s)
- Marta Bueno
- Department of Computational Biology, University of Pittsburgh, Pennsylvania, USA
| | | | | |
Collapse
|