Detailed knowledge on the formation of biomembrane domains, their structure, composition, and physical characteristics is scarce. Despite its frequently discussed importance in signaling, e.g., in obtaining localized non-homogeneous receptor compositions in the plasma membrane, the nanometer size as well as the dynamic and transient nature of domains impede their experimental characterization. In turn, atomistic molecular dynamics (MD) simulations combine both, high spatial and high temporal resolution. Here, using microsecond atomistic MD simulations, we characterize the spontaneous and unbiased formation of nano-domains in a plasma membrane model containing phosphatidylcholine (POPC), palmitoyl-sphingomyelin (PSM), and cholesterol (Chol) in the presence or absence of the neurotransmitter serotonin at different temperatures. In the ternary mixture, highly ordered and highly disordered domains of similar composition coexist at 303 K. The distinction of domains by lipid acyl chain order gets lost at lower temperatures of 298 and 294 K, suggesting a phase transition at ambient temperature. By comparison of domain ordering and composition, we demonstrate how the domain-specific binding of the neurotransmitter serotonin results in a modified domain lipid composition and a substantial downward shift of the phase transition temperature. Our simulations thus suggest a novel mode of action of neurotransmitters possibly of importance in neuronal signal transmission.