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.

Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics

The way in which basal tractions, associated with mantle convection, couples with the lithosphere is a fundamental problem in geodynamics. A successful lithosphere-mantle coupling model for the Earth will satisfy observations of plate motions, intraplate stresses, and the plate boundary zone deformation. We solve the depth integrated three-dimensional force balance equations in a global finite element model that takes into account effects of both topography and shallow lithosphere structure as well as tractions originating from deeper mantle convection. The contribution from topography and lithosphere structure is estimated by calculating gravitational potential energy differences. The basal tractions are derived from a fully dynamic flow model with both radial and lateral viscosity variations. We simultaneously fit stresses and plate motions in order to delineate a best-fit lithosphere-mantle coupling model. We use both the World Stress Map and the Global Strain Rate Model to constrain the models. We find that a strongly coupled model with a stiff lithosphere and 3-4 orders of lateral viscosity variations in the lithosphere are best able to match the observational constraints. Our predicted deviatoric stresses, which are dominated by contribution from mantle tractions, range between 20-70 MPa. The best-fitting coupled models predict strain rates that are consistent with observations. That is, the intraplate areas are nearly rigid whereas plate boundaries and some other continental deformation zones display high strain rates. Comparison of mantle tractions and surface velocities indicate that in most areas tractions are driving, although in a few regions, including western North America, tractions are resistive. Citation: Ghosh, A., W. E. Holt, and L. M. Wen (2013), Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics.

A peridynamic theory for linear elastic shells

A state-based peridynamic formulation for linear elastic shells is presented. The emphasis is on introducing, possibly for the first time, a general surface based peridynamic model to represent the deformation characteristics of structures that have one geometric dimension much smaller than the other two. A new notion of curved bonds is exploited to cater for force transfer between the peridynamic particles describing the shell. Starting with the three dimensional force and deformation states, appropriate surface based force, moment and several deformation states are arrived at. Upon application on the curved bonds, such states yield the necessary force and deformation vectors governing the motion of the shell. By incorporating a shear correction factor, the formulation also accommodates analysis of shells that have higher thickness. In order to attain this, a consistent second order approximation to the complementary energy density is considered and incorporated in peridynamics via constitutive correspondence. Unlike the uncoupled constitution for thin shells, a consequence of a first order approximation, constitutive relations for thick shells are fully coupled in that surface wryness influences the in-plane stress resultants and surface strain the moments. Our proposal on the peridynamic shell theory is numerically assessed against simulations on static deformation of spherical and cylindrical shells, that of flat plates and quasi-static fracture propagation in a cylindrical shell. (C) 2016 Elsevier Ltd. All rights reserved.

Note: Design and development of an integrated three-dimensional scanner for atomic force microscopy

A compact scanning head for the Atomic Force Microscope (AFM) greatly enhances the portability of AFM and facilitates easy integration with other tools. This paper reports the design and development of a three-dimensional (3D) scanner integrated into an AFM micro-probe. The scanner is realized by means of a novel design for the AFM probe along with a magnetic actuation system. The integrated scanner, the actuation system, and their associated mechanical mounts are fabricated and evaluated. The experimentally calibrated actuation ranges are shown to be over 1 mu m along all the three axes. (c) 2013 AIP Publishing LLC.

Direct Measurement of Three-Dimensional Forces in Atomic Force Microscopy

Direct measurement of three-dimensional (3-D) forces between an atomic force microscope (AFM) probe and the sample benefits diverse applications of AFM, including force spectroscopy, nanometrology, and manipulation. This paper presents the design and evaluation of a measurement system, wherein the deflection of the AFM probe is obtained at two points to enable direct measurement of all the three components of 3-D tip-sample forces in real time. The optimal locations for measurement of deflection on the probe are derived for a conventional AFM probe. Further, a new optimal geometry is proposed for the probe that enables measurement of 3-D forces with identical sensitivity and nearly identical resolution along all three axes. Subsequently, the designed measurement system and the optimized AFM probe are both fabricated and evaluated. The evaluation demonstrates accurate measurement of tip-sample forces with minimal cross-sensitivities. Finally, the real-time measurement system is employed as part of a feedback control system to regulate the normal component of the interaction force, and to perform force-controlled scribing of a groove on the surface of polymethyl methacrylate.

Combined Atomic Force Microscopy and Modeling Study of The Evolution of Octadecylamine Films on a Mica Surface

