Binding modes of diketo-acid inhibitors of HIV-1 integrase: A comparative molecular dynamics simulation study

HIV-1 integrase (IN) has become an attractive target since drug resistance against HIV-1 reverse transcriptase (RT) and protease (PR) has appeared. Diketo acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN: however the action mechanism is not well understood. Here, to study the inhibition mechanism of DKAs we performed 10 ns comparative molecular dynamics simulations on HIV-1 IN bound with three most representative DMA inhibitors: Shionogi inhibitor, S-1360 and two Merck inhibitors L-731,988 and L-708,906. Our simulations show that the acidic part of S-1360 formed salt bridge and cation-pi interactions with Lys159. In addition, the catalytic Glu152 in S-1360 was pushed away from the active site to form an ion-pair interaction with Arg199. The Merck inhibitors can maintain either one or both of these ion-pair interaction features. The difference in potencies of the DMA inhibitors is thus attributed to the different binding modes at the catalytic site. Such structural information at atomic level, not only demonstrates the action modes of DMA inhibitors but also provides a novel starting point for structural-based design of HIV-1 IN inhibitors.

Assessment of groundwater quality recovery strategies using nitrate transport modelling. Application to the Saône alluvial formations (Tournus, Saône-et-Loire)

To assess the efficiency of different agro-environmental strategies used to reduce groundwater pollution by nitrates, transport modelling in soils and groundwater has been carried out on two withdrawal areas in an alluvial plain. In a first time, the agro-environmental model AgriFlux allowed the simulation of water and nitrates fluxes flowing to groundwater. This model was calibrated for each agro-pedological unit of the studied territory. In a second time, the application of the hydrogeological model MODFLOW-MT3D allowed the simulation of nitrate transport in groundwater for the 1980-2004 period. This soil-groundwater coupled modelling has shown that soil nature is the first factor that conditions the vulnerability to nitrates. Thus, nitrate leaching occurs preferentially under sandy soils. Efficiency of different agro-environmental operations for groundwater quality recovery was quantified. The best results are obtained by combination of (1) grassland re-installation on sandy agricultural lots located in near well protection perimeter and (2) fertilization reduction on sandy agricultural lots located in the well alimentation area upstream the near protection perimeter. On other soils, the effect of grassland on groundwater quality improvement is more limited. Nevertheless, the control of nitrate fertilisation remains essential and is justified in both near and far well protection perimeters. Modelling thus allows optimising and priorizing agro-environmental actions in alluvial agricultural zones. [Comte J.-C., Banton O., Kockmann F., Villard A., Creuzot G. (2006), Assessment of groundwater quality recovery strategies using nitrate transport modelling. Application to the Saône alluvial formations (Tournus, Saône-et-Loire), Ingénieries Eau-Agriculture-Territoires, 45, 15-28]

A probabilistic modelling scheme for analysis of long-term failure of cemented femoral joint replacements

Reliable prediction of long-term medical device performance using computer simulation requires consideration of variability in surgical procedure, as well as patient-specific factors. However, even deterministic simulation of long-term failure processes for such devices is time and resource consuming so that including variability can lead to excessive time to achieve useful predictions. This study investigates the use of an accelerated probabilistic framework for predicting the likely performance envelope of a device and applies it to femoral prosthesis loosening in cemented hip arthroplasty.
A creep and fatigue damage failure model for bone cement, in conjunction with an interfacial fatigue model for the implant–cement interface, was used to simulate loosening of a prosthesis within a cement mantle. A deterministic set of trial simulations was used to account for variability of a set of surgical and patient factors, and a response surface method was used to perform and accelerate a Monte Carlo simulation to achieve an estimate of the likely range of prosthesis loosening. The proposed framework was used to conceptually investigate the influence of prosthesis selection and surgical placement on prosthesis migration.
Results demonstrate that the response surface method is capable of dramatically reducing the time to achieve convergence in mean and variance of predicted response variables. A critical requirement for realistic predictions is the size and quality of the initial training dataset used to generate the response surface and further work is required to determine the recommendations for a minimum number of initial trials. Results of this conceptual application predicted that loosening was sensitive to the implant size and femoral width. Furthermore, different rankings of implant performance were predicted when only individual simulations (e.g. an average condition) were used to rank implants, compared with when stochastic simulations were used. In conclusion, the proposed framework provides a viable approach to predicting realistic ranges of loosening behaviour for orthopaedic implants in reduced timeframes compared with conventional Monte Carlo simulations.

