In Situ Loading and Time-Resolved Synchrotron-Based Phase Contrast Tomography for the Mechanical Investigation of Connective Knee Tissues: A Proof-of-Concept Study.

Articular cartilage and meniscus transfer and distribute mechanical loads in the knee joint. Degeneration of these connective tissues occurs during the progression of knee osteoarthritis, which affects their composition, microstructure, and mechanical properties. A deeper understanding of disease progression can be obtained by studying them simultaneously. Time-resolved synchrotron-based X-ray phase-contrast tomography (SR-PhC-µCT) allows to capture the tissue dynamics. This proof-of-concept study presents a rheometer setup for simultaneous in situ unconfined compression and SR-PhC-µCT of connective knee tissues. The microstructural response of bovine cartilage (n = 16) and meniscus (n = 4) samples under axial continuously increased strain, or two steps of 15% strain (stress-relaxation) is studied. The chondrocyte distribution in cartilage and the collagen fiber orientation in the meniscus are assessed. Variations in chondrocyte density reveal an increase in the top 40% of the sample during loading, compared to the lower half. Meniscus collagen fibers reorient perpendicular to the loading direction during compression and partially redisperse during relaxation. Radiation damage, image repeatability, and image quality assessments show little to no effects on the results. In conclusion, this approach is highly promising for future studies of human knee tissues to understand their microstructure, mechanical response, and progression in degenerative diseases.

Multiscale In Silico Modeling of Cartilage Injuries.

Injurious loading of the joint can be accompanied by articular cartilage damage and trigger inflammation. However, it is not well-known which mechanism controls further cartilage degradation, ultimately leading to post-traumatic osteoarthritis. For personalized prognostics, there should also be a method that can predict tissue alterations following joint and cartilage injury. This chapter gives an overview of experimental and computational methods to characterize and predict cartilage degradation following joint injury. Two mechanisms for cartilage degradation are proposed. In (1) biomechanically driven cartilage degradation, it is assumed that excessive levels of strain or stress of the fibrillar or non-fibrillar matrix lead to proteoglycan loss or collagen damage and degradation. In (2) biochemically driven cartilage degradation, it is assumed that diffusion of inflammatory cytokines leads to degradation of the extracellular matrix. When implementing these two mechanisms in a computational in silico modeling workflow, supplemented by in vitro and in vivo experiments, it is shown that biomechanically driven cartilage degradation is concentrated on the damage environment, while inflammation via synovial fluid affects all free cartilage surfaces. It is also proposed how the presented in silico modeling methodology may be used in the future for personalized prognostics and treatment planning of patients with a joint injury.

Correction: A musculoskeletal finite element model of rat knee joint for evaluating cartilage biomechanics during gait.

[This corrects the article DOI: 10.1371/journal.pcbi.1009398.].

Comprehensive Automated Routine Implementation, Validation, and Benchmark of the Anisotropic Force Field (AUA4) Using Python and GROMACS.

Molecular simulation users are sometimes discouraged from using specific molecular models because of the inconvenience of finding the force field parameters and preparing and validating the topology files. To facilitate this process and make the accurate anisotropic force field AUA4 available to molecular dynamics users, we have created and validated an automated topology and coordinate file creation routine for the GROMACS molecular simulation software. In the present work, we describe the AUA4, explain its particularities and how it was implemented, thoroughly validating the implementation, and for the first time, perform a molecular dynamics benchmark for this transferable force field. Several properties were computed, namely, liquid density, vapor pressure, and vaporization enthalpy by conducting explicit vapor-liquid interface simulations. The results evidence the correct implementation showing slight deviations from the parametrization studies. The benchmark shows the superior predictive capability of the AUA4 in recreating liquid density (RMSD equal to 17.0 kg/m3) and vaporization enthalpy (RMSD equal to 1.3 kJ/mol) compared to other transferable force fields. In addition, its superior computational time performance doubles or even triples compared to an all-atom force field such as the OPLS, depending on whether the workstation counts with GPU.

Comparison of constitutive models for meniscus and their effect on the knee joint biomechanics during gait.