The time evolution of the film thickness and domain formation of octadecylamine molecules adsorbed oil a mica surface is investigated Using atomic force microscopy. The adsorbed Film thickness is determined by measuring the height profile across the mica-amine interface of a mica surface partially immersed in a 15 mM solution of octadecylamine in chloroform. Using this novel procedure, adsorption of amine on mica is found to occur in three distinct stages, with morphologically distinct domain Formation and growth occurring during each stage. In the first stage, where adsorption is primarily in the thin-film regime, all average Film thickness of 0.2 (+/- 0.3) nm is formed for exposure times below 30 s and 0.8 (+/- 0.2) nm for 60 s of immersion time. During this stage, large sample spanning domains are observed. The second stage, which occurs between 60-300 s, is associated with it regime of rapid film growth, and the film thickness increases from about 0.8 to 25 nm during this stage. Once the thick-film regime is established, further exposure to the amine solution results in all increase in the domain area, and it regime of lateral domain growth is observed. In this stage, the domain area coverage grows from 38 to 75%, and the FTIR spectra reveal an increased level of crystallinity in the film. Using it diffusion-controlled model and it two-step Langmuir isotherm, the time evolution of the film growth is quantitatively captured. The model predicts the time at which the thin to thick film transition occurs as well its the time required for complete film growth at longer times. The Ward-Tordai equation is also solved to determine the model parameters in the monolayer (thin-film) regime, which occurs during the initial stages of film growth.

Polymer drift in a solvent by force acting on one polymer end

We investigate the effect of hydrodynamic interactions on the non-equilibrium drift dynamics of an ideal flexible polymer pulled by a constant force applied at one polymer end using the perturbation theory and the renormalization group method. For moderate force, if the polymer elongation is small, the hydrodynamic interactions are not screened and the velocity and the longitudinal elongation of the polymer are computed using the renormalization group method. Both the velocity and elongation are nonlinear functions of the driving force in this regime. For large elongation we found two regimes. For large force but finite chain length L the hydrodynamic interactions are screened. For large chain lengths and a finite force the hydrodynamic interactions are only partially screened, which in three dimensions results in unusual logarithmic corrections to the velocity and the longitudinal elongation.

Miniature Compliant Grippers With Vision-Based Force Sensing

This paper is concerned with grasping biological cells in aqueous medium with miniature grippers that can also help estimate forces using vision-based displacement measurement and computation. We present the design, fabrication, and testing of three single-piece, compliant miniature grippers with parallel and angular jaw motions. Two grippers were designed using experience and intuition, while the third one was designed using topology optimization with implicit manufacturing constraints. These grippers were fabricated using different manufacturing techniques using spring steel and polydimethylsiloxane ( PDMS). The grippers also serve the purpose of a force sensor. Toward this, we present a vision-based force-sensing technique by solving Cauchy's problem in elasticity using an improved algorithm. We validated this technique at the macroscale, where there was an independent method to estimate the force. In this study, the gripper was used to hold a yeast ball and a zebrafish egg cell of less than 1 mm in diameter. The forces involved were estimated to be about 30 and 10 mN for the yeast ball and the zebrafish egg cell, respectively.

Analysis of continuum using boundary compatibility conditions of integrated force method

Static and vibration problems of an indeterminate continuum are traditionally analyzed by the stiffness method. The force method is more or less non-existent for such problems. This situation is primarily due to the incomplete state of development of the compatibility conditions which are essential for the analysis of indeterminate structures by the flexibility method. The understanding of the Compatibility Conditions (CC) has been substantially augmented. Based on the understanding of CC, a novel formulation termed the Integrated Force Method (IFM) has been established. In this paper IFM has been extended for the static and vibration analyses of a continuum. The IFM analysis is illustrated taking three examples: 1. (1) rectangular plate in flexure 2. (2) analysis of a cantilevered dam 3. (3) free vibration analysis of a beam. From the examples solved it is observed that the force response of an indeterminate continuum with mixed boundary conditions can be generated by IFM without any reference to displacements in the field or on the boundary. Displacements if required can be calculated by back substitution.

Measurement of Gibbs energy of formation of Ca2PbO4 using a solid-state cell with three electrodes

Phase relations in the system Ca-Pb-O at 1100 K have been determined by equilibrating 18 compositions in the ternary and identifying the phases present in quenched samples by X-ray diffraction and energy dispersive X-ray analysis (EDX). Only one ternary compound Ca2PbO4 was found to be present. The compound coexists with CaO and PbO. The intermetallic compounds Ca2Pb, Ca5Pb3 and CaPb and liquid alloys are in equilibrium with CaO. The standard Gibbs energies of formation of Ca2PbO4 (880 - 1100 K) and Pb3O4 (770 - 910 K) were determined using solid-state cells based on yttria-stabilized zirconia as the solid electrolyte. Pure oxygen gas at 0.1 MPa was used as the reference electrode. For measurements on Ca2PbO4, a novel cell design with three electrodes in series, separated by solid electrolyte membranes, was used to avoid polarization of the electrode containing three solid phases. Two three-phase electrodes were used. The first absorbs the electrochemical flux of oxygen from the reference electrode to the measuring electrode. The other three-phase electrode, which is unaffected by the oxygen flux through the solid electrolyte, is used for electromotive force (EMF) measurement. The results from EMF studies were cross-checked using thermogravimetry (TG) under controlled oxygen partial pressures. The stability of Pb3O4 was investigated using a conventional solid-state cell with RuO2 electrodes. The results can be summarized by the following equations: 2CaO + PbO +1/2O(2) --> Ca2PbO4 Delta(r)G degrees/J mol(-1) = (- 128340 + 93.21 T/K) +/- 200 3PbO + 1/2O(2) --> Pb3O4 Delta(r)G degrees/J mol(-1) = (- 70060 + 77.5 T/K) +/- 150

