The interaction of muscle moment arm,knee laxity,and torque in a multi-scale musculoskeletal model of the lower limb

IntroductionMusculoskeletal modeling allows insight into the interaction of muscle force and knee joint kinematics that cannot be measured in the laboratory. However, musculoskeletal models of the lower extremity commonly use simplified representations of the knee that may limit analyses of the interaction between muscle forces and joint kinematics. The goal of this research was to demonstrate how muscle forces alter knee kinematics and consequently muscle moment arms and joint torque in a musculoskeletal model of the lower limb that includes a deformable representation of the knee.MethodsTwo musculoskeletal models of the lower limb including specimen-specific articular geometries and ligament deformability at the knee were built in a finite element framework and calibrated to match mean isometric torque data collected from 12 healthy subjects. Muscle moment arms were compared between simulations of passive knee flexion and maximum isometric knee extension and flexion. In addition, isometric torque results were compared with predictions using simplified knee models in which the deformability of the knee was removed and the kinematics at the joint were prescribed for all degrees of freedom.ResultsPeak isometric torque estimated with a deformable knee representation occurred between 45° and 60° in extension, and 45° in flexion. The maximum isometric flexion torques generated by the models with deformable ligaments were 14.6% and 17.9% larger than those generated by the models with prescribed kinematics; by contrast, the maximum isometric extension torques generated by the models were similar. The change in hamstrings moment arms during isometric flexion was greater than that of the quadriceps during isometric extension (a mean RMS difference of 9.8 mm compared to 2.9 mm, respectively).DiscussionThe large changes in the moment arms of the hamstrings, when activated in a model with deformable ligaments, resulted in changes to flexion torque. When simulating human motion, the inclusion of a deformable joint in a multi-scale musculoskeletal finite element model of the lower limb may preserve the realistic interaction of muscle force with knee kinematics and torque.     

Effect of implementing magnetic resonance imaging for patient-specific OpenSim models on lower-body kinematics and knee ligament lengths

BackgroundOpenSim models are typically based on cadaver findings that are generalized to represent a wide range of populations, which curbs their validity. Patient-specific modelling through incorporating magnetic resonance imaging (MRI) improves the model’s biofidelity with respect to joint alignment and articulations, muscle wrapping, and ligament insertions. The purpose of this study was to determine if the inclusion of an MRI-based knee model would elicit differences in lower limb kinematics and resulting knee ligament lengths during a side cut task.MethodsEleven participants were analyzed with the popular Rajagopal OpenSim model, two variations of the same model to include three and six degrees of freedom knee (DOF), and a fourth version featuring a four DOF MRI-based knee model. These four models were used in an inverse kinematics analysis of a side cut task and the resulting lower limb kinematics and knee ligament lengths were analyzed.ResultsThe MRI-based model was more responsive to the movement task than the original Rajagopal model while less susceptible to soft tissue artifact than the unconstrained six DOF model. Ligament isometry was greatest in the original Rajagopal model and smallest in the six DOF model.ConclusionsWhen using musculoskeletal modelling software, one must acutely consider the model choice as the resulting kinematics and ligament lengths are dependent on this decision. The MRI-based knee model is responsive to the kinematics and ligament lengths of highly dynamic tasks and may prove to be the most valid option for continuing with late-stage modelling operations such as static optimization.     

Using Molecular Dynamics Simulations on Crambin to Evaluate the Suitability of Different Continuum Dielectric and Hydrogen Atom Models for Protein Simulations

Abstract

Molecular dynamics simulations of enzymes with enough explicit waters of solvation to realistically account for solute-solvent interactions can burden the computational resources required to perform the simulation by more than two orders of magnitude. Since enzyme simulations even with an implicit solvation model can be imposing for a supercomputer, it is important to assess the suitability of different continuum dielectric models for protein simulations. A series of 100-picosecond molecular dynamics simulations were performed on the X-ray crystal structure of the protein crambin to examine how well computed structures, obtained using seven continuum dielectric and two hydrogen atom models, agreed with the X-ray structure. The best level of agreement between computed and experimental structures was obtained using a constant dielectric of 2 and the all-hydrogen model. Continuum dielectric models of 1,1^*r, and 2^*r also led to computed structures in reasonably good agreement with the X-ray structure. In all cases, the all-hydrogen model gave better agreement than the united atom model, although, in one case, the difference was not significant. Dielectric models of 4, 80, and 4^*r with either hydrogen model yielded significantly poorer fits. It is especially noteworthy that the observed trends did not semiquantitatively converge until about 50 picoseconds into the simulations, suggesting that validation studies for protein calculations based on energy minimizations or short simulations should be viewed with caution.     

