Temperature dependent free energy surface of polymer folding from equilibrium and quench studies

Langevin dynamics simulation studies have been employed to calculate the temperature dependent free energy surface and folding characteristics of a 500 monomer long linear alkane (polyethylene) chain with a realistic interaction potential. Both equilibrium and temperature quench simulation studies have been carried out. Using the shape anisotropy parameter (S) of the folded molecule as the order parameter, we find a weakly first order phase transition between the high-temperature molten globule and low-temperature rodlike crystalline states separated by a small barrier of the order of k(B)T. Near the melting temperature (580 K), we observe an intriguing intermittent fluctuation with pronounced ``1/f noise characteristics'' between these two states with large difference in shape and structure. We have also studied the possibilities of different pathways of folding to states much below the melting point. At 300 K starting from the all-trans linear configuration, the chain folds stepwise into a very regular fourfold crystallite with very high shape anisotropy. Whereas, when quenched from a high temperature (900 K) random coil regime, we identify a two step transition from the random coiled state to a molten globulelike state and, further, to a anisotropic rodlike state. The trajectory reveals an interesting coupling between the two order parameters, namely, radius of gyration (R-g) and the shape anisotropy parameter (S). The rodlike final state of the quench trajectory is characterized by lower shape anisotropy parameter and significantly larger number of gauche defects as compared to the final state obtained through equilibrium simulation starting from all-trans linear chain. The quench study shows indication of a nucleationlike pathway from the molten globule to the rodlike state involving an underlying rugged energy landscape. (C) 2010 American Institute of Physics. doi:10.1063/1.3509398]

Out-of-equilibrium microrheology using optical tweezers to probe directional viscoelastic properties under shear

Many wormlike micellar systems exhibit appreciable shear thinning due to shear-induced alignment. As the micelles get aligned introducing directionality in the system, the viscoelastic properties are no longer expected to be isotropic. An optical-tweezers-based active microrheology technique enables us to probe the out-of-equilibrium rheological properties of a wormlike micellar system simultaneously along two orthogonal directions-parallel to the applied shear, as well as perpendicular to it. While the displacements of a trapped bead in response to active drag force carry signature of conventional shear thinning, its spontaneous position fluctuations along the perpendicular direction manifest an orthogonal shear thickening, an effect hitherto unobserved. Copyright (C) EPLA, 2010

Lattice embedding and equilibrium geometry of metal-halogen chains in the two-band extended Hubbard model

Two-band extended Hubbard model studies show that the shift in optical gap of the metal-halogen (MX) chain upon embedding in a crystalline environment depends upon alternation in the site-diagonal electron-lattice interaction parameter (epsilon(M)) and the strength of electron-electron interactions at the metal site (U(M)). The equilibrium geometry studies on isolated chains show that the MX chains tend to distort for alternating epsilon(M) and small U(M) values.

Coupled electrostatic-elastic analysis for topology optimization using material interpolation

