Structure, dynamics, and interactions of jacalin. Insights from molecular dynamics simulations examined in conjunction with results of X-ray studies

Molecular dynamics simulations have been carried out on all the jacalin-carbohydrate complexes of known structure, models of unliganded molecules derived from the complexes and also models of relevant complexes where X-ray structures are not available. Results of the simulations and the available crystal structures involving jacalin permit delineation of the relatively rigid and flexible regions of the molecule and the dynamical variability of the hydrogen bonds involved in stabilizing the structure. Local flexibility appears to be related to solvent accessibility. Hydrogen bonds involving side chains and water bridges involving buried water molecules appear to be important in the stabilization of loop structures. The lectin-carbohydrate interactions observed in crystal structures, the average parameters pertaining to them derived from simulations, energetic contribution of the stacking residue estimated from quantum mechanical calculations, and the scatter of the locations of carbohydrate and carbohydrate-binding residues are consistent with the known thermodynamic parameters of jacalin-carbohydrate interactions. The simulations, along with X-ray results, provide a fuller picture of carbohydrate binding by jacalin than provided by crystallographic analysis alone. The simulations confirm that in the unliganded structures water molecules tend to occupy the positions occupied by carbohydrate oxygens in the lectin-carbohydrate complexes. Population distributions in simulations of the free lectin, the ligands, and the complexes indicate a combination of conformational selection and induced fit. Proteins 2009; 77:760-777.

Soft matter lithium salt electrolytes: ion conduction and application to rechargeable batteries

Soft matter provides diverse opportunities for the development of electrolytes for all solid state lithium batteries. Here we review soft matter solid electrolytes for lithium batteriesthat are primarily obtained starting from liquid electrolytic systems. This concept of solid electrolyte synthesis from liquid is significantly different from prevalent approaches. The novelty of our approach is discussed in the light of various fundamental issues and in relation to its application to rechargeable lithium batteries.

Post-liquefaction undrained monotonic behaviour of sands: experiments and DEM simulations

This paper presents the results of a series of servo-controlled cyclic triaxial tests and numerical simulations using the three- dimensional discrete element method (DEM) on post-liquefaction undrained monotonic strength of granular materials. In a first test series,undrained monotonic tests were carried out after dissipating the excess pore water pressure developed during liquefaction. The influence of different parameters such as amplitude of axial strain,relative density and confining pressure prior to liquefaction on the post-liquefaction undrained response have been investigated.The results obtained highlight an insignificant influence of amplitude of axial strain, confining pressure and a significant influence of relative density on the post-liquefaction undrained monotonic stress-strain response.In the second series, undrained monotonic tests were carried out on similar triaxial samples without dissipating the excess pore water pressure developed during liquefaction. The results highlight that the amplitude of axial strain prior to liquefaction has a significant influence on the post-liquefaction undrained monotonic response.In addition,DEM simulations have been carried out on an assembly of spheres to simulate post-liquefaction behaviour.The simulations were very similar to the experiments with an objective to understand the behaviour of monotonic strength of liquefied samples from the grain scale. The numerical simulations using DEM have captured qualitatively all the features of the post-liquefaction undrained monotonic response in a manner similar to that of the experiments.In addition,a detailed study on the evolution of micromechanical parameters such as the average coordination number and induced anisotropic coefficients has been reported during the post-liquefaction undrained monotonic loading.

Precipitate growth with composition-dependent diffusivity: Comparison between theory and phase field simulations

We consider the growth of an isolated precipitate when the matrix diffusivity depends on the composition. We have simulated precipitate growth using the Cahn-Hilliard model, and find good agreement between our results and those from a sharp interface theory for systems with and without a dilatational misfit. With misfit, we report (and rationalize) an interesting difference between systems with a constant diffusivity and those with a variable diffusivity in the matrix.

Extending temporal simulations

Direct numerical simulations (DNS) of spatially growing turbulent shear layers may be performed as temporal simulations by solving the governing equations with some additional terms while imposing streamwise periodicity. These terms are functions of the means whose spatial growth is calculated easily and accurately from statistics of the temporal DNS. Equations for such simulations are derived.

Critical lattice size limit for synchronized chaotic state in one- and two-dimensional diffusively coupled map lattices

