Explicit solvation thermodynamics in ionic solution: extending grid inhomogeneous solvation theory to solvation free energy of salt–water mixtures

Hydration thermodynamics play a fundamental role in fields ranging from the pharmaceutical industry to environmental research. Numerous methods exist to predict solvation thermodynamics of compounds ranging from small molecules to large biomolecules. Arguably the most precise methods are those based on molecular dynamics (MD) simulations in explicit solvent. One theory that has seen increased use is inhomogeneous solvation theory (IST). However, while many applications require accurate description of salt–water mixtures, no implementation of IST is currently able to estimate solvation properties involving more than one solvent species. Here, we present an extension to grid inhomogeneous solvation theory (GIST) that can take salt contributions into account. At the example of carbazole in 1 M NaCl solution, we compute the solvation energy as well as first and second order entropies. While the effect of the first order ion entropy is small, both the water–water and water–ion entropies contribute strongly. We show that the water–ion entropies are efficiently approximated using the Kirkwood superposition approximation. However, this approach cannot be applied to the water–water entropy. Furthermore, we test the quantitative validity of our method by computing salting-out coefficients and comparing them to experimental data. We find a good correlation to experimental salting-out constants, while the absolute values are overpredicted due to the approximate second order entropy. Since ions are frequently used in MD, either to neutralize the system or as a part of the investigated process, our method greatly extends the applicability of GIST. The use-cases range from biopharmaceuticals, where many assays require high salt concentrations, to environmental research, where solubility in sea water is important to model the fate of organic substances.

Pyranose Ring Puckering Thermodynamics for Glycan Monosaccharides Associated with Vertebrate Proteins

The conformational properties of carbohydrates can contribute to protein structure directly through covalent conjugation in the cases of glycoproteins and proteoglycans and indirectly in the case of transmembrane proteins embedded in glycolipid-containing bilayers. However, there continue to be significant challenges associated with experimental structural biology of such carbohydrate-containing systems. All-atom explicit-solvent molecular dynamics simulations provide a direct atomic resolution view of biomolecular dynamics and thermodynamics, but the accuracy of the results depends on the quality of the force field parametrization used in the simulations. A key determinant of the conformational properties of carbohydrates is ring puckering. Here, we applied extended system adaptive biasing force (eABF) all-atom explicit-solvent molecular dynamics simulations to characterize the ring puckering thermodynamics of the ten common pyranose monosaccharides found in vertebrate biology (as represented by the CHARMM carbohydrate force field). The results, along with those for idose, demonstrate that the CHARMM force field reliably models ring puckering across this diverse set of molecules, including accurately capturing the subtle balance between 4C1 and 1C4 chair conformations in the cases of iduronate and of idose. This suggests the broad applicability of the force field for accurate modeling of carbohydrate-containing vertebrate biomolecules such as glycoproteins, proteoglycans, and glycolipids.

Ergodic Algorithmic Model (EAM), with Water as Implicit Solvent, in Chemical, Biochemical, and Biological Processes

For many years, we have devoted our research to the study of the thermodynamic properties of hydrophobic hydration processes in water, and we have proposed the Ergodic Algorithmic Model (EAM) for maintaining the thermodynamic properties of any hydrophobic hydration reaction at a constant pressure from the experimental determination of an equilibrium constant (or other potential functions) as a function of temperature. The model has been successfully validated by the statistical analysis of the information elements provided by the EAM model for about fifty compounds. The binding functions are convoluted functions, RlnKeq = {f(1/T)* g(T)} and RTlnKeq = {f(T)* g(lnT)}, where the primary linear functions f(1/T) and f(T) are modified and transformed into parabolic curves by the secondary functions g(T) and g(lnT), respectively. Convoluted functions are consistent with biphasic dual-structure partition function, {DS-PF} = {M-PF} ∙ {T-PF} ∙ {ζw}, composed by ({M-PF} (Density Entropy), {T-PF}) (Intensity Entropy), and {ζw} (implicit solvent). In the present paper, after recalling the essential aspects of the model, we outline the importance of considering the solvent as “implicit” in chemical and biochemical reactions. Moreover, we compare the information obtained by computer simulations using the models till now proposed with “explicit” solvent, showing the mess of information lost without considering the experimental approach of the EAM model.

Insights into the Ligand Binding to Bromodomain-Containing Protein 9 (BRD9): A Guide to the Selection of Potential Binders by Computational Methods

