Finger Muscle Attachments for an OpenSim Upper-Extremity Model

We determined muscle attachment points for the index, middle, ring and little fingers in an OpenSim upper-extremity model. Attachment points were selected to match both experimentally measured locations and mechanical function (moment arms). Although experimental measurements of finger muscle attachments have been made, models differ from specimens in many respects such as bone segment ratio, joint kinematics and coordinate system. Likewise, moment arms are not available for all intrinsic finger muscles. Therefore, it was necessary to scale and translate muscle attachments from one experimental or model environment to another while preserving mechanical function. We used a two-step process. First, we estimated muscle function by calculating moment arms for all intrinsic and extrinsic muscles using the partial velocity method. Second, optimization using Simulated Annealing and Hooke-Jeeves algorithms found muscle-tendon paths that minimized root mean square (RMS) differences between experimental and modeled moment arms. The partial velocity method resulted in variance accounted for (VAF) between measured and calculated moment arms of 75.5% on average (range from 48.5% to 99.5%) for intrinsic and extrinsic index finger muscles where measured data were available. RMS error between experimental and optimized values was within one standard deviation (S.D) of measured moment arm (mean RMS error = 1.5 mm < measured S.D = 2.5 mm). Validation of both steps of the technique allowed for estimation of muscle attachment points for muscles whose moment arms have not been measured. Differences between modeled and experimentally measured muscle attachments, averaged over all finger joints, were less than 4.9 mm (within 7.1% of the average length of the muscle-tendon paths). The resulting non-proprietary musculoskeletal model of the human fingers could be useful for many applications, including better understanding of complex multi-touch and gestural movements.


Introduction
Dexterous manipulation often involves complex movements of several fingers. For example, grasping or pinching can involve the coordinated activations of many hand muscles [1]. Movements similar to pinching and grasping are commonly used as inputs for human-computer interfaces (HCIs) such as a smart phone and tablet computer touch screen. The "multitouch" interfaces of these devices often require complex, multi-finger gestures or gesture sequences on the touch screen [2]. However, we have little understanding of whether the cumulative effects of long-term exposures can lead to injuries such as musculoskeletal disorders (MSDs). One way to identify exposure to forces that may lead to injury is to estimate muscle and tendon forces during repetitive activities. However, individual muscle tension or stress is difficult to measure in vivo [3,4].
Musculoskeletal models can be helpful for understanding many aspects of biomechanics and motor control. Biomechanical models can help to understand hand anatomy [5,6]. Models have also contributed to understanding upper-extremity biomechanics, including internal muscle loading during motor tasks [7][8][9][10]. Musculoskeletal models have also helped to understand neural control of hand and arm movements. For example, models helped to identify the function of the intrinsic hand [11,12] and thumb and index finger muscles [13,14], and to study movement coordination of the interphalangeal joints [15,16].
However, existing models have several limitations. First, many models are developed with proprietary or commercial software, and are not generally available for use. Second, many models are two-dimensional (2D), whereas movements are often in three-dimensional (3D) space. Third, many models focus on specific joints or sets of joints and do not include the entire upper extremity. Finally, even models that incorporate all of the major skeletal segments of the hand and arm do not currently include all of the muscles potentially responsible for actuating complex finger motions.
Muscle attachment points have been experimentally measured for both intrinsic and extrinsic muscles of the hand [5]. However, published attachments cannot be directly used in the OpenSim model [9]. The 50th percentile male used for the OpenSim model differs from the specimens used for experimental measurements. For example, the OpenSim model and experimental specimens have different sized segments, segment proportions, joint center positions, rotational axis orientations and anatomical coordinate systems. Therefore, directly adding experimentally measured muscle attachments to the OpenSim model results in different modeled moment arms than those experimentally measured for the same muscles. Muscle moment arm is an important functional measure because it determines the joint torques that result from muscle forces. Therefore, to create a musculoskeletal model capable of evaluating muscle function, it is necessary to determine sets of muscle attachments that result in moment arms that are representative of human subjects.
Several studies have measured moment arm in vivo [17][18][19] and in situ [20][21][22][23][24]. However, to our knowledge, moment arms for all muscles of the middle, ring and little fingers have not been reported. Moreover, even when moment arms are known, the specific muscle attachments are unknown and indeterminate: many potential muscle attachments can result in similar moment arms [25].
To overcome some of these limitations, we seek to add intrinsic and extrinsic finger muscles to a non-proprietary, 3D musculoskeletal model of the upper extremity [9,26]. The OpenSim upper-extremity model currently incorporates muscles at the shoulder and elbow. Adding finger muscles is a necessary component of using the model to better understand dexterous finger movements. The OpenSim model is particularly suited to help understand the role of segment interactions and the consequences of multi-joint articulation of finger muscles. Therefore, we sought to identify muscle attachments for the intrinsic and extrinsic finger muscles within the OpenSim upper-extremity model that resulted in the functionally important characteristic of matching experimentally measured moment arms.
The procedure we used to identify muscle attachments involved two steps. First, we determined muscle moment arms from experimentally measured tendon locations [5]. We used the "partial velocity" method calculation [27,28] to reproduce anatomical moment arms from experimentally measured muscle attachment points [5]. Successful prediction was considered to be calculated moment arms within one standard deviation (S.D) of experimental measurements (S.D = 2.5 mm; [6,20]). Second, given a moment arm curve, we identified muscle-tendon paths using a data-driven optimization that is capable of determining muscle attachments whose resulting moment arms matched experimental measurements [25]. We considered successful approximation to be optimized moment arms within 10% of experimentally derived values [29,30].

