Hydrogen-bonded Network Of Water Molecules Research Articles

Ni(II) triggered sensing and reversible uptake of ammonia by tetrazole based Ni(II) metallogel: A one-pot mechanistic approach

The release of ammonia gas and its removal through traditional approaches is a critical challenge inhibiting the sustainability of chemical industries. The quest to improve sustainability has triggered interest to develop chemoprobe and reusable adsorbent for toxic gases such as ammonia. The present study offers the first example of a novel one-pot synthesis of a stable Ni(II)metallogel (Ni(II)MG) from 1, 4 adiponitrile and sodium azide. The Ni(II)MG is composed of Ni(II) coordinated with four in-situ formed CN-Tz (5-(1H-tetrazol-5-yl)pentanenitrile) ligands and two water molecules forming predominant octahedral geometry. The Ni(II)MG has good strength, and viscoelastic, thixotropic, and self-healing properties. The Ni(II)MG and its xerogel (Ni(II)XG) have been thoroughly characterized and utilized for ammonia sensing and reversible adsorption at ambient condition. The mesoporous nature and presence of hydrogen-bonded water molecules in the pores make this material effective for reversible adsorption of ammonia. The lewis acidic Ni(II) increases the cloud network of hydrogen-bonded water molecules and triggers the coordinated water molecules to activate hydrogen-bonded water molecules for binding the ammonia molecules. The Ni(II) triggered extended hydrogen bonding through coordinated water molecule is also supported by DFT studies. The Ni(II)XG shows good ammonia adsorption capacity (8 mmol/gm) which is confirmed by TPD studies. The properties of Ni(II)XG are found to be intact after ammonia adsorption even after 20 cycles. Hence, it could be a promising low cost, sustainable and reusable material as a chemoprobe for removal of ammonia.

The Signature of Fluctuations of the Hydrogen Bond Network Formed by Water Molecules in the Interfacial Layer of Anionic Lipids

As the water molecules found at the interface of lipid bilayers exhibit distinct structural and reorientation dynamics compared to water molecules found in bulk, the fluctuations in their hydrogen bond (HB) network are expected to be different from those generated by the bulk water molecules. The research presented here aims to gain an insight into temperature-dependent fluctuations of a HB network of water molecules found in an interfacial layer of multilamellar liposomes (MLVs) composed of anionic 1,2-dimyristoyl-sn-glycero-3-phospho-L-serine (DMPS) lipids. Besides suspending DMPS lipids in phosphate buffer saline (PBS) of different pH values (6.0, 7.4, and 8.0), the changes in HB network fluctuations were altered by the incorporation of a non-polar flavonoid molecule myricetin (MCE) within the hydrocarbon chain region. By performing a multivariate analysis on the water combination band observed in temperature-dependent FTIR spectra, the results of which were further mathematically analyzed, the temperature-dependent fluctuations of interfacial water molecules were captured; the latter were the greatest for DMPS in PBS with a pH value of 7.4 and in general were greater for DMPS multibilayers in the absence of MCE. The presence of MCE made DMPS lipids more separated, allowing deeper penetration of water molecules towards the non-polar region and their restricted motion that resulted in decreased fluctuations. The experimentally observed results were supported by MD simulations of DMPS (+MCE) lipid bilayers.

Crystal structures and properties of Fe(H2O)2[Ni(CN)4]·nH2O hydrates

We report the crystal structures and properties of iron hydrate compounds. These hydrates have 2D topologies and crystallize with three different structural features depending on the number of weakly bonded water molecules. Coordinated water molecules occupy the axial sites of the iron ion and support a network of hydrogen-bonded water molecules in the interlayer region. Crystal structures of these hydrates are known for Mn, Co, and Ni but not for Fe. The structure for the solids considered herein was solved and refined from powder XRD patterns complemented by IR, UV/Vis-NIR, Mössbauer, and TG data. The relevance of this structural study is related to the application of Hofmann-like coordination polymers as spin-crossover materials when the coordinated water molecules are replaced by dipyridines, pyridine, and derivatives that form pillars between adjacent layers. Upon dehydration of these solids, a reactive material is formed by assembly of tetrahedral iron ions with square planar nickel ions.

Heme-copper and Heme O2-derived synthetic (bioinorganic) chemistry toward an understanding of cytochrome c oxidase dioxygen chemistry

