Enthalpy-Entropy Interplay in π-Stacking Interaction of Benzene Dimer in Water

Unveiling the impacts of microplastics on cadmium transfer in the soil-plant-human system: A review

The co-contamination of soils by microplastics (MPs) and cadmium (Cd), one of the most perilous heavy metals, is emerging as a significant global concern, posing risks to plant productivity and human health. However, there remains a gap in the literature concerning comprehensive evaluations of the combined effects of MPs and Cd on soil-plant-human systems. This review examines the interactions and co-impacts of MPs and Cd in soil-plant-human systems, elucidating their mechanisms and synergistic effects on plant development and health risks. We also review the origins and contamination levels of MPs and Cd, revealing that sewage, atmospheric deposition, and biosolid applications are contributors to the contamination of soil with MPs and Cd. Our meta-analysis demonstrates that MPs significantly (p<0.05) increase the bioavailability of soil Cd and the accumulation of Cd in plant shoots by 6.9 and 9.3 %, respectively. The MPs facilitate Cd desorption from soils through direct adsorption via surface complexation and physical adsorption, as well as indirectly by modifying soil physicochemical properties, such as pH and dissolved organic carbon, and altering soil microbial diversity. These interactions augment the bioavailability of Cd, along with MPs, adversely affect plant growth and its physiological functions. Moreover, the ingestion of MPs and Cd through the food chain significantly enhances the bioaccessibility of Cd and exacerbates histopathological alterations in human tissues, thereby amplifying the associated health risks. This review provides insights into the coexistence of MPs and Cd and their synergistic effects on soil-plant-human systems, emphasizing the need for further research in this critical subject area.

Unveiling the Role of Water on π-π Stacking Through Microwave Spectroscopy of (Thiophene)2-(Water)1-2 Clusters

There has been an ongoing debate about whether water enhances or hinders π-π stacking, a phenomenon crucial in various biological and chemical systems. In this study, the influence of water on π-π stacking is investigated by microwave spectroscopic observation of gas-phase hydrated clusters of thiophene dimers. Two isomers of (C4H4S)2-H2O and two isomers of (C4H4S)2-(H2O)2 have been unambiguously identified. These identifications are supported by quantum chemistry calculations and isotopic measurements. In each of these conformations, water molecules are situated between aromatic pairs, forming distinctive interactions. Water molecules engage with thiophene molecules either as hydrogen bond donors through OH···π interactions or as hydrogen bond acceptors through CH···O interactions. The energy decomposition analysis indicates that the bonding pattern of water molecules significantly affects the π···π interactions between aromatic rings. These findings offer valuable structural insights into the role of water in shaping π-π stacking.

Chitinase is a Potent Insecticidal Molecular Target of Camptothecin and Its Derivatives

Camptothecin (CPT) is a prominent molecule in natural product research because of its application prospects in medicine and agriculture. In this study, CPT and its derivatives were discovered to be competitive inhibitors of group II and group h insect chitinases, both of which are key components of insect chitinolytic systems. CPT and 7-ethyl-10-hydroxycamptothecin (SN-38) inhibited group II chitinase from Ostrinia furnacalis (OfChtII) with K_i values of 5.1 and 2.0 μM, respectively. Results from tryptophan fluorescence spectroscopy, molecular docking analysis, and molecular dynamics simulations revealed that both CPT and SN-38 inhibit OfChtII-C1 by interacting with solvent-exposed tryptophan residues in a substrate-binding cleft. CPT exhibited high insecticidal activity toward the orthopteran pest Locusta migratoria, possibly because of the midgut metabolism of CPT, with only moderate activities toward lepidopteran pests. Even though SN-38 exhibited much lower insecticidal activities than CPT, it still showed higher inhibitory activity toward chitinase. This study reports a new molecular target of CPT and provides insights into molecular design of CPT-based insecticides against different kinds of pests.

Cumulative Millisecond-Long Sampling for a Comprehensive Energetic Evaluation of Aqueous Ionic Liquid Effects on Amino Acid Interactions

The interactions of amino acid side-chains confer diverse energetic contributions and physical properties to a protein's stability and function. Various computational tools estimate the effect of changing a given amino acid on the protein's stability based on parametrized (free) energy functions. When parametrized for the prediction of protein stability in water, such energy functions can lead to suboptimal results for other solvents, such as ionic liquids (IL), aqueous ionic liquids (aIL), or salt solutions. However, to our knowledge, no comprehensive data are available describing the energetic effects of aIL on intramolecular protein interactions. Here, we present the most comprehensive set of potential of mean force (PMF) profiles of pairwise protein-residue interactions to date, covering 50 relevant interactions in water, the two biotechnologically relevant aIL [BMIM/Cl] and [BMIM/TfO], and [Na/Cl]. These results are based on a cumulated simulation time of >1 ms. aIL and salt ions can weaken, but also strengthen, specific residue interactions by more than 3 kcal mol^-1, depending on the residue pair, residue-residue configuration, participating ions, and concentration, necessitating considering such interactions specifically. These changes originate from a complex interplay of competitive or cooperative noncovalent ion-residue interactions, changes in solvent structural dynamics, or unspecific charge screening effects and occur at the contact distance but also at larger, solvent-separated distances. This data provide explanations at the atomistic and energetic levels for complex IL effects on protein stability and should help improve the prediction accuracies of computational tools that estimate protein stability based on (free) energy functions.

