Motivated by recent experimental results, we present a systematic theoretical analysis of dark-bright-soliton interactions and multiple-dark-bright-soliton complexes in atomic two-component Bose-Einstein condensates. We study analytically the interactions between two dark-bright solitons in a homogeneous condensate and then extend our considerations to the presence of the trap. We illustrate the existence of robust stationary dark-bright-soliton ``molecules,{''} composed of two or more solitons, which are formed due to the competition of the interaction forces between the dark-and bright-soliton components and the trap force. Our analysis is based on an effective equation of motion, derived for the distance between two dark-bright solitons. This equation provides equilibrium positions and characteristic oscillation frequencies of the solitons, which are found to be in good agreement with the eigenfrequencies of the anomalous modes of the system.