Simulating distal radius fracture strength using biomechanical tests: a modeling study examining the influence of boundary conditions

Distal radius fracture strength has been quantified using in vitro biomechanical testing. These tests are frequently performed using one of two methods: (1) load is applied directly to the embedded isolated radius or (2) load is applied through the hand with the wrist joint intact. Fracture loads established using the isolated radius method are consistently 1.5 to 3 times greater than those for the intact wrist method. To address this discrepancy, a validated finite element modeling procedure was used to predict distal radius fracture strength for 22 female forearms under boundary conditions simulating the isolated radius and intact wrist method. Predicted fracture strength was highly correlated between methods (r = 0.94; p < 0.001); however, intact wrist simulations were characterized by significantly reduced cortical shell load carriage and increased stress and strain concentrations. These changes resulted in fracture strength values less than half those predicted for the isolated radius simulations (2274 ± 824 N for isolated radius, 1124 ± 375 N for intact wrist; p < 0.001). The isolated radius method underestimated the mechanical importance of the trabecular compartment compared to the more physiologically relevant intact wrist scenario. These differences should be borne in mind when interpreting the physiologic importance of mechanical testing and simulation results.     

Biomechanical Comparison of Osteoporotic Distal Radius Fractures Fixed by Distal Locking Screws with Different Length

Objectives

To evaluate the postoperative stability of osteoporotic distal radius fractures fixed with distal locking screws with different length.

Methods

A comminuted extra-articular dorsally unstable distal radius fracture, treated with volar locking plate system, was created. The 18 specimens were randomized into 3 groups based on distal locked screws with different length: Group A had unicortical screws with 50% length to the dorsal cortex. Group B had unicortical screws with 75% length to the dorsal cortex. Group C had bicortical screws. Axial compression and bending loads were imposed on the models before and after cycling testing as well as load to clinical and catastrophic failure.

Results

Minimum change in stiffness was observed before and after fatigue for all groups. The final stiffness to bending forces was statistically similar in all groups, but stiffness to axial compression was statistically significant different: Group A approached significance with respect to groups B and C (P = 0.017, 0.009), whereas stiffness in group B and C was statistically similar (P = 0.93). Load to clinical failure was significantly less for group A (456.54±78.59 N) compared with groups B (580.24±73.85 N) and C (591.07±38.40 N). Load to catastrophic failure was statistically similar between groups, but mean values for Group A were 18% less than means for Group C.

Conclusions

The volar locking plate system fixed with unicortical locking screws with at least 75% length not only produced early stability for osteoporotic distal radius fractures, but also avoided extensor tendon complications due to dorsal screw protrusion.     

The effect of static muscle forces on the fracture strength of the intact distal radius in vitro in response to simulated forward fall impacts

The distal radius fracture (DRF) is a particularly dominant injury of the wrist, commonly resulting from a forward fall on an outstretched hand. In an attempt to reduce the prevalence, costs, and potential long-term pain/deformities associated with this injury, in vivo and in vitro investigations have sought to classify the kinematics and kinetics of DRFs. In vivo forward fall work has identified a preparatory muscle contraction that occurs in the upper extremity prior to peak impact force. The present investigation constitutes the first attempt to systematically determine the effect of static muscle forces on the fracture threshold of the distal radius in vitro. Paired human cadaveric forearm specimens were divided into two groups, one that had no muscle forces applied (i.e., right arms) and the other that had muscle forces applied to ECU, ECRL, FCU and FCR (i.e., left arms), with magnitudes based on peak muscle forces and in vivo lower bound forward fall activation patterns. The specimens were secured in a custom-built pneumatic impact loading device and subjected to incremental impacts at pre-fracture (25 J) and fracture (150  J) levels. Similar fracture forces (6565 (866) N and 8665 (5133) N), impulses (47 (6) Ns and 57 (30) Ns), and energies (152 (38) J and 144 (45) J) were observed for both groups of specimens (p>0.05). Accordingly, it is suggested that, at the magnitudes presently simulated, muscle forces have little effect on the way the distal radius responds to forward fall initiated impact loading.     

Prediction of Colles’ fracture load in human radius using cohesive finite element modeling

