Quantitative prediction of physical properties of liquids is important for many applications. Computational methods based on either explicit or implicit solvent models can be used to approximate thermodynamics properties of liquids. Here, we evaluate the predictive power of implicit solvent models for solvation free energy of organic molecules in organic solvents. We compared the results calculated with four generalized Born (GB) models (GB^Still, GB^HCT, GB^OBCI, and GB^OBCII), the Poisson–Boltzmann (PB) model, and the density-based solvent model SMD with previous solvation free energy calculations (Zhang et al. J. Chem. Inf. Model. 2015, 55, 1192–1201) and experimental data. The comparison indicates that both PB and GB give poor agreement with explicit solvent calculations and even worse agreement with experiments (root-mean-square deviation ≈ 15 kJ/mol). The main problem seems to be the prediction of the apolar contribution, which should include the solvent entropy. The quantum mechanical-based SMD model gives significantly better agreement with experimental data than do PB or GB, but it is not as good as explicit solvent calculation results. The dielectric constant ε of the solvent is found to be a powerful predictor for the polar contribution to the free energy in implicit models; however, the Onsager relation may not hold for realistic solvent, as suggested by explicit solvent and SMD calculations. From the comparison, we also find that with an optimization of the apolar contribution, the PB model gives slightly better agreement with experiments than the SMD model, whereas the correlation between the optimized GB models and experiments remains poor. Further optimization of the apolar contribution is needed for GB models to be able to treat solvents other than water.

The structure of ribonucleic acid (RNA) polymers is strongly dependent on the presence of, in particular Mg²⁺ cations to stabilize structural features. Only in high-resolution X-ray crystallography structures can ions be identified reliably. Here, we perform molecular dynamics simulations of 24 RNA structures with varying ion concentrations. Twelve of the structures were helical and the others complex folded. The aim of the study is to predict ion positions but also to evaluate the impact of different types of ions (Na⁺ or Mg²⁺) and the ionic strength on structural stability and variations of RNA. As a general conclusion Mg²⁺ is found to conserve the experimental structure better than Na⁺ and, where experimental ion positions are available, they can be reproduced with reasonable accuracy. If a large surplus of ions is present the added electrostatic screening makes prediction of binding-sites less reproducible. Distinct differences in ion-binding between helical and complex folded structures are found. The strength of binding (ΔG^‡ for breaking RNA atom-ion interactions) is found to differ between roughly 10 and 26 kJ/mol for the different RNA atoms. Differences in stability between helical and complex folded structures and of the influence of metal ions on either are discussed.

Data quality as well as library size are crucial issues for force field development. In order to predict molecular properties in a large chemical space, the foundation to build force fields on needs to encompass a large variety of chemical compounds. The tabulated molecular physicochemical properties also need to be accurate. Due to the limited transparency in data used for development of existing force fields it is hard to establish data quality and reusability is low. This paper presents the Alexandria library as an open and freely accessible database of optimized molecular geometries, frequencies, electrostatic moments up to the hexadecupole, electrostatic potential, polarizabilities, and thermochemistry, obtained from quantum chemistry calculations for 2704 compounds. Values are tabulated and where available compared to experimental data. This library can assist systematic development and training of empirical force fields for a broad range of molecules.