In general relativity without matter, the equation of motion of the metric field is described by the Hilbert's action or the Einstein tensor $G$. It's natural to lead to this conclusion once one realizes how natural the Riemannian scalar curvature is. In short,
scalar curvature ==> Einstein tensor $G$
Enter the matter field. By spacetime symmetries we have the energy-momentum tensor $T$, following Noether's theorem. This is also natural to me. In short,
symmetry ==> Energy-Momentum tensor $T$
What is not natural to me yet is: why we should couple metric and matter as Einstein did?
$$ G = T \mbox{ (Einstein equation)}$$
Of course, at the end we care about if the theory matches our observation, and as well-known it turned to be pretty good. But I'm still curious if there's some explanation that justifies, if any, why this method is canonical/natural/unique in any sense.