Molecular dynamics investigation on shearing between osteopontin and hydroxyapatite in biological materials

Bone, a hard biological material, possesses a combination of high stiffness and toughness, even though the main basic building blocks of bone are simply mineral platelets and protein molecules. Bone has a very complex microstructure with at least seven hierachical levels. This unique material characteristic attracts great attention, but the deformation mechanisms in bone have not been well understood. Simulation at nano-length scale such as molecular dynamics (MD) is proven to be a powerful tool to investigate bone nanomechanics for developing new artificial biological materials. This study focuses on the ultra large and thin layer of extrafibrillar protein matrix (thickness = ~ 1 nm) located between mineralized collagen fibrils (MCF). Non-collagenous proteins such as osteopontin (OPN) can be found in this protein matrix, while MCF consists mainly of hydroxyapatite (HA) nanoplatelets (thickness = 1.5 – 4.5 nm). By using molecular dynamics method, an OPN peptide was pulled between two HA mineral platelets with water in presence. Periodic boundary condition (PBC) was applied. The results indicate that the mechanical response of OPN peptide greatly depends on the attractive electrostatics interaction between the acidic residues in OPN peptide and HA mineral surfaces. These bonds restrict the movement of OPN peptide, leading to a high energy dissipation under shear loading.

Molecular dynamics simulation of mechanical behavior of osteopontin-hydroxyapatite interfaces

Bone is characterized with an optimized combination of high stiffness and toughness. The understanding of bone nanomechanics is critical to the development of new artificial biological materials with unique properties. In this work, the mechanical characteristics of the interfaces between osteopontin (OPN, a noncollagenous protein in extrafibrillar protein matrix) and hydroxyapatite (HA, a mineral nanoplatelet in mineralized collagen fibrils) were investigated using molecular dynamics method. We found that the interfacial mechanical behaviour is governed by the electrostatic attraction between acidic amino acid residues in OPN and calcium in HA. Higher energy dissipation is associated with the OPN peptides with a higher number of acidic amino acid residues. When loading in the interface direction, new bonds between some acidic residues and HA surface are formed, resulting in a stick-slip type motion of OPN peptide on the HA surface and high interfacial energy dissipation. The formation of new bonds during loading is considered to be a key mechanism responsible for high fracture resistance observed in bone and other biological materials.

Comprehensive contribution of filament thickness and crosslinker failure to the rheological property of F-actin cytoskeleton

Rheological property of F-actin cytoskeleton is significant to the restructuring of cytoskeleton under a variety of cell activities. This study numerically validates the rheological property of F-actin cytoskeleton is not only a result of kinetic energy dissipation of F-actin, but also greatly depends on the configuration remodeling of networks structure. Both filament geometry and crosslinker properties can affect the remodeling of F-actin cytoskeleton. The crosslinker unbinding is found to dissipate energy and induce prominent stress relaxation in the F-actin adjacent to cross-linkages. Coupled with F-actin elasticity, the energy dissipation and stress relaxation are more significant in bundled F-actin networks than in single F-actin networks.

Spatial B-jump at sudden channel enlargements with abrupt drop

An analytical and experimental study of the hydraulic jump in stilling basins with abrupt drop and sudden enlargement, called the spatial B-jump here, is carried out for finding the sequent depth ratio and resulting energy dissipation. The spatial B-jump studied has its toe downstream of the expansion section, and the stream lines at the toe are characterized by downward curvature. An expression is obtained for the sequent depth ratio based on the momentum equation with suitable assumptions for the extra pressure force term because of the abrupt drop in the bed and sudden enlargement in the basin width. Predictions compare favorably with experiments. It is shown that the spatial B-jump needs less tailwater depth, thereby enhancing the stability of the jump when compared either with spatial jump, which forms in sudden expanding channels, or with B-jump, which forms in a channel with an abrupt drop in bed. It is also shown that there is a significant increase in relative energy loss for the spatial B-jump compared to either the spatial jump or B-jump alone.

Coalescence of rigid droplets in a stirred dispersion-II. Band-limited force fluctuations