Cytochrome c oxidase (CcO), also widely known as mitochondrial electron-transport-chain complex IV, is a multi-subunit transmembrane protein responsible for catalyzing the last step of the electron transport chain, dioxygen reduction to water, which is essential to the establishment and maintenance of the membrane proton gradient that drives ATP synthesis. Although many intermediates in the CcO catalytic cycle have been spectroscopically and/or computationally authenticated, the specifics regarding the IP intermediate, hypothesized to be a heme-Cu (hydro)peroxo species whose O−O bond homolysis is supported by a hydrogen-bonding network of water molecules, are largely obscured by the fast kinetics of the A (FeIII–O2•–/CuI/Tyr) → PM (FeIV=O/CuII–OH/Tyr•) step. In this review, we have focused on the recent advancements in the design, development, and characterization of synthetic heme-peroxo‑copper model complexes, which can circumvent the abovementioned limitation, for the investigation of the formation of IP and its O−O cleavage chemistry. Novel findings regarding (a) proton and electron transfer (PT/ET) processes, together with their contributions to exogenous phenol induced O−O cleavage, (b) the stereo-electronic tunability of the secondary coordination sphere (especially hydrogen-bonding) on the geometric and spin state alteration of the heme-peroxo‑copper unit, and (c) a plausible mechanism for the Tyr-His cofactor biogenesis, are discussed in great detail. Additionally, since the ferric-superoxide and the ferryl-oxo (Compound II) species are critically involved in the CcO catalytic cycle, this review also highlights a few fundamental aspects of these heme-only (i.e., without copper) species, including the structural and reactivity influences of electron-donating trans-axial ligands and Lewis acid-promoted H-bonding.

Variations in the Protein Hydration and Hydrogen-Bond Network of Water Molecules Induced by the Changes in the Secondary Structures of Proteins Studied through Near-Infrared Spectroscopy.

This study investigated how the secondary structural changes of proteins in aqueous solutions affect their hydration and the hydrogen-bond network of water molecules using near-infrared (NIR) spectroscopy. The aqueous solutions of three types of proteins, i.e., ovalbumin, β-lactoglobulin, and bovine serum albumin, were denatured by heating, and changes in the NIR bands of water reflecting the states of hydrogen bonds induced via protein secondary structural changes were investigated. On heating, the intermolecular hydrogen bonds between water molecules as well as between water and protein molecules were broken, and protein molecules were no longer strongly bound by the surrounding water molecules. Consequently, the denaturation was observed to proceed depending on the thermodynamic properties of the proteins. When the aqueous solutions of proteins were cooled after denaturation, the hydrogen-bond network was reformed. However, the state of protein hydration was changed owing to the secondary structural changes of proteins, and the variation patterns were different depending on the protein species. These changes in protein hydration may be derived from the differences in the surface charges of proteins. The elucidation of the mechanism of protein hydration and the formation of the hydrogen-bond network of water molecules will afford a comprehensive understanding of the protein functioning and dysfunctioning derived from the structural changes in proteins.

Enhancing the Performance of Triboelectric Generator: A Novel Approach Using Solid-Liquid Interface-Treated Foam and Metal Contacts.

This work introduces a novel approach for enhancing the performance of a triboelectric generator (TEG) by using a solid-liquid interface-treated foam (SLITF) as its active layer, combined with two metal contacts of different work functions. SLITF is made by absorbing water into a cellulose foam, which enables charges generated by friction energy during the sliding motion to be separated and transferred through the conductive path formed by the hydrogen-bonded network of water molecules. Unlike traditional TEGs, the SLITF-TEG demonstrates an impressive current density of 3.57 A/m2 and can harvest electric power up to 0.174 W/m2 with an induced voltage of approximately 0.55 V. The device generates a direct current in the external circuit, eliminating the limitations of low current density and alternating current found in traditional TEGs. By connecting six-unit cells of SLITF-TEG in series and parallel, the peak voltage and current can be increased up to 3.2 V and 12.5 mA, respectively. Furthermore, the SLITF-TEG has the potential to serve as a self-powered vibration sensor with high accuracy (R2 = 0.99). The findings demonstrate the significant potential of the SLITF-TEG approach for efficiently harvesting low-frequency mechanical energy from the natural environment, with broad implications for a range of applications.

Microenvironment regulation of Ru(bda)L2 catalyst incorporated in metal-organic framework for effective photo-driven water oxidation

