Improving Small-Molecule Force Field Parameters in Ligand Binding Studies

Small molecules are major players of many chemical processes in diverse fields, from material science to biology. They are made by a combination of carbon and heteroatoms typically organized in system-specific structures of different complexity. This peculiarity hampers the application of standard force field parameters and their in silico study by means of atomistic simulations. Here, we combine quantum-mechanics and atomistic free-energy calculations to achieve an improved parametrization of the ligand torsion angles with respect to the state-of-the-art force fields in the paradigmatic molecular binding system benzamidine/trypsin. Funnel-Metadynamics calculations with the new parameters greatly reproduced the high-resolution crystallographic ligand binding mode and allowed a more accurate description of the binding mechanism, when the ligand might assume specific conformations to cross energy barriers. Our study impacts on future drug design investigations considering that the vast majority of marketed drugs are small-molecules.

Magnesium Force Fields for OPC Water with Accurate Solvation, Ion-Binding, and Water-Exchange Properties: Successful Transfer from SPC/E

Magnesium plays a vital role in a large variety of biological processes. To model such processes by molecular dynamics simulations, researchers rely on accurate force field parameters for Mg2+ and water. OPC is one of the most promising water models yielding an improved description of biomolecules in water. The aim of this work is to provide force field parameters for Mg2+ that lead to accurate simulation results in combination with OPC water. Using twelve different Mg2+ parameter sets, that were previously optimized with different water models, we systematically assess the transferability to OPC based on a large variety of experimental properties. The results show that the Mg2+ parameters for SPC/E are transferable to OPC and closely reproduce the experimental solvation free energy, radius of the first hydration shell, coordination number, activity derivative, and binding affinity toward the phosphate oxygen on RNA. Two optimal parameter sets are presented: MicroMg yields water exchange in OPC on the microsecond timescale in agreement with experiments. NanoMg yields accelerated exchange on the nanosecond timescale and facilitates the direct observation of ion binding events for enhanced sampling purposes.

Empirical optimization of molecular simulation force fields by Bayesian inference

The demands on the accuracy of force fields for classical molecular dynamics simulations are steadily growing as larger and more complex systems are studied over longer times. One way to meet these growing demands is to hand over the learning of force fields and their parameters to machines in a systematic (semi)automatic manner. Doing so, we can take full advantage of exascale computing, the increasing availability of experimental data, and advances in quantum mechanical computations and the calculation of experimental observables from molecular ensembles. Here, we discuss and illustrate the challenges one faces in this endeavor and explore a way forward by adapting the Bayesian inference of ensembles (BioEn) method [Hummer and Köfinger, J. Chem. Phys. (2015)] for force field parameterization. In the Bayesian inference of force fields (BioFF) method developed here, the optimization problem is regularized by a simplified prior on the force field parameters and an entropic prior acting on the ensemble. The latter compensates for the unavoidable over simplifications in the parameter prior. We determine optimal force field parameters using an iterative predictor–corrector approach, in which we run simulations, determine the reference ensemble using the weighted histogram analysis method (WHAM), and update the force field according to the BioFF posterior. We illustrate this approach for a simple polymer model, using the distance between two labeled sites as the experimental observable. By systematically resolving force field issues, instead of just reweighting a structural ensemble, the BioFF corrections extend to observables not included in ensemble reweighting. We envision future force field optimization as a formalized, systematic, and (semi)automatic machine-learning effort that incorporates a wide range of data from experiment and high-level quantum chemical calculations, and takes advantage of exascale computing resources. Graphic abstract

Predicting 1,9-Decadiene−Water Partition Coefficients Using the 3D-RISM-KH Molecular Solvation Theory

The Three-Dimensional Reference Interaction Site Model (3D-RISM) with Kovalenko−Hirata (KH) closure is applied to calculate the 1,9-Decadiene/Water partition coefficients for a diverse class of compounds. The liquid state of 1,9-Decadiene is represented with the united atom TraPPE force field parameters. The 3D-RISM-KH computed partition functions are in good agreement with the experimental results. Our computational scheme can be used for a quantitative structure partitioning prediction for decadiene-water system, which has been used in membrane-mimicking of the egg-lecithin/water permeability experiments.