In this paper, we present a novel analytical formulation for the coupled partial differential equations governing electrostatically actuated constrained elastic structures of inhomogeneous material composition. We also present a computationally efficient numerical framework for solving the coupled equations over a reference domain with a fixed finite-element mesh. This serves two purposes: (i) a series of problems with varying geometries and piece-wise homogeneous and/or inhomogeneous material distribution can be solved with a single pre-processing step, (ii) topology optimization methods can be easily implemented by interpolating the material at each point in the reference domain from a void to a dielectric or a conductor. This is attained by considering the steady-state electrical current conduction equation with a `leaky capacitor' model instead of the usual electrostatic equation. This formulation is amenable for both static and transient problems in the elastic domain coupled with the quasi-electrostatic electric field. The procedure is numerically implemented on the COMSOL Multiphysics (R) platform using the weak variational form of the governing equations. Examples have been presented to show the accuracy and versatility of the scheme. The accuracy of the scheme is validated for the special case of piece-wise homogeneous material in the limit of the leaky-capacitor model approaching the ideal case.

Topology of wolfgram proteins and 2?, 3?-cyclic nucleotide 3?-phosphodiesterase in CNS myelin: Studies with proteases

The topological disposition of Wolfgram proteins (WP) and their relationship with 2', 3'-cyclic nucleotide 3'-phosphodiesterase (CNPase) in human, rat, sheep, bovine, guinea pig and chicken CNS myelin was investigated. Controlled digestion of myelin with trypsin gave a 35KDa protein band (WP-t) when electrophoresed on dodecyl sulfate-polyacrylamide gel in all species. Western blot analysis showed that the WP-t was derived from WP. WP-t was also formed when rat myelin was treated with other proteases such as kallikrein, thermolysin and leucine aminopeptidase. Staining for CNPase activity on nitrocellulose blots showed that WP-t is enzymatically active. Much of the CNPase activity remained with the membrane fraction even after treatment with high concentrations of trypsin when WP were completely hydrolysed and no protein bands with M.W > 14KDa were detected on the gels. Therefore protein fragments of WP with M.W < 14KDa may contain CNPase activity. From these results, it is suggested that the topological disposition of all the various WP is such that a 35KDa fragment is embedded in the lipid bilayer and the remaining fragment exposed at the intraperiod line in the myelin structure which may play a role in the initiation of myelinogenesis.

Solid-state photobehaviour and crystal packing of o-chlorobenzylidene-DL-piperitone: influence of molecular topology on photobehaviour

C17H19ClO, M(r) = 274.7, triclinic, P1BAR, a = 11.154 (3), b = 12.685 (2), c = 12.713 (2) angstrom, alpha = 100.68 (1), beta = 113.58 (1), gamma = 104.50 (2)-degrees, V = 1511.1 (6) angstrom3, Z = 4, D(m) = 1.22, D(x) = 1.215 Mg m-3, Cu K-alpha, lambda = 1.5418 angstrom, mu = 2.16 mm-1, F(000) = 584, T = 293 K, R = 0.057 for 3481 observed reflections. The title compound is photostable in the crystalline state and lattice-energy calculations have been employed to rationalize the photobehaviour. The well-known beta-steering ability of the chloro group is not operative in this system as there are no Cl...Cl interactions in the crystal lattice. All five benzylidene-DL-piperitone structures so far studied are alpha-packed and the molecular topology appears to be a deciding factor even in the presence of steering groups.

The Equations of Equilibrium in Orthogonal Curvilinear Reference Coordinates

Analytical solutions to problems in finite elasticity are most often derived using the semi-inverse approach along with the spatial form of the equations of motion involving the Cauchy stress tensor. This procedure is somewhat indirect since the spatial equations involve derivatives with respect to spatial coordinates while the unknown functions are in terms of material coordinates, thus necessitating the use of the chain rule. In this classroom note, we derive compact expressions for the components of the divergence, with respect to orthogonal material coordinates, of the first Piola-Kirchhoff stress tensor. The spatial coordinate system is also assumed to be an orthogonal curvilinear one, although, not necessarily of the same type as the material coordinate system. We show by means of some example applications how analytical solutions can be derived more directly using the derived results.

Parallel implementation of alternate quadrant interlocking factorisation method on star topology

This paper discusses the parallel implementation of the solution of a set of linear equations using the Alternative Quadrant Interlocking Factorisation Methods (AQIF), on a star topology. Both the AQIF and LU decomposition methods are mapped onto star topology on an IBM SP2 system, with MPI as the internode communicator. Performance parameters such as speedup, efficiency have been obtained through experimental and theoretical means. The studies demonstrate (i) a mismatch of 15% between the theoretical and experimental results, (ii) scalability of the AQIF algorithm, and (iii) faster executing AQIF algorithm.

Synthesis of workspaces of planar manipulators with arbitrary topology using shape representation and simulated annealing

This paper presents a general methodology for the synthesis of the external boundary of the workspaces of a planar manipulator with arbitrary topology. Both the desired workspace and the manipulator workspaces are identified by their boundaries and are treated as simple closed polygons. The paper introduces the concept of best match configuration and shows that the corresponding transformation can be obtained by using the concept of shape normalization available in image processing literature. Introduction of the concept of shape in workspace synthesis allows highly accurate synthesis with fewer numbers of design variables. This paper uses a new global property based vector representation for the shape of the workspaces which is computationally efficient because six out of the seven elements of this vector are obtained as a by-product of the shape normalization procedure. The synthesis of workspaces is formulated as an optimization problem where the distance between the shape vector of the desired workspace and that of the workspace of the manipulator at hand are minimized by changing the dimensional parameters of the manipulator. In view of the irregular nature of the error manifold, the statistical optimization procedure of simulated annealing has been used. A number of worked-out examples illustrate the generality and efficiency of the present method. (C) 1998 Elsevier Science Ltd. All rights reserved.

Hyperbranched Polyacetals with Tunable Degradation Rates

We report the first synthesis of hyperbranched polyacetals via a melt transacetalization polymerization process. The process proceeds via the self-condensation of an AB(2) type monomer carrying a hydroxyl group and a dimethylacetal unit; the continuous removal of low boiling methanol drives the equilibrium toward polymer formation. Because of the susceptibility of the acetal linkage to hydrolysis, the polymer degrades readily under mildly acidic conditions to yield the corresponding hydroxyl aldehyde as the primary product. Furthermore, because of the unique topology of hyperbranched structures, the rate of polymer degradation was readily tuned by changing just the nature of the end-groups alone; instead of the dimethylacetal bearing monomer, longer chain dialkylacetals (dibutyl and dihexyl) monomers yielded hyperbranched polymers carrying longer alkyl groups at their molecular periphery. The highly branched topology and the relatively high volume fraction of the terminal alkyl groups resulted in a significant lowering of the ingress rates of the aqueous reagents to the loci of degradation, and consequently the degradation rates of the polymers were dramatically influenced by the hydrophobic nature of the terminal alkyl substituents. The simple synthesis and easy tunability of the degradation rates make these materials fairly attractive candidates for use as degradable scaffolds.