Methods
We expanded an existing musculoskeletal model of upper extremity on the OpenSim platform to include finger muscles (2.3.2, Simbios, Stanford, CA; [9,31]). Currently, the shoulder, elbow, forearm, wrist, thumb and index finger are modeled with 15 joint degrees of freedom (DOFs) and 50 muscle elements. The existing model is actuated by 50 muscle compartments including extrinsic muscles of the thumb and index finger. The rotational axes and centers for the thumb joints are based on measured values [32,33], and those for the index finger joints are modeled as the long axis of cylinders fit to the articular surfaces of the metacarpal and phalangeal bones.

Musculoskeletal Model
We added custom joints to the middle, ring and little fingers of the OpenSim model because only the thumb and index finger currently have active joints [9]. Each finger was modeled to have four DOFs linking four successive bones: metacarpal bone, proximal phalange, middle phalange and distal phalange. These four bone segments were linked with three joints: metacarpophalangeal (MCP), proximal interphalangeal (PIP) and distal interphalangeal (DIP). The MCP joint was capable of flex/extension and ab/adduction as a universal joint (2 DOFs). The PIP and DIP joints were limited to flex/extension as revolute joints (1 DOF). We constrained joint angles to be within their physiological range of motion (RoM) from experimental measurements [20,23,34]. The MCP joints were modeled with a RoM of 0°(extension:-) to 90°( flexion: +) as well as 0°(adduction:-) to 30°(abduction: +). The PIP and DIP joints were also modeled with a RoM of 0°to 90°(flexion: +) and 50°(flexion: +), respectively.

