Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Integration of local-scale hydrological and regional-scale geophysical based on a nonlinear Bayesian sequential simulation approach

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the regional scale represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed a downscaling procedure based on a non-linear Bayesian sequential simulation approach. The basic objective of this algorithm is to estimate the value of the sparsely sampled hydraulic conductivity at non-sampled locations based on its relation to the electrical conductivity, which is available throughout the model space. The in situ relationship between the hydraulic and electrical conductivities is described through a non-parametric multivariate kernel density function. This method is then applied to the stochastic integration of low-resolution, re- gional-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities. Finally, the overall viability of this downscaling approach is tested and verified by performing and comparing flow and transport simulation through the original and the downscaled hydraulic conductivity fields. Our results indicate that the proposed procedure does indeed allow for obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.

Nonlinear analysis of human physical activity patterns in health and disease.

The reliable and objective assessment of chronic disease state has been and still is a very significant challenge in clinical medicine. An essential feature of human behavior related to the health status, the functional capacity, and the quality of life is the physical activity during daily life. A common way to assess physical activity is to measure the quantity of body movement. Since human activity is controlled by various factors both extrinsic and intrinsic to the body, quantitative parameters only provide a partial assessment and do not allow for a clear distinction between normal and abnormal activity. In this paper, we propose a methodology for the analysis of human activity pattern based on the definition of different physical activity time series with the appropriate analysis methods. The temporal pattern of postures, movements, and transitions between postures was quantified using fractal analysis and symbolic dynamics statistics. The derived nonlinear metrics were able to discriminate patterns of daily activity generated from healthy and chronic pain states.

Stimulus statistics shape oscillations in nonlinear recurrent neural networks.

Rhythmic activity plays a central role in neural computations and brain functions ranging from homeostasis to attention, as well as in neurological and neuropsychiatric disorders. Despite this pervasiveness, little is known about the mechanisms whereby the frequency and power of oscillatory activity are modulated, and how they reflect the inputs received by neurons. Numerous studies have reported input-dependent fluctuations in peak frequency and power (as well as couplings across these features). However, it remains unresolved what mediates these spectral shifts among neural populations. Extending previous findings regarding stochastic nonlinear systems and experimental observations, we provide analytical insights regarding oscillatory responses of neural populations to stimulation from either endogenous or exogenous origins. Using a deceptively simple yet sparse and randomly connected network of neurons, we show how spiking inputs can reliably modulate the peak frequency and power expressed by synchronous neural populations without any changes in circuitry. Our results reveal that a generic, non-nonlinear and input-induced mechanism can robustly mediate these spectral fluctuations, and thus provide a framework in which inputs to the neurons bidirectionally regulate both the frequency and power expressed by synchronous populations. Theoretical and computational analysis of the ensuing spectral fluctuations was found to reflect the underlying dynamics of the input stimuli driving the neurons. Our results provide insights regarding a generic mechanism supporting spectral transitions observed across cortical networks and spanning multiple frequency bands.

Changes in running mechanics and spring-mass behaviour during a 5-km time trial.

Research into the biomechanical manifestation of fatigue during exhaustive runs is increasingly popular but additional understanding of the adaptation of the spring-mass behaviour during the course of strenuous, self-paced exercises continues to be a challenge in order to develop optimized training and injury prevention programs. This study investigated continuous changes in running mechanics and spring-mass behaviour during a 5-km run. 12 competitive triathletes performed a 5-km running time trial (mean performance: ̴17 min 30 s) on a 200 m indoor track. Vertical and anterior-posterior ground reaction forces were measured every 200 m by a 5-m long force platform system, and used to determine spring-mass model characteristics. After a fast start, running velocity progressively decreased (- 11.6%; P<0.001) in the middle part of the race before an end spurt in the final 400-600 m. Stride length (- 7.4%; P<0.001) and frequency (- 4.1%; P=0.001) decreased over the 25 laps, while contact time (+ 8.9%; P<0.001) and total stride duration (+ 4.1%; P<0.001) progressively lengthened. Peak vertical forces (- 2.0%; P<0.01) and leg compression (- 4.3%; P<0.05), but not centre of mass vertical displacement (+ 3.2%; P>0.05), decreased with time. As a result, vertical stiffness decreased (- 6.0%; P<0.001) during the run, whereas leg stiffness changes were not significant (+ 1.3%; P>0.05). Spring-mass behaviour progressively changes during a 5-km time trial towards deteriorated vertical stiffness, which alters impact and force production characteristics.

Robust Nonlinear Mapping of Soil Contamination Using Support Regression

