Thermodynamics of site-specific small molecular ion interactions with DNA duplex: a molecular dynamics study

Ionic Pairing and Selective Solvation of Butylmethylimidazolium Chloride Ion Pairs in DMSO-Water Mixtures: A Comprehensive Examination via Molecular Dynamics Simulations and Potentials of Mean Force Analysis

Ionic liquids (ILs) with dimethyl sulfoxide (DMSO) and water act as a promising solvent medium for the dissolution of cellulose in an efficient manner. To develop a proper solvent system, it is really important to understand the thermodynamics of the molecular solutions consisting of ILs, DMSO, and water. The ion-pairing propensity of the ILs in the presence of DMSO and water plays a crucial role in governing the property of the solvent mixtures. Employing all-atom molecular dynamics simulations, we estimate the potentials of mean force between BMIM+ and Cl- ions in DMSO-water mixtures. Analysis reveals a significant increase in the thermodynamic stability of both contact ion pair (CIP) and solvent-assisted ion pair (SAIP) states with a rising DMSO mole fraction. Thermodynamic assessments highlight the entropic stabilization of CIP states and SAIP states in pure water, in DMSO-water mixtures, and in pure DMSO. The structural analysis reveals that in comparison to the DMSO local density, the local water density is relatively very high around ion pairs, more specifically in the solvation shell of a chloride ion. Preferential binding coefficients also consistently indicate exclusion of DMSO from the ion pair in DMSO-water mixtures. To enhance our understanding regarding the solvent molecules kinetics around the ion pairs, the survival probabilities of DMSO and water are computed. The calculations reveal that the water molecules prefer a prolonged stay in the solvation shell of Cl- ions.

Role of Terminal Groups of cis-1,4-Polyisoprene Chains in the Formation of Physical Junction Points in Natural Rubber

The terminal structures of cis-1,4-polyisoprene (PI) chains play a vital role in the excellent comprehensive performance of Hevea natural rubber (NR) with properties such as high toughness, tear-resistance, and wet skid resistance. The cis-1,4-polyisoprene chain constituting NR exhibits a distinct composition of terminal groups comprising two distinct types, namely, the ω and α terminal groups. The structures of the ω terminal [dimethyl allyl (DMA)-(trans-1,4-isoprene)2] and six kinds of α end groups of the polymer chain of NR have been explored by utilizing a newly developed 2D NMR method. In the present work, we examine different kinds of PI melt systems, and we choose various combinations of terminal groups: Hydrogen, one DMA unit with two trans isoprene units as ω end groups and ester-terminated isopentene (α1), hydroxy-terminated isopentene (α2), ester-terminated isobutane (α3), hydroxy-terminated isobutane (α4), ester-terminated 1,4-cis-isoprene (α5), and hydroxy-terminated 1,4-cis-isoprene (α6), i.e., HPIH (PI0)-pure PI (Hydrogen terminal), ωPIα1 (PII), ωPIα2 (PIII), ωPIα3 (PIIII), ωPIα4 (PIIV), ωPIα5 (PIV), and ωPIα6 (PIVI). We evaluated dynamic and static properties of PI chains such as the end-to-end vector autocorrelation function (C(t)), its average relaxation time (τ), end-to-end distance (Ree), and radius of gyration (Rg). We also estimated the diffusion coefficients of polyisoprene chains and pair correlation functions [radial distribution functions (RDFs)], potentials of mean force (PMFs) in between end residues, and survival probability (P(τ)) of end groups around the end group by analyzing the equilibrated trajectories of full-atom MD simulations. As per the examination of C(t), rotational relaxation time τ, and RDFs, we discovered that the existence of a strong hydrogen bond in α2-α2, α4-α4, and α6-α6 residues makes the dynamics of hydroxy-terminated polyisoprene chains in ωPIα2,α4,α6 melt systems slower. From the analyses of RDFs and PMFs (W(r)), the association between [α2]-[α2], [α4]-[α4], and [α6]-[α6] terminals in ωPIα2,α4,α6 melt systems is significantly stronger than in [ISO]-[ISO] [Hydrogen terminated 1,4-cis-isoprene:(ISO)] in HPIH and ω-ω, [α1]-[α1], [α3]-[α3], and [α5]-[α5] in ωPIα1,α3,α5 systems. We quantified the fraction of cluster formation of terminal groups of a given size in the seven PI melt systems by employing the criteria of PMFs. It is revealed that no stable cluster exists in the HPIH, ωPIα1, ωPIα3, and ωPIα5 melt systems. Conversely, in the ωPIα2, ωPIα4, and ωPIα6 systems, we perceived stable clusters of [(α2)p] [(α4)p] and [(α6)p] end groups where p (2 ≤ x ≤ 6). These stable clusters validate the presence of physical junction points in between hydroxy-terminated polyisoprene chains through their α2, α4, and α6 terminals. These physical junction points might be crucial for superior properties of NR such as high toughness, crack growth resistance, and strain-induced crystallization.

In Silico Characterization of the Binding Modes of Surfactants with Bovine Serum Albumin

The binding interactions of the surfactants: anionic sodium dodecyl sulphate (SDS), cationic cetyltrimethylammonium bromide (CTAB), non-ionic octyl glucoside (OG), and zwitterionic 3-[Hexadecyl(dimethyl)ammonio]-1-propanesulfonate (HPS), with bovine serum albumin (BSA) were investigated by computer simulation. The results disclosed that the surfactants bound stably between hydrophobic subdomain IIA and IIIA where tryptophan-213 residue, an important intrinsic fluorophore in BSA is housed. The interactions of the surfactants with the BSA were electrostatic and hydrophobic interactions. The head-groups of SDS, HPS and OG formed hydrogen bonds with the BSA, while that of CTAB was shielded from intermolecular hydrogen-bonding due to intervening methyl groups. Subsequently, molecular dynamics (MD) simulation of the protein-surfactant complexes revealed that hydrogen bonds formed by OG were stronger than those of SDS and HPS. However, the decomposed force-field energies showed that OG had the least interaction energy with the BSA. In addition to MD simulation, it was found by density functional theory (DFT) that the differences in the coulomb interaction energies can be attributed to charge distribution in the surfactants. Overall, free energies calculated by linear interaction energy (LIE) proved that the binding of each surfactant was dominated by differences between van der Waals interactions in bound and free states.

Molecular dynamics simulation elucidates the preferential binding affinity of sodium and tetramethylammonium ions for tetrameric Nafion unit under aqueous conditions

Preferential binding interaction of sodium and tetramethylammonium cations with Nafion pendant and backbone sites respectively.