Scaling Model and Moment Arms
To normalize for differences among data sets and reproduce hand skin surface from the bony segments in the OpenSim model, we assumed that all linear dimensions scaled isometrically, as found for the ratio among the length of phalanx, the width and thickness of each joint [5,25,[36][37][38][39]. We scaled measured anthropometric data to the OpenSim model dimensions to describe muscle-tendon paths within OpenSim (Table 1), then normalized moment arms to the length of the middle phalanx [5,25].

Reproducing Moment Arms from Experimentally Measured Attachments
Because experimental moment arm curves are available only for the index finger at the MCP joint, it was necessary to calculate moment arm values from measured muscle attachment points for the middle, ring and little fingers. We calculated moment arms from tendon attachment locations [5] using the "partial velocity" method [28]. The partial velocity method provides a consistent technique to compute the moment arms of muscles crossing many types of joints [27].
To validate the moment arms calculated using the partial velocity method, we compared moment arms calculated from muscle attachments for the index finger [5] to measured moment arms [20]. Experimentally measured moment arms as a function of joint angle were only available for the MCP joint of the index finger, with the exception of the FDS muscle [40]. Therefore, we limited comparisons of moment arm curves to the index finger MCP joint and the FDS muscle for all of the fingers.

Muscle Attachment Determination
We used a data-driven optimization approach to identify muscle attachment sites that resulted in moment arms that most closely matched the experimentally derived curves over the joint range of motion. We performed optimizations using the same methods that were used to find attachment locations of the index finger [25]. The optimization consists of three parts: 1) an objective function, 2) boundary conditions and 3) an inequality constraint. We defined the objective function as the root mean square (RMS) error between the experimentally derived moment arms, r j ðq ! i Þ and the model-estimated moment arms,r j ðq ! i ; x ! Þ. The optimization searched for the minimum values of RMS error (f ðx ! Þ) over the domain of attachment points (x ! ) and joint angles (q ! ) that satisfy both flex/extension and ab/adduction moment arm relationships (g j (x)). Moment arm curves are different in magnitude and shape at every joint. However, we do not currently have criteria to weight errors at different joints. Therefore, we summed RMS We defined the optimization parameters and variables as: x ! Muscle attachment points (distal and proximal; x ! 2 R 6 ) q ! i Joint angle (resolution (i) of 100 increments (m) covering the RoM) r j ðq ! i Þ Experimentally measured moment armŝ r j ðq ! i ; x ! ÞModel-estimated moment arms i Joint motions (flex/extension or ab/adduction) j Individual muscles (j = 40) lb j Lower bound (bone surface) ub j Upper bound (skin surface) ε j Maximum standard deviation of experimental moment arms g j ðx ! Þ RMS error of ab/adduction moment arms during flexion/extension movements Second, the boundary conditions constrained muscle attachments to be between the bone (the lower bound, lb j [9]) and the skin surface (the upper bound, ub j [36,38,39]).
Third, we imposed an inequality constraint on ab/adduction (g j ðx ! Þ ε j ). Muscle-tendon paths have moment arms for both flex/extension and ab/adduction at the MCP joint. However, measured moment arms for ab/adduction are more uncertain than those for flex/extension. For example, ab/adduction moment arms have higher standard deviations than flex/extension moment arms, reflecting either measurement or anatomical variability [5,20,21,23]. Moreover, ab/adduction moment arm values depend on joint angles of the MCP, PIP and DIP joints [41]. However, the specific postures used to measure ab/adduction moment arm values are not available. Therefore, we included ab/adduction in the optimization by using an inequality constraint ðg j ðx ! Þ ε j Þ to find solutions where muscle attachments ðx ! Þ resulted in moment arms within the standard deviations about experimentally measured moment arms, (g j ðx ! Þ ε j ). The standard deviations for ab/adduction, ε j , were 2.5 mm for extrinsic muscles and 1.7 mm for intrinsic muscles [20]. Our optimization was implemented using Simulated Annealing [42] and Hooke-Jeeves [43] algorithms. Both algorithms were chosen because they have been used successfully for musculoskeletal modeling [25].
Initial parameter selection can influence optimization in complex search spaces [13,25,43]. Therefore, we performed sensitivity analyses for all attachments (x ! j ). We performed optimizations using 26 different starting points (x ! 0 ). The starting points included experimental estimates attachment locations [5] and where possible existing OpenSim model attachments [9]. Optimizations from different initial conditions can discover multiple solutions, which require a criterion to select a single set of attachment points [25]. We assumed that the optimized attachment points that were nearest to the experimentally measured locations and resulted in smooth tendon paths were the most anatomically reasonable.
To determine the locations of experimentally measured attachments in the OpenSim coordinate systems, we transformed attachments from the An et al. (1979) reference frames to OpenSim. We constructed a transformation matrix that mapped translation and rotation in three dimensional (3D) Cartesian coordinates (3D Euclidean space). We approximated kinematic parameters: segment rotation (θ), segment distance (d), joint twist (α) and joint offset (a) from the OpenSim model [9]. X coordinate axes were assigned to the joint axes (-: ulnar to +: radial deviation), Y coordinate axes were assigned to the phalangeal bone (-: distal to +: proximal phalanx), and Z coordinate axes were determined by the right hand rule (-: palmar to +: dorsal direction) in the OpenSim coordinate system. Primary coordinate systems of An et al. (1979), located at the approximate center of rotation of the distal, proximal phalangeal and metacarpal heads (2, 4 and 6 reference frames), were translated to OpenSim PIP, MCP joints and middle of metacarpal bone coordinate systems, respectively. Secondary coordinate systems of An et al. (1979), a translation of the proximal systems to the centers of the concave articular surface (1, 3 and 5 reference frames) were translated to OpenSim DIP, PIP and MCP joint coordinate systems, respectively.
Attachment points closest to experimentally measured values were chosen for analysis and presentation. If several candidate attachment points had Euclidean distances to the experimentally measured point within the standard deviation of experimental measurements (28% of middle phalanx length [5]), we chose the path with the largest radius of curvature (the "smoothest" path). Curvature is the rate of change of tendon path direction (i.e. deviation of the tendon path from a straight line). The radius of curvature was determined from three successive attachments, e.g., origin, via and insertion points [25]. The tendon path with the largest radius of curvature was considered the smoothest path and chosen for analysis and presentation. To measure the agreement between model-optimized and experiment-measured moment arm curves, we estimated the variance accounted for (VAF) using the equation: We implemented optimizations in Matlab (2012a, Mathworks, Natick, MA) using the OpenSim Application Programming Interface (API; OpenSim 2.0 Doxygen) to compute moment arms.

