Novel strategies in feedforward adaptation to a position-dependent perturbation

To investigate the control mechanisms used in adapting to position-dependent forces, subjects performed 150 horizontal reaching movements over 25 cm in the presence of a position-dependent parabolic force field (PF). The PF acted only over the first 10 cm of the movement. On every fifth trial, a virtual mechanical guide (double wall) constrained subjects to move along a straight-line path between the start and target positions. Its purpose was to register lateral force to track formation of an internal model of the force field, and to look for evidence of possible alternative adaptive strategies. The force field produced a force to the right, which initially caused subjects to deviate in that direction. They reacted by producing deviations to the left, into the force field, as early as the second trial. Further adaptation resulted in rapid exponential reduction of kinematic error in the latter portion of the movement, where the greatest perturbation to the handpath was initially observed, whereas there was little modification of the handpath in the region where the PF was active. Significant force directed to counteract the PF was measured on the first guided trial, and was modified during the first half of the learning set. The total force impulse in the region of the PF increased throughout the learning trials, but it always remained less than that produced by the PF. The force profile did not resemble a mirror image of the PF in that it tended to be more trapezoidal than parabolic in shape. As in previous studies of force-field adaptation, we found that changes in muscle activation involved a general increase in the activity of all muscles, which increased arm stiffness, and selectively-greater increases in the activation of muscles which counteracted the PF. With training, activation was exponentially reduced, albeit more slowly than kinematic error. Progressive changes in kinematics and EMG occurred predominantly in the region of the workspace beyond the force field. We suggest that constraints on muscle mechanics limit the ability of the central nervous system to employ an inverse dynamics model to nullify impulse-like forces by generating mirror-image forces. Consequently, subjects adopted a strategy of slightly overcompensating for the first half of the force field, then allowing the force field to push them in the opposite direction. Muscle activity patterns in the region beyond the boundary of the force field were subsequently adjusted because of the relatively-slow response of the second-order mechanics of muscle impedance to the force impulse.

Molecular dynamics simulations of the hydrophobin SC3 at a hydrophobic/hydrophilic interface