Mechanical behavior of meniscus can be modeled using constitutive material models of varying complexity, such as isotropic elastic or fibril reinforced poroelastic (FRPE). However, the FRPE material is complex to implement, computationally demanding in 3D geometries, and simulation is time-consuming. Hence, we aimed to quantify the most suitable and efficient constitutive model of meniscus for simulation of cartilage responses in the knee joint during walking. We showed that simpler constitutive material models can reproduce similar cartilage responses to a knee model with the FRPE meniscus, but only knee models that consider orthotropic elastic meniscus can also reproduce meniscus responses adequately.

Injury-related cell death and proteoglycan loss in articular cartilage: Numerical model combining necrosis, reactive oxygen species, and inflammatory cytokines.

Osteoarthritis (OA) is a common musculoskeletal disease that leads to deterioration of articular cartilage, joint pain, and decreased quality of life. When OA develops after a joint injury, it is designated as post-traumatic OA (PTOA). The etiology of PTOA remains poorly understood, but it is known that proteoglycan (PG) loss, cell dysfunction, and cell death in cartilage are among the first signs of the disease. These processes, influenced by biomechanical and inflammatory stimuli, disturb the normal cell-regulated balance between tissue synthesis and degeneration. Previous computational mechanobiological models have not explicitly incorporated the cell-mediated degradation mechanisms triggered by an injury that eventually can lead to tissue-level compositional changes. Here, we developed a 2-D mechanobiological finite element model to predict necrosis, apoptosis following excessive production of reactive oxygen species (ROS), and inflammatory cytokine (interleukin-1)-driven apoptosis in cartilage explant. The resulting PG loss over 30 days was simulated. Biomechanically triggered PG degeneration, associated with cell necrosis, excessive ROS production, and cell apoptosis, was predicted to be localized near a lesion, while interleukin-1 diffusion-driven PG degeneration was manifested more globally. Interestingly, the model also showed proteolytic activity and PG biosynthesis closer to the levels of healthy tissue when pro-inflammatory cytokines were rapidly inhibited or cleared from the culture medium, leading to partial recovery of PG content. The numerical predictions of cell death and PG loss were supported by previous experimental findings. Furthermore, the simulated ROS and inflammation mechanisms had longer-lasting effects (over 3 days) on the PG content than localized necrosis. The mechanobiological model presented here may serve as a numerical tool for assessing early cartilage degeneration mechanisms and the efficacy of interventions to mitigate PTOA progression.

Adaptation of Fibril-Reinforced Poroviscoelastic Properties in Rabbit Collateral Ligaments 8 Weeks After Anterior Cruciate Ligament Transection.

Ligaments of the knee provide stability and prevent excessive motions of the joint. Rupture of the anterior cruciate ligament (ACL), a common sports injury, results in an altered loading environment for other tissues in the joint, likely leading to their mechanical adaptation. In the collateral ligaments, the patterns and mechanisms of biomechanical adaptation following ACL transection (ACLT) remain unknown. We aimed to characterize the adaptation of elastic and viscoelastic properties of the lateral and medial collateral ligaments eight weeks after ACLT. Unilateral ACLT was performed in six rabbits, and collateral ligaments were harvested from transected and contralateral knee joints after eight weeks, and from an intact control group (eight knees from four animals). The cross-sectional areas were measured with micro-computed tomography. Stepwise tensile stress-relaxation testing was conducted up to 6% final strain, and the elastic and viscoelastic properties were characterized with a fibril-reinforced poroviscoelastic material model. We found that the cross-sectional area of the collateral ligaments in the ACL transected knees increased, the nonlinear elastic collagen network modulus of the LCL decreased, and the amount of fast relaxation in the MCL decreased. Our results indicate that rupture of the ACL leads to an early adaptation of the elastic and viscoelastic properties of the collagen fibrillar network in the collateral ligaments. These adaptations may be important to consider when evaluating whole knee joint mechanics after ACL rupture, and the results aid in understanding the consequences of ACL rupture on other tissues.

A novel knee joint model in FEBio with inhomogeneous fibril-reinforced biphasic cartilage simulating tissue mechanical responses during gait: data from the osteoarthritis initiative.

We developed a novel knee joint model in FEBio to simulate walking. Knee cartilage was modeled using a fibril-reinforced biphasic (FRB) formulation with depth-wise collagen architecture and split-lines to account for cartilage structure. Under axial compression, the knee model with FRB cartilage yielded contact pressures, similar to reported experimental data. Furthermore, gait analysis with FRB cartilage simulated spatial and temporal trends in cartilage fluid pressures, stresses, and strains, comparable to those of a fibril-reinforced poroviscoelastic (FRPVE) material in Abaqus. This knee joint model in FEBio could be used for further studies of knee disorders using physiologically relevant loading.

