Molecular solutes in nematic liquid crystals: Orientational order and electric field gradients