We consider diffusively coupled map lattices with P neighbors (where P is arbitrary) and study the stability of the synchronized state. We show that there exists a critical lattice size beyond which the synchronized state is unstable. This generalizes earlier results for nearest neighbor coupling. We confirm the analytical results by performing numerical simulations on coupled map lattices with logistic map at each node. The above analysis is also extended to two-dimensional P-neighbor diffusively coupled map lattices.

Generalized survival in step fluctuations

The properties of the generalized survival probability, that is, the probability of not crossing an arbitrary location R during relaxation, have been investigated experimentally (via scanning tunneling microscope observations) and numerically. The results confirm that the generalized survival probability decays exponentially with a time constant tau(s)(R). The distance dependence of the time constant is shown to be tau(s)(R)=tau(s0)exp[-R/w(T)], where w(2)(T) is the material-dependent mean-squared width of the step fluctuations. The result reveals the dependence on the physical parameters of the system inherent in the prior prediction of the time constant scaling with R/L-alpha, with L the system size and alpha the roughness exponent. The survival behavior is also analyzed using a contrasting concept, the generalized inside survival S-in(t,R), which involves fluctuations to an arbitrary location R further from the average. Numerical simulations of the inside survival probability also show an exponential time dependence, and the extracted time constant empirically shows (R/w)(lambda) behavior, with lambda varying over 0.6 to 0.8 as the sampling conditions are changed. The experimental data show similar behavior, and can be well fit with lambda=1.0 for T=300 K, and 0.5

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.

Melting of an Anchored Bilayer: Molecular Dynamics Simulations of the Structural Transition in (CnH2n+1NH3)(2)PbI4 (n=12, 14, 16, 18)

The thermally driven Structural phase transition in the organic-inorganic hybrid perovskite (CnH2n+1NH3)(2)PbI4 has been investigated using molecular dynamics (MD) simulations. This system consists of positively charged alkyl-amine chains anchored to a rigid negatively charged PbI4 sheet with the chains organized as bilayers with a herringbone arrangement. Atomistic simulations were performed using ail isothermal-isobaric ensemble over a wide temperature range from 65 to 665 K for different alkyl chain lengths, n = 12, 14, 16, and 18. The simulations are able to reproduce the essential Features of the experimental observations of this system, including the existence of a transition, the linear variation of the transition temperature with alkyl chain length, and the expansion of the bilayer thickness at the transition. By use of the distance fluctuation Criteria, it is Shown that the transition is associated With a Melting of the alkyl chains of the anchored bilayer. Ail analysis of the conformation of the alkyl chains shows increased disorder in the form of gauche defects above due melting transition. Simulations also show that the melting transition is characterized by the complete disappearance of all-trans alkyl chains in the anchored bilayer, in agreement with experimental observations. A conformationally disordered chain has a larger effective cross-sectional area, and above due transition a uniformly tilted arrangement of the anchored chains call no longer be Sustained. At the melt the angular distribution of the orientation of the chains are 110 longer uniform; the chains are splayed allowing for increased space for individual chains of the anchored bilayer. This is reflected in a sharp rise in the ratio of the mean head-to-head to tail-to-tail distance of the chains of the bilayer at the transition resulting in in expansion of the bilayer thickness. The present MD simulations provide a simple explanation as to how changes in conformation of individual alkyl-chains gives rise to the observed increase in the interlayer lattice spacing of (CnH2n+1NH3)(2)PbI4 at the melting transition.

Diffusion of Pure CH4 and Its Binary Mixture with CO2 in Faujasite NaY: A Combination of Neutron Scattering Experiments and Molecular Dynamics Simulations

The self-diffusion properties of pure CH4 and its binary mixture with CO2 within MY zeolite have been investigated by combining an experimental quasi-elastic neutron scattering (QENS) technique and classical Molecular dynamics simulations. The QENS measurements carried out at 200 K led to an unexpected self-diffusivity profile for Pure CH4 with the presence of a maximum for a loading of 32 CH4/unit cell, which was never observed before for the diffusion of apolar species in azeolite system With large windows. Molecular dynamics simulations were performed using two distinct microscopic models for representing the CH4/NaY interactions. Depending on the model, we are able to fairly reproduce either the magnitude or the profile of the self-diffusivity.Further analysis allowed LIS to provide some molecular insight into the diffusion mechanism in play. The QENS measurements report only a slight decrease of the self-diffusivity of CH4 in the presence of CO2 when the CO2 loading increases. Molecular dynamics simulations successfully capture this experimental trend and suggest a plausible microscopic diffusion mechanism in the case of this binary mixture.

