Among azobenzene derivatives, azothiophenes represent a relatively recent family of compounds that exhibit similar characteristics as dyes and photoreactive systems. Their technological applications are extensive thanks to the additional design flexibility conferred by the heteroaromatic ring. In this study, we present a comprehensive investigation of the structural and electronic properties of novel dyes derived from 3-thiophenamine, utilizing a multilevel approach. We thoroughly examined the potential energy surfaces of the E and Z isomers for three molecules, each bearing different substituents on the phenyl ring at the para position relative to the diazo group. This exploration was conducted through quantum chemistry calculations at various levels of theory, employing a continuum solvent model. Subsequently, we incorporated an explicit solvent (a dimethyl sulfoxide-water mixture) to simulate the most stable isomers using classical molecular dynamics, delivering a clear picture of the local solvation structure and intermolecular interactions. Finally, a hybrid quantum mechanics/molecular mechanics (QM/MM) approach was employed to accurately describe the evolution of the solutes' properties within their environment, accounting for finite temperature effects.