Modelling and Control of Shape Memory Alloy Actuators by using Preisach Model, Genetic Algorithm and Fuzzy Logic

Shapememoryalloy (SMA) actuators, which have the ability to return to a predetermined shape when heated, have many potential applications in aeronautics, surgical tools, robotics and so on. Nonlinearity hysteresis effects existing in SMA actuators present a problem in the motion control of these smart actuators. This paper investigates the control problem of SMA actuators in both simulation and experiment. In the simulation, the numerical Preisachmodel with geometrical interpretation is used for hysteresis modeling of SMA actuators. This model is then incorporated in a closed loop PID control strategy. The optimal values of PID parameters are determined by using geneticalgorithm to minimize the mean squared error between desired output displacement and simulated output. However, the control performance is not good compared with the simulation results when these parameters are applied to the real SMA control since the system is disturbed by unknown factors and changes in the surrounding environment of the system. A further automated readjustment of the PID parameters using fuzzylogic is proposed for compensating the limitation. To demonstrate the effectiveness of the proposed controller, real time control experiment results are presented.

Modelling and Instrumentation of Stretch Blow Moulding

The injection stretch blow moulding process is used to manufacture PET containers used in the soft drinks and carbonated soft drinks industry. The process consists of a test tube like specimen known as a preform which is heated, stretch and blown into a mould to form the container. This research is focused on developing a validated simulation of the process thus enabling manufacturers to design their products in a virtual environment without the need to waste time, material and energy. The simulation has been developed using the commercial FEA package Abaqus and has been validated using state of the art data acquisition system consisting of measurements for preform temperature (inner and outer wall) using a device known as THERMOscan (Figure 1), stretch rod force and velocity, internal pressure and air temperature inside the preform using an instrumented stretch rod and the?exact?timing of when the preform touches the mould wall using contact sensors.? In addition, validation studies have also been performed by blowing a perform without a mould and using high sped imaging technology in cooperation with an advanced digital image correlation system (VIC 3D) to provided new quantitative information on the behaviour of PET during blowing.? The approach has resulted in a realistic simulation in terms of accurate input parameters, preform shape evolution and prediction of final properties.

Application of a mechanobiological simulation technique to stents used clinically

Many cardiovascular diseases are characterised by the restriction of blood flow through arteries. Stents can be expanded within arteries to remove such restrictions; however, tissue in-growth into the stent can lead to restenosis. In order to predict the long-term efficacy of stenting, a mechanobiological model of the arterial tissue reaction to stress is required. In this study, a computational model of arterial tissue response to stenting is applied to three clinically relevant stent designs. We ask the question whether such a mechanobiological model can differentiate between stents used clinically, and we compare these predictions to a purely mechanical analysis. In doing so, we are testing the hypothesis that a mechanobiological model of arterial tissue response to injury could predict the long-term outcomes of stent design. Finite element analysis of the expansion of three different stent types was performed in an idealised, 3D artery. Injury was calculated in the arterial tissue using a remaining-life damage mechanics approach. The inflammatory response to this initial injury was modelled using equations governing variables which represented tissue-degrading species and growth factors. Three levels of inflammation response were modelled to account for inter-patient variability. A lattice-based model of smooth muscle cell behaviour was implemented, treating cells as discrete agents governed by local rules. The simulations predicted differences between stent designs similar to those found in vivo. It showed that the volume of neointima produced could be quantified, providing a quantitative comparison of stents. In contrast, the differences between stents based on stress alone were highly dependent on the choice of comparison criteria. These results show that the choice of stress criteria for stent comparisons is critical. This study shows that mechanobiological modelling may provide a valuable tool in stent design, allowing predictions of their long-term efficacy. The level of inflammation was shown to affect the sensitivity of the model to stent design. If this finding was verified in patients, this could suggest that high-inflammation patients may require alternative treatments to stenting.

Modelling of poly(ethylene terephthalate) in injection stretch-blow moulding

