Dynamics of loop formation in a semiflexible polymer

The dynamics of loop formation by linear polymer chains has been a topic of several theoretical and experimental studies. Formation of loops and their opening are key processes in many important biological processes. Loop formation in flexible chains has been extensively studied by many groups. However, in the more realistic case of semiflexible polymers, not much results are available. In a recent study [K. P. Santo and K. L. Sebastian, Phys. Rev. E 73, 031923 (2006)], we investigated opening dynamics of semiflexible loops in the short chain limit and presented results for opening rates as a function of the length of the chain. We presented an approximate model for a semiflexible polymer in the rod limit based on a semiclassical expansion of the bending energy of the chain. The model provided an easy way to describe the dynamics. In this paper, using this model, we investigate the reverse process, i.e., the loop formation dynamics of a semiflexible polymer chain by describing the process as a diffusion-controlled reaction. We make use of the ``closure approximation'' of Wilemski and Fixman [G. Wilemski and M. Fixman, J. Chem. Phys. 60, 878 (1974)], in which a sink function is used to represent the reaction. We perform a detailed multidimensional analysis of the problem and calculate closing times for a semiflexible chain. We show that for short chains, the loop formation time tau decreases with the contour length of the polymer. But for longer chains, it increases with length obeying a power law and so it has a minimum at an intermediate length. In terms of dimensionless variables, the closing time is found to be given by tau similar to L-n exp(const/L), where n=4.5-6. The minimum loop formation time occurs at a length L-m of about 2.2-2.4. These are, indeed, the results that are physically expected, but a multidimensional analysis leading to these results does not seem to exist in the literature so far.

Allostery and conformational free energy changes in human tryptophanyl-tRNA synthetase from essential dynamics and structure networks

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of llostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side-chains and then to compute the ligand-induced population shift. Finally, we obtain the free-energy landscape of the protein in equilibrium, characterizing the free-energy minima accessed by the protein complexes. We have chosen human tryptophanyl-tRNA synthetase (hTrpRS), a protein esponsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure.

The dynamics of magnetically trapped fluids. I - Implications for umbral dots and penumbral grains

A study of the magnetohydrodynamic system in which a nonmagnetized fluid in a gravitational field is surrounded by a fluid carrying a vertical magnetic field is presented. It is pointed out that this study can throw some light on the fine-structural features of a sunspot. The equilibrium configuration of the field-free fluid is a tapering column ending at an apex. The regions away form the apex can be studied by the slender flux tube approximation. A scheme developed to treat the apex indicates that, just below the apex, the radius of the tapering column opens up with a 3/2 power dependence on the depth below the apex. If the internal pressure of the field-free fluid is increased, the apex rises, and a static equilibrium may not be possible beyond a limit if the magnetic pressure drops quickly above a certain height. The nature of steady-flow solutions beyond this limit is investigated. Under conditions inside a sunspot, a column of field-free gas is found to rise with a velocity of about 100 km/hr. If umbral dots and penumbral grains are interpreted as regions where the field-free gas ultimately emerges, a very natural explanation of most of their observed properties is obtained.

A Model of Anomalous Chain Translocation Dynamics

A model of polymer translocation based on the stochastic dynamics of the number of monomers on one side of a pore-containing surface is formulated in terms of a one-dimensional generalized Langevin equation, in which the random force is assumed to be characterized by long-ranged temporal correlations. The model is introduced to rationalize anomalies in measured and simulated values of the average time of passage through the pore, which in general cannot be satisfactorily accounted for by simple Brownian diffusion mechanisms. Calculations are presented of the mean first passage time for barrier crossing and of the mean square displacement of a monomeric segment, in the limits of strong and weak diffusive bias. The calculations produce estimates of the exponents in various scaling relations that are in satisfactory agreement with available data.

An Insight to the Dynamics of Conserved Water Molecular Triad in IMPDH II (Human): Recognition of Cofactor and Substrate to Catalytic Arg 322

Inosine 5' monophosphate dehydrogenase (IMPDH II) is a key enzyme involved in the de novo biosynthesis pathway of purine nucleotides and is also considered to be an excellent target for cancer inhibitor design. The conserve R 322 residue (in human) is thought to play some role in the recognition of inhibitor and cofactor through the catalytic D 364 and N 303. The 15 ns simulation and the water dynamics of the three different PDB structures (1B3O, 1NF7, and 1NFB) of human IMPDH by CHARMM force field have clearly indicated the involvement of three conserved water molecules (W-L, W-M, and W-C) in the recognition of catalytic residues (R 322, D 364, and N 303) to inhibitor and cofactor. Both the guanidine nitrogen atoms (NH1 and NH 2) of the R 322 have anchored the di- and mono-nucleotide (cofactor and inhibitor) binding domains via the conserved W-C and W-L water molecules. Another conserved water molecule W-M seems to bridge the two domains including the R 322 and also the W-C and W-L through seven centers H-bonding coordination. The conserved water molecular triad (W-C - W-M - W-L) in the protein complex may thought to play some important role in the recognition of inhibitor and cofactor to the protein through R 322 residue.

