The fragment molecular orbital (FMO)-linear combination of molecular orbitals (LCMO) method incorporates as an efficient post-process calculation of one-electron orbitals of the whole system after the FMO total energy calculation. A straightforward way to increase the accuracy is inclusion of the trimer effect. Here, we derive a comprehensive formulation called the FMO3-LCMO method. To keep the computational costs of the trimer term low enough, we use a matrix-size reduction technique. We evaluated the accuracy and efficiency of the FMO3-LCMO scheme in model biological systems (alanine oligomer and chignolin). The results show that delocalized electronic orbitals with covalent and hydrogen bonds are better described at the trimer level, and the FMO3-LCMO method is applicable to quantitative evaluations of a wide range of frontier orbitals in large biosystems.