Simulations of the injection stretch-blow moulding process have been developed for the manufacture of poly(ethylene terephthalate) bottles using the commercial finite element package ABAQUS/standard. Initially a simulation of the manufacture of a 330 mL bottle was developed with three different material models (hyperelastic, creep, and a non-linear viscoelastic model (Buckley model)) to ascertain their suitability for modelling poly(ethylene terephthalate). The Buckley model was found to give results for the sidewall thickness that matched best with those measured from bottles off the production line. Following the investigation of the material models, the Buckley model was chosen to conduct a three-dimensional simulation of the manufacture of a 2 L bottle. It was found that the model was also capable of predicting the wall thickness distribution accurately for this bottle. In the development of the three-dimensional simulation a novel approach, which uses an axisymmetric model until the material reaches the petaloid base, was developed. This resulted in substantial savings in computing time. © 2000 IoM Communication Ltd.

Chemometric techniques in multivariate statistical modelling of process plant

The techniques of principal component analysis (PCA) and partial least squares (PLS) are introduced from the point of view of providing a multivariate statistical method for modelling process plants. The advantages and limitations of PCA and PLS are discussed from the perspective of the type of data and problems that might be encountered in this application area. These concepts are exemplified by two case studies dealing first with data from a continuous stirred tank reactor (CSTR) simulation and second a literature source describing a low-density polyethylene (LDPE) reactor simulation.

Flow due to multiple jets downstream of a barrage: Experiments, 3-D CFD and depth-averaged modelling

The flow through and downstream of a row of seven open draft tubes in a barrage has been investigated through laboratory experiments in a wide flume, a three-dimensional (3D) computational fluid dynamics simulation, and a two-dimensional depth-averaged computation. Agreement between the experiments and the 3D modeling is shown to be good, including the prediction of an asymmetric Coandă effect. One aim is to determine the distance downstream at which depth-averaged modeling provides a reasonable prediction; this is shown to be approximately 20 tube diameters downstream of the barrage. Upstream of this, the depth-averaged modeling inaccurately predicts water level, bed shear, and the 3D flow field. The 3D model shows that bed shear stress can be markedly magnified near the barrage, particularly where the jets become attached.

Multibody modelling of a TB31 and a TB32 crash test with vertical portable concrete barriers: model verification and sensitivity analysis.

In this article the multibody simulation software package MADYMO for analysing and optimizing occupant safety design was used to model crash tests for Normal Containment barriers in accordance with EN 1317. The verification process was carried out by simulating a TB31 and a TB32 crash test performed on vertical portable concrete barriers and by comparing the numerical results to those obtained experimentally. The same modelling approach was applied to both tests to evaluate the predictive capacity of the modelling at two different impact speeds. A sensitivity analysis of the vehicle stiffness was also carried out. The capacity to predict all of the principal EN1317 criteria was assessed for the first time: the acceleration severity index, the theoretical head impact velocity, the barrier working width and the vehicle exit box. Results showed a maximum error of 6% for the acceleration severity index and 21% for theoretical head impact velocity for the numerical simulation in comparison to the recorded data. The exit box position was predicted with a maximum error of 4°. For the working width, a large percentage difference was observed for test TB31 due to the small absolute value of the barrier deflection but the results were well within the limit value from the standard for both tests. The sensitivity analysis showed the robustness of the modelling with respect to contact stiffness increase of ±20% and ±40%. This is the first multibody model of portable concrete barriers that can reproduce not only the acceleration severity index but all the test criteria of EN 1317 and is therefore a valuable tool for new product development and for injury biomechanics research.

A scaling method for modelling the crashworthiness of novel roadside barrier designs

In this paper, a novel method for modelling a scaled vehicle–barrier crash test similar to the 20◦ angled barrier test specified in EN 1317 is reported. The intended application is for proof-of-concept evaluation of novel roadside barrier designs, and as a cost-effective precursor to full-scale testing or detailed computational modelling. The method is based on the combination of the conservation of energy law and the equation of motion of a spring mass system representing the impact, and shows, for the first time, the feasibility of applying classical scaling theories to evaluation of roadside barrier design. The scaling method is used to set the initial velocity of the vehicle in the scaled test and to provide scaling factors to convert the measured vehicle accelerations in the scaled test to predicted full-scale accelerations. These values can then be used to calculate the Acceleration Severity Index score of the barrier for a full-scale test. The theoretical validity of the method is demonstrated by comparison to numerical simulations of scaled and full-scale angled barrier impacts using multibody analysis implemented in the crash simulation software MADYMO. Results show a maximum error of 0.3% ascribable to the scaling method.