We present the theoretical and implementation details of the DFTB method with analytical gradients into the Gaussian package. SCF optimization performance with DIIS method is compared with modified Broyden method and Anderson methods, from which DIIS method is to be preferred. The number of geometry optimization steps can be significantly decreased with the Berny algorithm compared to the conjugated gradient method. From the comparison of geometry parameters as well as trends in the change of HOMO and LUMO energies, the results show clearly that DFTB can reproduce values obtained at the B3LYP/6-31G(d) level well. We also investigated the DFTB method as the low level quantum mechanical (QM) method in an ONIOM (QM:QM') description for all amino acid model system and found DFTB promising as a method to describe polarization and charge-transfer effects in large systems.