Conformational dynamics of HIV-1 protease: a comparative molecular dynamics simulation study with multiple amber force fields

Flap dynamics of HIV-1 protease (HIV-pr) controls the entry of inhibitors and substrates to the active site. Dynamical models from previous simulations are not all consistent with each other and not all are supported by the NMR results. In the present work, the er effect of force field on the dynamics of HIV-pr is investigated by MD simulations using three AMBER force fields ff99, ff99SB, and ff03. The generalized order parameters for amide backbone are calculated from the three force fields and compared with the NMR S2 values. We found that the ff99SB and ff03 force field calculated order parameters agree reasonably well with the NMR S2 values, whereas ff99 calculated values deviate most from the NMR order parameters. Stereochemical geometry of protein models from each force field also agrees well with the remarks from NMR S2 values. However, between ff99SB and ff03, there are several differences, most notably in the loop regions. It is found that these loops are, in general, more flexible in the ff03 force field. This results in a larger active site cavity in the simulation with the ff03 force field. The effect of this difference in computer-aided drug design against flexible receptors is discussed.

Force Biased Molecular Dynamics Simulation Study of Effect of Dendrimer Generation on Interaction with DNA

We have studied the effect of dendrimer generation on the interaction between dsDNA and the PAMAM dendrimer using force biased simulation of dsDNA with three generations of dendrimer: G3, G4, and G5. Our results for the potential of mean force (PMF) and the dendrimer asphericity along the binding pathway, combined with visualization of the simulations, demonstrate that dendrimer generation has a pronounced impact on the interaction. The PMF increases linearly with increasing generation of the dendrimer. While, in agreement with previous results, we see an increase in the extent to which the dendrimer bends the dsDNA with increasing dendrimer generation, we also see that the deformation of the dendrimer is greater with smaller generation of the dendrimer. The larger dendrimer forces the dsDNA to conform to its structure, while the smaller dendrimer is forced to conform to the structure of the dsDNA. Monitoring the number of bound cations at different values of force bias distance shows the expected effect of ions being expelled when the dendrimer binds dsDNA.

Dynamics of intermittent force fluctuations in vesicular nanotubulation

Irregular force fluctuations are seen in most nanotubulation experiments. The dynamics behind their presence has, however, been neither commented upon nor modeled. A simple estimate of the mean energy dissipated in force drops turns out to be several times the thermal energy. This coupled with the rate dependent nature of the deformation reported in several experiments point to a dynamical origin of the serrations. We simplify the whole process of tether formation through a three-stage model of successive deformations of sphere to ellipsoid, neck-formation, and tubule birth and extension. Based on this, we envisage a rate-softening frictional force at the neck that must be overcome before a nanotube can be pulled out. Our minimal model includes elastic and visco-elastic deformation of the vesicle, and has built-in dependence on pull velocity, vesicle radius, and other material parameters, enabling us to capture various kinds of serrated force-extension curves for different parameter choices. Serrations are predicted in the nanotubulation region. Other features of force-extension plots reported in the literature such as a plateauing serrated region beyond a force drop, serrated flow region with a small positive slope, an increase in the elastic threshold with pull velocity, force-extension curves for vesicles with larger radius lying lower than those for smaller radius, are all also predicted by the model. A toy model is introduced to demonstrate that the role of the friction law is limited to inducing stick-slip oscillations in the force, and all other qualitative and quantitative features emerging from the model can only be attributed to other physical mechanisms included in the deformation dynamics of the vesicle. (C) 2014 AIP Publishing LLC.

Mechanical force antagonizes the inhibitory effects of RecX on RecA filament formation in Mycobacterium tuberculosis

Efficient bacterial recombinational DNA repair involves rapid cycles of RecA filament assembly and disassembly. The RecX protein plays a crucial inhibitory role in RecA filament formation and stability. As the broken ends of DNA are tethered during homologous search, RecA filaments assembled at the ends are likely subject to force. In this work, we investigated the interplay between RecX and force on RecA filament formation and stability. Using magnetic tweezers, at single molecular level, we found that Mycobacterium tuberculosis (Mt) RecX could catalyze stepwise de-polymerization of preformed MtRecA filament in the presence of ATP hydrolysis at low forces (<7 pN). However, applying larger forces antagonized the inhibitory effects of MtRecX, and a partially de-polymerized MtRecA filament could repolymerize in the presence of MtRecX, which cannot be explained by previous models. Theoretical analysis of force-dependent conformational free energies of naked ssDNA and RecA nucleoprotein filament suggests that mechanical force stabilizes RecA filament, which provides a possible mechanism for the observation. As the antagonizing effect of force on the inhibitory function of RecX takes place in a physiological range; these findings broadly suggest a potential mechanosensitive regulation during homologous recombination.