A multistep transversal linearization (MTL) method in non-linear dynamics through a Magnus characterization

A new form of a multi-step transversal linearization (MTL) method is developed and numerically explored in this study for a numeric-analytical integration of non-linear dynamical systems under deterministic excitations. As with other transversal linearization methods, the present version also requires that the linearized solution manifold transversally intersects the non-linear solution manifold at a chosen set of points or cross-section in the state space. However, a major point of departure of the present method is that it has the flexibility of treating non-linear damping and stiffness terms of the original system as damping and stiffness terms in the transversally linearized system, even though these linearized terms become explicit functions of time. From this perspective, the present development is closely related to the popular practice of tangent-space linearization adopted in finite element (FE) based solutions of non-linear problems in structural dynamics. The only difference is that the MTL method would require construction of transversal system matrices in lieu of the tangent system matrices needed within an FE framework. The resulting time-varying linearized system matrix is then treated as a Lie element using Magnus’ characterization [W. Magnus, On the exponential solution of differential equations for a linear operator, Commun. Pure Appl. Math., VII (1954) 649–673] and the associated fundamental solution matrix (FSM) is obtained through repeated Lie-bracket operations (or nested commutators). An advantage of this approach is that the underlying exponential transformation could preserve certain intrinsic structural properties of the solution of the non-linear problem. Yet another advantage of the transversal linearization lies in the non-unique representation of the linearized vector field – an aspect that has been specifically exploited in this study to enhance the spectral stability of the proposed family of methods and thus contain the temporal propagation of local errors. A simple analysis of the formal orders of accuracy is provided within a finite dimensional framework. Only a limited numerical exploration of the method is presently provided for a couple of popularly known non-linear oscillators, viz. a hardening Duffing oscillator, which has a non-linear stiffness term, and the van der Pol oscillator, which is self-excited and has a non-linear damping term.

Buckled nano rod - a two state system: quantum effects on its dynamics

We consider a suspended elastic rod under longitudinal compression. The compression can be used to adjust potential energy for transverse displacements from the harmonic to the double well regime. The two minima in potential energy curve describe two possible buckled states. Using transition state theory (TST) we have calculated the rate of conversion from one state to other. If the strain epsilon = 4 epsilon c the simple TST rate diverges. We suggest a method to correct this divergence for quantum calculations. We also find that zero point energy contributions can be quite large so that single mode calculations can lead to large errors in the rate.

Dynamics and control of buckling type devices using SMA wire integrated beam

An analysis and design study using Shape Memory Alloy (SMA) wire integrated beam and its buckling shape control are reported. The dynamical system performance is analyzed with a mathematical set-up involving nonlocal and rate sensitive kinetics of phase transformation in the SMA wire. A standard phenomenological constitutive model reported by Brinson (1993) is modified by considering certain consistency conditions in the material property tensors and by eliminating spurious singularity. Considering the inhomogeneity effects, a finite element model of the SMA wire is developed. Simulations are carried out to study the buckling shape control of a beam integrated with SMA wire.

Dynamics of Barrierless and Activated Chemical Reactions in a Dispersive Medium within the Fractional Diffusion Equation Approach

Barrierless chemical reactions have often been modeled as a Brownian motion on a one-dimensional harmonic potential energy surface with a position-dependent reaction sink or window located near the minimum of the surface. This simple (but highly successful) description leads to a nonexponential survival probability only at small to intermediate times but exponential decay in the long-time limit. However, in several reactive events involving proteins and glasses, the reactions are found to exhibit a strongly nonexponential (power law) decay kinetics even in the long time. In order to address such reactions, here, we introduce a model of barrierless chemical reaction where the motion along the reaction coordinate sustains dispersive diffusion. A complete analytical solution of the model can be obtained only in the frequency domain, but an asymptotic solution is obtained in the limit of long time. In this case, the asymptotic long-time decay of the survival probability is a power law of the Mittag−Leffler functional form. When the barrier height is increased, the decay of the survival probability still remains nonexponential, in contrast to the ordinary Brownian motion case where the rate is given by the Smoluchowski limit of the well-known Kramers' expression. Interestingly, the reaction under dispersive diffusion is shown to exhibit strong dependence on the initial state of the system, thus predicting a strong dependence on the excitation wavelength for photoisomerization reactions in a dispersive medium. The theory also predicts a fractional viscosity dependence of the rate, which is often observed in the reactions occurring in complex environments.