Superplasticity in intermetallic NiAl nanowires via atomistic simulations

A novel superplastic deformation in an intermetallic B2-NiAl nanowire of cross-sectional dimensions of similar to 20 angstrom with failure strain as high as similar to 700% at 700 K temperature is reported. The minimum temperature under which the superplasticity has been observed is around 0.36 T-m, which is much lower than 0.5 T-m (T-m = melting temperature i.e. 1911 K for bulk B2-NiAl). Superplasticity is observed due to transformation from crystalline phase to amorphous phase after yielding of the nanowire. (C) 2010 Elsevier B.V. All rights reserved.

Stochastic simulations of the origins and implications of long-tailed distributions in gene expression

Gene expression noise results in protein number distributions ranging from long-tailed to Gaussian. We show how long-tailed distributions arise from a stochastic model of the constituent chemical reactions and suggest that, in conjunction with cooperative switches, they lead to more sensitive selection of a subpopulation of cells with high protein number than is possible with Gaussian distributions. Single-cell-tracking experiments are presented to validate some of the assumptions of the stochastic simulations. We also examine the effect of DNA looping on the shape of protein distributions. We further show that when switches are incorporated in the regulation of a gene via a feedback loop, the distributions can become bimodal. This might explain the bimodal distribution of certain morphogens during early embryogenesis.

Ion Transport in a Polymer-Plastic Solid Soft Matter Electrolyte in the Light of Solvent Dynamics and Ion Association

Ion transport in a recently demonstrated promising soft matter solid plastic-polymer electrolyte is discussed here in the context of solvent dynamics and ion association. The plastic-polymer composite electrolytes display liquid-like ionic conductivity in the solid state,compliable mechanical strength (similar to 1 MPa), and wide electrochemical voltage stability (>= 5 V). Polyacrylonitrile (PAN) dispersed in lithium perchlorate (LiClO4)-succinonitrile (SN) was chosen as the model system for the study (abbreviated LiClO4-SN:PAN). Systematic observation of various mid-infrared isomer and ion association bands as a function of temperature and polyme concentration shows an effective increase in trans conformer concentration along with free Li+ ion concentration. This strongly supports the view that enhancement in LiClO4-SN:PAN ionic conductivity over the neat plastic electrolyte (LiClO4-SN) is due to both increase in charge mobility and concentration. The ionic conductivity and infrared spectroscopy studies are supported by Brillouin light scattering. For the LiClO4-SN:PAN composites, a peak at 17 GHz was observed in addition to the normal trans-gauche isomerism (as in neat SN) at 12 GHz. The fast process is attributed to increased dynamics of those SN molecules whose energy barrier of transition from gauche to trans has reduced under influences induced by the changes in temperature and polymer concentration. The observations from ionic conductivity, spectroscopy, and light scattering studies were further supplemented by temperature dependent nuclear magnetic resonance H-1 and Li-7 line width measurements.

Comments on "Full-sphere simulations of circulation-dominated solar dynamo: Exploring the parity issue" - Reply

We respond to Dikpati et al.'s criticism of our recent solar dynamo model. A different treatment of the magnetic buoyancy is the most probable reason for their different results.

Elastic behaviour of matter under very high pressures - Uniform compression

In order to study the elastic behaviour of matter when subjected to very large pressures, such as occur for example in the interior of the earth, and to provide an explanation for phenomena like earthquakes, it is essential to be able to calculate the values of the elastic constants of a substance under a state of large initial stress in terms of the elastic constants of a natural or stress-free state. An attempt has been made in this paper to derive expressions for these quantities for a substance of cubic symmetry on the basis of non-linear theory of elasticity and including up to cubic powers of the strain components in the strain energy function. A simple method of deriving them directly from the energy function itself has been indicated for any general case and the same has been applied to the case of hydrostatic compression. The notion of an effective elastic energy-the energy require to effect an infinitesimal deformation over a state of finite strain-has been introduced, the coefficients in this expression being the effective elastic constants. A separation of this effective energy function into normal co-ordinates has been given for the particular case of cubic symmetry and it has been pointed out, that when any of such coefficients in this normal form becomes negative, elastic instability will set in, with associated release of energy.