Improved multi-scale modeling of uterine electrical activity

ObjectiveWe aim to simulate the uterine electrical activity at the myometrium level and at the skin level where it can be recorded non-invasively.Material and methodsWe use 2D models both at the myometrium scale and for the volume conductor. The multi-scale model has been implemented in Python, as a complete package including the needed tools. To speedup the slowest step, we integrated parallel execution.ResultsWe obtain realistic simulations of EHG signals as recorded by an electrode array placed on the woman abdomen.ConclusionThese simulations are still generic, the next problem to address will be the identification of the model's parameters to obtain patient-specific simulations.     

A practical solution to reduce soft tissue artifact error at the knee using adaptive kinematic constraints

Musculoskeletal modeling and simulations have vast potential in clinical and research fields, but face various challenges in representing the complexities of the human body. Soft tissue artifact from skin-mounted markers may lead to non-physiological representation of joint motions being used as inputs to models in simulations. To address this, we have developed adaptive joint constraints on five of the six degree of freedom of the knee joint based on in vivo tibiofemoral joint motions recorded during walking, hopping and cutting motions from subjects instrumented with intra-cortical pins inserted into their tibia and femur. The constraint boundaries vary as a function of knee flexion angle and were tested on four whole-body models including four to six knee degrees of freedom. A musculoskeletal model developed in OpenSim simulation software was constrained to these in vivo boundaries during level gait and inverse kinematics and dynamics were then resolved. Statistical parametric mapping indicated significant differences (p < 0.05) in kinematics between bone pin constrained and unconstrained model conditions, notably in knee translations, while hip and ankle flexion/extension angles were also affected, indicating the error at the knee propagates to surrounding joints. These changes to hip, knee, and ankle kinematics led to measurable changes in hip and knee transverse plane moments, and knee frontal plane moments and forces. Since knee flexion angle can be validly represented using skin mounted markers, our tool uses this reliable measure to guide the five other degrees of freedom at the knee and provide a more valid representation of the kinematics for these degrees of freedom.     

Coating of magnetic nanoparticles affects their interactions with model cell membranes

BackgroundThe use of functionalized iron oxide nanoparticles of various chemical properties and architectures offers a new promising direction in theranostic applications. The increasing applications of nanoparticles in medicine require that these engineered nanomaterials will contact human cells without damaging essential tissues. Thus, efficient delivery must be achieved, while minimizing cytotoxicity during passage through cell membranes to reach intracellular target compartments.MethodsDifferential Scanning Calorimetry (DSC), molecular modeling, and atomistic Molecular Dynamics (MD) simulations were performed for two magnetite nanoparticles coated with polyvinyl alcohol (PVA) and polyarabic acid (ARA) in order to assess their interactions with model DPPC membranes.ResultsDSC experiments showed that both nanoparticles interact strongly with DPPC lipid head groups, albeit to a different degree, which was further confirmed and quantified by MD simulations. The two systems were simulated, and dynamical and structural properties were monitored. A bimodal diffusion was observed for both nanoparticles, representing the diffusion in the water phase and in the proximity of the lipid bilayer. Nanoparticles did not enter the bilayer, but caused ordering of the head groups and reduced the area per lipid compared to the pure bilayer, with MAG-PVA interacting more strongly and being closer to the lipid bilayer.ConclusionsResults of DSC experiments and MD simulations were in excellent agreement. Our findings demonstrate that the external coating is a key factor that affects nanoparticle-membrane interactions. Magnetite nanoparticles coated with PVA and ARA did not destabilize the model membrane and can be considered promising platforms for biomedical applications.General significanceUnderstanding the physico-chemical interactions of different nanoparticle coatings in contact with model cell membranes is the first step for assessing toxic response and could lead to predictive models for estimating toxicity. DSC in combination with MD simulations is an effective strategy to assess physico-chemical interactions of coated nanoparticles with lipid bilayers.     

Development of an accurate three-dimensional finite element knee model

