Melanosomes, pigment granules in melanophores, play a principal role in physiological color adaptation of fish and frog. Melanophores regulate melanosome trafficking on cytoskeletal filaments to generate a range of spatiotemporal patterns. Here, we present the first comprehensive model of spatiotemporal evolution of melanosome patterns. The model encompasses both physical and biochemical aspects of melanosome dynamics. It consists of (i) a kinetic description of biochemical reactions involved in intracellular signaling, (ii) a system of macroscopic reaction-diffusion-convection equations for melanosome concentration, and (iii) a set of constitutive relationships for coupling transport with the biochemical network. The model relates molecular-level regulatory actions to cell-level melanosome distribution, allowing unification of existing experimental observations and qualitative hypotheses into an integrated, consistent framework. The model reproduces salient features of melanosome patterns, both during transient and steady state. It gives useful insights into how cells coordinate motor-assisted transport to maintain and adapt spatial organization of intracellular organelles. In particular, we calculate the optimal transition paths from aggregation to dispersion in fish melanophores. The calculations suggest that fish melanophores optimally control intracellular signaling to maximize the efficiency of motor-assisted transport during dispersion.