Clay minerals are implicated in the retention of biomolecules within organic matter in many soil environments. Spectroscopic studies have proposed several mechanisms for biomolecule adsorption on clays. Here, we employ molecular dynamics simulations to investigate these mechanisms in hydrated adsorbate conformations of montmorillonite, a smectite-type clay, with ten biomolecules of varying chemistry and structure, including sugars related to cellulose and hemicellulose, lignin-related phenolic acid, and amino acids with different functional groups. Our molecular modeling captures biomolecule–clay and biomolecule–biomolecule interactions that dictate selectivity and competition in adsorption retention and interlayer nanopore trapping, which we determine experimentally by NMR and X-ray diffraction, respectively. Specific adsorbate structures are important in facilitating the electrostatic attraction and Van der Waals energies underlying the hierarchy in biomolecule adsorption. Stabilized by a network of direct and water-bridged hydrogen bonds, favorable electrostatic interactions drive this hierarchy whereby amino acids with positively charged side chains are preferentially adsorbed on the negatively charged clay surface compared to the sugars and carboxylate-rich aromatics and amino acids. With divalent metal cations, our model adsorbate conformations illustrate hydrated metal cation bridging of carboxylate-containing biomolecules to the clay surface, thus explaining divalent cation-promoted adsorption from our experimental data. Adsorption experiments with a mixture of biomolecules reveal selective inhibition in biomolecule adsorption, which our molecular modeling attributes to electrostatic biomolecule–biomolecule pairing that is more energetically favorable than the biomolecule–clay complex. In sum, our findings highlight chemical and structural features that can inform hypotheses for predicting biomolecule adsorption at water–clay interfaces.