Hydrophobins are small (similar to 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/ hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily P-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in a-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.

Hydrotalcite intercalated siRNA:computational characterization of the interlayer environment

Using molecular dynamics (MD) simulations, we explore the structural and dynamical properties of siRNA within the intercalated environment of a Mg:Al 2:1 Layered Double Hydroxide (LDH) nanoparticle. An ab initio force field (Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies: COMPASS) is used for the MD simulations of the hybrid organic-inorganic systems. The structure, arrangement, mobility, close contacts and hydrogen bonds associated with the intercalated RNA are examined and contrasted with those of the isolated RNA. Computed powder X-ray diffraction patterns are also compared with related LDH-DNA experiments. As a method of probing whether the intercalated environment approximates the crystalline or rather the aqueous state, we explore the stability of the principle parameters (e.g., the major groove width) that differentiate both A- and A'- crystalline forms of siRNA and contrast this with recent findings for the same siRNA simulated in water. We find the crystalline forms remain structurally distinct when intercalated, whereas this is not the case in water. Implications for the stability of hybrid LDH-RNA systems are discussed.

Modelagem molecular voltada para a construção de um nanobiossensor para detecção do herbicida Imazaquin

In this study, our goal was develop and describe a molecular model of the enzyme-inhibiting interaction which can be used for an optimized projection of a Microscope Force Atomic nanobiosensor to detect pesticides molecules, used in agriculture, to evaluate its accordance with limit levels stipulated in valid legislation for its use. The studied herbicide (imazaquin) is a typical member of imidazolinone family and is an inhibitor of the enzymatic activity of Acetohydroxiacid Synthase (AHAS) enzyme that is responsible for the first step of pathway for the synthesis of side-chains in amino acids. The analysis of this enzyme property in the presence of its cofactors was made to obtain structural information and charge distribution of the molecular surface to evaluate its capacity of became immobilized on the Microscopy Atomic Force tip. The computational simulation of the system, using Molecular Dynamics, was possible with the force-field parameters for the cofactor and the herbicides obtained by the online tool SwissParam and it was implemented in force-field CHARMM27, used by software GROMACS; then appropriated simulations were made to validate the new parameters. The molecular orientation of the AHAS was defined based on electrostatic map and the availability of the herbicide in the active site. Steered Molecular Dynamics (SMD) Simulations, followed by quantum mechanics calculations for more representative frames, according to the sequential QM/MM methodology, in a specific direction of extraction of the herbicide from the active site. Therefore, external harmonic forces were applied with similar force constants of AFM cantilever for to simulate herbicide detection experiments by the proposed nanobiosensor. Force value of 1391 pN and binding energy of -14048.52 kJ mol-1 were calculated.

Computational and experimental study of the behavior of cyano-based ionic liquids in aqueous solution

The solvation of cyano- (CN-) based ionic liquids (ILs) and their capacity to establish hydrogen bonds (H-bonds) with water was studied by means of experimental and computational approaches. Experimentally, water activity data were measured for aqueous solutions of ILs based on 1-butyl-3-methylimidazolium ([BMIM](+)) cation combined with one of the following anions: thiocyanate ([SCN](-)), dicyanamide ([DCA](-)), or tricyanomethanide ([TCM](-)), and of 1-ethyl-3-methylimidazolium tetracyanoborate ([EMIM][TCB]). From the latter data, water activity coefficients were estimated showing that [BMIM][SCN] and [BMIM][DCA], unlike [BMIM][TCM] and [EMIM][TCB], are able to establish favorable interactions with water. Computationally, the conductor like screening model for real solvents (COSMO-RS) was used to estimate the water activity coefficients which compare well with the experimental ones. From the COSMO-RS results, it is suggested that the polarity of each ion composing the ILs has a strong effect on the solvation phenomena. Furthermore, classical molecular dynamics (MD) simulations were performed for obtaining an atomic level picture of the local molecular neighborhood of the different species. From the experimental and computational data it is showed that increasing the number of CN groups in the ILs' anions does not enhance their ability to establish H-bonds with water but decreases their polarities, being [BMIM][DCA] and [BMIM][SCN] the ones presenting higher propensity to interact.

A robust and reproducible procedure for cross-linking thermoset polymers using molecular simulation

Molecular simulation can provide valuable guidance in establishing clear links between structure and function to enable the design of new polymer-based materials. However, molecular simulation of thermoset polymers in particular, such as epoxies, present specific challenges, chiefly in the credible preparation of polymerised samples. Despite this need, a comprehensive, reproducible and robust process for accomplishing this using molecular simulation is still lacking. Here, we introduce a clear and reproducible cross-linking protocol to reliably generate three dimensional epoxy cross-linked polymer structures for use in molecular simulations. This protocol is sufficiently detailed to allow complete reproduction of our results, and is applicable to any general thermoset polymer. Amongst our developments, key features include a reproducible procedure for calculation of partial atomic charges, a reliable process for generating and validating an equilibrated liquid precursor mixture, and establishment of a novel, robust and reproducible protocol for generating the three-dimensional cross-linked solid polymer. We use these structures as input to subsequent molecular dynamics simulations to calculate a range thermo-mechanical properties, which compare favourably with experimental data. Our general protocol provides a benchmark for the process of simulating epoxy polymers, and can be readily translated to prepare and model epoxy samples that are dynamically cross-linked in the presence of surfaces and nanostructures.

A THEORETICAL STUDY OF INTERACTION OF NANOPARTICLES WITH BIOMOLECULE

Many types of materials at nanoscale are currently being used in everyday life. The production and use of such products based on engineered nanomaterials have raised concerns of the possible risks and hazards associated with these nanomaterials. In order to evaluate and gain a better understanding of their effects on living organisms, we have performed first-principles quantum mechanical calculations and molecular dynamics simulations. Specifically, we will investigate the interaction of nanomaterials including semiconducting quantum dots and metallic nanoparticles with various biological molecules, such as dopamine, DNA nucleobases and lipid membranes. Firstly, interactions of semiconducting CdSe/CdS quantum dots (QDs) with the dopamine and the DNA nucleobase molecules are investigated using similar quantum mechanical approach to the one used for the metallic nanoparticles. A variety of interaction sites are explored. Our results show that small-sized Cd4Se4 and Cd4S4 QDs interact strongly with the DNA nucleobase if a DNA nucleobase has the amide or hydroxyl chemical group. These results indicate that these QDs are suitable for detecting subcellular structures, as also reported by experiments. The next two chapters describe a preparation required for the simulation of nanoparticles interacting with membranes leading to accurate structure models for the membranes. We develop a method for the molecular crystalline structure prediction of 1,2-Dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC), 1,2-Dimyristoyl-sn-glycero-3-phosphorylethanolamine (DMPE) and cyclic di-amino acid peptide using first-principles methods. Since an accurate determination of the structure of an organic crystal is usually an extremely difficult task due to availability of the large number of its conformers, we propose a new computational scheme by applying knowledge of symmetry, structural chemistry and chemical bonding to reduce the sampling size of the conformation space. The interaction of metal nanoparticles with cell membranes is finally carried out by molecular dynamics simulations, and the results are reported in the last chapter. A new force field is developed which accurately describes the interaction forces between the clusters representing small-sized metal nanoparticles and the lipid bilayer molecules. The permeation of nanoparticles into the cell membrane is analyzed together with the RMSD values of the membrane modeled by a lipid bilayer. The simulation results suggest that the AgNPs could cause the same amount of deformation as the AuNPs for the dysfunction of the membrane.