A musculoskeletal finite element model of rat knee joint for evaluating cartilage biomechanics during gait.

Abnormal loading of the knee due to injuries or obesity is thought to contribute to the development of osteoarthritis (OA). Small animal models have been used for studying OA progression mechanisms. However, numerical models to study cartilage responses under dynamic loading in preclinical animal models have not been developed. Here we present a musculoskeletal finite element model of a rat knee joint to evaluate cartilage biomechanical responses during a gait cycle. The rat knee joint geometries were obtained from a 3-D MRI dataset and the boundary conditions regarding loading in the joint were extracted from a musculoskeletal model of the rat hindlimb. The fibril-reinforced poroelastic (FRPE) properties of the rat cartilage were derived from data of mechanical indentation tests. Our numerical results showed the relevance of simulating anatomical and locomotion characteristics in the rat knee joint for estimating tissue responses such as contact pressures, stresses, strains, and fluid pressures. We found that the contact pressure and maximum principal strain were virtually constant in the medial compartment whereas they showed the highest values at the beginning of the gait cycle in the lateral compartment. Furthermore, we found that the maximum principal stress increased during the stance phase of gait, with the greatest values at midstance. We anticipate that our approach serves as a first step towards investigating the effects of gait abnormalities on the adaptation and degeneration of rat knee joint tissues and could be used to evaluate biomechanically-driven mechanisms of the progression of OA as a consequence of joint injury or obesity.

Crack propagation in articular cartilage under cyclic loading using cohesive finite element modeling.

Severe joint injuries often involve cartilage defects that propagate after mechanical loading. The propagation of these lesions may contribute to the development of post-traumatic osteoarthritis (PTOA). However, the mechanisms behind their propagation remain unknown. Currently, no numerical predictive methods exist for estimating crack propagation in cartilage under cyclic loading, yet they would provide essential insights into crack growth in injured tissue after trauma. Here, we present a numerical approach to estimate crack propagation in articular cartilage under cyclic loading using a cohesive damage model. Four different material models for cartilage (hyperelastic, poro-hyperelastic, poro-hyper-viscoelastic, and fibril-reinforced poro-hyperelastic (FRPHE) with different collagen orientations) were implemented. Our numerical cohesive damage model was able to replicate the experimental crack length reported in the literature, showing greater crack length with an increasing number of loading cycles. Damage initiation stress (4.35-4.73 MPa) and fracture energy (0.97-1.55 N/mm) values obtained for the poro-hyperelastic, poro-hyper-viscoelastic, and parallel-FRPHE models were within the range of what has been reported previously. The crack growth predictions obtained by the FRPHE models showed the influence of anisotropy of the fibrillar matrix on the cartilage response. Our results indicate that our cohesive damage model could potentially be used to estimate the adverse conditions in injured soft tissue such as osteochondral lesions, menisci tears, or partial ligament ruptures under (ab)normal biomechanical scenarios.

Shear strain and inflammation-induced fixed charge density loss in the knee joint cartilage following ACL injury and reconstruction: A computational study.

Excessive tissue deformation near cartilage lesions and acute inflammation within the knee joint after anterior cruciate ligament (ACL) rupture and reconstruction surgery accelerate the loss of fixed charge density (FCD) and subsequent cartilage tissue degeneration. Here, we show how biomechanical and biochemical degradation pathways can predict FCD loss using a patient-specific finite element model of an ACL reconstructed knee joint exhibiting a chondral lesion. Biomechanical degradation was based on the excessive maximum shear strains that may result in cell apoptosis, while biochemical degradation was driven by the diffusion of pro-inflammatory cytokines. We found that the biomechanical model was able to predict substantial localized FCD loss near the lesion and on the medial areas of the lateral tibial cartilage. In turn, the biochemical model predicted FCD loss all around the lesion and at intact areas; the highest FCD loss was at the cartilage-synovial fluid-interface and decreased toward the deeper zones. Interestingly, simulating a downturn of an acute inflammatory response by reducing the cytokine concentration exponentially over time in synovial fluid led to a partial recovery of FCD content in the cartilage. Our novel numerical approach suggests that in vivo FCD loss can be estimated in injured cartilage following ACL injury and reconstruction. Our novel modeling platform can benefit the prediction of PTOA progression and the development of treatment interventions such as disease-modifying drug testing and rehabilitation strategies.

