Conformational study of valinomycin: a molecular dynamics approach

Valinomycin is a highly flexible cyclic dodecadepsipeptide that transports ions across membranes. Such a flexibility in the conformation is required for its biological function since it has to encounter a variety of environments and liganding state. Exploration of conformational space of this molecule is therefore important and is one of the objectives of the present study that has been carried out by means of high temperature Molecular Dynamics. Further, the stability of the known bracelet-like structure of the uncomplexed valinomycin and the inherent flexibility around this structure has been investigated. The uncomplexed form of valinomycin has been simulated at 75–100 K for 1 ns in order to elucidate the average conformational properties. An alanine-analog of valinomycin has been simulated under identical conditions in order to evaluate the effect of sidechain on the conformational properties, The studies confirm the effect of sidechain on conformational equilibrium.

Buckling, stiffening, and negative dissipation in the dynamics of a biopolymer in an active medium

We present a generic theory for the dynamics of a stiff filament under tension, in an active medium with orientational correlations, such as a microtubule in contractile actin. In sharp contrast to the case of a passive medium, we find the filament can stiffen, and possibly oscillate or buckle, depending on both the contractile or tensile nature of the activity and the filament-medium anchoring interaction. We also demonstrate a strong violation of the fluctuation-dissipation (FD) relation in the effective dynamics of the filament, including a negative FD ratio. Our approach is also of relevance to the dynamics of axons, and our model equations bear a remarkable formal similarity to those in recent work [Martin P, Hudspeth AJ, Juelicher F (2001) Proc Natl Acad Sci USA 98: 14380-14385] on auditory hair cells. Detailed tests of our predictions can be made by using a single filament in actomyosin extracts or bacterial suspensions.

A lower bound to the survival probability and an approximate first passage time distribution for Markovian and non-Markovian dynamics in phase space

We derive a very general expression of the survival probability and the first passage time distribution for a particle executing Brownian motion in full phase space with an absorbing boundary condition at a point in the position space, which is valid irrespective of the statistical nature of the dynamics. The expression, together with the Jensen's inequality, naturally leads to a lower bound to the actual survival probability and an approximate first passage time distribution. These are expressed in terms of the position-position, velocity-velocity, and position-velocity variances. Knowledge of these variances enables one to compute a lower bound to the survival probability and consequently the first passage distribution function. As examples, we compute these for a Gaussian Markovian process and, in the case of non-Markovian process, with an exponentially decaying friction kernel and also with a power law friction kernel. Our analysis shows that the survival probability decays exponentially at the long time irrespective of the nature of the dynamics with an exponent equal to the transition state rate constant.

Thermal anomalies in ternary Ge42-xPbxSe58 glasses near the charge carrier reversal threshold

The carrier type reversal (CTR) from p- to n-type in semiconducting chalcogenide glasses is an important and a long standing problem in glass science. Ge-Se glasses exhibit CTR when the metallic elements Bi and Pb are added. For example, bulk Ge42-xSe58Pbx glasses exhibit CTR around 8-9 at. % of Pb. These glasses have been prepared by melt quenching method. Glass transition temperature (T-g), Specific heat change between the liquid and the glassy states (Delta C-p) at T-g and the nonreversing heat flow (Delta H-nr) measured by modulated differential scanning calorimetry exhibit anomalies at 9 at. % of Pb. These observed anomalies are interpreted on the basis of the nano scale phase separation occurring in these glasses.

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.

Evolutionary Significance of Self-acylation Property in Acyl Carrier Proteins

Acyl carrier protein is an integral component of many cellular metabolic processes. A number of studies have reported self-acylation behavior in acyl carrier proteins. Although AM exhibit high levels of similarity in their primary and tertiary structures, self-acylation behavior is restricted to only some ACPs that can be classified into two major families based on their function. The first family of ACPs is involved in polyketide biosynthesis, whereas the second family participates in fatty acid synthesis. Facilitated by the growing number of genome sequences available for analyses, large-scale phylogenetic studies were used in these studies to uncover as to how self-acylation behavior of acyl carrier proteins is linked with the evolution of metabolic pathways in organisms. These studies show that self-acylation behavior in acyl carrier proteins was lost during the course of evolution, with certain organisms and organelles viz. plastids, retaining it for specified functions. (C) 2009 IUBMB IUBMB Life, 61(8): 853-859, 2009

Structural Insights into the acyl-intermediates of plasmodium falciparum fatty acid synthesis pathway; the mechanism of expansion of acyl carrier protein core

Acyl carrier protein (ACP) plays a central role in fatty acid biosynthesis. However, the molecular machinery that mediates its function is not yet fully understood. Therefore, structural studies were carried out on the acyl-ACP intermediates of Plasmodium falciparum using NMR as a spectroscopic probe. Chemical shift perturbation studies put forth a new picture of the interaction of ACP molecule with the acyl chain, namely, the hydrophobic core can protect up to 12 carbon units, and additional carbons protrude out from the top of the hydrophobic cavity. The latter hypothesis stems from chemical shift changes observed in C-alpha and C-beta of Ser-37 in tetradecanoyl-ACP. C-13, N-15-Double-filtered nuclear Overhauser effect (NOE) spectroscopy experiments further substantiate the concept; in octanoyl (C-8)- and dodecanoyl (C-12)-ACP, a long range NOE is observed within the phosphopantetheine arm, suggesting an arch-like conformation. This NOE is nearly invisible in tetradecanoyl (C-14)-ACP, indicating a change in conformation of the prosthetic group. Furthermore, the present study provides insights into the molecular mechanism of ACP expansion, as revealed from a unique side chain-to-backbone hydrogen bond between two fairly conserved residues, Ile-55 HN and Glu-48 O. The backbone amide of Ile-55 HN reports a pK(a) value for the carboxylate, similar to 1.9 pH units higher than model compound value, suggesting strong electrostatic repulsion between helix II and helix III. Charge-charge repulsion between the helices in combination with thrust from inside due to acyl chain would energetically favor the separation of the two helices. Helix III has fewer structural restraints and, hence, undergoes major conformational change without altering the overall-fold of P. falciparum ACP.