We study the adsorption process of model peptides, nucleobases,
and selected standard ligands on gold through the development of a
computational protocol based on fully atomistic classical molecular
dynamics (MD) simulations combined with umbrella sampling techniques.
The specific features of the interface components, namely, the molecule,
the metallic substrate, and the solvent, are taken into account through
different combinations of force fields (FFs), which are found to strongly
affect the results, especially changing absolute and relative adsorption
free energies and trends. Overall, noncovalent interactions drive
the process along the adsorption pathways. Our findings also show
that a suitable choice of the FF combinations can shed light on the
affinity, position, orientation, and dynamic fluctuations of the target
molecule with respect to the surface. The proposed protocol may help
the understanding of the adsorption process at the microscopic level
and may drive the in-silico design of biosensors for detection purposes.