Spatial data analysis mapping and visualization is of great importance in various fields: environment, pollution, natural hazards and risks, epidemiology, spatial econometrics, etc. A basic task of spatial mapping is to make predictions based on some empirical data (measurements). A number of state-of-the-art methods can be used for the task: deterministic interpolations, methods of geostatistics: the family of kriging estimators (Deutsch and Journel, 1997), machine learning algorithms such as artificial neural networks (ANN) of different architectures, hybrid ANN-geostatistics models (Kanevski and Maignan, 2004; Kanevski et al., 1996), etc. All the methods mentioned above can be used for solving the problem of spatial data mapping. Environmental empirical data are always contaminated/corrupted by noise, and often with noise of unknown nature. That's one of the reasons why deterministic models can be inconsistent, since they treat the measurements as values of some unknown function that should be interpolated. Kriging estimators treat the measurements as the realization of some spatial randomn process. To obtain the estimation with kriging one has to model the spatial structure of the data: spatial correlation function or (semi-)variogram. This task can be complicated if there is not sufficient number of measurements and variogram is sensitive to outliers and extremes. ANN is a powerful tool, but it also suffers from the number of reasons. of a special type ? multiplayer perceptrons ? are often used as a detrending tool in hybrid (ANN+geostatistics) models (Kanevski and Maignank, 2004). Therefore, development and adaptation of the method that would be nonlinear and robust to noise in measurements, would deal with the small empirical datasets and which has solid mathematical background is of great importance. The present paper deals with such model, based on Statistical Learning Theory (SLT) - Support Vector Regression. SLT is a general mathematical framework devoted to the problem of estimation of the dependencies from empirical data (Hastie et al, 2004; Vapnik, 1998). SLT models for classification - Support Vector Machines - have shown good results on different machine learning tasks. The results of SVM classification of spatial data are also promising (Kanevski et al, 2002). The properties of SVM for regression - Support Vector Regression (SVR) are less studied. First results of the application of SVR for spatial mapping of physical quantities were obtained by the authorsin for mapping of medium porosity (Kanevski et al, 1999), and for mapping of radioactively contaminated territories (Kanevski and Canu, 2000). The present paper is devoted to further understanding of the properties of SVR model for spatial data analysis and mapping. Detailed description of the SVR theory can be found in (Cristianini and Shawe-Taylor, 2000; Smola, 1996) and basic equations for the nonlinear modeling are given in section 2. Section 3 discusses the application of SVR for spatial data mapping on the real case study - soil pollution by Cs137 radionuclide. Section 4 discusses the properties of the modelapplied to noised data or data with outliers.

Influence of the world's most challenging mountain ultra-marathon on energy cost and running mechanics

PURPOSE: To examine the effects of the world's most challenging mountain ultra-marathon (Tor des Géants(®) 2012) on the energy cost of three types of locomotion (cycling, level and uphill running) and running kinematics. METHODS: Before (pre-) and immediately after (post-) the competition, a group of ten male experienced ultra-marathon runners performed in random order three submaximal 4-min exercise trials: cycling at a power of 1.5 W kg(-1) body mass; level running at 9 km h(-1) and uphill running at 6 km h(-1) at an inclination of +15 % on a motorized treadmill. Two video cameras recorded running mechanics at different sampling rates. RESULTS: Between pre- and post-, the uphill-running energy cost decreased by 13.8 % (P = 0.004); no change was noted in the energy cost of level running or cycling (NS). There was an increase in contact time (+10.3 %, P = 0.019) and duty factor (+8.1 %, P = 0.001) and a decrease in swing time (-6.4 %, P = 0.008) in the uphill-running condition. CONCLUSION: After this extreme mountain ultra-marathon, the subjects modified only their uphill-running patterns for a more economical step mechanics.

Outcome of unilateral ankle arthrodesis and total ankle replacement in terms of bilateral gait mechanics.

Previous studies assessed the outcome of ankle arthrodesis (AA) and total ankle replacement (TAR) surgeries; however, the extent of postoperative recovery towards bilateral gait mechanics (BGM) is unknown. We evaluated the outcome of the two surgeries at least 2 years post rehabilitation, focusing on BGM. 36 participants, including 12 AA patients, 12 TAR patients, and 12 controls were included. Gait assessment over 50 m distance was performed utilizing pressure insoles and 3D inertial sensors, following which an intraindividual comparison was performed. Most spatiotemporal and kinematic parameters in the TAR group were indicative of good gait symmetry, while the AA group presented significant differences. Plantar pressure symmetry among the AA group was also significantly distorted. Abnormality in biomechanical behavior of the AA unoperated, contralateral foot was observed. In summary, our results indicate an altered BGM in AA patients, whereas a relatively fully recovered BGM is observed in TAR patients, despite the quantitative differences in several parameters when compared to a healthy population. Our study supports a biomechanical assessment and rehabilitation of both operated and unoperated sides after major surgeries for ankle osteoarthrosis.

Modelling of iron-filled magneto-active polymers with a dispersed chain-like microstructure