This paper presents the development of a detailed articulating three-dimensional finite-element model of the human knee, derived from MRI scan images. The model utilises precise material models and many contact interfaces in order to produce a realistic kinematic response. The behaviour of the model was examined within two fields of biomechanical simulations: general life and car-crash. These simulations were performed with the non-linear explicit dynamic code PAM-SAFE trade mark. The knee model produced results that compared favourably with existing literature. Such a model (together with other joint models that could be constructed using the same techniques) would be a valuable tool for examining new designs of prosthesis and mechanisms of injury.     

The effects of knee support on the sagittal lower-body joint kinematics and kinetics of deep squats

Little work has been done to examine the deep squat position (>130° sagittal knee flexion). In baseball and softball, catchers perform this squat an average of 146 times per nine-inning game. To alleviate some of the stress on their knees caused by this repetitive loading, some catchers wear foam knee supports.ObjectivesThis work quantifies the effects of knee support on lower-body joint kinematics and kinetics in the deep squat position.MethodsSubjects in this study performed the deep squat with no support, foam support, and instrumented support. In order to measure the force through the knee support, instrumented knee supports were designed and fabricated. We then developed an inverse dynamic model to incorporate the support loads. From the model, joint angles and moments were calculated for the three conditions.ResultsWith support there is a significant reduction in the sagittal moment at the knee of 43% on the dominant side and 63% on the non-dominant side compared to without support. These reductions are a result of the foam supports carrying approximately 20% of body weight on each side.ConclusionKnee support reduces the moment necessary to generate the deep squat position common to baseball catchers. Given the short moment arm of the patella femoral tendon, even small changes in moment can have a large effect in the tibial-femoral contact forces, particularly at deep squat angles. Reducing knee forces may be effective in decreasing incidence of osteochondritis dissecans.     

基于X线图像的膝关节周围原发性骨肿瘤辅助诊断的机器学习模型研究

摘要 目的：开发机器学习模型,并评估其在膝关节周围原发性骨肿瘤诊断方面的准确性。方法：本文将深度卷积神经网络（DCNN）这一深度学习方法应用于膝关节X线图像的影像组学分析,探讨其辅助诊断膝关节周围原发性骨肿瘤的临床价值。结果：该深度学习模型在区分正常与肿瘤影像方面展现出优异的诊断准确性,使用DCNN模型进行5轮测试的总体准确性为（99.8±0.4）%,而阳性预测值和阴性预测值分别为（100.0±0.0）%和（99.6±0.8）%,各个数据集的曲线下面积（AUC）分别为0.99、1.00、1.00、1.0和1.0,平均AUC为（0.998±0.004）;进一步使用DCNN模型进行了10轮测试显示其在区分良性与恶性骨肿瘤方面的总体准确性为（71.2±1.6）%,且达到了强阳性预测值（91.9±8.5）%,各个数据集的AUC分别为0.63、0.63、0.58、0.69、0.55、0.63、0.54、0.57、0.73、0.63,平均AUC为（0.62±0.06）。结论：本文是首个将人工智能技术应用于骨肿瘤诊断的X线图像影像组学分析方面的研究,人工智能影像组学模型能够帮助医生自动地快速筛查骨肿瘤,确定良性或恶性肿瘤时,阳性预测值较高。     

Numerical Models for the Simulation of Flexible Artificial Heart Valves: Part I - Computational Methods

Abstract

A numerical model of the coupled motion of a flexing surface in a high Reynolds number flow is presented for the simulation of flexible polyurethane heart valves in the aortic position. This is achieved by matching a Lagrangian dynamic leaflet model with a panel method based flow solver. The two models are coupled via the time-dependent pressure field using the unsteady Bernoulli equation.

Incorporation of sub-cycling in the dynamic model equations and fast pre conditioning techniques in the panel method solver yields efficient convergence and near real-time simulations of valve motion. The generality of dynamic model allows different material properties and/or geometries to be studied easily and interactively. This interactivity is realized by embedding the models within a design environment created using the software IRIS Explorer ^TM.

Two flow domains are developed, an infinite domain and an internal domain using conformal mapping theory. In addition bending stress on the valve is computed using a simple stress model based on spline and circle equation techniques.     

Wide step width reduces knee abduction moment of obese adults during stair negotiation

