A consideration of the surrounding environment is necessary for a meaningful analysis of the reaction activity in large molecular systems. We propose an approach to perform unrestricted Hartree-Fock (UHF) calculations within the framework of the fragment molecular orbital (FMO) method (FMO-UHF) to study large systems with unpaired electrons. Prior to an energy analysis one has to optimize geometry, which requires an accurate analytic energy gradient. We derive the FMO-UHF energy and its analytic gradient and implement them into GAMESS. The performance of FMO-UHF is evaluated for a solvated organic molecule and a solvated metal complex, as well as for the active part of a protein, in terms of energy, gradient, and geometry optimization. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4737860]