Osteoporotic and age-related fractures are a significant public health problem. One of the most common osteoporotic fracture sites in the aging population is distal radius. There is evidence in the literature that distal radius fractures (Colles’ fracture) are an indicative of increased risk of future spine and hip fractures. In this study, a nonlinear fracture mechanics-based finite element method is applied to human radius to assess its fracture load as a function of cortical bone geometry and material properties. Seven three-dimensional finite element models of radius were created and the fracture loads were determined by using cohesive finite element modeling which explicitly represents the crack and the fracture process zone behavior. The fracture loads found in the simulations (731–6793 N) were in the range of experimental values reported in the literature. The fracture loads predicted by the simulations decreased by 4–5% per decade based only on material level changes and by 6–20% per decade when geometrical changes were also included. Cortical polar moment of inertia at 15% distal radius showed the highest correlation to fracture load (r²=0.97). These findings demonstrate the strength of fracture mechanics-based finite element modeling and show that combining geometrical and material properties provides a better assessment of fracture risk in human radius.     

Off-axis loads cause failure of the distal radius at lower magnitudes than axial loads: a finite element analysis

Distal radius (Colles') fractures are a common fall-related injury in older adults and frequently result in long-term pain and reduced ability to perform activities of daily living. Because the occurrence of a fracture during a fall depends on both the strength of the bone and upon the kinematics and kinetics of the impact itself, we sought to understand how changes in bone mineral density (BMD) and loading direction affect the fracture strength and fracture initiation location in the distal radius. A three-dimensional finite element model of the radius, scaphoid, and lunate was used to examine changes of +/-2% and +/-4% BMD, and both axial and physiologically relevant off-axis loads on the radius. Changes in BMD resulted in similar percent changes in fracture strength. However, modifying the applied load to include dorsal and lateral components (assuming a dorsal view of the wrist, rather than an anatomic view) resulted in a 47% decrease in fracture strength (axial failure load: 2752N, off-axis: 1448N). Loading direction also influenced the fracture initiation site. Axially loaded radii failed on the medial surface immediately proximal to the styloid process. In contrast, off-axis loads, containing dorsal and lateral components, caused failure on the dorsal-lateral surface. Because the radius appears to be very sensitive to loading direction, the results suggest that much of the variability in fracture strength seen in cadaver studies may be attributed to varying boundary conditions. The results further suggest that interventions focused on reducing the incidence of Colles' fractures when falls onto the upper extremities are unavoidable may benefit from increasing the extent to which the radius is loaded along its axis.     

Experimental production of extra- and intra-articular fractures of the os calcis

Although studies have been conducted in the past to duplicate traumatic fractures of the os calcis, biomechanical force data as a function of extra- and intra-articular fractures are not available. Consequently, in this study, a dynamic single impact model was used to provide such information. Using intact human cadaver lower extremities, impact loading was applied to the plantar surface of the foot using a mini-sled pendulum equipment. The proximal tibia was fixed in polymethylmethacrylate. Following impact, pathology to the os calcis was classified into intact (no injury; 14 cases), and extra-articular (6 cases) and intra-articular (6 cases) fractures. Peak dynamic forces were used to conduct statistical analysis. Mean forces for the intact and (both) fracture groups were 4144 N (standard error, SE: 689) and 7802 N (SE: 597). Mean forces for the extra- and intra-articular fracture groups were 7445 N (SE: 711) and 8159 N (SE: 1006). The peak force influenced injury outcome (ANOVA, p<0.005). Differences in the forces were found between intact and injured specimens (p<0.01); intact specimens and specimens with extra-articular pathology (p<0.001); intact specimens and specimens with intra-articular pathology (p<0.005). The present experimental protocol, which successfully reproduced clinically relevant os calcis pathology, can be extended to accommodate other variables such as the simulation of Achilles tendon force, the inclusion of other angles of force application, and the application of the impact force to limited regions of the plantar force of the foot in order to study other injury mechanisms.     

Simulation of elbow and forearm motion in vitro using a load controlled testing apparatus

The purpose of this study was to compare passive to active testing on the kinematics of the elbow and forearm using a load-controlled testing apparatus that simulates muscle loading. Ten fresh-frozen upper extremities were tested. Active control was achieved by employing computer-controlled pneumatic actuators attached to the tendons of the brachialis, biceps, triceps, brachioradialis and pronator teres. Motion of the radius and ulna relative to the humerus was measured with an electromagnetic tracking system. Active elbow flexion produced more repeatable motion of the radius and ulna than when tested passively (p<0.05). The decrease in variability, as determined from the standard deviation of five successive trials in each specimen, was 76.5 and 58.0% for the varus-valgus and internal-external motions respectively (of the ulna relative to the humerus). The variability in flexion during simulated active forearm supination was 30.6% less than during passive testing. Thus under passive control, in the absence of stability provided by muscular loading across the joint, these uncontrolled motions produce increased variability amongst trials. The smooth and repeatable motions resulting from active control, that probably model more closely the physiologic state, appear to be beneficial in the evaluation of unconstrained kinematics of the intact elbow and forearm.     

A new methodology to measure load transfer through the forearm using multiple universal force sensors

Previous approaches to measuring forces in the forearm have made the assumption that forces acting in the radius and ulna are uniaxial near the wrist and elbow. To accurately describe forces in the forearm and the forces in the interosseous ligament, we have developed a new methodology to quantitatively determine the 3-D force vectors acting in forearm structures when a compressive load is applied to the hand. A materials testing machine equipped with a six degree-of-freedom universal force–moment sensor (UFS) was employed to apply a uniaxial compressive force to cadaveric forearms gripped at the hand and humerus. Miniature UFSs were implanted into the distal radius and proximal ulna to measure force vectors there. A 3-D digitizing device was used to measure transformations between UFS coordinate systems, utilized for calculating the force vectors in the distal ulna, proximal radius, and the interosseous ligament (IOL). This method was found to be repeatable to within 3 N, and accurate to within 2 N for force magnitudes. Computer models of the forearm, generated from CT scans, were used to visualize the force vectors in 3-D. Application of this methodology to eight forearm specimens showed that the radius carries most of the load at the wrist while force in the IOL relieves load acting in the radius at the mid-forearm. For a 136 N applied hand force, the force in the IOL was 36±21 N. Advantages of this methodology include the determination of 3-D force vectors, especially those in the IOL, as well as computer generated 3-D visualization of results.     

A dynamic simulator to evaluate distal radio-ulnar joint kinematics

In order to perform cadaveric biomechanical studies of the human forearm and distal radio-ulnar joint, a dynamic simulator has been constructed. The device is based upon a Plexiglas frame, to which the ulna is secured in a vertical orientation and the humerus in a horizontal orientation. The hand is secured in a sliding bar linkage to a stepper-motor that is used to rotate the forearm. The tendons to be loaded are connected to pneumatic actuators that provide agonist and antagonist muscle loading resulting in torque along the forearm axis. The muscle loading profiles and magnitudes are programmable as a function of the pronation-supination position and direction. A magnetic tracking system is used to collect three-dimensional kinematics data of up to four segments, in conjunction with the muscle tendon loads, forearm torque and other prescribed experimental measures. All functions are under PC control using custom software written with LabVIEW (National Instruments, Austin, TX). For the DRUJ testing, the validity of the tendon loading protocol to produce physiologic torque/rotation patterns was verified using in vivo data. The relationship of individual muscle forces to forearm torque was determined by a cadaveric study.     

Biomechanical simulations of forward fall arrests: effects of upper extremity arrest strategy,gender and aging-related declines in muscle strength

Computer simulation was used to predict the extent to which age-related muscle atrophy may adversely affect the safe arrest of a forward fall onto the arms. The biomechanical factors affecting the separate risks for wrist fracture or head impact were examined using a two-dimensional, 5-link, forward dynamic model. The hypothesis was tested in older females that age-related loss in muscular strength renders the use of the arms ineffective in arresting a forward fall without either a torso impact exceeding 0.5m/s or distal forearm loads sufficient to fracture the wrist. The results demonstrate that typical age-related decline in arm muscle strength substantially reduces the ability to arrest a forward fall without the elbows buckling and, therefore, a risk of torso and/or head impact. The model predicted that older women with below-average bone strength risk a Colles fracture when arresting typical falls, particularly with an extended arm.     

Equine cortical bone exhibits rising R-curve fracture mechanics

Previous studies of the fracture properties of cortical bone have suggested that the fracture toughness increases with crack length, which is indicative of rising R-curve behavior. Based on this indirect evidence and the similarity of bone to ceramic matrix composites, we hypothesized that bone would exhibit rising R-curve behavior in the transverse orientation and that the characteristics of the R-curves would be regionally dependent within the cortex due to variations in bone microstructure and toughening mechanisms. To test these hypotheses, we conducted R-curve experiments on specimens from equine third metacarpal bones using standard fracture mechanics testing methods. Compact type specimens from the dorsal and lateral regions in the middle of the diaphysis were oriented for crack propagation transverse to the longitudinal axis of the bone.The test results demonstrate that equine cortical bone exhibits rising R-curve behavior during transverse crack propagation as hypothesized. Statistical analyses of the crack growth initiation toughness, K0, the peak toughness, Kpeak, and the crack extension at peak toughness, deltaa, revealed significant regional differences in these characteristics. Specifically, the lateral cortex displayed higher crack growth initiation and peak toughnesses. The dorsal cortex exhibited greater crack extension at the peak of crack growth resistance. Scanning electron microscopy revealed osteon pullout on fracture surfaces from the dorsal cortex and but not in the lateral cortex. Taken together, the significant differences in R-curves and the SEM fractography indicate that the fracture mechanisms acting in equine cortical bone are regionally dependent.     

Systematic assessment of training-induced changes in corticospinal output to hand using frameless stereotaxic transcranial magnetic stimulation

Measuring changes in the characteristics of corticospinal output has become a critical part of assessing the impact of motor experience on cortical organization in both the intact and injured human brain. In this protocol we describe a method for systematically assessing training-induced changes in corticospinal output that integrates volumetric anatomical MRI with transcranial magnetic stimulation (TMS). A TMS coil is sited to a target grid superimposed onto a 3D MRI of cortex using a stereotaxic neuronavigation system. Subjects are then required to exercise the first dorsal interosseus (FDI) muscle on two different tasks for a total of 30 min. The protocol allows for reliably and repeatedly detecting changes in corticospinal output to FDI muscle in response to brief periods of motor training.     

Static and dynamic error of a biplanar videoradiography system using marker-based and markerless tracking techniques

The use of biplanar videoradiography technology has become increasingly popular for evaluating joint function in vivo. Two fundamentally different methods are currently employed to reconstruct 3D bone motions captured using this technology. Marker-based tracking requires at least three radio-opaque markers to be implanted in the bone of interest. Markerless tracking makes use of algorithms designed to match 3D bone shapes to biplanar videoradiography data. In order to reliably quantify in vivo bone motion, the systematic error of these tracking techniques should be evaluated. Herein, we present new markerless tracking software that makes use of modern GPU technology, describe a versatile method for quantifying the systematic error of a biplanar videoradiography motion capture system using independent gold standard instrumentation, and evaluate the systematic error of the W.M. Keck XROMM Facility's biplanar videoradiography system using both marker-based and markerless tracking algorithms under static and dynamic motion conditions. A polycarbonate flag embedded with 12 radio-opaque markers was used to evaluate the systematic error of the marker-based tracking algorithm. Three human cadaveric bones (distal femur, distal radius, and distal ulna) were used to evaluate the systematic error of the markerless tracking algorithm. The systematic error was evaluated by comparing motions to independent gold standard instrumentation. Static motions were compared to high accuracy linear and rotary stages while dynamic motions were compared to a high accuracy angular displacement transducer. Marker-based tracking was shown to effectively track motion to within 0.1?mm and 0.1 deg under static and dynamic conditions. Furthermore, the presented results indicate that markerless tracking can be used to effectively track rapid bone motions to within 0.15 deg for the distal aspects of the femur, radius, and ulna. Both marker-based and markerless tracking techniques were in excellent agreement with the gold standard instrumentation for both static and dynamic testing protocols. Future research will employ these techniques to quantify in vivo joint motion for high-speed upper and lower extremity impacts such as jumping, landing, and hammering.     

Radius of Paranthropus robustus from member 1, Swartkrans formation, South Africa

Recently recovered hominid postcrania from Member 1, Swartkrans Formation include the proximal and distal ends of a right radius attributed to a single individual of Paranthropus robustus. These fossils are essentially similar to Australopithecus afarensis, A. africanus, and P. boisei homologues. The head manifests an ape-like circumferentia articularis, and the distal end has prominent medial, dorsal, and lateral tubercles and a well developed brachioradialis crest, features also commonly exhibited by extant great apes. The volar set of the P. robustus radiocarpal joint, like that of Australopithecus homologues, more closely resembles the neutral condition exhibited by Homo than the greater flexion evinced by living apes. Compared with fossil and recent specimens of Homo, the configuration of the P. robustus radial head suggests enhanced stability against medial displacement during pronation and supination; the strong crest for the attachment of brachioradialis may attest to enhanced forearm flexor capability. In addition, this crest and the prominent dorsal tubercles may indicate enhanced hand extensor and, therefore, hand flexor capabilities. The differences in radial morphology between Paranthropus and Homo may relate to significant behavioral differences between these two synchronic taxa.     

Neutral zone and range of motion in the spine are greater with stepwise loading than with a continuous loading protocol. An in vitro porcine investigation

Biomechanical testing of the spine has traditionally been performed to help understand the normal function of the spine as well as to evaluate the effects of injury and surgical procedures on spinal behaviour. The overall objective of this investigation was to compare traditional stepwise loading with the recently introduced continuous loading protocol, determining the effect of loading protocol on the mechanical behaviour of the spine. For all tests, a custom spine testing machine was used to apply pure moments of flexion extension, axial rotation, and lateral bending to a maximum of 2 Nm, using six porcine cervical spine specimens (C2-C4). Motions of C2 with respect to C4 were measured with an optoelectronic camera system. Motion parameters calculated were range of motion (ROM), neutral zone (NZ), and the ratio of NZ and ROM. The continuous loading protocol had smaller values for all motion parameters in each loading direction (p<0.05). ROM for the continuous test ranged between 88% and 93% of that of stepwise for the three loading directions. The continuous protocol NZ was 56-75% of that of the stepwise test. The findings of the study demonstrate that the two loading protocols provide differing spinal behaviours.     

Computed Tomographic Imaging of Subchondral Fatigue Cracks in the Distal End of the Third Metacarpal Bone in the Thoroughbred Racehorse Can Predict Crack Micromotion in an Ex-Vivo Model

Articular stress fracture arising from the distal end of the third metacarpal bone (MC3) is a common serious injury in Thoroughbred racehorses. Currently, there is no method for predicting fracture risk clinically. We describe an ex-vivo biomechanical model in which we measured subchondral crack micromotion under compressive loading that modeled high speed running. Using this model, we determined the relationship between subchondral crack dimensions measured using computed tomography (CT) and crack micromotion. Thoracic limbs from 40 Thoroughbred racehorses that had sustained a catastrophic injury were studied. Limbs were radiographed and examined using CT. Parasagittal subchondral fatigue crack dimensions were measured on CT images using image analysis software. MC3 bones with fatigue cracks were tested using five cycles of compressive loading at -7,500N (38 condyles, 18 horses). Crack motion was recorded using an extensometer. Mechanical testing was validated using bones with 3 mm and 5 mm deep parasagittal subchondral slots that modeled naturally occurring fatigue cracks. After testing, subchondral crack density was determined histologically. Creation of parasagittal subchondral slots induced significant micromotion during loading (p<0.001). In our biomechanical model, we found a significant positive correlation between extensometer micromotion and parasagittal crack area derived from reconstructed CT images (S_R = 0.32, p<0.05). Correlations with transverse and frontal plane crack lengths were not significant. Histologic fatigue damage was not significantly correlated with crack dimensions determined by CT or extensometer micromotion. Bones with parasagittal crack area measurements above 30 mm² may have a high risk of crack propagation and condylar fracture in vivo because of crack micromotion. In conclusion, our results suggest that CT could be used to quantify subchondral fatigue crack dimensions in racing Thoroughbred horses in-vivo to assess risk of condylar fracture. Horses with parasagittal crack arrays that exceed 30 mm² may have a high risk for development of condylar fracture.     

In situ Compressive Loading and Correlative Noninvasive Imaging of the Bone-periodontal Ligament-tooth Fibrous Joint

This study demonstrates a novel biomechanics testing protocol. The advantage of this protocol includes the use of an in situ loading device coupled to a high resolution X-ray microscope, thus enabling visualization of internal structural elements under simulated physiological loads and wet conditions. Experimental specimens will include intact bone-periodontal ligament (PDL)-tooth fibrous joints. Results will illustrate three important features of the protocol as they can be applied to organ level biomechanics: 1) reactionary force vs. displacement: tooth displacement within the alveolar socket and its reactionary response to loading, 2) three-dimensional (3D) spatial configuration and morphometrics: geometric relationship of the tooth with the alveolar socket, and 3) changes in readouts 1 and 2 due to a change in loading axis, i.e. from concentric to eccentric loads. Efficacy of the proposed protocol will be evaluated by coupling mechanical testing readouts to 3D morphometrics and overall biomechanics of the joint. In addition, this technique will emphasize on the need to equilibrate experimental conditions, specifically reactionary loads prior to acquiring tomograms of fibrous joints. It should be noted that the proposed protocol is limited to testing specimens under ex vivo conditions, and that use of contrast agents to visualize soft tissue mechanical response could lead to erroneous conclusions about tissue and organ-level biomechanics.     