Moment Arms Calculated from Muscle Attachments Fitted Experimentally Measured Values for the Index Finger
Moment arms calculated using the partial velocity method fitted those measured at the MCP joint of the index finger (Fig 1; Table 2). Over the entire range of motion, calculated moment arms derived from anatomical muscle attachments [5] were within one standard deviation of experimentally measured moment arms [20]. Variance accounted for (VAF) averaged 75.5% across all index finger muscles, ranging from minimum (min.) 48.2% (UI) to maximum (max.) 99.5% (FDS). VAF for the UI muscle was low because of the small value for the moment arms of this muscle. However, RMS error for UI muscle was within one standard deviation (S.D) of experimentally measured moment arms (RMS error = 2.0 mm < measured S.D = 2.1 mm), implying that the calculated moment arm was reasonable. VAF were 96.5% for extrinsic tendons and 86.9% for intrinsic muscles (flex/extension).

Moment Arms Calculated from Muscle Attachments were Reasonable for All Finger Muscles
Calculated moment arms for the middle, ring, and little finger muscles were reasonable based on moment arm magnitude and variation with joint angle. Specifically, all calculated flex/extension moment arms had an ordering of magnitude that was comparable to the ordering of finger dimensions. Moment arms for the middle finger were largest, ring or index second and little smallest (Figs 2-4). Moreover, calculated moment arms were consistent functions of joint angle. Calculated moment arms of FDP, FDS, RI, UI and LU increased with finger flexion, while those of TE, ES and EDC decreased with flexion for the index, middle, ring and little fingers. Moment arm variation with flex/extension was consistent across all finger joints.
Unlike flex/extension moment arms, ab/adduction moment arms did not show ordering that paralleled finger dimensions. For example, calculated FDS and EDC ab/adduction moment arms for the little finger were larger than those for the middle finger (Fig 3). However,  All are expressed in millimeters (mm) with flexion (+)/ extension (-) and abduction (+)/ adduction (-) at the MCP joint (degree) of the index finger. (VAF) represents Variance Accounted For, RMS represents root mean square error (Δ) between derived moment arms [5] and experimentally measured values [20], and (Error) represents percentage error, respectively.   little finger moment arms were consistent with FDS and EDC muscle locations, which were longer than measured middle finger muscle locations from the MCP joint (Euclidean distance, d little = 0.52 > d middle = 0.40 × middle phalanx length). Therefore, the partial velocity calculation predicted moment arm curves consistent with expectations based on finger anatomy.

Optimized Moment Arms Matched Calculated Values Derived from Anatomical Attachments
Optimization identified multiple muscle attachments that locally minimized RMS error (Fig 5). The average Euclidean distance from experimentally measured attachment points to the closest, most smooth optimized muscle attachments (i.e., those selected for analysis and presentation) was 4.9 mm, or 7.1% of the average length along the muscle-tendon pathways.
Model-optimized moment arms agreed with experimentally derived values for the index, middle, ring and little fingers at the MCP joint (Figs 6-9). VAF averaged 84.3, 80.4, 86.6 and 84.1% for the index, middle, ring and little fingers, respectively ( Table 3, 4, 5, 6). At the middle finger (min. VAF = 80.4%), RMS error ranged from min. 0.3 mm (UI) to max. 2.4 mm (EDC). Because S.D of moment arms has not been reported for any finger but the index finger, we compared the RMS errors to the S.D of index finger moment arms. Max. RMS error between computationally optimized and experimentally derived moment arm values was within one S. D of the experimentally measured index finger moment arm for all muscles (mean RMS error = 1.5 mm < measured S.D = 2.5 mm; [20]).
For the PIP and DIP joints, model-optimized moment arms also matched experimentally derived data (Figs 6-9). VAF averaged 67.4, 81.8, 92.2 and 62.6% for muscles at the index, Finger Muscle Attachments middle, ring and little finger PIP joints, respectively. For the DIP joints, VAF averaged 73.5, 81.4, 87.1 and 89.6% across all muscles for the index, middle, ring and little fingers, respectively. RMS error between computationally optimized and experimentally derived moment arm values was within one S.D of the experimentally measured index finger moment arm for all muscles (mean RMS error = 0.7 mm < measured S.D = 1.3 mm; [20]).