A Ru(bda)L2 molecular catalyst (H2bda = 2,2′-bipyridine-6,6′-dicarboxylic acid, L represents ligands) was incorporated into a UiO-66 metal-organic framework (MOF) as isolating sites for chemical- and photo-driven water oxidation with the water nucleophilic attack mechanism. An impressive turn-over number of 566 was obtained for photo-driven water oxidation in the phosphate buffer saline, which was nearly 20 times that of the homogenous counterpart in the presence of ruthenium tris(bipyridine) and sodium persulfate. Infrared spectroscopy, X-ray diffraction and X-ray spectroscopy techniques were used to characterize the composition and structure of the catalysts. Kinetic isotope effect (KIE), proton inventory experiment and proton nuclear magnetic resonance were conducted to understand the photo-driven O2 evolution mechanism. It was observed that the protons participated in O2 evolution with a normal KIE value, and the uptake of phosphates by UiO-66 matrix improved the affinity towards H2O via hydrogen bonding. The high O2 output and the limited photosensitizer oxidative decomposition in PBS suggest that the phosphate acted as a proton mediator to assist proton and/or proton-couple electron transfer through a hydrogen-bonded network of water molecules. This study uses a heterogenous matrix to emulate an enzyme-like microenvironment capable of achieving efficient artificial photosynthesis.

Probing collective terahertz vibrations of a hydrogen-bonded water network at buried electrochemical interfaces.

The exceptional properties of liquid water such as thermodynamic, physical, and dielectric anomalies originate mostly from the hydrogen-bond networks of water molecules. The structural and dynamic properties of the hydrogen-bond networks have a significant impact on many biological and chemical processes in aqueous systems. In particular, the properties of interfacial water molecules with termination of the network at a solid surface are crucial to understanding the role of water in heterogeneous reactions. However, direct monitoring of the dynamics of hydrogen-bonded interfacial water molecules has been limited because of the lack of a suitable surface-selective spectroscopic means in the terahertz (THz) frequency range where collective vibrations of water exist. Here we show that hydrogen-bond vibrations below 9 THz can be measured in situ at an electrochemical interface, which is buried between two THz-opaque media, by using a density of states format of surface-enhanced inelastic light scattering spectra. The interpretation of the obtained spectra over the range 0.3-6 THz indicates that the negatively charged surface accelerates collective translational motions of water molecules in the lateral direction with the increase of hydrogen-bond defects. Alternatively, the positively charged surface results in suppression of lateral mobility. This work gives a new perspective on in situ spectroscopic investigations in heterogeneous reactions.

Investigated hydrogen-bond network kinetics of acetone-water solutions by spontaneous and stimulated Raman spectroscopy.

The hydrogen bond (HB) network structure and kinetics of the acetone-water mixed solutions were investigated by the spontaneous Raman and stimulated Raman scattering (SRS) spectra. The HB network of water molecules was enhanced when the volume fraction of acetone ranged from 0 to 0.25. Two new SRS peaks of water at 3272 and 3380 cm-1 were obtained, resulting from the cooperation of the polar carbonyl (C = O)-enhanced HB and the ice-like structure formed around the methyl groups. However, when the volume fraction went beyond 0.25, the spontaneous Raman main peak at 3445 cm-1 showed a significant blue-shift, and the corresponding SRS signal disappeared, indicating that the HB of water was weakened, which originated from the self-association of acetone. In the meantime, the fully tetrahedral HB structure among water molecules was destroyed at the higher volume fraction (≥ 0.8). Hopefully, our study here would advance the study of HB network structures and kinetics in other aqueous solutions.

Structural basis of leukotriene B4 receptor 1 activation

Leukotriene B4 receptor 1 (BLT1) plays crucial roles in the acute inflammatory responses and is a valuable target for anti-inflammation treatment, however, the mechanism by which leukotriene B4 (LTB4) activates receptor remains unclear. Here, we report the cryo-electron microscopy (cryo-EM) structure of the LTB4 -bound human BLT1 in complex with a Gi protein in an active conformation at resolution of 2.91 Å. In combination of molecule dynamics (MD) simulation, docking and site-directed mutagenesis, our structure reveals that a hydrogen-bond network of water molecules and key polar residues is the key molecular determinant for LTB4 binding. We also find that the displacement of residues M1013.36 and I2717.39 to the center of receptor, which unlock the ion lock of the lower part of pocket, is the key mechanism of receptor activation. In addition, we reveal a binding site of phosphatidylinositol (PI) and discover that the widely open ligand binding pocket may contribute the lack of specificity and efficacy for current BLT1-targeting drug design. Taken together, our structural analysis provides a scaffold for understanding BLT1 activation and a rational basis for designing anti-leukotriene drugs.

Structure of an aqueous RbCl solution in the gigapascal pressure range by neutron diffraction combined with empirical potential structure refinement modeling