PARAMETRIC STUDY OF REAXFF SIMULATION PARAMETERS FOR MOLECULAR DYNAMICS MODELING OF REACTIVE CARBON GASES

Abstract The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although the Reax Force Field (ReaxFF) can be used to simulate the chemical behavior of carbon-based systems, the simulation settings required for accurate predictions have not been fully explored. Using the ReaxFF, molecular dynamics (MD) simulations are used to simulate the chemical behavior of pure carbon and hydrocarbon reactive gases that are involved in the formation of carbon structures such as graphite, buckyballs, amorphous carbon, and carbon nanotubes. It is determined that the maximum simulation time step that can be used in MD simulations with the ReaxFF is dependent on the simulated temperature and selected parameter set, as are the predicted reaction rates. It is also determined that different carbon-based reactive gases react at different rates, and that the predicted equilibrium structures are generally the same for the different ReaxFF parameter sets, except in the case of the predicted formation of large graphitic structures with the Chenoweth parameter set under specific conditions.

Molecular mechanism of stabilization of thin films for improved water evaporation protection

All-atom molecular dynamics simulations and experimental characterization have been used to examine the structure and dynamics of novel evaporation-suppressing films where the addition of a water-soluble polymer to an ethylene glycol monooctadecyl ether monolayer leads to improved water evaporation resistance. Simulations and Langmuir trough experiments demonstrate the surface activity of poly(vinyl pyrrolidone) (PVP). Subsequent MD simulations performed on the thin films supported by the PVP sublayer show that, at low surface pressures, the polymer tends to concentrate at the film/water interface. The simulated atomic concentration profiles, hydrogen bonding patterns, and mobility analyses of the water-polymer-monolayer interfaces reveal that the presence of PVP increases the atomic density near the monolayer film, improves the film stability, and reduces the mobility of interfacial waters. These observations explain the molecular basis of the improved efficacy of these monolayer/polymer systems for evaporation protection of water and can be used to guide future development of organic thin films for other applications.

Effet de la douleur expérimentale tonique sur l’acquisition et la rétention d’une adaptation locomotrice