Load transfer analysis of the distal radius from in-vivo high-resolution CT-imaging.

Prevention of osteoporotic bone fractures requires accurate diagnostic methods to detect the increase in bone fragility at an early stage of osteoporosis. However, today's bone fracture risk prediction, primarily based on bone density measurement, is not sufficiently precise. There is increasing evidence that, in addition to bone density, also the bone microarchitecture and its mechanical loading conditions are important factors determining the fracture risk. Recently, it has been shown that new high-resolution imaging techniques in combination with new computer modeling techniques based on the finite-element (FE) method can account for these additional factors. These techniques might provide information that is more relevant for the prediction of bone fracture risk. So far, however, these new imaged-based FE techniques have not been feasible in-vivo. The objectives of this study were to quantify the load transfer through the trabecular network in a distal radius using a computer model based on in-vivo high-resolution images and to determine if common regions of fractures can be explained as a result of high tissue loading in these regions. The left distal radius and the two adjacent carpal bones of a healthy volunteer were imaged using a high-resolution three-dimensional CT system providing an isotropic resolution of 165 microm. The bone representation was converted into a FE-model that was used to calculate stresses and strains in the trabecular network. The two carpal bones were loaded using different load ratios (for each load case 1000 N in total) representing impact forces on the hand either in near-neutral position or ulnar/radial deviation. The load transfer through the trabecular network of the radius was characterized by the tissue strain energy density (SED) distribution for all load cases. It was found that the distribution of the tissue loading depends on the ratio of the forces acting on the carpal bones. For all load cases the higher SED values (on average: 0.02 +/- 0.08 (S.D.) N mm(-2)) are found in a 10 mm region adjacent to the articular surface which corresponds well with the region where Colles- or Chauffeur-fractures occur. We expect that, eventually, this new approach can lead to a better prediction of the fracture risk than methods based on bone density alone since it accounts for the bone microstructure as well as its loading conditions.     