A novel free floating accelerometer force balance system for shock tunnel applications

In order to overcome the interference of the model mounting system with the external aerodynamics of the body during shock tunnel testing, a new free floating internally mountable balance system that ensures unrestrained model motion during testing has been designed, fabricated and tested. Minimal friction ball bearings are used for ensuring the free floating condition of the model during tunnel testing. The drag force acting on a blunt leading edge flat plate at hypersonic Mach number has been measured using the new balance system. Finite element model (FEM) and CFD are exhaustively used in the design as well as for calibrating the new balance system. The experimentally measured drag force on the blunt leading edge flat plate at stagnation enthalpy of 0.7 and 1.2 MJ/kg and nominal Mach number of 5.75 matches well with FEM results. The concept can also be extended for measuring all the three fundamental aerodynamic forces in short duration test facilities like free piston driven shock tunnels.

Three-Component Self-Assembly of a Series of Triply Interlocked Pd12 Coordination Prisms and Their Non-Interlocked Pd6 Analogues

Template-assisted formation of multicomponent Pd6 coordination prisms and formation of their self-templated triply interlocked Pd12 analogues in the absence of an external template have been established in a single step through Pd?N/Pd?O coordination. Treatment of cis-[Pd(en)(NO3)2] with K3tma and linear pillar 4,4'-bpy (en=ethylenediamine, H3tma=benzene-1,3,5-tricarboxylic acid, 4,4'-bpy=4,4'-bipyridine) gave intercalated coordination cage [{Pd(en)}6(bpy)3(tma)2]2[NO3]12 (1) exclusively, whereas the same reaction in the presence of H3tma as an aromatic guest gave a H3tma-encapsulating non-interlocked discrete Pd6 molecular prism [{Pd(en)}6(bpy)3(tma)2(H3tma)2][NO3]6 (2). Though the same reaction using cis-[Pd(NO3)2(pn)] (pn=propane-1,2-diamine) instead of cis-[Pd(en)(NO3)2] gave triply interlocked coordination cage [{Pd(pn)}6(bpy)3(tma)2]2[NO3]12 (3) along with non-interlocked Pd6 analogue [{Pd(pn)}6(bpy)3(tma)2](NO3)6 (3'), and the presence of H3tma as a guest gave H3tma-encapsulating molecular prism [{Pd(pn)}6(bpy)3(tma)2(H3tma)2][NO3]6 (4) exclusively. In solution, the amount of 3' decreases as the temperature is decreased, and in the solid state 3 is the sole product. Notably, an analogous reaction using the relatively short pillar pz (pz=pyrazine) instead of 4,4'-bpy gave triply interlocked coordination cage [{Pd(pn)}6(pz)3(tma)2]2[NO3]12 (5) as the single product. Interestingly, the same reaction using slightly more bulky cis-[Pd(NO3)2(tmen)] (tmen=N,N,N',N'-tetramethylethylene diamine) instead of cis-[Pd(NO3)2(pn)] gave non-interlocked [{Pd(tmen)}6(pz)3(tma)2][NO3]6 (6) exclusively. Complexes 1, 3, and 5 represent the first examples of template-free triply interlocked molecular prisms obtained through multicomponent self-assembly. Formation of the complexes was supported by IR and multinuclear NMR (1H and 13C) spectroscopy. Formation of guest-encapsulating complexes (2 and 4) was confirmed by 2D DOSY and ROESY NMR spectroscopic analyses, whereas for complexes 1, 3, 5, and 6 single-crystal X-ray diffraction techniques unambiguously confirmed their formation. The gross geometries of H3tma-encapsulating complexes 2 and 4 were obtained by universal force field (UFF) simulations.

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.

Two-axis force sensing and control of a Reorientable scanning probe

This paper presents the design and implementation of a reorientable scanning probe that is capable of two-axis force sensing and control in the 2-D scanning (X-Z) plane. The probe is comprised of three major components, namely a compliant manipulator, laser measurement system, and magnetic actuation system. Control of the position and orientation of the probe tip is realized by means of magnetic actuation combined with a novel structural design. The design of the manipulator's compliance and that of the optical path of the laser measurement system together enable achieving sensitivity to lateral (X) forces that is nearly identical to that of normal (Z) forces. The achieved sensitivity ratio, of about 0.6, is significantly higher than that of conventional scanning probe systems. The theoretical bases for the structural design and the sensitivity of the two-axis force sensing system are presented. Subsequently, fabrication of the manipulator is described and the result of experimental evaluation of the scanning probe's features is discussed. The scanning probe is used to access the vertical and re-entrant features on the two sides of a cylindrical micropipette, which are subsequently scanned by regulating the lateral force of tip-sample interaction.

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.