Discussion
The goal of this study was to determine attachment points for muscles of the index, middle, ring and little fingers in the OpenSim upper-extremity model. Moment arms calculated from experimentally measured muscle attachments using the partial velocity method matched experimentally measured moment arms for the index finger, where both attachments and moment arm data were available. Optimization found OpenSim muscle attachments resulting in moment arms that matched moment arms derived from experimental measurements.

Limitations
Our approach had several limitations. First, the muscle wrapping objects that we used during optimization were simple shapes that in some circumstances can result in discontinuous moment arm curves [9]. Muscle wrapping can be used to better describe muscle-tendon geometry [44,45]. However, to our knowledge quantitative, experimental data for wrapping surfaces are not available for finger muscles. Therefore, we chose to use minimal wrapping surfaces that prevented the muscle-tendon paths from penetrating the bone, but did not otherwise change the muscle-tendon paths.
Second, compiling data drawn from different reference frames and sources into a common model could introduce errors. We attempted to mitigate these potential sources of error by primarily basing our calculations on data collected from all fingers from the same specimens [5,20]. Moreover, we normalized not only by middle phalanx length for flex/extension moment arms at the PIP and DIP joints, but also normalized by MCP thickness for flex/extension moment arms and MCP width for ab/adduction values. Third, although it was possible to estimate experimentally derived muscle attachment locations in OpenSim, the coordinate transformations involved several assumptions or approximations. Attachment estimations required both scaling and coordinate transformations from 2D to 3D. Coordinate transformations are complicated by the absence of quantitative information for both 3D rotation and translation in the 2D model. We were therefore required to make assumptions about both 3D reference frame orientations and origins. However, the discovery of muscle attachments that were similar to finger anatomy suggests that our assumptions were reasonable (i.e., low ΔR = 4.9 mm, which was less than the 6.9 mm interspecimen phalanx variation observed experimentally; Tables 7-10). Moreover, optimization resulted in muscle moment arms that matched experimentally measured values. Therefore, even though the attachment sites are for a hypothetical 50 th percentile male and not for a specific specimen, muscle function is as close to experimentally measured values as feasible.
Finally, we did not conduct experiments to validate the model at a functional level [13]. Full model validation was outside of the scope of the present study because it would require the determination of muscle physiological parameters and measurements of muscle activation during Finger Muscle Attachments diverse tasks. Moreover, modeling potentially complex structures such as the extensor mechanism may require additional data, assumptions, and validation [46]. We therefore limited model validation to the anatomical level of muscle attachment locations. Future experiments will seek to functionally validate the model and resolve discrepancies among experimental measurements [20,21,40].

