We present an order N method for calculating electrostatic interactions that has been integrated into the molecular dynamics portion of the TINKER Molecular Modeling package. This method, termed the Image-Charge Solvation Model (ICSM), and introduced previously by Dr. Lin et al. (1) in 2009, is a hybrid electrostatic approach that combines the strengths of both explicit and implicit representations of the solvent. In this model, a multiple-image method is used to calculate reaction fields due to the implicit solvent while the Fast Multipole Method (FMM) is used to calculate the Coulomb interactions for all charges, including the charges in the explicit solvent part.The integrated package is validated through simulations of liquid water. The results are compared with those obtained by the Particle Mesh Ewald (PME) method that is built in the TINKER package. Timing performance of TINKER with the integrated ICSM is benchmarked on bulk water as a function of the size of the system. In particular, timing analysis results show that the ICSM outperforms the PME for sufficiently large systems with the break-even point at around 30,000 particles in the simulated system. To demonstrate the capability of the package on large macromolecules, the model is used to simulate the potassium channel KcsA.