Relative solvation free energies calculated using an ab initio QM/MM-based free energy perturbation method: dependence of results on simulation length

Dopamine-Decorated TiO2 Nanoparticles in Water: A QM/MM vs an MM Description

Nanoparticle functionalization is a modern strategy in nanotechnology to build up devices for several applications. Modeling fully decorated metal oxide nanoparticles of realistic size (few nanometers) in an aqueous environment is a challenging task. In this work, we present a case study relevant for solar-light exploitation and for biomedical applications, i.e., a dopamine-functionalized TiO₂ nanoparticle (1700 atoms) in bulk water, for which we have performed an extensive comparative investigation with both MM and QM/MM approaches of the structural properties and of the conformational dynamics. We have used a combined multiscale protocol for a more efficient exploration of the complex conformational space. On the basis of the results of this study and of some QM and experimental data, we have defined strengths and limitations of the existing force field parameters. Our findings will be useful for an improved modeling and simulation of many other similar hybrid bioinorganic nanosystems in an aqueous environment that are pivotal in a broad range of nanotechnological applications.

λ-Meta Dynamics Approach To Compute Absolute Solvation Free Energy

We present a new approach to combine λ dynamics with meta-dynamics (named λ-meta dynamics) to compute free energy surface with respect to λ. Particularly, the λ-meta dynamics method extends meta-dynamics to a single virtual variable λ, i.e., the coupling parameter between solute and solvent, to compute absolute solvation free energy as an exemplary application. We demonstrate that λ-meta dynamics simulations can recover the accurate potential of mean force surface with respect to λ compared to the benchmark results from traditional λ-dynamics with umbrella sampling. The solvation free energy results for five small organic molecules from λ-meta dynamics simulations using the same filling scheme show that the statistical errors are within ±0.5 kcal/mol. The new λ-meta dynamics method is general and other variables such as order parameters to describe conformational changes can be easily combined with λ-meta dynamics. This should allow for efficient samplings on high-dimension free energy landscapes.

Minimum MD simulation length required to achieve reliable results in free energy perturbation calculations: case study of relative binding free energies of fructose-1,6-bisphosphatase inhibitors

In an attempt to establish the criteria for the length of simulation to achieve the desired convergence of free energy calculations, two studies were carried out on chosen complexes of FBPase-AMP mimics. Calculations were performed for varied length of simulations and for different starting configurations using both conventional- and QM/MM-FEP methods. The results demonstrate that for small perturbations, 1248 ps simulation time could be regarded a reasonable yardstick to achieve convergence of the results. As the simulation time is extended, the errors associated with free energy calculations also gradually tapers off. Moreover, when starting the simulation from different initial configurations of the systems, the results are not changed significantly, when performed for 1248 ps. This study carried on FBPase-AMP mimics corroborates well with our previous successful demonstration of requirement of simulation time for solvation studies, both by conventional and ab initio FEP. The establishment of aforementioned criteria of simulation length serves a useful benchmark in drug design efforts using FEP methodologies, to draw a meaningful and unequivocal conclusion.