Introduction : Une proportion importante des individus ayant recours à des services de réadaptation physique vit avec de la douleur et des incapacités locomotrices. Plusieurs interventions proposées par les professionnels de la réadaptation afin de cibler leurs difficultés locomotrices nécessitent des apprentissages moteurs. Toutefois, très peu d’études ont évalué l’influence de la douleur sur l’apprentissage moteur et aucune n’a ciblé l’apprentissage d’une nouvelle tâche locomotrice. L’objectif de la thèse était d’évaluer l’influence de stimulations nociceptives cutanée et musculaire sur l’acquisition et la rétention d’une adaptation locomotrice. Méthodologie : Des individus en santé ont participé à des séances de laboratoire lors de deux journées consécutives. Lors de chaque séance, les participants devaient apprendre à marcher le plus normalement possible en présence d’un champ de force perturbant les mouvements de leur cheville, produit par une orthèse robotisée. La première journée permettait d’évaluer le comportement des participants lors de la phase d’acquisition de l’apprentissage. La seconde journée permettait d’évaluer leur rétention. Selon le groupe expérimental, l’apprentissage se faisait en présence d’une stimulation nociceptive cutanée, musculaire ou d’aucune stimulation (groupe contrôle). Initialement, l’application du champ de force provoquait d’importantes déviations des mouvements de la cheville (i.e. erreurs de mouvement), que les participants apprenaient graduellement à réduire en compensant activement la perturbation. L’erreur de mouvement moyenne durant la phase d’oscillation (en valeur absolue) a été quantifiée comme indicateur de performance. Une analyse plus approfondie des erreurs de mouvement et de l’activité musculaire a permis d’évaluer les stratégies motrices employées par les participants. Résultats : Les stimulations nociceptives n’ont pas affecté la performance lors de la phase d’acquisition de l’apprentissage moteur. Cependant, en présence de douleur, les erreurs de mouvement résiduelles se trouvaient plus tard dans la phase d’oscillation, suggérant l’utilisation d’une stratégie motrice moins anticipatoire que pour le groupe contrôle. Pour le groupe douleur musculaire, cette stratégie était associée à une activation précoce du muscle tibial antérieur réduite. La présence de douleur cutanée au Jour 1 interférait avec la performance des participants au Jour 2, lorsque le test de rétention était effectué en absence de douleur. Cet effet n’était pas observé lorsque la stimulation nociceptive cutanée était appliquée les deux jours, ou lorsque la douleur au Jour 1 était d’origine musculaire. Conclusion : Les résultats de cette thèse démontrent que dans certaines circonstances la douleur peut influencer de façon importante la performance lors d’un test de rétention d’une adaptation locomotrice, malgré une performance normale lors de la phase d’acquisition. Cet effet, observé uniquement avec la douleur cutanée, semble cependant plus lié au changement de contexte entre l’acquisition des habiletés motrices et le test de rétention (avec vs. sans douleur) qu’à une interférence directe avec la consolidation des habiletés motrices. Par ailleurs, malgré l’absence d’influence de la douleur sur la performance des participants lors de la phase d’acquisition de l’apprentissage, les stratégies motrices utilisées par ceux-ci étaient différentes de celles employées par le groupe contrôle.

Symmetry of in-shoe force signatures in athletes under field conditions : an application for injury prevention

