Analysis and approximation of mixed-dimensional PDEs on 3D-1D domains coupled with Lagrange multipliers