PurposeAn increased likelihood of developing obesity-related knee osteoarthritis may be associated with increased peak internal knee abduction moments (KAbM). Increases in step width (SW) may act to reduce this moment. The purpose of this study was to determine the effects of increased SW on knee biomechanics during stair negotiation of healthy-weight and obese participants.MethodsParticipants (24: 10 obese and 14 healthy-weight) used stairs and walked over level ground while walking at their preferred speed in two different SW conditions – preferred and wide (200% preferred). A 2 × 2 (group × condition) mixed model analysis of variance was performed to analyze differences between groups and conditions (p < 0.05).ResultsIncreased SW increased the loading-response peak knee extension moment during descent and level gait, decreased loading-response KAbMs, knee extension and abduction range of motion (ROM) during ascent, and knee adduction ROM during descent. Increased SW increased loading-response peak mediolateral ground reaction force (GRF), increased peak knee abduction angle during ascent, and decreased peak knee adduction angle during descent and level gait. Obese participants experienced disproportionate changes in loading-response mediolateral GRF, KAbM and peak adduction angle during level walking, and peak knee abduction angle and ROM during ascent.ConclusionIncreased SW successfully decreased loading-response peak KAbM. Implications of this finding are that increased SW may decrease medial compartment knee joint loading, decreasing pain and reducing joint deterioration. Increased SW influenced obese and healthy-weight participants differently and should be investigated further.     

Gait Biomechanics and Patient-Reported Function as Predictors of Response to a Hip Strengthening Exercise Intervention in Patients with Knee Osteoarthritis

ObjectiveMuscle strengthening exercises have been shown to improve pain and function in adults with mild-to-moderate knee osteoarthritis, but individual response rates can vary greatly. Predicting individuals who respond and those who do not is important in developing a more efficient and effective model of care for knee osteoarthritis (OA). Therefore, the purpose of this study was to use pre-intervention gait kinematics and patient-reported outcome measures to predict post-intervention response to a 6-week hip strengthening exercise intervention in patients with mild-to-moderate knee OA.MethodsThirty-nine patients with mild-to-moderate knee osteoarthritis completed a 6-week hip-strengthening program and were subgrouped as Non-Responders, Low-Responders, or High-Responders following the intervention based on their change in Knee injury Osteoarthritis Outcome Score (KOOS). Predictors of responder subgroups were retrospectively determined from baseline patient-reported outcome measures and kinematic gait parameters in a discriminant analysis of principal components. A 3–4 year follow-up on 16 of the patients with knee OA was also done to examine long-term changes in these parameters.ResultsA unique combination of patient-reported outcome measures and kinematic factors was able to successfully subgroup patients with knee osteoarthritis with a cross-validated classification accuracy of 85.4%. Lower patient-reported function in daily living (ADL) scores and hip frontal plane kinematics during the loading response were most important in classifying High-Responders from other sub-groups, while a combination of hip, knee, ankle kinematics were used to classify Non-Responders from Low-Responders.ConclusionPatient-reported outcome measures and objective biomechanical gait data can be an effective method of predicting individual treatment success to an exercise intervention. Measuring gait kinematics, along with patient-reported outcome measures in a clinical setting can be useful in helping make evidence-based decisions regarding optimal treatment for patients with knee OA.     

Validation of model-predicted tibial tray-synthetic bone relative motion in cementless total knee replacement during activities of daily living

As fixation of cementless total knee replacement components during the first 4–6 weeks after surgery is crucial to establish bony ingrowth into the porous surface, several studies have quantified implant-bone micromotion. Relative motion between the tray and bone can be measured in vitro, but the full micromotion contour map cannot typically be accessed experimentally. Finite element models have been employed to estimate the full micromotion map, but have not been directly validated over a range of loading conditions. The goal of this study was to develop and validate computational models for the prediction of tray-bone micromotion under simulated activities of daily living. Gait, stair descent and deep knee bend were experimentally evaluated on four samples of a cementless tibial tray implanted into proximal tibial Sawbones™ constructs. Measurements of the relative motion between the tray and the anterior cortical shell were collected with digital image correlation and used to validate a finite element model that replicated the experiment. Additionally, a probabilistic analysis was performed to account for experimental uncertainty and determine model sensitivity to alignment and frictional parameters. The finite element models were able to distinguish between activities and capture the experimental trends. Best-matching simulations from the probabilistic analysis matched measured displacement with an average root mean square (RMS) difference of 14.3 µm and Pearson-product correlation of 0.93, while the mean model presented an average RMS difference of 27.1 µm and a correlation of 0.8. Maximum deviations from average experimental measurements were 40.5 and 87.1 µm for the best-matching and average simulations, respectively. The computational pipeline developed in this study can facilitate and enhance pre-clinical assessment of novel implant components.     