The estimation of the binding of a set of molecules against BRD9 protein was carried out through an in silico molecular dynamics-driven exhaustive analysis to guide the identification of potential novel ligands. Starting from eight crystal structures of this protein co-complexed with known binders and one apo form, we conducted an exhaustive molecular docking/molecular dynamics (MD) investigation. To balance accuracy and an affordable calculation time, the systems were simulated for 100 ns in explicit solvent. Moreover, one complex was simulated for 1 µs to assess the influence of simulation time on the results. A set of MD-derived parameters was computed and compared with molecular docking-derived and experimental data. MM-GBSA and the per-residue interaction energy emerged as the main indicators for the good interaction between the specific binder and the protein counterpart. To assess the performance of the proposed analysis workflow, we tested six molecules featuring different binding affinities for BRD9, obtaining promising outcomes. Further insights were reported to highlight the influence of the starting structure on the molecular dynamics simulations evolution. The data confirmed that a ranking of BRD9 binders using key parameters arising from molecular dynamics is advisable to discard poor ligands before moving on with the synthesis and the biological tests.

Perspectives for the reconstruction of 3D chromatin conformation using single cell Hi-C data

Construction of chromosomes 3D models based on single cell Hi-C data constitute an important challenge. We present a reconstruction approach, DPDchrom, that incorporates basic knowledge whether the reconstructed conformation should be coil-like or globular and spring relaxation at contact sites. In contrast to previously published protocols, DPDchrom can naturally form globular conformation due to the presence of explicit solvent. Benchmarking of this and several other methods on artificial polymer models reveals similar reconstruction accuracy at high contact density and DPDchrom advantage at low contact density. To compare 3D structures insensitively to spatial orientation and scale, we propose the Modified Jaccard Index. We analyzed two sources of the contact dropout: contact radius change and random contact sampling. We found that the reconstruction accuracy exponentially depends on the number of contacts per genomic bin allowing to estimate the reconstruction accuracy in advance. We applied DPDchrom to model chromosome configurations based on single-cell Hi-C data of mouse oocytes and found that these configurations differ significantly from a random one, that is consistent with other studies.

Beyond Continuum Solvent Models in Computational Homogeneous Catalysis

In homogeneous catalysis solvent is an inherent part of the catalytic system. As such, it must be considered in the computational modeling. The most common approach to include solvent effects in quantum mechanical calculations is by means of continuum solvent models. When they are properly used, average solvent effects are efficiently captured, mainly those related with solvent polarity. However, neglecting atomistic description of solvent molecules has its limitations, and continuum solvent models all alone cannot be applied to whatever situation. In many cases, inclusion of explicit solvent molecules in the quantum mechanical description of the system is mandatory. The purpose of this article is to highlight through selected examples what are the reasons that urge to go beyond the continuum models to the employment of micro-solvated (cluster-continuum) of fully explicit solvent models, in this way setting the limits of continuum solvent models in computational homogeneous catalysis. These examples showcase that inclusion of solvent molecules in the calculation not only can improve the description of already known mechanisms but can yield new mechanistic views of a reaction. With the aim of systematizing the use of explicit solvent models, after discussing the success and limitations of continuum solvent models, issues related with solvent coordination and solvent dynamics, solvent effects in reactions involving small, charged species, as well as reactions in protic solvents and the role of solvent as reagent itself are successively considered.

Charge Transfer Landscape Manifesting Structure-Rate Relationship in the Condensed Phase via Machine Learning

In this work, we develop a machine learning (ML) strategy to map molecular structure to condensed-phase charge transfer (CT) properties including CT rate constants, energy levels, electronic couplings, energy gaps, reorganization energies, and reaction free energies, which are called CT fingerprints. The CT fingerprints of selected landmark structures covering the conformation space of an organic photovoltaic molecule dissolved in explicit solvent are computed and used to train ML models using kernel ridge regression. The ML models show high predictive power with R2>0.97, and both mean absolute error and root mean square error within chemical accuracy. The CT landscape for millions of molecular dynamics sampled structures is thus constructed, which allows for instant prediction of CT rate properties given any molecular structure. The unprecedented CT landscape will shed light on real-time CT dynamics in nanoscale and mesoscale condensed-phase systems, and the optimal fabrication design for homogeneous and heterogeneous optoelectronic devices.