The rotation of a magnetised accreting neutron star (NS) in a binary system is described by its spin period and two angles: spin inclination 𝛼 with respect to the orbital momentum and magnetic angle 𝜒 between the spin and the magnetic moment. Magnetospheric accretion spins the NS up and adjusts its rotation axis, decreasing 𝛼 to nearly perfect alignment. Its effect upon the magnetic angle is more subtle and relatively unstudied. In this work, we model the magnetic angle evolution of a rigid spherical accreting NS. We find that the torque spinning the NS up may affect the magnetic angle while both 𝛼 and 𝜒 significantly deviate from zero, and the spin-up torque varies with the phase of the spin period. As the rotation axis of the NS is being aligned with the spin-up torque, the magnetic axis becomes misaligned with the rotation axis. Under favourable conditions, magnetic angle may increase by Δ𝜒 ∼ 15 • − 20 • . This orthogonalisation may be an important factor in the evolution of millisecond pulsars, as it partially compensates the 𝜒 decrease potentially caused by pulsar torques. If the direction of the spin-up torque changes randomly with time, as in wind-fed high-mass X-ray binaries, both the rotation axis of the NS and its magnetic axis become involved in a non-linear random-walk evolution. The ultimate attractor of this process is a bimodal distribution in 𝜒 peaking at 𝜒 = 0 and 𝜒 = 90 • .