Prediction of local fixed charge density loss in cartilage following ACL injury and reconstruction: A computational proof-of-concept study with MRI follow-up.

The purpose of this proof-of-concept study was to develop three-dimensional patient-specific mechanobiological knee joint models to simulate alterations in the fixed charged density (FCD) around cartilage lesions during the stance phase of the walking gait. Two patients with anterior cruciate ligament (ACL) reconstructed knees were imaged at 1 and 3 years after surgery. The magnetic resonance imaging (MRI) data were used for segmenting the knee geometries, including the cartilage lesions. Based on these geometries, finite element (FE) models were developed. The gait of the patients was obtained using a motion capture system. Musculoskeletal modeling was utilized to calculate knee joint contact and lower extremity muscle forces for the FE models. Finally, a cartilage adaptation algorithm was implemented in both FE models. In the algorithm, it was assumed that excessive maximum shear and deviatoric strains (calculated as the combination of principal strains), and fluid velocity, are responsible for the FCD loss. Changes in the longitudinal T1ρ and T2 relaxation times were postulated to be related to changes in the cartilage composition and were compared with the numerical predictions. In patient 1 model, both the excessive fluid velocity and strain caused the FCD loss primarily near the cartilage lesion. T1ρ and T2 relaxation times increased during the follow-up in the same location. In contrast, in patient 2 model, only the excessive fluid velocity led to a slight FCD loss near the lesion, where MRI parameters did not show evidence of alterations. Significance: This novel proof-of-concept study suggests mechanisms through which a local FCD loss might occur near cartilage lesions. In order to obtain statistical evidence for these findings, the method should be investigated with a larger cohort of subjects.

Three-dimensional nonlinear finite element model to estimate backflow during flow-controlled infusions into the brain.

Convection-enhanced delivery is a technique to bypass the blood-brain barrier and deliver therapeutic drugs into the brain tissue. However, animal investigations and preliminary clinical trials have reported reduced efficacy to transport the infused drug in specific zones, attributed mainly to backflow, in which an annular gap is formed outside the catheter and the fluid preferentially flows toward the surface of the brain rather than through the tissue in front of the cannula tip. In this study, a three-dimensional human brain finite element model of backflow was developed to study the influence of anatomical structures during flow-controlled infusions. Predictions of backflow length were compared under the influence of ventricular pressure and the distance between the cannula and the ventricles. Simulations with zero relative ventricle pressure displayed similar backflow length predictions for larger cannula-ventricle distances. In addition, infusions near the ventricles revealed smaller backflow length and the liquid was observed to escape to the longitudinal fissure and ventricular cavities. Simulations with larger cannula-ventricle distances and nonzero relative ventricular pressure showed an increase of fluid flow through the tissue and away from the ventricles. These results reveal the importance of considering both the subject-specific anatomical details and the nonlinear effects in models focused on analyzing current and potential treatment options associated with convection-enhanced delivery optimization for future clinical trials.

Mechanobiological model for simulation of injured cartilage degradation via pro-inflammatory cytokines and mechanical stimulus.

Post-traumatic osteoarthritis (PTOA) is associated with cartilage degradation, ultimately leading to disability and decrease of quality of life. Two key mechanisms have been suggested to occur in PTOA: tissue inflammation and abnormal biomechanical loading. Both mechanisms have been suggested to result in loss of cartilage proteoglycans, the source of tissue fixed charge density (FCD). In order to predict the simultaneous effect of these degrading mechanisms on FCD content, a computational model has been developed. We simulated spatial and temporal changes of FCD content in injured cartilage using a novel finite element model that incorporates (1) diffusion of the pro-inflammatory cytokine interleukin-1 into tissue, and (2) the effect of excessive levels of shear strain near chondral defects during physiologically relevant loading. Cytokine-induced biochemical cartilage explant degradation occurs near the sides, top, and lesion, consistent with the literature. In turn, biomechanically-driven FCD loss is predicted near the lesion, in accordance with experimental findings: regions near lesions showed significantly more FCD depletion compared to regions away from lesions (p<0.01). Combined biochemical and biomechanical degradation is found near the free surfaces and especially near the lesion, and the corresponding bulk FCD loss agrees with experiments. We suggest that the presence of lesions plays a role in cytokine diffusion-driven degradation, and also predisposes cartilage for further biomechanical degradation. Models considering both these cartilage degradation pathways concomitantly are promising in silico tools for predicting disease progression, recognizing lesions at high risk, simulating treatments, and ultimately optimizing treatments to postpone the development of PTOA.