The coalescence of nearly rigid liquid droplets in a turbulent flow field is viewed as the drainage of a thin film of liquid under the action of a stochastic force representing the effect of turbulence. The force squeezing the drop pair is modelled as a correlated random function of time. The drops are assumed to coalesce once the film thickness becomes smaller than a critical thickness while they are regarded as separated if their distance of separation is larger than a prescribed distance. A semi-analytical solution is derived to determine the coalescence efficiency. The veracity of the solution procedure is established via a Monte-Carlo solution scheme. The model predicts a reversing trend of the dependence of the coalescence efficiency on the drop radii, the film liquid viscosity and the turbulence energy dissipation per unit mass, as the relative fluctuation increases. However, the dependence on physical parameters is weak (especially at high relative fluctuation) so that for the smallest droplets (which are nearly rigid) the coalescence efficiency may be treated as an empirical constant. The predictions of this model are compared with those of a white-noise force model. The results of this paper and those in Muralidhar and Ramkrishna (1986, Ind. Engng Chem. Fundam. 25, 554-56) suggest that dynamic drop deformation is the key factor that influences the coalescence efficiency.

Studying the diurnal and seasonal acclimation of photosystem Ii using chlorophyll-a fluorescence

A small fraction of the energy absorbed in the light reactions of photosynthesis is re-emitted as chlorophyll-a fluorescence. Chlorophyll-a fluorescence and photochemistry compete for excitation energy in photosystem II (PSII). Therefore, changes in the photochemical capacity can be detected through analysis of chlorophyll fluorescence. Chlorophyll fluorescence techniques have been widely used to follow the diurnal (fast), and the seasonal (slow) acclimation in the energy partitioning between photochemical and non-photochemical processes in PSII. Energy partitioning in PSII estimated through chlorophyll fluorescence can be used as a proxy of the plant physiological status, and measured at different spatial and temporal scales. However, a number of technical and theoretical limitations still limit the use of chlorophyll fluorescence data for the study of the acclimation of PSII. The aim of this Thesis was to study the diurnal and seasonal acclimation of PSII in field conditions through the development and testing of new chlorophyll fluorescence-based tools, overcoming these limitations. A new model capable of following the fast acclimation of PSII to rapid fluctuations in light intensity was developed. The model was used to study the rapid acclimation in the electron transport rate under fluctuating light. Additionally, new chlorophyll fluorescence parameters were developed for estimating the seasonal acclimation in the sustained rate constant of thermal energy dissipation and photochemistry. The parameters were used to quantitatively evaluate the effect of light and temperature on the seasonal acclimation of PSII. The results indicated that light environment not only affected the degree but also the kinetics of response of the acclimation to temperature, which was attributed to differences in the structural organization of PSII during seasonal acclimation. Furthermore, zeaxanthin-facilitated thermal dissipation appeared to be the main mechanisms modulating the fraction of absorbed energy being dissipated thermally during winter in field Scots pine. Finally, the integration between diurnal and seasonal acclimation mechanisms was studied using a recently developed instrument MONI-PAM (Walz GmbH, Germany) capable of continuously monitoring the energy partitioning in PSII.

Dynamics of dense sheared granular flows. Part II. The relative velocity distributions

The distribution of relative velocities between colliding particles in shear flows of inelastic spheres is analysed in the Volume fraction range 0.4-0.64. Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to line joining the centres). The distribution or pre-collisional normal relative velocities (along the line Joining the centres of the particles) is Found to be an exponential distribution for particles with low normal coefficient of restitution in the range 0.6-0.7. This is in contrast to the Gaussian distribution for the normal relative velocity in all elastic fluid in the absence of shear. A composite distribution function, which consists of an exponential and a Gaussian component, is proposed to span the range of inelasticities considered here. In the case of roughd particles, the relative velocity tangential to the surfaces at contact is also evaluated, and it is found to be close to a Gaussian distribution even for highly inelastic particles.Empirical relations are formulated for the relative velocity distribution. These are used to calculate the collisional contributions to the pressure, shear stress and the energy dissipation rate in a shear flow. The results of the calculation were round to be in quantitative agreement with simulation results, even for low coefficients of restitution for which the predictions obtained using the Enskog approximation are in error by an order of magnitude. The results are also applied to the flow down an inclined plane, to predict the angle of repose and the variation of the volume fraction with angle of inclination. These results are also found to be in quantitative agreement with previous simulations.

On the coalescence of twisted flux tubes

We study the problem of the coalescence of twisted flux tubes by assuming that the azimuthal field lines reconnect at a current sheet during the coalescence process and everywhere else the magnetic field is frozen in the fluid. We derive relations connecting the topology of the coalesced flux tube with the topologies of the initial flux tubes, and then obtain a structure equation for calculating the field configuration of the coalesced flux tube from the given topology. Some solutions for the two extreme cases of low-β plasma and high-β plasma are discussed. The coalesced flux tube has less twist than the initial flux tube. Magnetic helicity is found to be exactly conserved during the coalescence, but the assumptions in the model put a constraint on the energy dissipation so that we do not get a relaxation to the minimum-energy Taylor state in the low-β case. It is pointed out that the structure equation connecting the topology and the equilibrium configuration is quite general and can be of use in many two-dimensional flux tube problems.