An aqueous 3 m (=mol/kg) RbCl solution in D2O is measured at 298 K/0.1 MPa, 298 K/1 GPa, 523 K/1 GPa, and 523 K/4 GPa by neutron diffraction. The interference functions are analyzed by an empirical potential structure refinement (EPSR) modeling. The ion hydration and association and hydrogen-bonded water structure as a function of temperature and pressure are revealed in the pair correlation function, coordination number, triangular distribution, and spatial density function (3D structure). The results show that structure parameters are temperature and pressure-dependent. The second hydration layer of Rb+ and Cl− becomes evident under elevated pressure and temperature conditions. The average oxygen coordination number of Rb+ in the first hydration layer increases from 6.3 ± 1.7 at ambient pressure to 8.9 ± 2.3 at 4 GPa with decreasing coordination distance of 0.290 nm to 0.288 nm. The average oxygen coordination number of Cl− in the first hydration shell increases from 5.9 ± 1.5 at ambient pressure to 9.1 ± 2.1 at 4 GPa, and the corresponding coordination distance decreases from 0.322 nm to 0.314 nm with increasing pressure. The orientation of the water dipole in the first solvation shell of Rb+ and a central water molecule is sensitive to pressure, but that in the first solvation shell of Cl- does not change very much regardless of the thermodynamic condition. The number of contact-ion pairs Rb+-Cl− decreases with elevated temperature and increases with elevated pressure. Water molecules are closely packed, and the tetrahedral hydrogen-bonded network of water molecules no longer exists in extreme conditions.

Soft X-ray Emission Studies on Hydrate-Melt Electrolytes.

Highly salt-concentrated aqueous solutions are a new class of electrolytes, which provide a wide potential window exceeding 3 V and, hence, realize possibly inexpensive, safe, and high-energy-density storage devices. Herein, we investigate the evolution of the coordination structure and electronic state depending on the salt concentration through soft X-ray emission spectroscopy and first-principles molecular dynamics calculations. Close to the concentration limit, categorized as a "hydrate melt," a long-range hydrogen-bond network of water molecules disappears with emerging localized electronic states that resemble those in the gas phase. Such localized electronic states are attributed not only to their geometrically isolated nature but also to their dominant electrostatic interaction with Li+ cations. Therefore, the electrical properties of water in the hydrate melt can be more gaslike than liquidlike.

Rotational Dynamics of Water at the Phospholipid Bilayer Depending on the Head Groups Studied by Molecular Dynamics Simulations

Hydration states are a crucial factor that affect the self-assembly and properties of soft materials and biomolecules. Although previous experiments have revealed that the hydration state strongly depends on the chemical structure of lipid molecules, the mechanisms at the molecular level remain unknown. Classical and density-functional tight-binding (DFTB) molecular dynamics (MD) simulations are employed to determine the mechanisms underlying dissimilar water dynamics between lipid membranes with phosphatidylcholine (PC) and phosphatidylethanolamine (PE) head groups. The classical MD simulation shows that rotational relaxations of water are faster on the PE lipid than on the PC lipid, which is consistent with a previous experimental study using terahertz spectroscopy. Furthermore, DFTB-MD simulation of N(CH3)4+ and NH4+ ions, which correspond to the different head groups in PC and PE, respectively, shows qualitative agreement with the classical MD simulation. Remarkably, the PE lipids and the NH4+ ions break hydrogen bonds between water molecules in the secondary hydration shell. In contrast, the PC lipids and the N(CH3)4+ ions bind water molecules weakly in both the primary and secondary hydration shells and increase hydrogen bonds between water. Our atomistic simulations show that these changes in the hydrogen-bond network of water molecules cause the different rotational relaxation of water molecules between the two lipids.

Potential-Induced Adsorption and Structuring of Water at the Pt(111) Electrode Surface in Contact with an Ionic Liquid

Water adsorption is important in many fields from surface electrochemistry to electrocatalysis where molecular level information is much needed in order to gain a detailed understand on the role of interfacial water. Here we report on water at Pt(111) surfaces in contact with an [EIMIM][BF4] ionic liquid, which was spectroscopically resolved by using in situ sum-frequency generation (SFG). O-H modes are used to study water adsorption and water structure as a function of electrode potential, while analysis of C-H modes is used to infer orientational changes of [EMIM] cations at the interface. Different from the bulk where free water molecules are found, SFG spectra provide evidence an interfacial layer with an extended network of hydrogen-bonded water molecules exists and grows with increasing absolute potential which is used to identify the potential of zero charge at +0.1 V SHE where a pronounced minimum in O-H intensity is found.

