skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Computational predictions of metal–macrocycle stability constants require accurate treatments of local solvent and pH effects
Title: Computational predictions of metal–macrocycle stability constants require accurate treatments of local solvent and pH effects
Rational design of molecular chelating agents requires a detailed understanding of physicochemical ligand–metal interactions in solvent phase. Computational quantum chemistry methods should be able to provide this, but computational reports have shown poor accuracy when determining absolute binding constants for many chelating molecules. To understand why, we compare and benchmark static- and dynamics-based computational procedures for a range of monovalent and divalent cations binding to a conventional cryptand molecule: 2.2.2-cryptand ([2.2.2]). The benchmarking comparison shows that dynamics simulations using standard OPLS-AA classical potentials can reasonably predict binding constants for monovalent cations, but these procedures fail for divalent cations. We also consider computationally efficient static procedure using Kohn–Sham density functional theory (DFT) and cluster-continuum modeling that accounts for local microsolvation and pH effects. This approach accurately predicts binding energies for monovalent and divalent cations with an average error of 3.2 kcal mol −1 compared to experiment. This static procedure thus should be useful for future molecular screening efforts, and high absolute errors in the literature may be due to inadequate modeling of local solvent and pH effects.  more » « less
Award ID(s):
1705592
PAR ID:
10310191
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
23
Issue:
15
ISSN:
1463-9076
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, Proceedings of the National Academy of Sciences)
    RNA molecules cannot fold in the absence of counterions. Experiments are typically performed in the presence of monovalent and divalent cations. How to treat the impact of a solution containing a mixture of both ion types on RNA folding has remained a challenging problem for decades. By exploiting the large concentration difference between divalent and monovalent ions used in experiments, we develop a theory based on the reference interaction site model (RISM), which allows us to treat divalent cations explicitly while keeping the implicit screening effect due to monovalent ions. Our theory captures both the inner shell and outer shell coordination of divalent cations to phosphate groups, which we demonstrate is crucial for an accurate calculation of RNA folding thermodynamics. The RISM theory for ion–phosphate interactions when combined with simulations based on a transferable coarse-grained model allows us to predict accurately the folding of several RNA molecules in a mixture containing monovalent and divalent ions. The calculated folding free energies and ion-preferential coefficients for RNA molecules (pseudoknots, a fragment of the rRNA, and the aptamer domain of the adenine riboswitch) are in excellent agreement with experiments over a wide range of monovalent and divalent ion concentrations. Because the theory is general, it can be readily used to investigate ion and sequence effects on DNA properties. 
    more » « less
  2. ; ; ; ; (, Journal of Membrane Science Letters)
    Understanding the mechanisms of molecular transport in polyamide membranes is imperative to improve their solute-specific selectivity. We explored the partitioning behaviors of water and salts in polyamide membranes to elucidate the role of ion-membrane interactions in the transport. Quartz crystal microbalance (QCM) was employed to quantify the mass uptake at different temperatures and determine partition energies (Ek) for water and salts under two different pH values. Zeta potential and permeability tests were conducted to support the ionmembrane affinity trends observed with QCM and link these trends to ion-ion selectivity. Our results demonstrate a high affinity of water to the polyamide membrane (Ek < 0), with a significant swelling effect attributed to dipole interactions and hydrogen bonding. Ion partitioning revealed distinct differences between monovalent and divalent cations, as well as between kosmotropic and chaotropic anions. Specifically, divalent cations (Ca2+ and Mg2+) exhibited considerably lower partition energies (-0.99 and 0.29 kcal mol-1, respectively) and more efficient charge neutralization, indicating stronger interactions with the membrane compared to monovalent cations (~2.2 kcal mol-1). The partition energies of the chaotropic iodide and kosmotropic sulphate anions were substantially different (-5.5 and 4.0 kcal mol-1, respectively), likely due to the different tendency of these anions to shed their hydration shell and stick to the polymer. Last, our permeability tests indicate the potential existence of an intrinsic tradeoff between ion partitioning and intrapore diffusion, presumably due to the opposite effects that ion-membrane interactions have on these transport steps. Overall, our work underscores the role of ionspecific interactions in membrane transport and selectivity. 
    more » « less
  3. ; ; ; ; ; ; (, Chemistry – An Asian Journal)
    Abstract A reverse‐binding‐selectivity between monovalent and divalent cations was observed for two different self‐assembly G16‐hexadecamer and G8‐octamer systems. The dissociation constant between G4‐quadruplex and monomer was calculated via VT‐1H NMR experiments. Quantitative energy profiles revealed entropy as the key factor for the weaker binding toward Ba2+compared with K+in the G8‐octamer system despite stronger ion‐dipole interactions. This study is the first direct comparison of the G4‐quartet binding affinity between mono and divalent cations and will benefit future applications of G‐quadruplex‐related research. Further competition experiments between the G8‐octamer and 18‐crown‐6 with K+demonstrated the potential of this G8system as a new potassium receptor. 
    more » « less
  4. ; ; (, Physical Chemistry Chemical Physics)
    null (Ed.)
    The sodium anion (Na − ) was once thought to behave like a ‘genuine’ anion, with both the [Ne] core and the 3s valence shell interacting very weakly with their environments. In the present work, following a recent study of the surprisingly small quadrupolar line widths of Na − , NMR shielding calculations were carried out for the Na − /Na + [2.2.2]cryptand system solvated in methylamine, based on ab initio molecular dynamics simulations, followed by detailed analyses of the shielding constants. The results confirm that Na − does not act like a quasi-free ion that interacts only weakly with its surroundings. Rather, the filled 3s shell of Na − interacts strongly with its chemical environment, but only weakly with the ion's own core and the nucleus, and it isolates the core from the chemical environment. As a consequence, the Na − ion appears in NMR experiments like a free ion. 
    more » « less
  5. ; ; ; ; ; ; (, Soft Matter)
    Electric double layers are crucial to energy storage and electrocatalytic device performance. While double layer formation originates in electrostatic interactions, electric double layer properties are governed by a balance of both electrostatic and entropic driving forces; favorable ion-surface electrostatic interactions attract counterions to charged surfaces to compensate, or "screen," potentials, but the confinement of these same ions from a bulk reservoir to the interface incurs an entropic penalty. Here, we use a dicationic imidazolium ionic liquid and its monovalent analogue to explore how cation valence and entropy influence double layer formation and electrochemical reactivity using CO2 electroreduction as a model reaction. We find that divalent and monovalent cations display similar CO2 reduction kinetics but differ vastly in steady-state reactivity due to rapid electrochemically induced precipitation of insulating dicationic (bi)carbonate films. Using in situ surface-enhanced Raman scattering spectroscopy, we find that potential-dependent reorientation occurs at similar potentials between the two ionic liquids, but the introduction of a covalent link in the divalent cation imparts a more ordered double layer structure that favors (bi)carbonate precipitation. In mixed monovalent-divalent electrolytes, we find that the divalent cations dominate interfacial properties by preferentially accumulating at surfaces even at very low relative concentrations. Our findings confirm that ion entropy plays a key role in modulating local electrochemical environments and highlight how double layer properties are very sensitive to the properties of counterions that pay the lowest entropic penalty to accumulate at interfaces. Overall, we illustrate that ion entropy provides a new knob to tune reaction microenvironments and unveil how entropy plays a major role in modulating electrochemical reactivity in mixed ion electrolytes. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email