Pressure and thermally induced stages of wear in dry sliding of a steel ball against an aluminium-silicon alloy flat

Wear of etched near-eutectic aluminium silicon alloy slid against a steel ball under ambient is explored. The sliding velocity is kept low (0.01 m/s) and the nominal contact pressure is varied in a 15-40 MPa range. Four stages of wear are identified; ultra mild wear, mild wear, severe wear and post severe oxidative wear. The first transition is controlled by the protrusions of silicon particles, projecting out of the aluminium alloy matrix. Once these protrusions disappear under pressure and sliding, oxidation and bulk energy dissipation mechanisms take over to institute transitions to other stages of wear. The phenomenological characteristics of wear stages are explored using a variety of techniques including nanoindentation, focused ion beam milling, electron microscopy, X-ray photoelectron spectroscopy (XPS), energy dispersive X-ray spectroscopy (EDS) and optical interferometry. (c) 2010 Elsevier B.V. All rights reserved.

Structural studies of decaying fluid turbulence: Effect of initial conditions

We present results from a systematic numerical study of structural properties of an unforced, incompressible, homogeneous, and isotropic three-dimensional turbulent fluid with an initial energy spectrum that develops a cascade of kinetic energy to large wave numbers. The results are compared with those from a recently studied set of power-law initial energy spectra [C. Kalelkar and R. Pandit, Phys. Rev. E 69, 046304 (2004)] which do not exhibit such a cascade. Differences are exhibited in plots of vorticity isosurfaces, the temporal evolution of the kinetic energy-dissipation rate, and the rates of production of the mean enstrophy along the principal axes of the strain-rate tensor. A crossover between "non-cascade-type" and "cascade-type" behavior is shown numerically for a specific set of initial energy spectra.

The hard-particle model for a dense granular flow down an inclined plane