Mechanism of Inward Proton Transport in an Antarctic Microbial Rhodopsin.

Although the outward-directed proton transport across biological membranes is well studied and its importance for bioenergetics is clearly understood, inward-directed light-driven proton pumping by microbial rhodopsins has remained a mystery both physiologically and mechanistically. A new family of Antarctic rhodopsins, which is a subgroup within a novel class of schizorhodopsins reported recently, includes a member, denoted as AntR, which proved amenable to extensive characterization with experiments and computation. Phylogenetic analyses identify AntR as distinct from the well-studied microbial rhodopsins that function as outward-directed ion pumps, and bioinformatics sequence analyses reveal amino acid substitutions at conserved sites essential for outward proton pumping. Modeling and numerical simulations of AntR, combined with advanced analyses using the graph theory and centrality measures from social sciences, identify the dynamic three-dimensional network of hydrogen-bonded water molecules and amino acid residues that function as communication hubs in AntR. This network undergoes major rearrangement upon retinal isomerization, showing important changes in the connectivity of the active center, retinal Schiff base, to the opposing sides of the membrane, as required for proton transport. Numerical simulations and experimental studies of the photochemical cycle of AntR by spectroscopy and site-directed mutagenesis allowed us to identify pathways that could conduct protons in the direction opposite to that commonly known for outward-directed pumps.

Synthesis and solid state self-assembly of a 1,4-diazepine derivative: Water cluster as molecular glue and conformational isomerism

Diazepine (a seven membered heterocycle) derivatives have emerged as “privileged structures” in chemical biology and medicinal chemistry. In this article, synthesis of a 1, 4-diazepine derivative 4 from flavone 3 in excellent yield is described. The solid state self-assembly pertaining to diazepine derivative 4 was investigated with X-ray crystallography. Fascinatingly, the solid state structure of 4 represents a remarkable example of water cluster acting as structural glue to instigate conformational isomerism and hence facilitate the crystallization process by compensating the imbalance of hydrogen-bond donors and acceptors. NBO and QTAIM based calculations show that the hydrogen-bonding network of water molecules and hyper conjugative interactions play a vital role for stabilization of structure in the solid state. To our knowledge, presented here is the first illustration of any solvent cluster stabilizing the solid-state structure in a diazepine family. In addition, the crystal structure of flavone 3 hitherto not reported is also reported.

Unraveling a New Chemical Mechanism of Missing Sulfate Formation in Aerosol Haze: Gaseous NO2 with Aqueous HSO3-/SO32.

A source of missing sulfate production associated with high-level fine-particle pollution in the megacities of China is believed to stem from the oxidation of a notable fraction of sulfur dioxide (SO2) by nitrogen dioxide (NO2) in aqueous aerosol environments, suggesting that an unknown reaction pathway exists for aqueous sulfur oxidation. At weakly acidic aerosols, the dissolved SO2 mainly exists in the form of HSO3-, whereas at neutral aerosols, SO32- becomes the main form. Herein, by using both ab initio molecular metadynamics simulations and high-level quantum mechanical calculations, we show a hitherto unreported chemical mechanism for the formation of sulfate through the reaction between HSO3-/SO32- anions at the surface/in the interior of a water nanodroplet and gas-phase NO2 molecules. For weakly acidic aerosols, contrary to the conventional high-barrier electron-transfer pathway in the gas phase, HSO3- at the water nanodroplet surface can transfer an electron to NO2 with a low free-energy barrier of 4.7 kcal/mol through a water bridge. For neutral aerosols, the electron-transfer pathway between SO32- in the interior of the water nanodroplet and NO2 needs to overcome a lower free-energy barrier of 3.6 kcal/mol to form SO3-, with the assistance of the hydrogen-bonding network of water molecules. This new reaction pathway for the sulfate formation from HSO3-/SO32- via water nanodroplets and gaseous NO2 provides a new perspective on the growth of haze particles from pre-existing aqueous aerosols and suggests that new control strategies are needed to address haze pollution.

Water Dynamics in Polyacrylamide Hydrogels.

