Computer simulations of alkali-acetate solutions: Accuracy of the forcefields in difference concentrations

Estimating the Gibbs Hydration Energies of Actinium and Trans-Plutonium Actinides

The use of actinides for medical, scientific and technological purposes has gained momentum in the recent years. This creates a need to understand their interactions with biomolecules, both at the interface and as they become complexed. Calculation of the Gibbs binding energies of the ions to biomolecules, i. e., the Gibbs energy change associated with a transfer of an ion from the water phase to its binding site, could help to understand the actinides' toxicities and to design agents that bind them with high affinities. To this end, there is a need to obtain accurate reference values for actinide hydration, that for most actinides are not available from experiment. In this study, a set of ionic radii is developed that enables future calculations of binding energies for Pu³⁺ and five actinides with renewed scientific and technological interest: Ac³⁺ , Am³⁺ , Cm³⁺ , Bk³⁺ and Cf³⁺ . Reference hydration energies were calculated using quantum chemistry and ion solvation theory and agree well for all ions except Ac³⁺ , where ion solvation theory seems to underestimate the magnitude of the Gibbs hydration energy. The set of radii and reference energies that are presented here provide means to calculate binding energies for actinides and biomolecules.

Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19

We present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.

Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19

We present a supercomputer-driven pipeline for in-silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. We also describe preliminary results obtained for 23 systems involving eight protein targets of the proteome of SARS CoV-2. THe MD performed is temperature replica-exchange enhanced sampling, making use of the massively parallel supercomputing on the SUMMIT supercomputer at Oak Ridge National Laboratory, with which more than 1ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to ten configurations of each of the 23 SARS CoV-2 systems using AutoDock Vina. We also demonstrate that using Autodock-GPU on SUMMIT, it is possible to perform exhaustive docking of one billion compounds in under 24 hours. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and AI methods to cluster MD trajectories and rescore docking poses.

Structure of aqueous sodium acetate solutions by X-Ray scattering and density functional theory

The structure of aq. sodium acetate solution (CH3COONa, NaOAc) was studied by X-ray scattering and density function theory (DFT). For the first hydrated layer of Na+, coordination number (CN) between Na+ and O(W, I) decreases from 5.02 ± 0.85 at 0.976 mol/L to 3.62 ± 1.21 at 4.453 mol/L. The hydration of carbonyl oxygen (OC) and hydroxyl oxygen (OOC) of CH3COO− were investigated separately and the OC shows a stronger hydration bonds comparing with OOC. With concentrations increasing, the hydration shell structures of CH3COO− are not affected by the presence of large number of ions, each CH3COO− group binds about 6.23 ± 2.01 to 7.35 ± 1.73 water molecules, which indicates a relatively strong interaction between CH3COO− and water molecules. The larger uncertainty of the CN of Na+ and OC(OOC) reflects the relative looseness of Na-OC and Na-OOC ion pairs in aq. NaOAc solutions, even at the highest concentration (4.453 mol/L), suggesting the lack of contact ion pair (CIP) formation. In aq. NaOAc solutions, the so called “structure breaking” property of Na+ and CH3COO− become effective only for the second hydration sphere of bulk water. The DFT calculations of CH3COONa (H2O)n=5–7 clusters suggest that the solvent-shared ion pair (SIP) structures appear at n = 6 and become dominant at n = 7, which is well consistent with the result from X-ray scattering.

Specific Ion and Concentration Effects in Acetate Solutions with Na+ , K+ and Cs

How salt ions affect solutes and the water beyond the solvation shell is not well understood. Molecular dynamics simulations of alkali-acetate solutions were analysed here in order to examine if, and how, different cations and solute concentrations affect the water structure and the interactions between water and acetates. The results revealed that water structure is perturbed to more than 1 nm away from the acetates and that this effect is more pronounced in physiological than in molar electrolyte concentrations. Analysis of simulations of a soluble protein revealed that the water orientation is perturbed to at least 1.5 nm from the protein structure. Furthermore, modifications to the orientation of water around carboxylate side chains were shown to depend on the local environment on the protein surface, and could extend to well over 1 nm, which may have an effect on protein dynamics during MD simulations in small water boxes.

Polarizable Force Fields for Biomolecular Simulations: Recent Advances and Applications

Realistic modeling of biomolecular systems requires an accurate treatment of electrostatics, including electronic polarization. Due to recent advances in physical models, simulation algorithms, and computing hardware, biomolecular simulations with advanced force fields at biologically relevant timescales are becoming increasingly promising. These advancements have not only led to new biophysical insights but also afforded opportunities to advance our understanding of fundamental intermolecular forces. This article describes the recent advances and applications, as well as future directions, of polarizable force fields in biomolecular simulations.