This paper entails the application of the energy shaping methodology to control a flexible, elastic Cosserat rod model of a single octopus arm. The principal focus and novel contribution of this work is two-fold: (i) reduced order control oriented modeling of the realistic internal muscular architecture in an octopus arm; and (ii) incorporation of such models into the energy shaping methodology, extending our prior work by formally accounting for muscle constraints. Extension of the control scheme to the under-actuated muscle control case involves two steps: (i) design of a desired potential energy function whose static minimizer solves a given control task; and (ii) implementing the resulting energy shaping control input into the dynamic model. Due to the muscle actuator constraints, the desired potential energy function may not be arbitrarily chosen. Indeed, the desired energy must now satisfy a partial differential equation, known as the matching condition, which is derived for the infinite dimensional Hamiltonian control system. A particular solution to those matching conditions is described, paving the way to the application of energy shaping methodology. The overall control design methodology including muscle models is implemented and demonstrated in a dynamic simulation environment. Results of several bio-inspired numerical experiments involving the control of octopus arms are reported.