A novel mechanobiological model can predict how physiologically relevant dynamic loading causes proteoglycan loss in mechanically injured articular cartilage.

Cartilage provides low-friction properties and plays an essential role in diarthrodial joints. A hydrated ground substance composed mainly of proteoglycans (PGs) and a fibrillar collagen network are the main constituents of cartilage. Unfortunately, traumatic joint loading can destroy this complex structure and produce lesions in tissue, leading later to changes in tissue composition and, ultimately, to post-traumatic osteoarthritis (PTOA). Consequently, the fixed charge density (FCD) of PGs may decrease near the lesion. However, the underlying mechanisms leading to these tissue changes are unknown. Here, knee cartilage disks from bovine calves were injuriously compressed, followed by a physiologically relevant dynamic compression for twelve days. FCD content at different follow-up time points was assessed using digital densitometry. A novel cartilage degeneration model was developed by implementing deviatoric and maximum shear strain, as well as fluid velocity controlled algorithms to simulate the FCD loss as a function of time. Predicted loss of FCD was quite uniform around the cartilage lesions when the degeneration algorithm was driven by the fluid velocity, while the deviatoric and shear strain driven mechanisms exhibited slightly discontinuous FCD loss around cracks. Our degeneration algorithm predictions fitted well with the FCD content measured from the experiments. The developed model could subsequently be applied for prediction of FCD depletion around different cartilage lesions and for suggesting optimal rehabilitation protocols.

The effect of constitutive representations and structural constituents of ligaments on knee joint mechanics.

Ligaments provide stability to the human knee joint and play an essential role in restraining motion during daily activities. Compression-tension nonlinearity is a well-known characteristic of ligaments. Moreover, simpler material representations without this feature might give reasonable results because ligaments are primarily in tension during loading. However, the biomechanical role of different constitutive representations and their fibril-reinforced poroelastic properties is unknown. A numerical knee model which considers geometric and material nonlinearities of meniscus and cartilages was applied. Five different constitutive models for the ligaments (spring, elastic, hyperelastic, porohyperelastic, and fibril-reinforced porohyperelastic (FRPHE)) were implemented. Knee joint forces for the models with elastic, hyperelastic and porohyperelastic properties showed similar behavior throughout the stance, while the model with FRPHE properties exhibited lower joint forces during the last 50% of the stance phase. The model with ligaments as springs produced the lowest joint forces at this same stance phase. The results also showed that the fibril network contributed substantially to the knee joint forces, while the nonfibrillar matrix and fluid had small effects. Our results indicate that simpler material models of ligaments with similar properties in compression and tension can be used when the loading is directed primarily along the ligament axis in tension.

Dynamic analysis of forces in the lumbar spine during bag carrying.

OBJECTIVE: The intervertebral disc supports axial and shear forces generated during tasks such as lifting and carrying weights. The objective of this study was to determine the forces in the lumbar spine of workers carrying a bag on the head, on the shoulder and on the anterior part of the trunk. METHODS: Kinematic measurements were recorded for 10 subjects carrying bags of 10, 20 and 25 kg on each of the three aforementioned positions. A simple dynamic model implemented in a custom program was then developed to determine the lumbar forces using the accelerations and positions obtained from the kinematic analysis. RESULTS: The analyses yielded a maximum compressive force of 2338.4 ± 422 N when a 25-kg bag was carried on the anterior part of the trunk. CONCLUSION: Carrying bags on the anterior part of the trunk generated higher lumbar forces compared to those developed by carrying the bag on the head or on the shoulder. Force levels suggest that this activity represents a moderate risk for the subjects. However, future biomechanical models should be developed to analyze the cumulative effect in the discs when longer periods of time are spent in this activity.

Physical Absorption of Green House Gases in Amines: The Influence of Functionality, Structure, and Cross-Interactions.