Perfectly hard particles are those which experience an infinite repulsive force when they overlap, and no force when they do not overlap. In the hard-particle model, the only static state is the isostatic state where the forces between particles are statically determinate. In the flowing state, the interactions between particles are instantaneous because the time of contact approaches zero in the limit of infinite particle stiffness. Here, we discuss the development of a hard particle model for a realistic granular flow down an inclined plane, and examine its utility for predicting the salient features both qualitatively and quantitatively. We first discuss Discrete Element simulations, that even very dense flows of sand or glass beads with volume fraction between 0.5 and 0.58 are in the rapid flow regime, due to the very high particle stiffness. An important length scale in the shear flow of inelastic particles is the `conduction length' delta = (d/(1 - e(2))(1/2)), where d is the particle diameter and e is the coefficient of restitution. When the macroscopic scale h (height of the flowing layer) is larger than the conduction length, the rates of shear production and inelastic dissipation are nearly equal in the bulk of the flow, while the rate of conduction of energy is O((delta/h)(2)) smaller than the rate of dissipation of energy. Energy conduction is important in boundary layers of thickness delta at the top and bottom. The flow in the boundary layer at the top and bottom is examined using asymptotic analysis. We derive an exact relationship showing that the a boundary layer solution exists only if the volume fraction in the bulk decreases as the angle of inclination is increased. In the opposite case, where the volume fraction increases as the angle of inclination is increased, there is no boundary layer solution. The boundary layer theory also provides us with a way of understanding the cessation of flow when at a given angle of inclination when the height of the layer is decreased below a value h(stop), which is a function of the angle of inclination. There is dissipation of energy due to particle collisions in the flow as well as due to particle collisions with the base, and the fraction of energy dissipation in the base increases as the thickness decreases. When the shear production in the flow cannot compensate for the additional energy drawn out of the flow due to the wall collisions, the temperature decreases to zero and the flow stops. Scaling relations can be derived for h(stop) as a function of angle of inclination.

Dynamics of sheared inelastic dumbbells

We study the dynamical properties of the homogeneous shear flow of inelastic dumbbells in two dimensions as a first step towards examining the effect of shape on the properties of flowing granular materials. The dumbbells are modelled as smooth fused disks characterized by the ratio of the distance between centres (L) and the disk diameter (D), with an aspect ratio (L/D) varying between 0 and 1 in our simulations. Area fractions studied are in the range 0.1-0.7, while coefficients of normal restitution (e(n)) from 0.99 to 0.7 are considered. The simulations use a modified form of the event-driven methodology for circular disks. The average orientation is characterized by an order parameter S, which varies between 0 (for a perfectly disordered fluid) and 1 (for a fluid with the axes of all dumbbells in the same direction). We investigate power-law fits of S as a function of (L D) and (1 - e(n)(2)) There is a gradual increase in ordering as the area fraction is increased, as the aspect ratio is increased or as the coefficient of restitution is decreased. The order parameter has a maximum value of about 0.5 for the highest area fraction and lowest coefficient of restitution considered here. The mean energy of the velocity fluctuations in the flow direction is higher than that in the gradient direction and the rotational energy, though the difference decreases as the area fraction increases, due to the efficient collisional transfer of energy between the three directions. The distributions of the translational and rotational velocities are Gaussian to a very good approximation. The pressure is found to be remarkably independent of the coefficient of restitution. The pressure and dissipation rate show relatively little variation when scaled by the collision frequency for all the area fractions studied here, indicating that the collision frequency determines the momentum transport and energy dissipation, even at the lowest area fractions studied here. The mean angular velocity of the particles is equal to half the vorticity at low area fractions, but the magnitude systematically decreases to less than half the vorticity as the area fraction is increased, even though the stress tensor is symmetric.

Rocking response of rectangular rigid blocks under random noise base excitations

The response of a rigid rectangular block resting on a rigid foundation and acted upon simultaneously by a horizontal and a vertical random white-noise excitation is considered. In the equation of motion, the energy dissipation is modeled through a viscous damping term. Under the assumption that the body does not topple, the steady-state joint probability density function of the rotation and the rotational velocity is obtained using the Fokker-Planck equation approach. Closed form solution is obtained for a specific combination of system parameters. A more general but approximate solution to the joint probability density function based on the method of equivalent non-linearization is also presented. Further, the problem of overturning of the block is approached in the framework of the diffusion methods for first passage failure studies. The overturning of the block is deemed incipient when the response trajectories in the phase plane cross the separatrix of the conservative unforced system. Expressions for the moments of first passage time are obtained via a series solution to the governing generalized Pontriagin-Vitt equations. Numerical results illustra- tive of the theoretical solutions are presented and their validity is examined through limited amount of digital simulations.

Microsomal redox systems in brown adipose tissue: high lipid peroxidation, low cholesterol biosynthesis and no detectable cytochrome P-450

The presence of redox systems in microsomes of brown adipose tissue (BAT) in cold exposed rats was investigated and compared with liver. BAT microsomes showed high activity of lipid peroxidation measured both by the formation of malondialdehyde (MDA) and by oxygen uptake. NADH and NADPH dependent cytochrome c reductase activity were present in both BAT and liver microsomes. Aminopyrine demethylase and aniline hydroxylase activities, the characteristic detoxification enzymes in liver microsomes could not be detected in BAT microsomes. BAT minces showed very poor incorporation of [1-14C]acetate and [2-14C]-mevalonate in unsaponifiable lipid fraction compared to liver. Biosynthesis of cholesterol and ubiquinone, but not fatty acids, and the activity of 3-hydroxy-3-methyl glutaryl CoA reductase appear to be very low in BAT. Examination of difference spectra showed the presence of only cytochrome b 5 in BAT microsomes. In addition to the inability to detect the enzyme activities dependent on cytochrome P-450, a protein with the characteristic spectrum, molecular size in SDS-PAGE and interaction with antibodies in double diffusion test, also could not be detected in BAT microsomes. The high activity of lipid peroxidation in microsomes, being associated with large oxygen uptake and oxidation of NADPH, will also contribute to the energy dissipation as heat in BAT, considered important in thermogenesis.

Direct numerical simulations of statistically steady, homogeneous, isotropic fluid turbulence with polymer additives

We carry out a direct numerical simulation (DNS) study that reveals the effects of polymers on statistically steady, forced, homogeneous, and isotropic fluid turbulence. We find clear manifestations of dissipation-reduction phenomena: on the addition of polymers to the turbulent fluid, we obtain a reduction in the energy dissipation rate; a significant modification of the fluid-energy spectrum, especially in the deep-dissipation range; and signatures of the suppression of small-scale structures, including a decrease in small-scale vorticity filaments. We also compare our results with recent experiments and earlier DNS studies of decaying fluid turbulence with polymer additives.