This paper presents a novel neuromusculoskeletal (NMS) model of the human lower limb that uses the electromyographic (EMG) signals from 16 muscles to estimate forces generated by 34 musculotendon actuators and the resulting joint moments at the hip, knee and ankle joints during varied contractile conditions. Our proposed methodology allows overcoming limitations on force computation shown by currently available NMS models, which constrain the operation of muscles to satisfy joint moments about one single degree of freedom (DOF) only (i.e. knee flexion-extension). The design of advanced human machine interfaces can benefit from the application of our proposed multi-DOF NMS model. The better estimates of the human internal state it provides with respect to single-DOF NMS models, will allow designing more intuitive human-machine interfaces for the simultaneous EMG-driven actuation of multiple joints in lower limb powered orthoses.