Ion Pairing of the Neurotransmitters Acetylcholine and Glutamate in Aqueous Solutions

Cation-Binding of Glutamate in Aqueous Solution

Interactions of the cations Li+, Na+, Mg2+, and Ca2+ with L-glutamate (Glu-) in aqueous solution were studied at room temperature with dielectric relaxation spectroscopy in the gigahertz region. Spectra of ∼0.4 M NaGlu with added LiCl, NaCl, MgCl2, or CaCl2 (c(MCln) ≤ 1.5 M) were evaluated and experiments supplemented by density functional theory and 3D reference interaction site model (3D-RISM) calculations. In addition to the modes found for aqueous NaGlu, namely, the reorientation of free Glu- ions (peaking at ∼1.6 GHz), of moderately retarded H2O molecules hydrating the carboxylate moieties of Glu- (∼8.4 GHz), of the cooperative resettling of the H-bond network of bulk water (∼20 GHz), and its preceding fast H-bond flip (∼400 GHz), an additional low-frequency relaxation at ∼0.4 GHz was detected upon the addition of the four salts. In the case of NaGlu + MgCl2(aq) and NaGlu + CaCl2(aq), this mode could be unequivocally assigned to an ion pair formed by the cation and the side-chain carboxylate moiety of Glu-. For NaGlu + LiCl(aq), either this species or a backbone-[Li+-H2O-Cl--Glu-] triple ion is formed. Binding constants increase in the order Li+ Collapse

Key Words Collapse

MESH Headings Collapse

Grants

• Fonds der Chemischen Industrie
• G.A. Krestov Institute of Solution Chemistry, Russian Academy of Sciences

Collapse

Affiliation(s)

Sergej Friesen Institut für Physikalische und Theoretische Chemie, Universität Regensburg, Regensburg D-93040, Germany
Sergey E. Kruchinin G. A. Krestov Institute of Solution Chemistry, Russian Academy of Sciences, Akademicheskaya st. 1, Ivanovo 153045, Russian Federation
Marina V. Fedotova G. A. Krestov Institute of Solution Chemistry, Russian Academy of Sciences, Akademicheskaya st. 1, Ivanovo 153045, Russian Federation
Richard Buchner Institut für Physikalische und Theoretische Chemie, Universität Regensburg, Regensburg D-93040, Germany

Collapse

Molecular-scale understanding of diluent effects on ligand assembly for metal ion separations

Use of metal-selective ligands in solvent extraction is instrumental in extraction of critical materials and recycling, yet, diluent effects on extraction performance are not well understood. Experimental and empirical solvent parameters have been proposed to correlate with extraction performance, but are often inadequate predictors. We follow the hypothesis that the diluents' primary influence on extraction efficiency is whether or not it hinders assembly of the bulky extracting ligands into a geometry necessary for metal complexation. This behavior is readily accessible with molecular dynamics (MD), where the atomistic description of molecules can be applied to arbitrary extractant-solvent molecules and their mixtures. Several simulated quantities are considered, from both pairwise and graph theoretical analyses, and compared to experimental distribution ratio data for americium extraction by TODGA in a series of inert, non-interacting diluents. These simple properties, especially the formation of closed triplets corresponding to the 3 : 1 ligand : metal stoichiometric solvate, suggest a potential predictive power of this approach. This methodology provides a path forward to comprehensively understand and predict diluent effects in more complex systems involving different extracting ligands and multi-component diluent mixtures.

The Three-Dimensional Reference Interaction Site Model Approach as a Promising Tool for Studying Hydrated Viruses and Their Complexes with Ligands

Viruses are the most numerous biological form living in any ecosystem. Viral diseases affect not only people but also representatives of fauna and flora. The latest pandemic has shown how important it is for the scientific community to respond quickly to the challenge, including critically assessing the viral threat and developing appropriate measures to counter this threat. Scientists around the world are making enormous efforts to solve these problems. In silico methods, which allow quite rapid obtention of, in many cases, accurate information in this field, are effective tools for the description of various aspects of virus activity, including virus-host cell interactions, and, thus, can provide a molecular insight into the mechanism of virus functioning. The three-dimensional reference interaction site model (3D-RISM) seems to be one of the most effective and inexpensive methods to compute hydrated viruses, since the method allows us to provide efficient calculations of hydrated viruses, remaining all molecular details of the liquid environment and virus structure. The pandemic challenge has resulted in a fast increase in the number of 3D-RISM calculations devoted to hydrated viruses. To provide readers with a summary of this literature, we present a systematic overview of the 3D-RISM calculations, covering the period since 2010. We discuss various biophysical aspects of the 3D-RISM results and demonstrate capabilities, limitations, achievements, and prospects of the method using examples of viruses such as influenza, hepatitis, and SARS-CoV-2 viruses.

Specific Proton-Donor Properties of Glycine Betaine. Metric Parameters and Enthalpy of Noncovalent Interactions in its Dimer, Water Complexes and Crystalline Hydrate