Beer and wine consumption and risk of knee or hip osteoarthritis: a case control study

IntroductionThe aim of this study was to investigate the association between alcoholic and non-alcoholic beverages and knee or hip osteoarthritis (OA).MethodsWe conducted a case–control study of Caucasian men and women aged 45 to 86 years of age from Nottingham, UK. Cases had clinically severe symptoms and radiographic knee or hip OA; controls had no symptoms and no radiographic knee or hip OA. Exposure information was sought using interview-based questionnaires and a semi-quantitative food frequency questionnaire to assess beverage consumption at ages 21 to 50 years. Odds ratios (ORs), adjusted ORs (aORs), 95% confidence intervals (CI) and P values were estimated using logistic regression models.ResultsA total of 1,001 knee OA, 993 hip OA and 933 control participants were included in the study. Increasing beer consumption was associated with an increasing risk of OA (P for trend ≤0.001). Compared to those who did not consume beer, aORs for people who consumed 20 or more servings of beer were 1.93 (95% CI 1.26 to 2.94) and 2.15 (95% CI 1.45 to 3.19) for knee OA and hip OA, respectively. In contrast, increasing levels of wine consumption were associated with decreased likelihood of knee OA (P for trend <0.001). Compared to those who did not consume wine, aOR for knee OA among those who consumed 4 to 6 glasses of wine per week and ≥7 glasses of wine per week was 0.55 (95% CI 0.34 to 0.87) and 0.48 (95% CI 0.29 to 0.80), respectively. No association was identified between non-alcoholic beverages and knee or hip OA.ConclusionsBeer consumption appears to be a risk factor for knee and hip OA whereas consumption of wine has a negative association with knee OA. The mechanism behind these findings is speculative but warrants further study.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0534-4) contains supplementary material, which is available to authorized users.     

Proposing a Radial Basis Function and CSDM Indices to Predict the Traumatic Brain Injury Risk

BackgroundMany head injury indices and finite element (FE) head models have been proposed to predict traumatic brain injury (TBI). Although FE head models are suitable methods with high accuracy, they are computationally intensive. Head motion-based brain injury criteria are usually fast tools with lower accuracy. So, the objective of this study is to propose new criteria along with an artificial neural network model to predict TBI risks, which can be fast and accurate.MethodsFor this purpose, 250 FE head simulations have been carried out at 5 magnitudes and 50 rotational impact directions using the SIMon model. The effects of directions and magnitudes of rotational impacts were assessed for cumulative strain damage measure (CSDM) values. Next, statistical analysis and neural network were applied to predict CSDM values.ResultsThe results of the present research showed that the direction of rotation in the sagittal and frontal planes had a considerable effect on the CSDM values. Furthermore, new brain injury indices and a radial basis function neural network have been proposed to predict CSDM values which having high correlation coefficients with SIMon responses.ConclusionsThe results of this research demonstrated that rotational impact directions should be used to develop new head injury criteria being able to predict CSDM values. However, findings of present research proved that head motion-based brain injury criteria and RBF network can be used to predict FE head model responses with high speed and accuracy.     

Testing Models of Insulin Binding in Rat Adipocytes using Network Thermodynamic Computer Simulations

Abstract

Many different models have been proposed to explain the complex binding behavior of insulin to its receptor, but a systematic comparison of models with experimental results is lacking. We have used network thermodynamic computer simulations to compare models of insulin binding against the results of several experimental tests designed to differentiate between the models. Six models of insulin binding were tested (simple, diffusion-reaction, conversion, dissociation, heterogeneous site, and two-step intramembrane) against results reported in the literature for isolated rat adipocytes. Although still a matter of experimental controversy, the criteria selected for modeling were curvilinear Scatchard plots, bi-or multi-exponential dissociation, insulin-accelerated dissociation, lack of dependence of the overall dissociation constant on receptor number, and receptor reserve. Using a given set of parameter values most appropriate for each model, none was able to account for all of the observed experimental results. This indicates both the complexity of the binding reaction and the need for further model development. The approach of using computer simulations to systematically test models against experimental results affords not only insight into the critical features of a model enabling it to pass a test, but also indicates potential experiments which might differentiate between models.     