Dynamic Long-Range Interactions Influence Substrate Binding and Catalysis by Human Histidine Triad Nucleotide-Binding Proteins (HINTs), Key Regulators of Multiple Cellular Processes and Activators of Antiviral ProTides

Human histidine triad nucleotide-binding (hHINT) proteins catalyze nucleotide phosphoramidase and acyl-phosphatase reactions that are essential for the activation of antiviral proTides, such as Sofosbuvir and Remdesivir. hHINT1 and hHINT2 are highly homologous but exhibit disparate roles as regulators of opioid tolerance (hHINT1) and mitochondrial activity (hHINT2). NMR studies of hHINT1 reveal a pair of dynamic surface residues (Q62, E100), which gate a conserved water channel leading to the active site 13 Å away. hHINT2 crystal structures identify analogous residues (R99, D137) and water channel. hHINT1 Q62 variants significantly alter the steady-state k_cat and K_m for turnover of the fluorescent substrate (TpAd), while stopped-flow kinetics indicate that K_D also changes. hHINT2, like hHINT1, exhibits a burst phase of adenylation, monitored by fluorescent tryptamine release, prior to rate-limiting hydrolysis and nucleotide release. hHINT2 exhibits a much smaller burst-phase amplitude than hHINT1, which is further diminished in hHINT2 R99Q. Kinetic simulations suggest that amplitude variations can be accounted for by a variable fluorescent yield of the E·S complex from changes in the environment of bound TpAd. Isothermal titration calorimetry measurements of inhibitor binding show that these hHINT variants also alter the thermodynamic binding profile. We propose that these altered surface residues engender long-range dynamic changes that affect the orientation of bound ligands, altering the thermodynamic and kinetic characteristics of hHINT active site function. Thus, studies of the cellular roles and proTide activation potential by hHINTs should consider the importance of long-range interactions and possible protein binding surfaces far from the active site.

Thermodynamics of π–π Interactions of Benzene and Phenol in Water

The π–π interaction is a major driving force that stabilizes protein assemblies during protein folding. Recent studies have additionally demonstrated its involvement in the liquid–liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs). As the participating residues in IDPs are exposed to water, π–π interactions for LLPS must be modeled in water, as opposed to the interactions that are often established at the hydrophobic domains of folded proteins. Thus, we investigated the association of free energies of benzene and phenol dimers in water by integrating van der Waals (vdW)-corrected density functional theory (DFT) and DFT in classical explicit solvents (DFT-CES). By comparing the vdW-corrected DFT and DFT-CES results with high-level wavefunction calculations and experimental solvation free energies, respectively, we established the quantitative credibility of these approaches, enabling a reliable prediction of the benzene and phenol dimer association free energies in water. We discovered that solvation influences dimer association free energies, but not significantly when no direct hydrogen-bond-type interaction exists between two monomeric units, which can be explained by the enthalpy–entropy compensation. Our comprehensive computational study of the solvation effect on π–π interactions in water could help us understand the molecular-level driving mechanism underlying the IDP phase behaviors.

Molecular photoswitches in aqueous environments

Molecular photoswitches enable dynamic control of processes with high spatiotemporal precision, using light as external stimulus, and hence are ideal tools for different research areas spanning from chemical biology to smart materials. Photoswitches are typically organic molecules that feature extended aromatic systems to make them responsive to (visible) light. However, this renders them inherently lipophilic, while water-solubility is of crucial importance to apply photoswitchable organic molecules in biological systems, like in the rapidly emerging field of photopharmacology. Several strategies for solubilizing organic molecules in water are known, but there are not yet clear rules for applying them to photoswitchable molecules. Importantly, rendering photoswitches water-soluble has a serious impact on both their photophysical and biological properties, which must be taken into consideration when designing new systems. Altogether, these aspects pose considerable challenges for successfully applying molecular photoswitches in aqueous systems, and in particular in biologically relevant media. In this review, we focus on fully water-soluble photoswitches, such as those used in biological environments, in both in vitro and in vivo studies. We discuss the design principles and prospects for water-soluble photoswitches to inspire and enable their future applications.

Conformational Shifts of Stacked Heteroaromatics: Vacuum vs. Water Studied by Machine Learning

