Metamaterial resonant structures made from arrays of superconducting lumped-circuit elements can exhibit microwave-mode spectra with left-handed dispersion for the standing-wave resonances, resulting in a high density of modes in the same frequency range where superconducting qubits are typically operated, as well as a bandgap at lower frequencies that extends down to dc. Using this regime for multimode circuit quantum electrodynamics, we perform a series of measurements of such a superconducting metamaterial resonator coupled to a flux-tunable transmon qubit. Through microwave measurements of the metamaterial, we observe the coupling of the qubit to each of the modes that it passes through. Using a separate readout resonator, we probe the qubit dispersively and characterize the qubit energy relaxation as a function of frequency, which is strongly affected by the Purcell effect in the presence of the dense mode spectrum. Additionally, we investigate the ac Stark shift of the qubit as the photon number in the various metamaterial modes is varied. Through numerical simulations, we explore designs based on this scheme with enhanced coupling so that the coupling energy between the qubit and the metamaterial modes can exceed the intermode spacing in future devices with achievable circuit parameters. The ability to tailor the dense mode spectrum through the choice of circuit parameters and manipulate the photonic state of the metamaterial through interactions with qubits makes this a promising platform for analog quantum simulation with microwave photons and quantum memories.