The influence of fluctuation fields on the force-free field .2. The local expansion

In Paper I (Hu, 1982), we discussed the the influence of fluctuation fields on the force-free field for the case of conventional turbulence and demonstrated the general relationships. In the present paper, by using the approach of local expansion, the equation of average force-free field is obtained as (1+b)×B 0=(#x002B;a)B 0#x002B;a (1)×B 0#x002B;K. The average coefficientsa,a (1),b, andK show the influence of the fluctuation fields in small scale on the configurations of magnetic field in large scale. As the average magnetic field is no longer parallel to the average electric current, the average configurations of force-free fields are more general and complex than the usual ones. From the view point of physics, the energy and momentum of the turbulent structures should have influence on the equilibrium of the average fields. Several examples are discussed, and they show the basic features of the fluctuation fields and the influence of fluctuation fields on the average configurations of magnetic fields. The astrophysical environments are often in the turbulent state, the results of the present paper may be applied to the turbulent plasma where the magnetic field is strong.

Stabilizing to disruptive transition of focal adhesion response to mechanical forces

Strong mechanical forces can, obviously, disrupt cell-cell and cell-matrix adhesions, e.g., cyclic uniaxial stretch induces instability of cell adhesion, which then causes the reorientation of cells away from the stretching direction. However, recent experiments also demonstrated the existence of force dependent adhesion growth (rather than dissociation). To provide a quantitative explanation for the two seemingly contradictory phenomena, a microscopic model that includes both integrin-integrin interaction and integrin-ligand interaction is developed at molecular level by treating the focal adhesion as an adhesion cluster. The integrin clustering dynamics and integrin-ligand binding dynamics are then simulated within one unified theoretical frame with Monte Carlo simulation. We find that the focal adhesion will grow when the traction force is higher than a relative small threshold value, and the growth is dominated by the reduction of local chemical potential energy by the traction force. In contrast, the focal adhesion will rupture when the traction force exceeds a second threshold value, and the rupture is dominated by the breaking of integrin-ligand bonds. Consistent with the experiments, these results suggest a force map for various responses of cell adhesion to different scales of mechanical force. PMID: 20542514

The Synthesis of Stable Force-Closure Grasps

This thesis addresses the problem of synthesizing grasps that are force-closure and stable. The synthesis of force-closure grasps constructs independent regions of contact for the fingertips, such that the motion of the grasped object is totally constrained. The synthesis of stable grasps constructs virtual springs at the contacts, such that the grasped object is stable, and has a desired stiffness matrix about its stable equilibrium. A grasp on an object is force-closure if and only if we can exert, through the set of contacts, arbitrary forces and moments on the object. So force-closure implies equilibrium exists because zero forces and moment is spanned. In the reverse direction, we prove that a non-marginal equilibrium grasp is also a force-closure grasp, if it has at least two point contacts with friction in 2D, or two soft-finger contacts or three hard-finger contacts in 3D. Next, we prove that all force-closure grasps can be made stable, by using either active or passive springs at the contacts. The thesis develops a simple relation between the stability and stiffness of the grasp and the spatial configuration of the virtual springs at the contacts. The stiffness of the grasp depends also on whether the points of contact stick, or slide without friction on straight or curved surfaces of the object. The thesis presents fast and simple algorithms for directly constructing stable fore-closure grasps based on the shape of the grasped object. The formal framework of force-closure and stable grasps provides a partial explanation to why we stably grasp objects to easily, and to why our fingers are better soft than hard.

Evaporation-induced evolution of the capillary force between two grains

© 2014, Springer-Verlag Berlin Heidelberg.The evolution of capillary forces during evaporation and the corresponding changes in the geometrical characteristics of liquid (water) bridges between two glass spheres with constant separation are examined experimentally. For comparison, the liquid bridges were also tested for mechanical extension (at constant volume). The obtained results reveal substantial differences between the evolution of capillary force due to evaporation and the evolution due to extension of the liquid bridges. During both evaporation and extension, the change of interparticle capillary forces consists in a force decrease to zero either gradually or via rupture of the bridge. At small separations between the grains (short & wide bridges) during evaporation and at large volumes during extension, there is a slight initial increase of force. During evaporation, the capillary force decreases slowly at the beginning of the process and quickly at the end of the process; during extension, the capillary force decreases quickly at the beginning and slowly at the end of the process. Rupture during evaporation of the bridges occurs most abruptly for bridges with wider separations (tall and thin), sometimes occurring after only 25% of the water volume was evaporated. The evolution (pinning/depinning) of two geometrical characteristics of the bridge, the diameter of the three-phase contact line and the “apparent” contact angle at the solid/liquid/gas interface, seem to control the capillary force evolution. The findings are of relevance to the mechanics of unsaturated granular media in the final phase of drying.

Transfer of resistance training to enhance rapid coordinated force production by older adults

The purpose of this study was to examine the capacity of resistance training to enhance the rapid and coordinated production of force by older people. Thirty adults (greater than or equal to 60 years) completed a visually guided aiming task that required the generation of isometric torque in 2 df about the elbow prior to and following a 4-week training period. Groups of six participants were allocated to two progressive ( 40 - 100% maximal voluntary contraction (MVC)) resistance-training (PRT) groups, to two constant low-load (10% MVC) training groups (CLO) and to one no-training control group. Training movements required the generation of either combined flexion and supination (FLESUP), or combined extension and supination (EXTSUP). In response to training, target acquisition times in the aiming task decreased for all groups; however, both the nature of the training load and the training movement influenced the pattern and magnitude of improvements (EXTSUP_ CLO: 36%, FLESUP_ PRT 26%, EXTSUP_ PRT 22%, FLESUP_ CLO 20%, CONTROL 15%). For one group that trained with progressively increasing loads, there arose a subsequent decrease in performance in one condition of the transfer task. For each group, these adaptations were accompanied by systematic changes in the coordination of muscles about the elbow joint, particularly the biceps brachii.

An investigation into the acceleration response of a damaged beam-type structure to a moving force

In recent years there have been a growing number of publications on procedures for damage detection in beams from analysing their dynamic response to the passage of a moving force. Most of this research demonstrates their effectiveness by showing that a singularity that did not appear in the healthy structure is present in the response of the damaged structure. This paper elucidates from first principles how the acceleration response can be assumed to consist of ‘static’ and ‘dynamic’ components, and where the beam has experienced a localised loss in stiffness, an additional ‘damage’ component. The combination of these components establishes how the damage singularity will appear in the total response. For a given damage severity, the amplitude of the ‘damage’ component will depend on how close the damage location is to the sensor, and its frequency content will increase with higher velocities of the moving force. The latter has implications for damage detection because if the frequency content of the ‘damage’ component includes bridge and/or vehicle frequencies, it becomes more difficult to identify damage. The paper illustrates how a thorough understanding of the relationship between the ‘static‘ and ‘damage’ components contributes to establish if damage has occurred and to provide an estimation of its location and severity. The findings are corroborated using accelerations from a planar finite element simulation model where the effects of force velocity and bridge span are examined.

Force potentiation as a modulator of contractile performance: Implications for control of skeletal muscle force and energetics of work

This thesis investigated the modulation of dynamic contractile function and energetics of work by posttetanic potentiation (PTP). Mechanical experiments were conducted in vitro using software-controlled protocols to stimulate/determine contractile function during ramp shortening, and muscles were frozen during parallel incubations for biochemical analysis. The central feature of this research was the comparison of fast hindlimb muscles from wildtype and skeletal myosin light chain kinase knockout (skMLCK-/-) mice that does not express the primary mechanism for PTP: myosin regulatory light chain (RLC) phosphorylation. In contrast to smooth/cardiac muscles where RLC phosphorylation is indispensable, its precise physiological role in skeletal muscle is unclear. It was initially determined that tetanic potentiation was shortening speed dependent, and this sensitivity of the PTP mechanism to muscle shortening extended the stimulation frequency domain over which PTP was manifest. Thus, the physiological utility of RLC phosphorylation to augment contractile function in vivo may be more extensive than previously considered. Subsequent experiments studied the contraction-type dependence for PTP and demonstrated that the enhancement of contractile function was dependent on force level. Surprisingly, in the absence of RLC phosphorylation, skMLCK-/- muscles exhibited significant concentric PTP; consequently, up to ~50% of the dynamic PTP response in wildtype muscle may be attributed to an alternate mechanism. When the interaction of PTP and the catchlike property (CLP) was examined, we determined that unlike the acute augmentation of peak force by the CLP, RLC phosphorylation produced a longer-lasting enhancement of force and work in the potentiated state. Nevertheless, despite the apparent interference between these mechanisms, both offer physiological utility and may be complementary in achieving optimal contractile function in vivo. Finally, when the energetic implications of PTP were explored, we determined that during a brief period of repetitive concentric activation, total work performed was ~60% greater in wildtype vs. skMLCK-/- muscles but there was no genotype difference in High-Energy Phosphate Consumption or Economy (i.e. HEPC: work). In summary, this thesis provides novel insight into the modulatory effects of PTP and RLC phosphorylation, and through the observation of alternative mechanisms for PTP we further develop our understanding of the history-dependence of fast skeletal muscle function.

Trajectory and Force Control of a Direct Drive Arm

Using the MIT Serial Link Direct Drive Arm as the main experimental device, various issues in trajectory and force control of manipulators were studied in this thesis. Since accurate modeling is important for any controller, issues of estimating the dynamic model of a manipulator and its load were addressed first. Practical and effective algorithms were developed fro the Newton-Euler equations to estimate the inertial parameters of manipulator rigid-body loads and links. Load estimation was implemented both on PUMA 600 robot and on the MIT Serial Link Direct Drive Arm. With the link estimation algorithm, the inertial parameters of the direct drive arm were obtained. For both load and link estimation results, the estimated parameters are good models of the actual system for control purposes since torques and forces can be predicted accurately from these estimated parameters. The estimated model of the direct drive arm was them used to evaluate trajectory following performance by feedforward and computed torque control algorithms. The experimental evaluations showed that the dynamic compensation can greatly improve trajectory following accuracy. Various stability issues of force control were studied next. It was determined that there are two types of instability in force control. Dynamic instability, present in all of the previous force control algorithms discussed in this thesis, is caused by the interaction of a manipulator with a stiff environment. Kinematics instability is present only in the hybrid control algorithm of Raibert and Craig, and is caused by the interaction of the inertia matrix with the Jacobian inverse coordinate transformation in the feedback path. Several methods were suggested and demonstrated experimentally to solve these stability problems. The result of the stability analyses were then incorporated in implementing a stable force/position controller on the direct drive arm by the modified resolved acceleration method using both joint torque and wrist force sensor feedbacks.

Force constants and dipole‐moment derivatives of molecules from perturbed Hartree–Fock calculations

General expressions for the force constants and dipole‐moment derivatives of molecules are derived, and the problems arising in their practical application are reviewed. Great emphasis is placed on the use of the Hartree–Fock function as an approximate wavefunction, and a number of its properties are discussed and re‐emphasised. The main content of this paper is the development of a perturbed Hartree–Fock theory that makes possible the direct calculation of force constants and dipole‐moment derivatives from SCF–MO wavefunctions. Essentially the theory yields ∂ϕi / ∂RJα, the derivative of an MO with respect to a nuclear coordinate.

Highly efficient separation of actinides from lanthanides by a phenanthroline-derived bis-triazine ligand

The synthesis, lanthanide complexation, and solvent ex- traction of actinide(III) and lanthanide(III) radiotracers from nitric acid solutions by a phenanthroline-derived quadridentate bis-triazine ligand are described. The ligand separates Am(III) and Cm(III) from the lanthanides with remarkably high efficiency, high selectivity, and fast extraction kinetics compared to its 2,2'-bipyridine counterpart. Structures of the 1:2 bis-complexes of the ligand with Eu(III) and Yb(III) were elucidated by X-ray crystallography and force field calculations, respec-tively. The Eu(III) bis-complex is the first 1:2 bis-complex of a quadridentate bis-triazine ligand to be characterized by crystallography. The faster rates of extraction were verified by kinetics measurements using the rotating membrane cell technique in several diluents. The improved kinetics of metal ion extraction are related to the higher surface activity of the ligand at the phase interface. The improvement in the ligand's properties on replacing the bipyridine unit with a phenanthroline unit far exceeds what was anticipated based on ligand design alone.

First-principles based force-field for the interaction of proteins with Au(100)(5 x 1): an extension of GolP-CHARMM

Noncovalent recognition between peptides and inorganic materials is an established phenomenon. Key to exploiting these interactions in a wide range of materials self-assembly applications would be to harness the facet-selective control of peptide binding onto these materials. Fundamental understanding of what drives facet-selectivity in peptide binding is developing, but as yet is not sufficient to enable design of predictable facet-specific sequences. Computational simulation of the aqueous peptide-gold interface, commonly used to understand the mechanisms driving adsorption at an atomic level, has thus far neglected the role that surface reconstruction might play in facet specificity. Here the polarizable GolP-CHARMM suite of force fields is extended to include the reconstructed Au(100) surface. The force field, compatible with the bio-organic force field CHARMM, is parametrized using first-principles data. Our extended force field is tailored to reproduce the heterogeneity of weak chemisorbing N and S species to specific locations in the Au(100)(5 × 1) surface identified from the first-principles calculations. We apply our new model to predict and compare the three-dimensional structure of liquid water at Au(111), Au(100)(1 × 1), and Au(100)(5 × 1) interfaces. Using molecular dynamics simulations, we predict an increased likelihood for water-mediated peptide adsorption at the aqueous-Au(100)(1 × 1) interface compared with the Au(100)(5 × 1) interface. Therefore, our findings suggest that peptide binding can discriminate between the native and reconstructed Au(100) interfaces and that the role of reconstruction on binding at the Au(100) interface should not be neglected. © 2013 American Chemical Society.