The Rheology of 6-Propyl Duodecane and 5-Dibutyl Nonane by Non-Equilibrium Molecular Dynamics Simulations

Abstract

In recent papers, we reported non-equilibrium molecular dynamics (NEMD) simulations of planar Couette flow for liquid n- and i-butane, and liquid n-decane and 4-propyl heptane, using two collapsed atom models and an atomistically detailed model. It was found that the collapsed atom models predict the viscosities of the n-butane and n-decane quite well, and that the atomistically detailed model does not yield quantitative agreement with the viscosity of the n-alkanes or the branched alkanes, but it does have the one positive feature that the calculated viscosities of the branched alkanes are higher than these of the n-alkanes. In the present paper, we report results of NEMD simulations of planar Couette flow for liquid 6-propyl duodecane and 5-dibutyl nonane at 296 K and 0.782 g/cc, using an expanded collapsed atom model for simplicity. The strain rate dependent viscosity shows three different regions—two shear thinning ones and a Newtonian one. The slopes of the log-log plot for the branched-chain alkanes at high strain rate are different from those at intermediate strain rate, which is characterized as a rheological behavior of branched-chain alkanes. The Newtonian viscosity of the branched-chain alkanes can be extrapolated from the plateau value of the strain rate dependent viscosity at low strain rate as for straight-chain alkanes [J. Chem. Phys., 105, 1214 (1996)]. The results indicate that more-branched alkanes have a larger viscosity than less-branched C₁₇ alkanes.     

Effect of body weight support variation on muscle activities during robot assisted gait: a dynamic simulation study

Background and Objectives: While body weight support (BWS) intonation is vital during conventional gait training of neurologically challenged subjects, it is important to evaluate its effect during robot assisted gait training. In the present research we have studied the effect of BWS intonation on muscle activities during robotic gait training using dynamic simulations. Methods: Two dimensional (2-D) musculoskeletal model of human gait was developed conjointly with another 2-D model of a robotic orthosis capable of actuating hip, knee and ankle joints simultaneously. The musculoskeletal model consists of eight major muscle groups namely; soleus (SOL), gastrocnemius (GAS), tibialis anterior (TA), hamstrings (HAM), vasti (VAS), gluteus maximus (GLU), uniarticular hip flexors (iliopsoas, IP), and Rectus Femoris (RF). BWS was provided at levels of 0, 20, 40 and 60% during the simulations. In order to obtain a feasible set of muscle activities during subsequent gait cycles, an inverse dynamics algorithm along with a quadratic minimization algorithm was implemented. Results: The dynamic parameters of the robot assisted human gait such as joint angle trajectories, ground contact force (GCF), human limb joint torques and robot induced torques at different levels of BWS were derived. The patterns of muscle activities at variable BWS were derived and analysed. For most part of the gait cycle (GC) the muscle activation patterns are quite similar for all levels of BWS as is apparent from the mean of muscle activities for the complete GC. Conclusions: Effect of BWS variation during robot assisted gait on muscle activities was studied by developing dynamic simulation. It is expected that the proposed dynamic simulation approach will provide important inferences and information about the muscle function variations consequent upon a change in BWS during robot assisted gait. This information shall be quite important while investigating the influence of BWS intonation on neuromuscular parameters of interest during robotic gait training.     

A Two-Ellipsoid Model Based Upon A Gaussian Overlap Potential

Abstract

The two-ellipsoid model (TEM) is proposed as a versatile single-site model which can be used in the study of liquid crystal phases. This TEM uses two ellipsoids to describe a molecule, one ellipsoid for the geometry and the other for the interaction strengths of the molecule. The present TEM can mimic asymmetric interactions of a liquid crystal molecule by separating the center of the interaction ellipsoid from that of the geometry ellipsoid. The potential energy surfaces of the present TEMs compare favorably with those of the corresponding Gay-Berne and the site–site models.

Monte Carlo simulations with 320 particles are performed for a symmetric interaction TEM and an asymmetric interaction TEM. The asymmetric interaction TEM displays a slightly higher transition temperature than the symmetric interaction TEM indicating that asymmetric interactions can be a driving force in a phase transition. Radial and cylindrical distribution functions of two models in the isotropic phase are similar, but those in the nematic phase are quite different.