Magneto-active polymers are a class of smart materials commonly manufactured by mixing micron-sized iron particles in a rubber-like matrix. When cured in the presence of an externally applied magnetic field, the iron particles arrange themselves into chain-like structures that lend an overall anisotropy to the material. It has been observed through electron micrographs and X-ray tomographs that these chains are not always perfect in structure, and may have dispersion due to the conditions present during manufacturing or some undesirable material properties. We model the response of these materials to coupled magneto-mechanical loading in this paper using a probability based structure tensor that accounts for this imperfect anisotropy. The response of the matrix material is decoupled from the chain phase, though still being connected through kinematic constraints. The latter is based on the definition of a 'chain deformation gradient' and a 'chain magnetic field'. We conclude with numerical examples that demonstrate the effect of chain dispersion on the response of the material to magnetoelastic loading.

Energetics and mechanics of walking in patients with chronic low back pain and healthy matched controls.

PURPOSE: Walking in patients with chronic low back pain (cLBP) is characterized by motor control adaptations as a protective strategy against further injury or pain. The purpose of this study was to compare the preferred walking speed, the biomechanical and the energetic parameters of walking at different speeds between patients with cLBP and healthy men individually matched for age, body mass and height. METHODS: Energy cost of walking was assessed with a breath-by-breath gas analyser; mechanical and spatiotemporal parameters of walking were computed using two inertial sensors equipped with a triaxial accelerometer and gyroscope and compared in 13 men with cLBP and 13 control men (CTR) during treadmill walking at standard (0.83, 1.11, 1.38, 1.67 m s(-1)) and preferred (PWS) speeds. Low back pain intensity (visual analogue scale, cLBP only) and perceived exertion (Borg scale) were assessed at each walking speed. RESULTS: PWS was slower in cLBP [1.17 (SD = 0.13) m s(-1)] than in CTR group [1.33 (SD = 0.11) m s(-1); P = 0.002]. No significant difference was observed between groups in mechanical work (P ≥ 0.44), spatiotemporal parameters (P ≥ 0.16) and energy cost of walking (P ≥ 0.36). At the end of the treadmill protocol, perceived exertion was significantly higher in cLBP [11.7 (SD = 2.4)] than in CTR group [9.9 (SD = 1.1); P = 0.01]. Pain intensity did not significantly increase over time (P = 0.21). CONCLUSIONS: These results do not support the hypothesis of a less efficient walking pattern in patients with cLBP and imply that high walking speeds are well tolerated by patients with moderately disabling cLBP.

New insight in the relationship between regional patterns of knee cartilage thickness, osteoarthritis disease severity, and gait mechanics.

To test if the relationship between knee kinetics during walking and regional patterns of cartilage thickness is influenced by disease severity we tested the following hypotheses in a cross-sectional study of medial compartment osteoarthritis (OA) subjects: (1) the peak knee flexion (KFM) and adduction moments (KAM) during walking are associated with regional cartilage thickness and medial-to-lateral cartilage thickness ratios, and (2) the associations between knee moments and cartilage thickness data are dependent on disease severity. Seventy individuals with medial compartment knee OA were studied. Gait analysis was used to determine the knee moments and cartilage thickness was measured from magnetic resonance imaging. Multiple linear regression analyses tested for associations between cartilage thickness and knee kinetics. Medial cartilage thickness and medial-to-lateral cartilage thickness ratios were lower in subjects with greater KAM for specific regions of the femoral condyle and tibial plateau with no associations for KFM in patients of all disease severities. When separated by severity, the association between KAM and cartilage thickness was found only in patients with more severe OA, and KFM was significantly associated with cartilage thickness only for the less severe OA subjects for specific tibial plateau regions. The results support the idea that the KAM is larger in patients with more severe disease and the KFM has greater influence early in the disease process, which may lessen as pain increases with disease severity. Each component influences different regions of cartilage. Thus the relative contributions of both KAM and KFM should be considered when evaluating gait mechanics and the influence of any intervention for knee OA.

Comparison of Four Sections for Analyzing Running Mechanics Alterations During Repeated Treadmill Sprints.

We compared different approaches to analyze running mechanics alterations during repeated treadmill sprints. Thirteen active male athletes performed five 5-second sprints with 25 seconds of recovery on an instrumented treadmill. This approach allowed continuous measurement of running kinetics/kinematics and calculation of vertical and leg stiffness variables that were subsequently averaged over 3 distinct sections of the 5-second sprint (steps 2-5, 7-10, and 12-15) and for all steps (steps 2-15). Independently from the analyzed section, propulsive power and step frequency decreased with fatigue, while contact time and step length increased (P < .05). Except for step frequency, all mechanical variables varied (P < .05) across sprint sections. The only parameters that highly depend on running velocity (propulsive power and vertical stiffness) showed a significant interaction (P < .05) between the analyzed sections, with smaller magnitude of fatigue-induced change observed for steps 2-5. Considering all steps or only a few steps during early, middle, or late phases of 5-second sprints provides similar mechanical outcomes during repeated treadmill sprinting, although acceleration induces noticeable differences between the sections studied. Furthermore, quantifying mechanical alterations from the early acceleration phase may not be readily detectable, and is not recommended.