Monte Carlo simulations were performed in the isothermal-isobaric ensemble (NPT) to calculate the Henry constants of methane (CH4), nitrous oxide (N2O), and carbon dioxide (CO2) in pure H2O, amines, and alkanolamines using the classical Lorentz-Berthelot combining rules (L-B). The Henry constants of N2O and CO2 in water are highly overestimated and motivated us to propose a new set of unlike interactions. Contrarily, the Henry constant of N2O in MEA is underestimated by around 40%, and again, a new reoptimized cross unlike parameter is able to reproduce the constant to within 10%. An analysis is given of the relationship between the physical absorption of these gases and the chemical structure or functionality of 12 molecules including amines and alkanolamines using the anisotropic united atom intermolecular potential (AUA4). Finally, the solubility of N2O in an aqueous solution of monoethanolamine (MEA) at 30% (wt) was also studied. A Henry constant within 7% of the experimental value was found by using the reoptimized parameters along with L-B to account for the MEA + H2O unlike interactions. This very good agreement without additional adjustments for the MEA + H2O system may be attributed to the good excess properties predictions found in previous works for the binary mixture (MEA + H2O). However, further work, including additional alkanolamines in aqueous solutions at several concentrations, is required to verify this particular point.

Evaluation of the friction coefficient, the radial stress, and the damage work during needle insertions into agarose gels.

Agarose hydrogels have been extensively used as a phantom material to mimic the mechanical behavior of soft biological tissues, e.g. in studies aimed to analyze needle insertions into the organs producing tissue damage. To better predict the radial stress and damage during needle insertions, this study was aimed to determine the friction coefficient between the material of commercial catheters and hydrogels. The friction coefficient, the tissue damage and the radial stress were evaluated at 0.2, 1.8, and 10mm/s velocities for 28, 30, and 32 gauge needles of outer diameters equal to 0.36, 0.31, and 0.23mm, respectively. Force measurements during needle insertions and retractions on agarose gel samples were used to analyze damage and radial stress. The static friction coefficient (0.295±0.056) was significantly higher than the dynamic (0.255±0.086). The static and dynamic friction coefficients were significantly smaller for the 0.2mm/s velocity compared to those for the other two velocities, and there was no significant difference between the friction coefficients for 1.8 and 10mm/s. Radial stress averages were 131.2±54.1, 248.3±64.2, and 804.9±164.3Pa for the insertion velocity of 0.2, 1.8, and 10mm/s, respectively. The radial stress presented a tendency to increase at higher insertion velocities and needle size, which is consistent with other studies. However, the damage work did not show to be a good predictor of tissue damage, which appears to be due to simplifications in the analytical model. Differently to other approaches, the method proposed here based on radial stress may be extended in future studies to quantity tissue damage in vivo along the entire needle track.

Molecular simulation of thermodynamic and transport properties for the H2O+NaCl system.

Molecular dynamics and Monte Carlo simulations have been carried out to obtain thermodynamic and transport properties of the binary mixture H2O+NaCl at temperatures from T = 298 to 473 K. In particular, vapor pressures, liquid densities, viscosities, and vapor-liquid interfacial tensions have been obtained as functions of pressure and salt concentration. Several previously proposed fixed-point-charge models that include either Lennard-Jones (LJ) 12-6 or exponential-6 (Exp6) functional forms to describe non-Coulombic interactions were studied. In particular, for water we used the SPC and SPC/E (LJ) models in their rigid forms, a semiflexible version of the SPC/E (LJ) model, and the Errington-Panagiotopoulos Exp6 model; for NaCl, we used the Smith-Dang and Joung-Cheatham (LJ) parameterizations as well as the Tosi-Fumi (Exp6) model. While none of the model combinations are able to reproduce simultaneously all target properties, vapor pressures are well represented using the SPC plus Joung-Cheathem model combination, and all LJ models do well for the liquid density, with the semiflexible SPC/E plus Joung-Cheatham combination being the most accurate. For viscosities, the combination of rigid SPC/E plus Smith-Dang is the best alternative. For interfacial tensions, the combination of the semiflexible SPC/E plus Smith-Dang or Joung-Cheatham gives the best results. Inclusion of water flexibility improves the mixture densities and interfacial tensions, at the cost of larger deviations for the vapor pressures and viscosities. The Exp6 water plus Tosi-Fumi salt model combination was found to perform poorly for most of the properties of interest, in particular being unable to describe the experimental trend for the vapor pressure as a function of salt concentration.