You can approximate a dirac delta function in that region, Implement point source in the variational formulation
Note that this will not be 100% accurate, and will work better on a finer mesh.
The optimal solution is to ensure vertices at the points you want to apply a boundary condition on.
You can Also find the Closest degree of freedom to that point, see: How can I find the closest node to an arbitrary coordinates?