A system of magnetite single-domain magnetic nanoparticles in an aqueous colloidal suspension at different temperatures is simulated. In this framework, we study the response of the magnetization of the system to the presence of a time-dependent magnetic field at certain frequency. To do that, a Hamiltonian that includes the Zeeman interaction and the uniaxial magneto-crystalline anisotropy energy is considered. The dynamics of the system is driven according to the solution of the stochastic Landau–Lifshitz–Gilbert differential equation, in combination with a torque term acting on the nanoparticle in order to consider the Brown rotation influenced by the viscosity of the solvent. Thus, both Néel and Brown rotation mechanisms are contemplated. Our results allow us to conclude that Brownian relaxation mediates the alignment of the anisotropy axes with the external field. Thus, depending on the solid or liquid state of the solvent, it is possible to determine the conditions under which a reinforcement of the magnetic anisotropy can take place in order to increase the remanence and squareness of the hysteresis loops.