Polymeric hydrogels have wide applications including electrophoresis, biocompatible materials, water superadsorbents, and contact lenses. The properties of hydrogels involve the poorly characterized molecular dynamics of water and solutes trapped within the three-dimensional cross-linked polymer networks. Here we apply ultrafast two-dimensional infrared (2D IR) vibrational echo and polarization-selective pump-probe (PSPP) spectroscopies to investigate the ultrafast molecular dynamics of water and a small molecular anion solute, selenocyanate (SeCN-), in polyacrylamide hydrogels. For all mass concentrations of polymer studied (5% and above), the hydrogen-bonding network reorganization (spectral diffusion) dynamics and reorientation dynamics reported by both water and SeCN- solvated by water are significantly slower than in bulk water. As the polymer mass concentration increases, molecular dynamics in the hydrogels slow further. The magnitudes of the slowing, measured with both water and SeCN-, are similar. However, the entire hydrogen-bonding network of water molecules appears to slow down as a single ensemble, without a difference between the core water population and the interface water population at the polymer-water surface. In contrast, the dissolved SeCN- do exhibit two-component dynamics, where the major component is assigned to the anions fully solvated in the confined water nanopools. The slower component has a small amplitude which is correlated with the polymer mass concentration and is assigned to adsorbed anions strongly interacting with the polymer fiber networks.

Substitution of the D1-Asn87 site in photosystem II of cyanobacteria mimics the chloride-binding characteristics of spinach photosystem II

Photoinduced water oxidation at the O2-evolving complex (OEC) of photosystem II (PSII) is a complex process involving a tetramanganese-calcium cluster that is surrounded by a hydrogen-bonded network of water molecules, chloride ions, and amino acid residues. Although the structure of the OEC has remained conserved over eons of evolution, significant differences in the chloride-binding characteristics exist between cyanobacteria and higher plants. An analysis of amino acid residues in and around the OEC has identified residue 87 in the D1 subunit as the only significant difference between PSII in cyanobacteria and higher plants. We substituted the D1-Asn87 residue in the cyanobacterium Synechocystis sp. PCC 6803 (wildtype) with alanine, present in higher plants, or with aspartic acid. We studied PSII core complexes purified from D1-N87A and D1-N87D variant strains to probe the function of the D1-Asn87 residue in the water-oxidation mechanism. EPR spectra of the S2 state and flash-induced FTIR spectra of both D1-N87A and D1-N87D PSII core complexes exhibited characteristics similar to those of wildtype Synechocystis PSII core complexes. However, flash-induced O2-evolution studies revealed a decreased cycling efficiency of the D1-N87D variant, whereas the cycling efficiency of the D1-N87A PSII variant was similar to that of wildtype PSII. Steady-state O2-evolution activity assays revealed that substitution of the D1 residue at position 87 with alanine perturbs the chloride-binding site in the proton-exit channel. These findings provide new insight into the role of the D1-Asn87 site in the water-oxidation mechanism and explain the difference in the chloride-binding properties of cyanobacterial and higher-plant PSII.

Binding Thermodynamics and Kinetics Calculations Using Chemical Host and Guest: A Comprehensive Picture of Molecular Recognition.

Understanding the fine balance between changes of entropy and enthalpy and the competition between a guest and water molecules in molecular binding is crucial in fundamental studies and practical applications. Experiments provide measurements. However, illustrating the binding/unbinding processes gives a complete picture of molecular recognition not directly available from experiments, and computational methods bridge the gaps. Here, we investigated guest association/dissociation with β-cyclodextrin (β-CD) by using microsecond-time-scale molecular dynamics (MD) simulations, postanalysis and numerical calculations. We computed association and dissociation rate constants, enthalpy, and solvent and solute entropy of binding. All the computed values of kon, koff, ΔH, ΔS, and ΔG using GAFF-CD and q4MD-CD force fields for β-CD could be compared with experimental data directly and agreed reasonably with experiment findings. In addition, our study further interprets experiments. Both force fields resulted in similar computed ΔG from independently computed kinetics rates, ΔG = -RT ln(kon·C0/koff), and thermodynamics properties, ΔG = ΔH - TΔS. The water entropy calculations show that the entropy gain of desolvating water molecules are a major driving force, and both force fields have the same strength of nonpolar attractions between solutes and β-CD as well. Water molecules play a crucial role in guest binding to β-CD. However, collective water/β-CD motions could contribute to different computed kon and ΔH values by different force fields, mainly because the parameters of β-CD provide different motions of β-CD, hydrogen-bond networks of water molecules in the cavity of free β-CD, and strength of desolvation penalty. As a result, q4MD-CD suggests that guest binding is mostly driven by enthalpy, while GAFF-CD shows that gaining entropy is the major driving force of binding. The study deepens our understanding of ligand-receptor recognition and suggests strategies for force field parametrization for accurately modeling molecular systems.