Promoting sun safety in the workplace – what works?

Objective Australia has one of the highest skin cancer incidence and mortality rates in the world. Outdoor workers are a high risk group. Australian workplaces are undergoing large scale safety related changes, yet the mandate to provide specific sun safe practices remains absent. With much of the previous research aiming to improve sun safety in the workplace being quantitative in nature, relatively little is known about why certain sun safe strategies will or will not be successful in workplaces. Methods This qualitative article explores the enablers and barriers identified during an 18-month mixed methods project conducted in Queensland, Australia which aimed to improve workplace sun safe interventions. Results A variety of key enablers and barriers to implementing sun safe interventions in the workplace were identified, including presence of an engaged workplace champion, ownership and innovation by the workers. Conclusions These findings were part of a broader integration of interlinked qualitative and quantitative methods to yield a more complete picture of the determinants of the issue, implementation process and likelihood of changes at the workplace. Implications The paper provides guidance for public and occupational health practitioners on the selection of the most promising strategies when assisting workplaces to become sun safe.

Truckies and health promotion: A “hard-to-reach” group without a “proper” workplace.

This literature review was undertaken to inform a settings based health promotion research project, conducted by a public health research team at the Queensland University of Technology (QUT). The aim of this project is to identify how transport workplaces can support their truck drivers to access healthy food options and increase physical activity. Truck drivers in Australia are at increased risk of numerous chronic diseases, in part due to the restrictions placed upon them by the environment in which they work. Barriers to good health through adequate nutrition and physical activity are the result of a complex interaction between government regulations, corporate policies, the built environment and individual factors. Few interventions target this population in a meaningful and sustainable way, though evidence exists for interventions which can be translated into truck drivers working environment.

Workers are doing it for themselves? Local labor solidarity in precarious times

Historically, organized labor has played a fundamental role in guaranteeing basic rights and privileges for screen media workers and defending union and guild members (however unevenly) from egregious abuses of power. Yet, despite the recent turn to labor in media and cultural studies, organized labor today has received only scant attention, even less so in locations outside Hollywood. This presentation thus intervenes in two significant ways: first, it acknowledges the ongoing global ‘undoing’ of organized labor as a consequence of footloose production and conglomeration within the screen industries, and second, it examines a case example of worker solidarity and political praxis taking shape outside formal labor institutions in response to those structural shifts. Accordingly, it links an empirical study of individual agency to broader debates associated with the spatial dynamics of screen media production, including local capacity, regional competition, and precariousness. Drawing from ethnographic interviews with local film and television workers in Glasgow, Scotland, I consider the political alliance among three nascent labor organizations in the city: one for below-the-line crew, one for facility operators, and (oddly enough) one for producers. Collectively, the groups share a desire to transform Glasgow into a global production hub, following the infrastructure developments in nearby cities like Belfast, Prague, and Budapest. They furthermore frame their objectives in political terms: establishing global scale is considered a necessary maneuver to improve local working conditions like workplace safety, income disparity, skills training, and job access. Ultimately, I argue these groups are a product of an inadequate union structure and outdated policy vision for the screen sector , once-supportive institutions currently out of sync with the global realities of media production. Furthermore, the groups’ advocacy efforts reveal the extent to which workers themselves (in additional to capital) can seek “spatial fixes” to suture their prospects to specific political and economic goals. Of course, such activities manifest under conditions outside of the workers’ control but nevertheless point to an important tension within capitalist social relations, namely that the agency to reshape the spatial relationships in their own lives recasts the geography of labor in terms that aren’t inherent or exclusive to the interests of global capital.

Insight into the early stages of thermal unfolding of peanut agglutinin by molecular dynamics simulations

Peanut agglutinin is a homotetrameric nonglycosylated protein. The protein has a unique open quaternary structure. Molecular dynamics simulations have been employed follow the atomistic details of its unfolding at different temperatures. The early events of the deoligomerization of the protein have been elucidated in the present study. Simulation trajectories of the monomer as well as those of the tetramer have been compared and the tetramer is found to be substantially more stable than its monomeric counterpart. The tetramer shows retention of most of its.. secondary structure but considerable loss of the tertiary structure at high temperature. e generation of a This observation impies the molten globule-like intermediate in the later stages of deoligomerization. The quaternary structure of the protein has weakened to a large extent, but none of the subunits are separated. In addition, the importance of the metal-binding to the stability of the protein structure has also been investigated. Binding of the metal ions not only enhances the local stability of the metal-ion binding loop, but also imparts a global stability to the overall structure. The dynamics of different interfaces vary significantly as probed through interface clusters. The differences are substantially enhanced at higher temperatures. The dynamics and the stability of the interfaces have been captured mainly by cluster analysis, which has provided detailed information on the thermal deoligomerization of the protein.