The aetiology behind overuse injuries such as stress fractures is complex and multi-factorial. In sporting events where the loading is likely to be uneven (e.g. hurdling and jumps), research has suggested that the frequency of stress fractures seems to favour the athlete’s dominant limb. The tendency for an individual to have a preferred limb for voluntary motor acts makes limb selection a possible factor behind the development of unilateral overuse injuries, particularly when repeatedly used during high loading activities. The event of sprint hurdling is well suited for the study of loading asymmetry as the hurdling technique is repetitive and the limb movement asymmetrical. Of relevance to this study is the high incidence of Navicular Stress Fractures (NSF) in hurdlers, with suggestions there is a tendency for the fracture to develop in the trail leg foot, although this is not fully accepted. The Ground Reaction Force (GRF) with each foot contact is influenced by the hurdle action, with research finding step-to-step loading variations. However, it is unknown if this loading asymmetry extends to individual forefoot joints, thereby influencing stress fracture development. The first part of the study involved a series of investigations using a commercially available matrix style in-shoe sensor system (FscanTM, Tekscan Inc.). The suitability of insole sensor systems and custom made discrete sensors for use in hurdling-related training activities was assessed. The methodology used to analyse foot loading with each technology was investigated. The insole and discrete sensors systems tested proved to be unsuitable for use during full pace hurdling. Instead, a running barrier task designed to replicate the four repetitive foot contacts present during hurdling was assessed. This involved the clearance of a series of 6 barriers (low training hurdles), place in a straight line, using 4 strides between each. The second part of the study involved the analysis of "inter-limb" and "within foot loading asymmetries" using stance duration as well as vertical GRF under the Hallux (T1), the first metatarsal head (M1) and the central forefoot peak pressure site (M2), during walking, running, and running with barrier clearances. The contribution to loading asymmetry that each of the four repetitive foot contacts made during a series of barrier clearances was also assessed. Inter-limb asymmetry, in forefoot loading, occurred at discrete forefoot sites in a non-uniform manner across the three gait conditions. When the individual barrier foot contacts were compared, the stance duration was asymmetrical and the proportion of total forefoot load at M2 was asymmetrical. There were no significant differences between the proportion of forefoot load at M1, compared to M2; for any of the steps involved in the barrier clearance. A case study testing experimental (discrete) sensors during full pace sprinting and hurdling found that during both gait conditions, the trail limb experienced the greater vertical GRF at M1 and M2. During full pace hurdling, increased stance duration and vertical loading was a characteristic of the trail limb hurdle foot contacts. Commercially available in-shoe systems are not suitable for on field assessment of full pace hurdling. For the use of discrete sensor technology to become commonplace in the field, more robust sensors need to be developed.

Field measurement of wheel-rail impact force at insulated rail joint

When wheels pass over insulated rail joints (IRJs) a vertical impact force is generated. The ability to measure the impact force is valuable as the force signature helps understand the behaviour of the IRJs, in particular their potential for failure. The impact forces are thought to be one of the main factors that cause damage to the IRJ and track components. Study of the deterioration mechanism helps finding new methods to improve the service life of IRJs in track. In this research, the strain-gage-based wheel load detector, for the first time, is employed to measure the wheel–rail contact-impact force at an IRJ in a heavy haul rail line. In this technique, the strain gages are installed within the IRJ assembly without disturbing the structural integrity of IRJ and arranged in a full wheatstone bridge to form a wheel load detector. The instrumented IRJ is first tested and calibrated in the lab and then installed in the field. For comparison purposes, a reference rail section is also instrumented with the same strain gage pattern as the IRJ. In this paper the measurement technique, the process of instrumentation, and tests as well as some typical data obtained from the field and the inferences are presented.

Studies on Field Dependent Domain Structures in Multi-Grained 0.85PbMg(1/3)Nb(2/3)O3-0.15PbTiO(3) Thin Films by Scanning Force Microscopy

0.85PbMg(1/3)Nb(2/3)O(3)-0.15PbTiO(3) (0.85PMN-0.15PT) ferroelectric relaxor thin films have been deposited on La0.5Sr0.5CoO3/(111) Pt/TiO2/SiO2/Si by pulsed laser ablation by varying the oxygen partial pressures from 50 mTorr to 400 mTorr. The X-ray diffraction pattern reveals a pyrochlore free polycrystalline film. The grain morphology of the deposited films was studied using scanning electron microscopy and was found to be affected by oxygen pressure. By employing dynamic contact-electrostatic force microscopy we found that the distribution of polar nanoregions is majorly affected by oxygen pressure. Finally, the electric field induced switching in these films is discussed in terms of domain wall pinning.