1
|
|
2
|
Utiramerur S, Paulaitis M. Analysis of Cooperativity and Group Additivity in the Hydration of 1,2-Dimethoxyethane. J Phys Chem B 2021; 125:1660-1666. [DOI: 10.1021/acs.jpcb.0c10729] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Sowmi Utiramerur
- William G. Lowrie Department of Chemical and Biomolecular Engineering, Ohio State University, Columbus, Ohio 43210, United States
| | - Michael Paulaitis
- William G. Lowrie Department of Chemical and Biomolecular Engineering, Ohio State University, Columbus, Ohio 43210, United States
- The Center for Nanomedicine at the Wilmer Eye Institute, Johns Hopkins University School of Medicine, Baltimore, Maryland 21231, United States
| |
Collapse
|
3
|
The Molecular Interaction Process. J Pharm Sci 2020; 109:154-160. [DOI: 10.1016/j.xphs.2019.10.045] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2019] [Revised: 10/16/2019] [Accepted: 10/24/2019] [Indexed: 01/14/2023]
|
4
|
Kuo TC, Huang YC, Matulis D, Chen WY. Molecular self-interactions of ribonuclease A revealed by isothermal titration calorimetry and self-interaction chromatography – Effects of anisotropy of protein surface charges. J Taiwan Inst Chem Eng 2019. [DOI: 10.1016/j.jtice.2018.12.008] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
5
|
In Silico Prediction of Diffusion Interaction Parameter (kD), a Key Indicator of Antibody Solution Behaviors. Pharm Res 2018; 35:193. [DOI: 10.1007/s11095-018-2466-6] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2018] [Accepted: 07/24/2018] [Indexed: 12/11/2022]
|
6
|
Meric G, Robinson AS, Roberts CJ. Driving Forces for Nonnative Protein Aggregation and Approaches to Predict Aggregation-Prone Regions. Annu Rev Chem Biomol Eng 2017; 8:139-159. [DOI: 10.1146/annurev-chembioeng-060816-101404] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Gulsum Meric
- Department of Chemical and Biomolecular Engineering, University of Delaware, Newark, Delaware 19716
| | - Anne S. Robinson
- Department of Chemical and Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118
| | - Christopher J. Roberts
- Department of Chemical and Biomolecular Engineering, University of Delaware, Newark, Delaware 19716
| |
Collapse
|
7
|
Ma K, Zhao L. The Opposite Effect of Metal Ions on Short-/Long-Range Water Structure: A Multiple Characterization Study. Int J Mol Sci 2016; 17:ijms17050602. [PMID: 27120598 PMCID: PMC4881438 DOI: 10.3390/ijms17050602] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2016] [Revised: 04/14/2016] [Accepted: 04/18/2016] [Indexed: 12/21/2022] Open
Abstract
Inorganic electrolyte solutions are very important in our society as they dominate many biochemical and geochemical processes. Herein, an in-depth study was performed to illustrate the ion-induced effect on water structure by coupling NMR, viscometer, Raman and Molecular Dynamic (MD) simulations. The NMR coefficient (BNMR) and diffusion coefficient (D) from NMR, and viscosity coefficient (Bvis) from a viscometer all proved that dissolved metal ions are capable of enhancing the association degree of adjacent water molecules, and the impact on water structure decreased in the order of Cr3+ > Fe3+ > Cu2+ > Zn2+. This regularity was further evidenced by Raman analysis; however, the deconvoluted Raman spectrum indicated the decrease in high association water with salt concentration and the increase in low association water before 200 mmol·L−1. By virtue of MD simulations, the opposite changing manner proved to be the result of the opposite effect on short-/long-range water structure induced by metal ions. Our results may help to explain specific protein denaturation induced by metal ions.
Collapse
Affiliation(s)
- Kai Ma
- School of Environmental Science and Engineering, Tianjin University, Tianjin 300072, China.
| | - Lin Zhao
- School of Environmental Science and Engineering, Tianjin University, Tianjin 300072, China.
| |
Collapse
|
8
|
Kokkinos P, Syngouna VI, Tselepi MA, Bellou M, Chrysikopoulos CV, Vantarakis A. Transport of Human Adenoviruses in Water Saturated Laboratory Columns. FOOD AND ENVIRONMENTAL VIROLOGY 2015; 7:122-131. [PMID: 25578176 DOI: 10.1007/s12560-014-9179-8] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/30/2014] [Accepted: 12/23/2014] [Indexed: 06/04/2023]
Abstract
Groundwater may be contaminated with infective human enteric viruses from various wastewater discharges, sanitary landfills, septic tanks, agricultural practices, and artificial groundwater recharge. Coliphages have been widely used as surrogates of enteric viruses, because they share many fundamental properties and features. Although a large number of studies focusing on various factors (i.e. pore water solution chemistry, fluid velocity, moisture content, temperature, and grain size) that affect biocolloid (bacteria, viruses) transport have been published over the past two decades, little attention has been given toward human adenoviruses (hAdVs). The main objective of this study was to evaluate the effect of pore water velocity on hAdV transport in water saturated laboratory-scale columns packed with glass beads. The effects of pore water velocity on virus transport and retention in porous media was examined at three pore water velocities (0.39, 0.75, and 1.22 cm/min). The results indicated that all estimated average mass recovery values for hAdV were lower than those of coliphages, which were previously reported in the literature by others for experiments conducted under similar experimental conditions. However, no obvious relationship between hAdV mass recovery and water velocity could be established from the experimental results. The collision efficiencies were quantified using the classical colloid filtration theory. Average collision efficiency, α, values decreased with decreasing flow rate, Q, and pore water velocity, U, but no significant effect of U on α was observed. Furthermore, the surface properties of viruses and glass beads were used to construct classical DLVO potential energy profiles. The results revealed that the experimental conditions of this study were unfavorable to deposition and that no aggregation between virus particles is expected to occur. A thorough understanding of the key processes governing virus transport is pivotal for public health protection.
Collapse
Affiliation(s)
- P Kokkinos
- Environmental Microbiology Unit, Department of Public Health, School of Medicine, University of Patras, 26500, Patras, Greece,
| | | | | | | | | | | |
Collapse
|
9
|
Wakamatsu T. Method and apparatus for characterization of electric field-induced aggregation in pre-crystalline protein solutions. THE REVIEW OF SCIENTIFIC INSTRUMENTS 2015; 86:015112. [PMID: 25638127 DOI: 10.1063/1.4906328] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The article presents a method and an apparatus for the characterization of protein aggregation under an applied internal electric field. The method is based on a forward light scattering technique that is highly sensitive to aggregates in pre-crystalline protein solutions. Transparent conductive films are used as electrodes for a planar thin sample cell, which enables precise measurement of the forward light scattering at small angles through the electrodes. Evaluation of the protein aggregation under applied electric fields was demonstrated for a model lysozyme protein. In situ measurements of crystallizing lysozyme solutions under a low applied voltage revealed that the forward static light scattering profiles changed with time into power law profiles. This indicates the formation of lysozyme fractal clusters under applied electric fields in the pre-crystalline state. The method and the apparatus presented here can sensitively evaluate the promotion process in protein crystallization under an applied electric field.
Collapse
Affiliation(s)
- Takashi Wakamatsu
- Department of Electrical and Electronic System Engineering, National Institute of Technology, Ibaraki College, Hitachinaka, Ibaraki 312-8508, Japan
| |
Collapse
|
10
|
Abstract
We investigate the role of water molecules in 89 protein–RNA complexes taken from the Protein Data Bank. Those with tRNA and single-stranded RNA are less hydrated than with duplex or ribosomal proteins. Protein–RNA interfaces are hydrated less than protein–DNA interfaces, but more than protein–protein interfaces. Majority of the waters at protein–RNA interfaces makes multiple H-bonds; however, a fraction do not make any. Those making H-bonds have preferences for the polar groups of RNA than its partner protein. The spatial distribution of waters makes interfaces with ribosomal proteins and single-stranded RNA relatively ‘dry’ than interfaces with tRNA and duplex RNA. In contrast to protein–DNA interfaces, mainly due to the presence of the 2′OH, the ribose in protein–RNA interfaces is hydrated more than the phosphate or the bases. The minor groove in protein–RNA interfaces is hydrated more than the major groove, while in protein–DNA interfaces it is reverse. The strands make the highest number of water-mediated H-bonds per unit interface area followed by the helices and the non-regular structures. The preserved waters at protein–RNA interfaces make higher number of H-bonds than the other waters. Preserved waters contribute toward the affinity in protein–RNA recognition and should be carefully treated while engineering protein–RNA interfaces.
Collapse
Affiliation(s)
- Amita Barik
- Department of Biotechnology, Indian Institute of Technology Kharagpur, Kharagpur-721302, India
| | - Ranjit Prasad Bahadur
- Department of Biotechnology, Indian Institute of Technology Kharagpur, Kharagpur-721302, India
| |
Collapse
|
11
|
Bhirde AA, Hassan SA, Harr E, Chen X. Role of Albumin in the Formation and Stabilization of Nanoparticle Aggregates in Serum Studied by Continuous Photon Correlation Spectroscopy and Multiscale Computer Simulations. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2014; 118:16199-16208. [PMID: 25221633 PMCID: PMC4159775 DOI: 10.1021/jp5034068] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/07/2014] [Revised: 06/25/2014] [Indexed: 05/20/2023]
Abstract
Recently, small (<5 nm diameter) nanoparticles (NPs) have shown improved in vivo biocompatibility compared to that of larger (>10 nm) NPs. However, the fate of small NPs under physiological conditions is poorly understood and remains unexplored. Here, the long-term aggregation behavior of gold nanoparticles (AuNPs) exposed to serum proteins in a near-physiological setup is studied using continuous photon correlation spectroscopy and computer simulations. It is found that the medium, temperature, and NP concentration affect the aggregation of AuNPs, but the observed aggregates are much smaller than previously reported. Simulations show that a single layer of albumin is deposited on the NP surface, but the properties of the aggregates (size, shape, and internal structure) depend critically on the charge distribution on the proteins, which changes with the conditions of the solution. These results explain the seemingly conflicting data reported in the literature regarding the size of aggregates and the morphology of the albumin corona. The simulations suggest that controlling the concentration of NPs as well as the pH and ionic strength of the solution prior to intravenous administration may help to preserve properties of the functionalized NPs in the bloodstream.
Collapse
Affiliation(s)
- Ashwinkumar A. Bhirde
- Laboratory
of Molecular Imaging and Nanomedicine, National Institute of Biomedical
Imaging and Bioengineering, and Center for Molecular Modeling, Division of
Computational Bioscience, CIT, National
Institutes of Health, Bethesda, Maryland 20892, United States
| | - Sergio A. Hassan
- Laboratory
of Molecular Imaging and Nanomedicine, National Institute of Biomedical
Imaging and Bioengineering, and Center for Molecular Modeling, Division of
Computational Bioscience, CIT, National
Institutes of Health, Bethesda, Maryland 20892, United States
- (S.A.H.) E-mail: ; Phone: 301-402-1382
| | - Erick Harr
- Laboratory
of Molecular Imaging and Nanomedicine, National Institute of Biomedical
Imaging and Bioengineering, and Center for Molecular Modeling, Division of
Computational Bioscience, CIT, National
Institutes of Health, Bethesda, Maryland 20892, United States
| | - Xiaoyuan Chen
- Laboratory
of Molecular Imaging and Nanomedicine, National Institute of Biomedical
Imaging and Bioengineering, and Center for Molecular Modeling, Division of
Computational Bioscience, CIT, National
Institutes of Health, Bethesda, Maryland 20892, United States
- (X.C.) E-mail: ;
Phone: 301-451-4246
| |
Collapse
|
12
|
Quang LJ, Sandler SI, Lenhoff AM. Anisotropic Contributions to Protein–Protein Interactions. J Chem Theory Comput 2014; 10:835-45. [DOI: 10.1021/ct4006695] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Leigh J. Quang
- Department of Chemical and
Biomolecular Engineering, University of Delaware, Newark, Delaware 19716, United States of America
| | - Stanley I. Sandler
- Department of Chemical and
Biomolecular Engineering, University of Delaware, Newark, Delaware 19716, United States of America
| | - Abraham M. Lenhoff
- Department of Chemical and
Biomolecular Engineering, University of Delaware, Newark, Delaware 19716, United States of America
| |
Collapse
|
13
|
Gillespie CM, Asthagiri D, Lenhoff AM. Polymorphic Protein Crystal Growth: Influence of Hydration and Ions in Glucose Isomerase. CRYSTAL GROWTH & DESIGN 2014; 14:46-57. [PMID: 24955067 PMCID: PMC4061714 DOI: 10.1021/cg401063b] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
Crystal polymorphs of glucose isomerase were examined to characterize the properties and to quantify the energetics of protein crystal growth. Transitions of polymorph stability were measured in poly(ethylene glycol)/NaCl solutions, and one transition point was singled out for more detailed quantitative analysis. Single crystal x-ray diffraction was used to confirm space groups and identify complementary crystal structures. Crystal polymorph stability was found to depend on the NaCl concentration, with stability transitions requiring > 1 M NaCl combined with a low concentration of PEG. Both salting-in and salting-out behavior was observed and was found to differ for the two polymorphs. For NaCl concentrations above the observed polymorph transition, the increase in solubility of the less stable polymorph together with an increase in the osmotic second virial coefficient suggests that changes in protein hydration upon addition of salt may explain the experimental trends. A combination of atomistic and continuum models was employed to dissect this behavior. Molecular dynamics simulations of the solvent environment were interpreted using quasi-chemical theory to understand changes in protein hydration as a function of NaCl concentration. The results suggest that protein surface hydration and Na+ binding may introduce steric barriers to contact formation, resulting in polymorph selection.
Collapse
|
14
|
Blanco MA, Sahin E, Robinson AS, Roberts CJ. Coarse-grained model for colloidal protein interactions, B(22), and protein cluster formation. J Phys Chem B 2013; 117:16013-28. [PMID: 24289039 DOI: 10.1021/jp409300j] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Reversible protein cluster formation is an important initial step in the processes of native and non-native protein aggregation, but involves relatively long time and length scales for detailed atomistic simulations and extensive mapping of free energy landscapes. A coarse-grained (CG) model is presented to semiquantitatively characterize the thermodynamics and key configurations involved in the landscape for protein oligomerization, as well as experimental measures of interactions such as the osmotic second virial coefficient (B22). Based on earlier work (Grüenberger et al., J. Phys. Chem. B 2013, 117, 763), this CG model treats proteins as rigid bodies composed of one bead per amino acid, with each amino acid having specific parameters for its size, hydrophobicity, and charge. The net interactions are a combination of steric repulsions, short-range attractions, and screened long-range charge-charge interactions. Model parametrization was done by fitting simulation results against experimental value of B22 as a function of solution ionic strength for α-chymotrypsinogen A and γD-Crystallin (gD-Crys). The CG model is applied to characterize the pairwise interactions and dimerization of gD-Crys and the dependence on temperature, protein concentration, and ionic strength. The results illustrate that at experimentally relevant conditions where stable dimers do not form, the entropic contributions are predominant in the free-energy of protein cluster formation and colloidal protein interactions, arguing against interpretations that treat B22 primarily from energetic considerations alone. Additionally, the results suggest that electrostatic interactions help to modulate the population of the different stable configurations for protein nearest-neighbor pairs, while short-range attractions determine the relative orientations of proteins within these configurations. Finally, simulation results are combined with Principal Component Analysis to identify those amino-acids/surface patches that form interprotein contacts at conditions that favor dimerization of gD-Crys. The resulting regions agree with previously found aggregation-prone sites, as well as suggesting new ones that may be important.
Collapse
Affiliation(s)
- Marco A Blanco
- Department of Chemical and Biomolecular Engineering and Center for Molecular and Engineering Thermodynamics, University of Delaware , Newark, Delaware 19176, United States
| | | | | | | |
Collapse
|
15
|
|
16
|
Pellicane G, Cavero M. Theoretical study of interactions of BSA protein in a NaCl aqueous solution. J Chem Phys 2013; 138:115103. [DOI: 10.1063/1.4794919] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
17
|
Grünberger A, Lai PK, Blanco MA, Roberts CJ. Coarse-grained modeling of protein second osmotic virial coefficients: sterics and short-ranged attractions. J Phys Chem B 2013; 117:763-70. [PMID: 23245189 DOI: 10.1021/jp308234j] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
A series of coarse-grained models, with different levels of structural resolution, were tested to calculate the steric contributions to protein osmotic second virial coefficients (B(22,S)) for proteins ranging from small single-domain molecules to large multidomain molecules, using the recently developed Mayer sampling method. B(22,S) was compared for different levels of coarse-graining: four-beads-per-amino-acid (4bAA), one-bead-per-amino-acid (1bAA), one-sphere-per-domain (1sD), and one-sphere-per-protein (1sP). Values for the 1bAA and 4bAA models were quantitatively indistinguishable for both spherical and nonspherical proteins, and the agreement with values from all-atom models improved with increasing protein size, making the CG approach attractive for large proteins of biotechnological interest. Interestingly, in the absence of detailed structural information, the hydrodynamic radius (R(h)) along with a simple 1sP approximation provided reasonably accurate values for B(22,S) for both globular and highly asymmetric protein structures, while other 1sP approximations gave poorer agreement; this helps to justify the currently empirical practice of estimating B(22,S) from R(h) for large proteins such as antibodies. The results also indicate that either 1bAA or 4bAA CG models may be good starting points for incorporating short-range attractions. Comparison of gD-crystallin B(22) values including both sterics and short-range attractions shows that 1bAA and 4bAA models give equivalent results when properly scaled to account for differences in the number of surface beads in the two CG descriptions. This provides a basis for future work that will also incorporate long-ranged electrostatic attractions and repulsions.
Collapse
Affiliation(s)
- Alexander Grünberger
- Department of Chemical and Biomolecular Engineering, University of Delaware, Newark, Delaware 19716, USA
| | | | | | | |
Collapse
|
18
|
Vekilov PG. Phase diagrams and kinetics of phase transitions in protein solutions. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2012; 24:193101. [PMID: 22495288 DOI: 10.1088/0953-8984/24/19/193101] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
The phase behavior of proteins is of interest for fundamental and practical reasons. The nucleation of new phases is one of the last major unresolved problems of nature. The formation of protein condensed phases (crystals, polymers, and other solid aggregates, as well as dense liquids and gels) underlies pathological conditions, plays a crucial role in the biological function of the respective protein, or is an essential part of laboratory and industrial processes. In this review, we focus on phase transitions of proteins in their properly folded state. We first summarize the recently acquired understanding of physical processes underlying the phase diagrams of the protein solutions and the thermodynamics of protein phase transitions. Then we review recent findings on the kinetics of nucleation of dense liquid droplets and crystals. We explore the transition from nucleation to spinodal decomposition for liquid-liquid separation and introduce the new concept of solution-to-crystal spinodal. We review the two-step mechanism of protein crystal nucleation, in which mesoscopic metastable protein clusters serve as precursors to the ordered crystal nuclei. The concepts and mechanisms reviewed here provide powerful tools for control of the nucleation process by varying the solution thermodynamic parameters.
Collapse
Affiliation(s)
- Peter G Vekilov
- Department of Chemical and Biomolecular Engineering and Department of Chemistry, University of Houston, Houston, TX 77204-4004, USA.
| |
Collapse
|
19
|
Zhang F, Roosen-Runge F, Skoda MWA, Jacobs RMJ, Wolf M, Callow P, Frielinghaus H, Pipich V, Prévost S, Schreiber F. Hydration and interactions in protein solutions containing concentrated electrolytes studied by small-angle scattering. Phys Chem Chem Phys 2012; 14:2483-93. [DOI: 10.1039/c2cp23460b] [Citation(s) in RCA: 75] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
20
|
|
21
|
Priya MH, Pratt LR, Paulaitis ME. Effect of PEG end-group hydrophobicity on lysozyme interactions in solution characterized by light scattering. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2011; 27:13713-13718. [PMID: 21958073 DOI: 10.1021/la203138k] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
We compare protein-protein and protein-polymer osmotic virial coefficients measured by static light scattering for aqueous solutions of lysozyme with low-molecular-weight, hydroxy-terminated (hPEG) and methyl-terminated (mPEG) poly(ethylene glycol) at two solution conditions: pH 7.0 and 0.01 M ionic strength, and pH 6.2 and 0.8 M ionic strength. We find that adding PEG to aqueous lysozyme solutions makes a net repulsive contribution to lysozyme-lysozyme interactions, independent of ionic strength and PEG end-group hydrophobicity. PEG end-group hydrophobicity has a profound effect on the magnitude of this contribution, however, at low ionic strength where mPEG-lysozyme attractive interactions become significant. The enhanced attractions promote mPEG-lysozyme preferential interactions at the expense of lysozyme self-interactions, which leads to lysozyme-lysozyme interactions that are more repulsive in the presence of mPEG. These preferential interactions also lead to the preferential exclusion of diffusable ions locally around the protein, which results in a pronounced ionic strength dependence of mPEG-mediated lysozyme-lysozyme interactions.
Collapse
Affiliation(s)
- M Hamsa Priya
- William G. Lowrie Department of Chemical and Biomolecular Engineering, Ohio State University, Columbus, Ohio 43210, United States
| | | | | |
Collapse
|
22
|
Mason BD, Zhang L, Remmele RL, Zhang J. Opalescence of an IgG2 monoclonal antibody solution as it relates to liquid–liquid phase separation. J Pharm Sci 2011; 100:4587-96. [DOI: 10.1002/jps.22650] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2011] [Revised: 05/08/2011] [Accepted: 05/13/2011] [Indexed: 11/07/2022]
|
23
|
Brubaker WD, Freites JA, Golchert KJ, Shapiro RA, Morikis V, Tobias DJ, Martin RW. Separating instability from aggregation propensity in γS-crystallin variants. Biophys J 2011; 100:498-506. [PMID: 21244846 DOI: 10.1016/j.bpj.2010.12.3691] [Citation(s) in RCA: 60] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2010] [Revised: 12/01/2010] [Accepted: 12/06/2010] [Indexed: 02/02/2023] Open
Abstract
Molecular dynamics (MD) simulations, circular dichroism (CD), and dynamic light scattering (DLS) measurements were used to investigate the aggregation propensity of the eye-lens protein γS-crystallin. The wild-type protein was investigated along with the cataract-related G18V variant and the symmetry-related G106V variant. The MD simulations suggest that local sequence differences result in dramatic differences in dynamics and hydration between these two apparently similar point mutations. This finding is supported by the experimental measurements, which show that although both variants appear to be mostly folded at room temperature, both display increased aggregation propensity. Although the disease-related G18V variant is not the most strongly destabilized, it aggregates more readily than either the wild-type or the G106V variant. These results indicate that γS-crystallin provides an excellent model system for investigating the role of dynamics and hydration in aggregation by locally unfolded proteins.
Collapse
Affiliation(s)
- William D Brubaker
- Department of Molecular Biology and Biochemistry, University of California, Irvine, USA
| | | | | | | | | | | | | |
Collapse
|
24
|
Viscosity Analysis of High Concentration Bovine Serum Albumin Aqueous Solutions. Pharm Res 2011; 28:1973-83. [DOI: 10.1007/s11095-011-0424-7] [Citation(s) in RCA: 88] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2010] [Accepted: 03/08/2011] [Indexed: 11/26/2022]
|
25
|
|
26
|
|
27
|
Garber Cohen IP, Castello PR, González Flecha FL. Ice-induced partial unfolding and aggregation of an integral membrane protein. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2010; 1798:2040-7. [PMID: 20691147 DOI: 10.1016/j.bbamem.2010.07.035] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/10/2010] [Revised: 07/10/2010] [Accepted: 07/28/2010] [Indexed: 11/27/2022]
Abstract
Although the deleterious effects of ice on water-soluble proteins are well established, little is known about the freeze stability of membrane proteins. Here we explore this issue through a combined kinetic and spectroscopic approach using micellar-purified plasma membrane calcium pump as a model. The ATPase activity of this protein significantly diminished after freezing using a slow-cooling procedure, with the decrease in the activity being an exponential function of the storage time at 253K, with t(½)=3.9±0.6h. On the contrary, no significant changes on enzyme activity were detected when a fast cooling procedure was performed. Regardless of the cooling rate, successive freeze-thaw cycles produced an exponential decrease in the Ca(2+)-ATPase activity, with the number of cycles at which the activity was reduced to half being 9.2±0.3 (fast cooling) and 3.7±0.2 (slow cooling). PAGE analysis showed that neither degradation nor formation of SDS-stable aggregates of the protein takes place during protein inactivation. Instead, the inactivation process was found to be associated with the irreversible partial unfolding of the polypeptide chain, as assessed by Trp fluorescence, far UV circular dichroism, and 1-anilino-naphtalene-8-sulfonate binding. This inactive protein undergoes, in a later stage, a further irreversible transformation leading to large aggregates.
Collapse
|
28
|
Utiramerur S, Paulaitis ME. Cooperative hydrophobic/hydrophilic interactions in the hydration of dimethyl ether. J Chem Phys 2010; 132:155102. [DOI: 10.1063/1.3367977] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
29
|
Salinas BA, Sathish HA, Bishop SM, Harn N, Carpenter JF, Randolph TW. Understanding and modulating opalescence and viscosity in a monoclonal antibody formulation. J Pharm Sci 2010; 99:82-93. [PMID: 19475558 PMCID: PMC3927449 DOI: 10.1002/jps.21797] [Citation(s) in RCA: 138] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Opalescence and high viscosities can pose challenges for high concentration formulation of antibodies. Both phenomena result from protein-protein intermolecular interactions that can be modulated with solution ionic strength. We studied a therapeutic monoclonal antibody (mAb) that exhibits high viscosity in solutions at low ionic strength ( approximately 20 cP at 90 mg/mL and 23 degrees C) and significant opalescence at isotonic ionic strength (approximately 100 nephelometric turbidity units at 90 mg/mL and 23 degrees C). The intermolecular interactions responsible for these effects were characterized using membrane osmometry, static light scattering, and zeta potential measurements. The net protein-protein interactions were repulsive at low ionic strength ( approximately 4 mM) and attractive at isotonic ionic strengths. The high viscosities are attributed to electroviscous forces at low ionic strength and the significant opalescence at isotonic ionic strength is correlated with attractive antibody interactions. Furthermore, there appears to be a connection to critical phenomena and it is suggested that the extent of opalescence is dependent on the proximity to the critical point. We demonstrate that by balancing the repulsive and attractive forces via intermediate ionic strengths and by increasing the mAb concentration above the apparent critical concentration both opalescence and viscosity can be simultaneously minimized.
Collapse
Affiliation(s)
- Branden A Salinas
- Department of Chemical and Biological Engineering, University of Colorado, Boulder, Colorado, USA.
| | | | | | | | | | | |
Collapse
|
30
|
Le Brun V, Friess W, Schultz-Fademrecht T, Muehlau S, Garidel P. Lysozyme-lysozyme self-interactions as assessed by the osmotic second virial coefficient: Impact for physical protein stabilization. Biotechnol J 2009; 4:1305-19. [DOI: 10.1002/biot.200800274] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
31
|
Thakur SS, Deepalakshmi P, Gayathri P, Banerjee M, Murthy M, Balaram P. Detection of the protein dimers, multiple monomeric states and hydrated forms of Plasmodium falciparum triosephosphate isomerase in the gas phase. Protein Eng Des Sel 2009; 22:289-304. [DOI: 10.1093/protein/gzp005] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
32
|
Distinguishing thermodynamic and kinetic views of the preferential hydration of protein surfaces. Biophys J 2008; 95:2219-25. [PMID: 18515399 DOI: 10.1529/biophysj.108.133553] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Motivated by a quasi-chemical view of protein hydration, we define specific hydration sites on the surface of globular proteins in terms of the local water density at each site relative to bulk water density. The corresponding kinetic definition invokes the average residence time for a water molecule at each site and the average time that site remains unoccupied. Bound waters are identified by high site occupancies using either definition. In agreement with previous molecular dynamics simulation studies, we find only a weak correlation between local water densities and water residence times for hydration sites on the surface of two globular proteins, lysozyme and staphylococcal nuclease. However, a strong correlation is obtained when both the average residence and vacancy times are appropriately taken into account. In addition, two distinct kinetic regimes are observed for hydration sites with high occupancies: long residence times relative to vacancy times for a single water molecule, and short residence times with high turnover involving multiple water molecules. We also correlate water dynamics, characterized by average occupancy and vacancy times, with local heterogeneities in surface charge and surface roughness, and show that both features are necessary to obtain sites corresponding to kinetically bound waters.
Collapse
|
33
|
Protein-protein interaction on lysozyme crystallization revealed by rotational diffusion analysis. Biophys J 2008; 94:4484-92. [PMID: 18310245 DOI: 10.1529/biophysj.107.111872] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Intermolecular interactions between protein molecules diffusing in various environments underlie many biological processes as well as control protein crystallization, which is a crucial step in x-ray protein structure determinations. Protein interactions were investigated through protein rotational diffusion analysis. First, it was confirmed that tetragonal lysozyme crystals containing fluorescein-tagged lysozyme were successfully formed with the same morphology as that of native protein. Using this nondisruptive fluorescent tracer system, we characterized the effects of sodium chloride and ammonium sulfate concentrations on lysozyme-lysozyme interactions by steady-state and time-resolved fluorescence anisotropy measurements and the introduction of a novel interaction parameter, k(rot). The results suggested that the specific attractive interaction, which was reflected in the retardation of the protein rotational diffusion, was induced depending on the salt type and its concentration. The change in the attractive interactions also correlated with the crystallization/precipitation behavior of lysozyme. Moreover, we discuss the validity of our rotational diffusion analysis through comparison with the osmotic second virial coefficient, B(22), previously reported for lysozyme and those estimated from k(rot).
Collapse
|
34
|
Abstract
Normal mode analysis (NMA) has received much attention as a direct approach to extract the collective motions of macromolecules. However, the stringent requirement of computational resources by classical all-atom NMA limits the size of the macromolecules to which the method is normally applied. We implemented a novel coarse-grained normal mode approach based on partitioning the all-atom Hessian matrix into relevant and nonrelevant parts. It is interesting to note that, using classical all-atom NMA results as a reference, we found that this method generates more accurate results than do other coarse-grained approaches, including elastic network model and block normal mode approaches. Moreover, this new method is effective in incorporating the energetic contributions from the nonrelevant atoms, including surface water molecules, into the coarse-grained protein motions. The importance of such improvements is demonstrated by the effect of surface water to shift vibrational modes to higher frequencies and by an increase in overlap of the coarse-grained eigenvector space (the motion directions) with that obtained from molecular dynamics simulations of solvated protein in a water box. These results not only confirm the quality of our method but also point out the importance of incorporating surface structural water in studying protein dynamics.
Collapse
|
35
|
Reichmann D, Phillip Y, Carmi A, Schreiber G. On the Contribution of Water-Mediated Interactions to Protein-Complex Stability. Biochemistry 2007; 47:1051-60. [PMID: 18161993 DOI: 10.1021/bi7019639] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Dana Reichmann
- Department of Biological Chemistry, Weizmann Institute of Science, Rehovot, 76100, Israel
| | - Yael Phillip
- Department of Biological Chemistry, Weizmann Institute of Science, Rehovot, 76100, Israel
| | - Asaf Carmi
- Department of Biological Chemistry, Weizmann Institute of Science, Rehovot, 76100, Israel
| | - Gideon Schreiber
- Department of Biological Chemistry, Weizmann Institute of Science, Rehovot, 76100, Israel
| |
Collapse
|
36
|
Kim YC, Hummer G. Coarse-grained models for simulations of multiprotein complexes: application to ubiquitin binding. J Mol Biol 2007; 375:1416-33. [PMID: 18083189 DOI: 10.1016/j.jmb.2007.11.063] [Citation(s) in RCA: 207] [Impact Index Per Article: 12.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2007] [Revised: 11/19/2007] [Accepted: 11/19/2007] [Indexed: 10/22/2022]
Abstract
We develop coarse-grained models and effective energy functions for simulating thermodynamic and structural properties of multiprotein complexes with relatively low binding affinity (K(d) >1 microM) and apply them to binding of Vps27 to membrane-tethered ubiquitin. Folded protein domains are represented as rigid bodies. The interactions between the domains are treated at the residue level with amino-acid-dependent pair potentials and Debye-Hückel-type electrostatic interactions. Flexible linker peptides connecting rigid protein domains are represented as amino acid beads on a polymer with appropriate stretching, bending, and torsion-angle potentials. In simulations of membrane-attached protein complexes, interactions between amino acids and the membrane are described by residue-dependent short-range potentials and long-range electrostatics. We parameterize the energy functions by fitting the osmotic second virial coefficient of lysozyme and the binding affinity of the ubiquitin-CUE complex. For validation, extensive replica-exchange Monte Carlo simulations are performed of various protein complexes. Binding affinities for these complexes are in good agreement with the experimental data. The simulated structures are clustered on the basis of distance matrices between two proteins and ranked according to cluster population. In approximately 70% of the complexes, the distance root-mean-square is less than 5 A from the experimental structures. In approximately 90% of the complexes, the binding interfaces on both proteins are predicted correctly, and in all other cases at least one interface is correct. Transient and nonspecifically bound structures are also observed. With the validated model, we simulate the interaction between the Vps27 multiprotein complex and a membrane-tethered ubiquitin. Ubiquitin is found to bind preferentially to the two UIM domains of Vps27, but transient interactions between ubiquitin and the VHS and FYVE domains are observed as well. These specific and nonspecific interactions are found to be positively cooperative, resulting in a substantial enhancement of the overall binding affinity beyond the approximately 300 microM of the specific domains. We also find that the interactions between ubiquitin and Vps27 are highly dynamic, with conformational rearrangements enabling binding of Vps27 to diverse targets as part of the multivesicular-body protein-sorting pathway.
Collapse
Affiliation(s)
- Young C Kim
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892-0520, USA
| | | |
Collapse
|
37
|
Dumetz AC, Snellinger-O'brien AM, Kaler EW, Lenhoff AM. Patterns of protein protein interactions in salt solutions and implications for protein crystallization. Protein Sci 2007; 16:1867-77. [PMID: 17766383 PMCID: PMC2206983 DOI: 10.1110/ps.072957907] [Citation(s) in RCA: 119] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2007] [Revised: 06/19/2007] [Accepted: 06/20/2007] [Indexed: 10/22/2022]
Abstract
The second osmotic virial coefficients of seven proteins-ovalbumin, ribonuclease A, bovine serum albumin, alpha-lactalbumin, myoglobin, cytochrome c, and catalase-were measured in salt solutions. Comparison of the interaction trends in terms of the dimensionless second virial coefficient b(2) shows that, at low salt concentrations, protein-protein interactions can be either attractive or repulsive, possibly due to the anisotropy of the protein charge distribution. At high salt concentrations, the behavior depends on the salt: In sodium chloride, protein interactions generally show little salt dependence up to very high salt concentrations, whereas in ammonium sulfate, proteins show a sharp drop in b(2) with increasing salt concentration beyond a particular threshold. The experimental phase behavior of the proteins corroborates these observations in that precipitation always follows the drop in b(2). When the proteins crystallize, they do so at slightly lower salt concentrations than seen for precipitation. The b(2) measurements were extended to other salts for ovalbumin and catalase. The trends follow the Hofmeister series, and the effect of the salt can be interpreted as a water-mediated effect between the protein and salt molecules. The b(2) trends quantify protein-protein interactions and provide some understanding of the corresponding phase behavior. The results explain both why ammonium sulfate is among the best crystallization agents, as well as some of the difficulties that can be encountered in protein crystallization.
Collapse
Affiliation(s)
- André C Dumetz
- Center for Molecular and Engineering Thermodynamics, Department of Chemical Engineering, University of Delaware, Newark, Delaware 19716, USA
| | | | | | | |
Collapse
|
38
|
McGuffee SR, Elcock AH. Atomically detailed simulations of concentrated protein solutions: the effects of salt, pH, point mutations, and protein concentration in simulations of 1000-molecule systems. J Am Chem Soc 2007; 128:12098-110. [PMID: 16967959 DOI: 10.1021/ja0614058] [Citation(s) in RCA: 89] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
An ability to accurately simulate the dynamic behavior of concentrated macromolecular solutions would be of considerable utility in studies of a wide range of biological systems. With this goal in mind, a Brownian dynamics (BD) simulation method is reported here that allows systems to be modeled that comprise in excess of 1000 protein molecules, all of which are treated in atomic detail. Intermolecular forces are described in the method using an energy function that incorporates electrostatic and hydrophobic interactions and that is calibrated to reproduce experimental thermodynamic information with a single adjustable parameter. Using the method, BD simulations have been performed over a wide range of pH and ionic strengths for three proteins: hen egg white lysozyme (HEWL), chymotrypsinogen, and T4 lysozyme. The simulations reproduce experimental trends in second virial coefficients (B(22)) and translational diffusion coefficients, correctly capture changes in B(22) values due to single amino acid substitutions, and reveal a new explanation for the difficulties reported previously in the literature in reproducing B(22) values for protein solutions of very low ionic strength. In addition, a strong correlation is found between a residue's probability of being involved in a protein-protein contact in the simulations and its probability of being involved in an experimental crystal contact. Finally, exploratory simulations of HEWL indicate that the simulation model also gives a promising description of behavior at very high protein concentrations (approximately 250 g/L), suggesting that it may provide a suitable computational framework for modeling the complex behavior exhibited by macromolecules in cellular conditions.
Collapse
Affiliation(s)
- Sean R McGuffee
- Department of Biochemistry, University of Iowa, Iowa City, Iowa 52242, USA
| | | |
Collapse
|
39
|
Reichmann D, Rahat O, Cohen M, Neuvirth H, Schreiber G. The molecular architecture of protein-protein binding sites. Curr Opin Struct Biol 2007; 17:67-76. [PMID: 17239579 DOI: 10.1016/j.sbi.2007.01.004] [Citation(s) in RCA: 147] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2006] [Revised: 12/13/2006] [Accepted: 01/10/2007] [Indexed: 11/16/2022]
Abstract
The formation of specific protein interactions plays a crucial role in most, if not all, biological processes, including signal transduction, cell regulation, the immune response and others. Recent advances in our understanding of the molecular architecture of protein-protein binding sites, which facilitates such diversity in binding affinity and specificity, are enabling us to address key questions. What is the amino acid composition of binding sites? What are interface hotspots? How are binding sites organized? What are the differences between tight and weak interacting complexes? How does water contribute to binding? Can the knowledge gained be translated into protein design? And does a universal code for binding exist, or is it the architecture and chemistry of the interface that enable diverse but specific binding solutions?
Collapse
Affiliation(s)
- Dana Reichmann
- Department of Biological Chemistry, Weizmann Institute of Science, Rehovot 76100, Israel
| | | | | | | | | |
Collapse
|
40
|
Trilisky EI, Lenhoff AM. Sorption processes in ion-exchange chromatography of viruses. J Chromatogr A 2007; 1142:2-12. [PMID: 17240385 DOI: 10.1016/j.chroma.2006.12.094] [Citation(s) in RCA: 68] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2006] [Revised: 12/19/2006] [Accepted: 12/28/2006] [Indexed: 11/19/2022]
Abstract
Purified viruses are used in gene therapy and vaccine production. Ion-exchange chromatography (IEC) is the most common method for large-scale downstream purification of viruses and proteins. Published IEC protocols provide details for specific separations but not general methods for selecting operating parameters. To make the selection more systematic, we study adenovirus type 5 (Ad5) as a model virus and develop batch uptake and light scattering methods for optimizing the ionic strength and pH of adsorption, as well as providing heuristics for resin geometry. The static capacity for Ad5 was found to go through a maximum with increasing ionic strength. Comparison to a protein-resin system shows that resin capacity for the virus is at least an order of magnitude lower, even on a wide-pore resin. Virus penetration into the wide-pore resin is only partial and the uptake rate is an order of magnitude slower than the uptake onto a narrow-pore resin.
Collapse
Affiliation(s)
- E I Trilisky
- Department of Chemical Engineering, University of Delaware, Newark, DE 19716, USA
| | | |
Collapse
|
41
|
Roy S, Kar S, Chaudhuri S, Dasgupta AK. Potential of cadmium sulphide nanorods as an optical microscopic probe to the folding state of cytochrome C. Biophys Chem 2006; 124:52-61. [PMID: 16814922 DOI: 10.1016/j.bpc.2006.05.011] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2006] [Revised: 05/12/2006] [Accepted: 05/12/2006] [Indexed: 11/22/2022]
Abstract
The folding behavior of cytochrome C (Cyt-C) conjugated with CdS nanorods (CdSnr) is amenable to monitoring by bright field microscopy, the porosity and percolating behavior of such protein conjugated nanoclusters depending on the folding history prior to the conjugation. The method has been used to predict the thermal melting behavior as well as guanidine hydrochloride induced unfolding of Cyt-C. Dynamic light scattering studies indicate that the size distribution of the nanoforms widens in presence of the protein. Furthermore, there is emergence of clusters with higher conductivity and altered zeta potential. Increase of second virial coefficient of CdS nanoforms in the presence of Cyt-C (obtained from static light scattering experiments) implies presence of protein coat over the hydrophobic nanosurface. The results are supported by morphological changes observed through scanning electron microscopy (SEM). Accordingly, the X-ray diffraction pattern shows a change of crystallographic orientations of CdSnr in presence of Cyt-C.
Collapse
Affiliation(s)
- Shibsekhar Roy
- Department of Biochemistry, University of Calcutta, 35-Ballygunge Circular Road, Kolkata-700019, India
| | | | | | | |
Collapse
|
42
|
Paliwal A, Asthagiri D, Pratt LR, Ashbaugh HS, Paulaitis ME. An analysis of molecular packing and chemical association in liquid water using quasichemical theory. J Chem Phys 2006; 124:224502. [PMID: 16784293 DOI: 10.1063/1.2202350] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We calculate the hydration free energy of liquid TIP3P water at 298 K and 1 bar using a quasi-chemical theory framework in which interactions between a distinguished water molecule and the surrounding water molecules are partitioned into chemical associations with proximal (inner-shell) waters and classical electrostatic-dispersion interactions with the remaining (outer-shell) waters. The calculated free energy is found to be independent of this partitioning, as expected, and in excellent agreement with values derived from the literature. An analysis of the spatial distribution of inner-shell water molecules as a function of the inner-shell volume reveals that water molecules are preferentially excluded from the interior of large volumes as the occupancy number decreases. The driving force for water exclusion is formulated in terms of a free energy for rearranging inner-shell water molecules under the influence of the field exerted by outer-shell waters in order to accommodate one water molecule at the center. The results indicate a balance between chemical association and molecular packing in liquid water that becomes increasingly important as the inner-shell volume grows in size.
Collapse
Affiliation(s)
- A Paliwal
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | | | | | | | | |
Collapse
|