Trimethylglycine (glycine betaine, GB) is an important organic osmolyte that accumulates in various plant species in response to environmental stresses and has significant potential as a bioactive agent with low environmental impact. It is assumed that the hydration of GB is playing an important role in the protective mechanism. The hydration and aggregation properties of GB have not yet been studied in detail at the atomistic level. In this work, noncovalent interactions in the GB dimer and its complexes with water and crystalline monohydrate are studied. Depending on the object, periodic and non-periodic DFT calculations are used. Particular attention is paid to the metric parameters and enthalpies of intermolecular hydrogen bonds. The identification of noncovalent interactions is carried out by means of the Bader analysis of periodic or non-periodic electron density. The enthalpy of hydrogen bonds is estimated using the Rosenberg formula (PCCP 2 (2000) 2699). The specific proton donor properties of glycine betaine are due to its ability to form intermolecular C-H∙∙∙O bonds with the oxygen atom of a water molecule or the carboxylate group of a neighboring GB. The enthalpy of these bonds can be significantly greater than 10 kJ/mol. The water molecule that forms a hydrogen bond with the carboxylate group of GB also interacts with its CH groups through lone pairs of electrons. The C-H∙∙∙O bonds contribute up to 40% of the total entropy of the GB-water interaction, which is about 45 kJ/mol. The possibility of identifying C-H∙∙∙O bonds by the proton nuclear magnetic resonance method is discussed.

Tebuconazole mediates cognitive impairment via the microbe-gut-brain axis (MGBA) in mice

Tebuconazole, one of the most widely used triazole fungicides, is reported to potentially pose a risk of inducing neurological disorders in human beings. Considering the increasing exposure, whether it could influence cognitive function remains to be elucidated. Herein, we used a mouse model to evaluate the potential cognitive risks and possible mechanisms from the continuous edible application of tebuconazole at low concentrations. Our study revealed that tebuconazole deteriorated spatial learning and memory and downregulated the expression of glutamate receptor subunits. Importantly, metagenomic analysis indicated that tebuconazole not only led to significant shifts in the composition and diversity of the gut microbiota but also changed intestinal homeostasis. Specifically, after exposure, tebuconazole circulated in the bloodstream and largely entered the gut-brain axis for disruption, including disturbing the Firmicutes/Bacteroidetes ratio, interrelated neurotransmitters and systemic immune factors. Moreover, pretreatment with probiotics improved immune factor expression and restored the deterioration of synaptic function and spatial learning and memory. The current study provides novel insights concerning perturbations of the gut microbiome and its functions as a potential new mechanism by which tebuconazole exposes cognitive function-related human health.

Computational Advances of Protein/Neurotransmitter-membrane Interactions Involved in Vesicle Fusion and Neurotransmitter Release

Vesicle fusion is of crucial importance to neuronal communication at neuron terminals. The exquisite but complex fusion machinery for neurotransmitter release is tightly controlled and regulated by protein/neurotransmitter-membrane interactions. Computational 'microscopies', in particular molecular dynamics simulations and related techniques, have provided notable insight into the physiological process over the past decades, and have made enormous contributions to fields such as neurology, pharmacology and pathophysiology. Here we review the computational advances of protein/neurotransmitter-membrane interactions related to presynaptic vesicle-membrane fusion and neurotransmitter release, and outline the in silico challenges ahead for understanding this important physiological process.

Molecular Dynamics Simulation of Association Processes in Aqueous Solutions of Maleate Salts of Drug-like Compounds: The Role of Counterion

The study of the formation of microstructures during the interaction of a protonated drug-like compound (API) with a maleic acid monoanion sheds light on the assembly processes in an aqueous solution at the molecular level. Molecular dynamics (MD) simulations coupled with density functional theory (DFT) calculations made it possible to find initial hydrogen bonding motifs during the assembly process, leading to the formation of heterodimers and trimers. The process of trimer formation [protonated API—maleic acid monoanion—protonated API] proceeds through the formation of three intermolecular H-bonds by the CO₂⁻ group of the maleic acid monoanion in both systems. The total enthalpy/energy of these H-bonds is more than 70 kJ/mol. Thus, the maleic acid monoanion plays a key role in the processes of association in aqueous solution, and the interaction of the maleic acid monoanion with API is more preferable than the interaction of API molecules with each other. DFT computations in the discrete continuum approximation reveal the spectral features of heterodimers and trimers, and the ATR-IR spectra confirmed these findings. MD simulations followed by DFT calculations made it possible to describe the initial stages of the formation of pharmaceutical cocrystals in an aqueous solution.

Diclofenac Ion Hydration: Experimental and Theoretical Search for Anion Pairs

Self-assembly of organic ions in aqueous solutions is a hot topic at the present time, and substances that are well-soluble in water are usually studied. In this work, aqueous solutions of sodium diclofenac are investigated, which, like most medicinal compounds, is poorly soluble in water. Classical MD modeling of an aqueous solution of diclofenac sodium showed equilibrium between the hydrated anion and the hydrated dimer of the diclofenac anion. The assignment and interpretation of the bands in the UV, NIR, and IR spectra are based on DFT calculations in the discrete-continuum approximation. It has been shown that the combined use of spectroscopic methods in various frequency ranges with classical MD simulations and DFT calculations provides valuable information on the association processes of medical compounds in aqueous solutions. Additionally, such a combined application of experimental and calculation methods allowed us to put forward a hypothesis about the mechanism of the effect of diclofenac sodium in high dilutions on a solution of diclofenac sodium.

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.