Improved Sampling in Ab Initio Free Energy Calculations of Biomolecules at Solid–Liquid Interfaces: Tight-Binding Assessment of Charged Amino Acids on TiO2 Anatase (101)

Estimating Free-Energy Surfaces and Their Convergence from Multiple, Independent Static and History-Dependent Biased Molecular-Dynamics Simulations with Mean Force Integration

Addressing the sampling problem is central to obtaining quantitative insight from molecular dynamics simulations. Adaptive biased sampling methods, such as metadynamics, tackle this issue by perturbing the Hamiltonian of a system with a history-dependent bias potential, enhancing the exploration of the ensemble of configurations and estimating the corresponding free energy surface (FES). Nevertheless, efficiently assessing and systematically improving their convergence remains an open problem. Here, building on mean force integration (MFI), we develop and test a metric for estimating the convergence of FESs obtained by combining asynchronous, independent simulations subject to diverse biasing protocols, including static biases, different variants of metadynamics, and various combinations of static and history-dependent biases. The developed metric and the ability to combine independent simulations granted by MFI enable us to devise strategies to systematically improve the quality of FES estimates. We demonstrate our approach by computing FES and assessing the convergence of a range of systems of increasing complexity, including one- and two-dimensional analytical FESs, alanine dipeptide, a Lennard-Jones supersaturated vapor undergoing liquid droplet nucleation, and the model of a colloidal system crystallizing via a two-step mechanism. The methods presented here can be generally applied to biased simulations and are implemented in pyMFI, a publicly accessible, open-source Python library.

Adsorption of Glycine on TiO2 in Water from On-the-fly Free-Energy Calculations and In Situ Electrochemical Impedance Spectroscopy

We report here an experimental-computational study of hydrated TiO2 anatase nanoparticles interacting with glycine, where we obtain quantitative agreement of the measured adsorption free energies. Ab initio simulations are performed within the tight binding and density functional theory in combination with enhanced free-energy sampling techniques, which exploit the thermodynamic integration of the unbiased mean forces collected on-the-fly along the molecular dynamics trajectories. The experiments adopt a new and efficient setup for electrochemical impedance spectroscopy measurements based on portable screen-printed gold electrodes, which allows fast and in situ signal assessment. The measured adsorption free energy is -30 kJ/mol (both from experiment and calculation), with preferential interaction of the charged NH3+ group which strongly adsorbs on the TiO2 bridging oxygens. This highlights the importance of the terminal amino groups in the adsorption mechanism of amino acids on hydrated metal oxides. The excellent agreement between computation and experiment for this amino acid opens the doors to the exploration of the interaction free energies for other moderately complex bionano systems.

Biomolecular Adsorprion at ZnS Nanomaterials: A Molecular Dynamics Simulation Study of the Adsorption Preferences, Effects of the Surface Curvature and Coating

The understanding of interactions between nanomaterials and biological molecules is of primary importance for biomedical applications of nanomaterials, as well as for the evaluation of their possible toxic effects. Here, we carried out extensive molecular dynamics simulations of the adsorption properties of about 30 small molecules representing biomolecular fragments at ZnS surfaces in aqueous media. We computed adsorption free energies and potentials of mean force of amino acid side chain analogs, lipids, and sugar fragments to ZnS (110) crystal surface and to a spherical ZnS nanoparticle. Furthermore, we investigated the effect of poly-methylmethacrylate (PMMA) coating on the adsorption preferences of biomolecules to ZnS. We found that only a few anionic molecules: aspartic and glutamic acids side chains, as well as the anionic form of cysteine show significant binding to pristine ZnS surface, while other molecules show weak or no binding. Spherical ZnS nanoparticles show stronger binding of these molecules due to binding at the edges between different surface facets. Coating of ZnS by PMMA changes binding preferences drastically: the molecules that adsorb to a pristine ZnS surface do not adsorb on PMMA-coated surfaces, while some others, particularly hydrophobic or aromatic amino-acids, show high binding affinity due to binding to the coating. We investigate further the hydration properties of the ZnS surface and relate them to the binding preferences of biomolecules.

Supercooled liquid-like dynamics in water near a fully hydrated titania surface: Decoupling of rotational and translational diffusion

We report an ab initio molecular dynamics (MD) simulation investigating the effect of a fully hydrated surface of TiO₂ on the water dynamics. It is found that the universal relation between the rotational and translational diffusion characteristics of bulk water is broken in the water layers near the surface with the rotational diffusion demonstrating progressive retardation relative to the translational diffusion when approaching the surface. This kind of rotation-translation decoupling has so far only been observed in the supercooled liquids approaching glass transition, and its observation in water at a normal liquid temperature is of conceptual interest. This finding is also of interest for the application-significant studies of the water interaction with fully hydrated nanoparticles. We note that this is the first observation of rotation-translation decoupling in an ab initio MD simulation of water.

Atomistic Perspective on Biomolecular Adsorption on Functionalized Carbon Nanomaterials under Ambient Conditions

The use of carbon-based nanomaterials is tremendously increasing in various areas of technological, bioengineering, and biomedical applications. The functionality of carbon-based nanomaterials can be further broadened via chemical functionalization of carbon nanomaterial surfaces. On the other hand, concern is rising on possible adverse effects when nanomaterials are taken up by biological organisms. In order to contribute into understanding of interactions of carbon-based nanomaterials with biological matter, we have investigated adsorption of small biomolecules on nanomaterials using enhanced sampling molecular dynamics. The biomolecules included amino acid side chain analogues, fragments of lipids, and sugar monomers. The adsorption behavior on unstructured amorphous carbon, pristine graphene and its derivatives (such as few-layer graphene, graphene oxide, and reduced graphene oxide) as well as pristine carbon nanotubes, and those functionalized with OH-, COOH-, COO-, NH2-, and NH3+ groups was investigated with respect to surface concentration. An adsorption profile, that is, the free energy as a function of distance from the nanomaterial surfaces, was determined for each molecule and surface using the Metadynamics approach. The results were analyzed in terms of chemical specificity, surface charge, and surface concentration. It was shown that although morphology of the nanomaterial has a limited effect on the adsorption properties, functionalization of the surface by various molecular groups can drastically change the adsorption behavior that can be used in the design of nanosurfaces with highly selective adsorption properties and safe for human health and environment.