Stacking interactions play a crucial role in drug design, as we can find aromatic cores or scaffolds in almost any available small molecule drug. To predict optimal binding geometries and enhance stacking interactions, usually high-level quantum mechanical calculations are performed. These calculations have two major drawbacks: they are very time consuming, and solvation can only be considered using implicit solvation. Therefore, most calculations are performed in vacuum. However, recent studies have revealed a direct correlation between the desolvation penalty, vacuum stacking interactions and binding affinity, making predictions even more difficult. To overcome the drawbacks of quantum mechanical calculations, in this study we use neural networks to perform fast geometry optimizations and molecular dynamics simulations of heteroaromatics stacked with toluene in vacuum and in explicit solvation. We show that the resulting energies in vacuum are in good agreement with high-level quantum mechanical calculations. Furthermore, we show that using explicit solvation substantially influences the favored orientations of heteroaromatic rings thereby emphasizing the necessity to include solvation properties starting from the earliest phases of drug design.

Prediction of octanol-air partition coefficients for PCBs at different ambient temperatures based on the solvation free energy and the dimer ratio

Temperature-dependent octanol-air partition coefficients (K_OA) are of great importance in assessing the environmental behavior and fate of persistent organic pollutants including polychlorinated biphenyls (PCBs). Due to the tremendous amounts of time, effort and cost needed for the experimental determination of K_OA, it is desirable to develop a rapid and precise predictive method to estimate K_OA just based on molecular structure. In the present study, a predictive model for log K_OA of PCBs at ambient temperatures was developed based on the thermodynamic relationship between K_OA and the solvation free energy from air to octanol (ΔG_OA). For the calculation of ΔG_OA of PCBs, the optimal combination of theoretical method and basis-set was identified to be HF/MIDI!6D for both geometry optimization and energy calculation. Dimer formation could affect the partition behavior and promote the apparent K_OA values of PCBs. After taking the effect of dimer formation into account, the goodness-of-fit, predictive ability, and robustness of the predictive model were significantly improved. Apparent log K_OA values of PCBs at different ambient temperatures ranging from 283.15 to 303.15 K were predicted. Compared with other reported models, the model developed in the present study had not only comparable goodness-of-fit and predictive ability, but also a universal application domain and the relative independency of experimental data. Therefore, the solvation free energy method could be a promising method for the prediction of K_OA.

Progressive Hydrophobicity of Fluorobenzenes

The potentials of mean force for the dimers of fluorobenzenes sample both π-stacked and T-shaped structures for partially fluorinated benzenes, namely, 1,4-difluorobenzene, 1,3,5-trifluorobenzene, and 1,2,4,5-tetrafluorobenzene, and sample only the T-shaped structures for benzene and hexafluorobenzene. While the free energy for the dimerization in water is very weakly dependent on the number of fluorine atoms, the formation of π-stacked structures is entropy-driven and the T-shaped structures appear due to an enthalpic minimum. Interestingly, the solvation behavior suggests that the accumulation of water around the contact and solvent-separated pairs decreases with the increase in the number of fluorine atoms, which signifies progressive hydrophobicity of fluorobenzenes.

Hydration of Aromatic Heterocycles as an Adversary of π-Stacking

Hydration is one of the key players in the protein-ligand binding process. It not only influences the binding process per se, but also the drug's absorption, distribution, metabolism, and excretion properties. To gain insights into the hydration of aromatic cores, the solvation thermodynamics of 40 aromatic mono- and bicyclic systems, frequently occurring in medicinal chemistry, are investigated. Thermodynamics is analyzed with two different methods: grid inhomogeneous solvation theory (GIST) and thermodynamic integration (TI). Our results agree well with previously published experimental and computational solvation free energies. The influence of adding heteroatoms to aromatic systems and how the position of these heteroatoms impacts the compound's interactions with water is studied. The solvation free energies of these heteroaromatics are highly correlated to their gas phase interaction energies with benzene: compounds showing a high interaction energy also have a high solvation free energy value. Therefore, replacing a compound with one having a higher gas phase interaction energy might not result in the expected improvement in affinity. The desolvation costs counteract the higher stacking interactions, hence weakening or even inverting the expected gain in binding free energy.

Structure, Dynamics, and Wettability of Water at Metal Interfaces

The water/metal interface often governs important chemophysical processes in various technologies. Therefore, from scientific and engineering perspectives, the detailed molecular-level elucidation of the water/metal interface is of high priority, but the related research is limited. In experiments, the surface-science techniques, which can provide full structural details of the surface, are not easy to directly apply to the interfacial systems under ambient conditions, and the well-defined facets cannot be entirely free from contamination at the contact with water. To answer long-standing debates regarding the wettability, structure, and dynamics of water at metal interfaces, we here develop reliable first-principles-based multiscale simulations. Using the state-of-the-art simulations, we find that the clean metal surfaces are actually superhydrophilic and yield zero contact angles. Furthermore, we disclose an inadequacy of widespread ice-like bilayer model of the water adlayers on metal surfaces from both averaged structural and dynamic points of view. Our findings on the nature of water on metal surfaces provide new molecular level perspectives on the tuning and design of water/metal interfaces that are at the heart of many energy applications.