Moment Arms and Muscle Attachments were Reasonable
Moment arms calculated using the partial velocity method for the middle, ring and little finger muscles (whose experimentally measured moment arms are not available) were reasonable based on moment arm magnitude and variation with joint angle. Flex/extension moment arms had the same ordering of magnitude as the finger dimensions, e.g., moment arms for the middle finger largest, ring or index second, and little smallest [40]. Joint thickness or phalanx lengths influenced moment arms, e.g., moment arm magnitudes increased with phalanx lengths [20,21,35,40]. Moreover, calculated moment arm variation with joint angle was consistent among muscles at all finger joints. For example, calculated flexor moment arms increased with finger flexion; calculated extensor values decreased. Experimentally measured moment arms for extrinsic (FDP/FDS) and intrinsic (RI/UI/LU) flexors increased with flexion for the index, middle, ring and little fingers, whereas extensors (TE/ES/EDC) increased with extension as observed in previous studies [5,6,20,23,24,36]. Therefore, the partial velocity method predicted moment arm values whose relationships to joint angle agreed with experimental data.
Muscle attachments found by optimization were also reasonable. Optimization identified multiple muscle attachments that locally minimized RMS error (Fig 5). However, at least one set of muscle attachments were close to experimentally measured attachments. The average distance between modeled and measured muscle attachments fell within 7.1% of the average length of the muscle paths and 71.0% of one standard deviation of specimens' phalanx length. Moreover, the agreement of the OpenSim moment arms with both directly measured moment arms [20] and moment arms calculated from muscle attachments [5] suggested that the Open-Sim muscle attachments are reasonable within the OpenSim coordinate system.
Although optimization found moment arm curves that closely matched experimental values for most muscles, some muscles such as UI showed poorer fits. Differences in VAF may in part be related to variability in experimentally measured data. An et al. (1983) reported more variability in ab/adduction moment arm (S.D = 2.5 mm) than that of flex/extension (S.D = 2.1 mm). Ab/adduction moment arms could also be influenced by finger postures to a greater extent than flex/extension moment arms [41]. Moreover, although RMS error was the same for ab/adduction and flex/extension, VAF was lower for ab/adduction because of smaller overall moment arms than flex/extension. Ab/adduction moment arms (mean MA = 3.4 mm) were less than half the magnitude of flex/extension values (mean MA = 7.2 mm). Therefore, low VAF values for some muscles may reflect variability, experimental uncertainty, and small overall moment arm values. The optimization procedure we proposed could potentially help to translate data from one musculoskeletal model to another. We demonstrated that reasonable muscle attachments can be determined using muscle function (e.g., moment arms) as a constraint. Scaling or individualizing musculoskeletal models can be difficult. Often, scaling is based on body segment lengths or limb circumferences. For example, scaling has been based on finger middle phalanx length [5], scaling elbow moment arm by cross-sectional area [47], or other linear dimensions [48]. However, scaling alone does not account for other differences among models such as different axes of rotation [34]. Using both scaling and data-driven optimization of functional parameters can identify muscle-tendon paths that both maintain anatomical proportionality and muscle function at the same time.

Muscle Parameters for a Dynamic Model
We added only anatomical information to the OpenSim model. Musculoskeletal anatomy alone can be helpful for understanding some aspects of function such as joint stability [49]. However, to dynamically simulate movement, modeling additional aspects of muscle anatomy, physiology, and motor control will be necessary. Muscle anatomical properties can be drawn from literature values, including physiological cross-sectional area and optimal fiber length [22], pennation angle [50][51][52][53], and tendon slack length [54,55]. Likewise, physiological properties such as passive and active force-length, series elasticity, force-velocity, force-activation, maximum isometric muscle force, and optimal muscle fiber length can be modeled based on experimental measurements [9,56,57]. Muscle activity from experimental measurements can be used to predict joint loading and to better understand joint function during specific tasks [58,59]. The muscle attachments we report will allow for dynamic models of tasks involving many fingers or the entire hand.
Addition of finger muscles to the OpenSim upper-extremity model will facilitate efforts to use the model to better characterize dexterous finger movements. Coupled with experimental data, the model will allow estimates of internal musculoskeletal loading during multi-touch tasks involving many fingers. The model will therefore aid understanding of complex multitouch and gestural movements, and potentially guide the design of technologies that reduce injury risk.       The coordinate system of the OpenSim model is attached to metacarpal (secondmc), proximal (proxph2), middle (midph2) and distal (distph2) phalanges.
x, y and z components indicate radioulnar (+ points out, perpendicular to the palm plan), axial (+ points from distal to proximal side) and dorsolar (+ points up, from palm to hand side), respectively. Optimized muscle-tendon locations of OpenSim are compared with those of experimental measurements (numbers in parentheses). ΔR is the Euclidean distance between OpenSim model and experimentally measured attachment sites [5]. doi:10.1371/journal.pone.0121712.t007 Finger Muscle Attachments Parentheses () values indicate measured muscle-tendon locations [5], transformed to OpenSim coordinate. ΔR is the Euclidean distance between OpenSim model and experimentally measured attachment sites [5]. doi:10.1371/journal.pone.0121712.t008 Finger Muscle Attachments