Nitroaromatic compounds are important intermediates in organic synthesis. The classic method used to synthesize them is chemical nitration, which involves the use of nitric acid diluted in water or acetic acid, both harmful to the environment. With the development of green chemistry, environmental friendly enzyme catalysis is increasingly employed in chemical processes. In this work, we adopted a non-aqueous horseradish peroxidase (HRP)/NaNO2/H2O2 reaction system to study the structural characteristics of aromatic compounds potentially nitrated by enzyme catalysis, as well as the relationship between the charges on carbon atoms in benzene ring and the nitro product distribution. Investigation of various reaction parameters showed that mild reaction conditions (ambient temperature and neutral pH), plus appropriate use of H2O2 and NaNO2 could prevent inactivation of HRP and polymerization of the substrates. Compared to aqueous-organic co-solvent reaction media, the aqueous-organic two-liquid phase system had great advantages in increasing the dissolved concentration of substrate and alleviating substrate inhibition. Analysis of the aromatic compounds' structural characteristics indicated that substrates containing substituents of-NH2 or-OH were readily catalyzed. Furthermore, analysis of the relationship between natural bond orbital (NBO) charges on carbon atoms in benzene ring, as calculated by the density functional method, and the nitro product distribution characteristics, demonstrated that the favored nitration sites were the ortho and para positions of substituents in benzene ring, similar to the selectivity of chemical nitration. (c) 2015 Elsevier Inc. All rights reserved.