Wrist impact velocities are smaller in forward falls than backward falls from standing

The wrist is a common fracture site for both young and older adults. The purpose of this study was to compare wrist kinematics in backward and forward falls with different fall protective responses. We carried out within-subject comparison of impact velocities and maximum velocities during descent of the distal radius among three different fall configurations: (a) backward falls with knees flexed, (b) backward falls with knees extended and (c) forward falls with knees flexed. We also examined the effect of fall configuration on fall durations, elbow flexion, trunk flexion and forearm angles at impact. Forward falls resulted in smaller impact velocities of the distal radius, longer fall duration, longer braking duration, greater elbow flexion and more horizontal landing position of the forearm compared to backward falls. The distal radius impact velocity during forward falls (1.33 m/s) was significantly lower than in backward falls, and among the backward falls the impact velocity of the flexed knee strategy (2.01 m/s) was significantly lower than the extended knee strategy (2.27 m/s). These impact velocities were significantly reduced from the maximum velocities observed during descent (forward falls=3.57 m/s, backward falls with knee flexed=3.16 m/s, backward falls with knees extended=3.52 m/s). We conclude that (1) smaller impact velocities of the wrists in forward falls could imply a lower fracture risk compared to backward falls, and (2) fall protective responses reduced wrist impact velocities in all fall directions.     

Estimating the prevalence of two or more diseases using outcomes from multiplex group testing

When screening a population for infectious diseases, pooling individual specimens (e.g., blood, swabs, urine, etc.) can provide enormous cost savings when compared to testing specimens individually. In the biostatistics literature, testing pools of specimens is commonly known as group testing or pooled testing. Although estimating a population-level prevalence with group testing data has received a large amount of attention, most of this work has focused on applications involving a single disease, such as human immunodeficiency virus. Modern methods of screening now involve testing pools and individuals for multiple diseases simultaneously through the use of multiplex assays. Hou et al. (2017, Biometrics, 73, 656–665) and Hou et al. (2020, Biostatistics, 21, 417–431) recently proposed group testing protocols for multiplex assays and derived relevant case identification characteristics, including the expected number of tests and those which quantify classification accuracy. In this article, we describe Bayesian methods to estimate population-level disease probabilities from implementing these protocols or any other multiplex group testing protocol which might be carried out in practice. Our estimation methods can be used with multiplex assays for two or more diseases while incorporating the possibility of test misclassification for each disease. We use chlamydia and gonorrhea testing data collected at the State Hygienic Laboratory at the University of Iowa to illustrate our work. We also provide an online R resource practitioners can use to implement the methods in this article.