Adsorption of dissolved molecules onto solid surfaces can be extremely sensitive to the atomic-scale properties of the solute and surface, causing difficulties for the design of fluidic systems in industrial, medical and technological applications. In this communication, we show that the Langmuir isotherm for adsorption of a small molecule to a realistic, heterogeneous surface can be predicted from atomic structures of the molecule and surface through molecular dynamics (MD) simulations. We highlight the method by studying the adsorption of dimethyl-methylphosphonate (DMMP) to amorphous silica substrates and show that subtle differences in the atomic-scale surface properties can have drastic effects on the Langmuir isotherm. The sensitivity of the method presented is sufficient to permit the optimization of fluidic devices and to determine fundamental design rules for controlling adsorption at the nanoscale.