Thermodynamics of ligand (substrate end product) binding to endoxylanase from Chainia sp. (NCL-82-5-1): isothermal calorimetry and fluorescence titration studies

The binding of xylo-oligosaccharides to Chainia endoxylanase resulted in a decrease in fluorescence intensity of the enzyme with the formation of 1:1 complex. Equilibrium and thermodynamic parameters of ligand binding were determined by fluorescence titrations and titration calorimetry. The affinity of xylanase for the oligosaccharides increases in the order X-2 < X-3 < X-4 less than or equal to X-5. Contributions from the enthalpy towards the free energy change decreased with increasing chain length from X-2 to X-4, whereas an increase in entropy was observed, the change in enthalpy and entropy of binding being compensatory. The entropically driven binding process suggested that hydrophobic interactions as well as hydrogen bonds play a predominant role in ligand binding.

Multiscale modeling of lamellar mesophases

The mesoscale simulation of a lamellar mesophase based on a free energy functional is examined with the objective of determining the relationship between the parameters in the model and molecular parameters. Attention is restricted to a symmetric lamellar phase with equal volumes of hydrophilic and hydrophobic components. Apart from the lamellar spacing, there are two parameters in the free energy functional. One of the parameters, r, determines the sharpness of the interface, and it is shown how this parameter can be obtained from the interface profile in a molecular simulation. The other parameter, A, provides an energy scale. Analytical expressions are derived to relate these parameters to r and A to the bending and compression moduli and the permeation constant in the macroscopic equation to the Onsager coefficient in the concentration diffusion equation. The linear hydrodynamic response predicted by the theory is verified by carrying out a mesoscale simulation using the lattice-Boltzmann technique and verifying that the analytical predictions are in agreement with simulation results. A macroscale model based on the layer thickness field and the layer normal field is proposed, and the relationship between the parameters in the macroscale model from the parameters in the mesoscale free energy functional is obtained.

The nature of non-equilibrium solidification of undercooled Agsingle bondCu alloys in contact with primary copper dendrites

Liquids of silver-copper alloys with near eutectic compositions embedded in a copper matrix were undercooled. The structural and microstructural investigations of these alloys solidified from undercooled temperature indicated the absence of both the eutectic reaction and diffusionless transformation below the equal free energy curve (T0). Instead the liquid maintained local equilibrium with the copper dendrites continuously until it intersected the extended solidus of the silver rich solid solution.

Molecular conformation of N,N -diarylthioureas: An assessment by 1H NMR and infrared spectroscopy

Several N,N -dipyridyl- and N-phenyl-N -pyridyl-thioureas were examined in different solvents at various temperatures by 1H NMR in order to study their conformational properties. The influence of concentration and the methyl substituent in the pyridine ring on the chemical shifts of the NH and pyridine groups was investigated. The observed chemical shifts are analysed in terms of the conformational properties of the molecules. Free energy barriers to the internal rotation about the C N bonds have been determined. Infrared spectra have been measured to supplement the NMR studies. Intramolecular hydrogen bonding played a major role in the preferred conformation of pyridylthioureas. The data further revealed an interesting dynamic exchange phenomenon occurring in symmetric N,N -dipyridylthioureas between two intramolecularly hydrogen bonded conformers.

Quantum-Ohmic Resistance Fluctuation In Disordered Conductors - An Invariant Imbedding Approach

It is now well known that in extreme quantum limit, dominated by the elastic impurity scattering and the concomitant quantum interference, the zero-temperature d.c. resistance of a strictly one-dimensional disordered system is non-additive and non-self-averaging. While these statistical fluctuations may persist in the case of a physically thin wire, they are implicitly and questionably ignored in higher dimensions. In this work, we have re-examined this question. Following an invariant imbedding formulation, we first derive a stochastic differential equation for the complex amplitude reflection coefficient and hence obtain a Fokker-Planck equation for the full probability distribution of resistance for a one-dimensional continuum with a Gaussian white-noise random potential. We then employ the Migdal-Kadanoff type bond moving procedure and derive the d-dimensional generalization of the above probability distribution, or rather the associated cumulant function –‘the free energy’. For d=3, our analysis shows that the dispersion dominates the mobilitly edge phenomena in that (i) a one-parameter B-function depending on the mean conductance only does not exist, (ii) an approximate treatment gives a diffusion-correction involving the second cumulant. It is, however, not clear whether the fluctuations can render the transition at the mobility edge ‘first-order’. We also report some analytical results for the case of the one dimensional system in the presence of a finite electric fiekl. We find a cross-over from the exponential to the power-low length dependence of resistance as the field increases from zero. Also, the distribution of resistance saturates asymptotically to a poissonian form. Most of our analytical results are supported by the recent numerical simulation work reported by some authors.

An absolute method for the determination of surface tension of liquids using pendent drop profiles

Drop formation at conical tips which is of relevance to metallurgists is investigated based on the principle of minimization of free energy using the variational approach. The dimensionless governing equations for drop profiles are computer solved using the fourth order Runge-Kutta method. For different cone angles, the theoretical plots of XT and ZT vs their ratio, are statistically analyzed, where XT and ZT are the dimensionless x and z coordinates of the drop profile at a plane at the conical tip, perpendicular to the axis of symmetry. Based on the mathematical description of these curves, an absolute method has been proposed for the determination of surface tension of liquids, which is shown to be preferable in comparison with the earlier pendent-drop profile methods.

Allostery in tRNA Synthetases Elucidated from MD Simulations and Protein Structure Networks

tRNA synthetases (aaRS) are enzymes crucial in the translation of genetic code. The enzyme accylates the acceptor stem of tRNA by the congnate amino acid bound at the active site, when the anti-codon is recognized by the anti-codon site of aaRS. In a typical aaRS, the distance between the anti-codon region and the amino accylation site is approximately 70 Å. We have investigated this allosteric phenomenon at molecular level by MD simulations followed by the analysis of protein structure networks (PSN) of non-covalent interactions. Specifically, we have generated conformational ensembles by performing MD simulations on different liganded states of methionyl tRNA synthetase (MetRS) from Escherichia coli and tryptophenyl tRNA synthetase (TrpRS) from Human. The correlated residues during the MD simulations are identified by cross correlation maps. We have identified the amino acids connecting the correlated residues by the shortest path between the two selected members of the PSN. The frequencies of paths have been evaluated from the MD snapshots[1]. The conformational populations in different liganded states of the protein have been beautifully captured in terms of network parameters such as hubs, cliques and communities[2]. These parameters have been associated with the rigidity and plasticity of the protein conformations and can be associated with free energy landscape. A comparison of allosteric communication in MetRS and TrpRS [3] elucidated in this study highlights diverse means adopted by different enzymes to perform a similar function. The computational method described for these two enzymes can be applied to the investigation of allostery in other systems.

Photo-induced double-strand DNA and site-specific protein cleavage activity of L-histidine (mu-oxo)diiron(III) complexes of heterocyclic bases

Three oxo-bridged diiron(III) complexes of L-histidine and heterocyclic bases [Fe-2(mu-O)(L-his)(2)(B)(2)](ClO4)(2) (1-3), where B is 2,2'-bipyridine (bpy),1,10-phenanthroline (phen), dipyrido[3,2-d:2',3'-f]quinoxaline (dpq), were prepared and characterized. The bpy complex 1 was structurally characterized by X-ray crystallography. The molecular structure showed a {Fe-2(mu-O)} core in which iron(III) in a FeN4O2 coordination is bound to tridentate monoanionic L-histidine and bidentate bpy ligands. The Fe center dot center dot center dot Fe distance is similar to 3.5 angstrom. The Fe-O-Fe unit is essentially linear, giving a bond angle of similar to 172 degrees. The complexes showed irreversible cyclic voltammetric cathodic response near -0.1 V vs. SCE in H2O-0.1 M KCl. The binuclear units displayed antiferromagnetic interaction between two high-spin (S = 5/2) iron(III) centers giving a -J value of -110 cm(-1). The complexes showed good DNA binding propensity giving a binding constant value of similar to 10(5) M-1. Isothermal titration calorimetric data indicated single binding mode to the DNA. The binding was found to be driven by negative free energy change and enthalpy. The dpq complex 3 showed oxidative double-strand DNA cleavage on exposure to UV-A and visible light. The phen complex 2 displayed single-strand photocleavage of DNA. The DNA double-strand breaks were rationalized from theoretical molecular docking calculations. Mechanistic investigations showed formation of hydroxyl radicals as the reactive species through photodecarboxylation of the L-histidine ligand. The complexes exhibited good binding propensity to bovine serum albumin (BSA) protein in Tris-HCl/NaCl buffer medium. The dpq complex 3 showed UV-A light-induced site-specific oxidative BSA cleavage forming fragments of similar to 45 kDa and similar to 20 kDa molecular weights via SOH pathway.

Mesoscale description of an asymmetric lamellar phase

The relationship between the parameters in a description based on a mesoscale free energy functional for the concentration field and the macroscopic properties, such as the bending and compression moduli and the permeation constant, are examined for an asymmetric lamellar phase where the mass fractions of the hydrophobic and hydrophilic parts are not equal. The difference in the mass fractions is incorporated using a cubic term in the free energy functional, in addition to the usual quadratic and quartic terms in the Landau–Ginsburg formulation. The relationship between the coefficient of the cubic term and the difference in the mass fractions of the hydrophilic and hydrophobic parts is obtained. For a lamellar phase, it is important to ensure that the surface tension is zero due to symmetry considerations. The relationship between the parameters in the free energy functional for zero surface tension is derived. When the interface between the hydrophilic and hydrophobic parts is diffuse, it is found that the bending and compression moduli, scaled by the parameters in the free energy functional, do increase as the asymmetry in the bilayer increases. When the interface between the hydrophilic and hydrophobic parts is sharp, the scaled bending and compression moduli show no dependence on the asymmetry in the bilayer. The ratio of the permeation constant in between the water and bilayer in a molecular description and the Onsager coefficient in the mesoscale description is O(1) for both sharp and diffuse interfaces and it increases as the difference in the mass fractions is increased.

Self-Organization of n-Alkane Chains in Water: Length Dependent Crossover from Helix and Toroid to Molten Globule

We demonstrate a chain length dependent crossover in the structural properties of linear hydrocarbon (n-alkane) chains using detailed atomistic simulations in explicit water. We identify a number of exotic structures of the polymer chain through energy minimization of representative snapshots collected from molecular dynamics trajectory. While the collapsed state is ring-like (circular) for small chains (CnH2n+2; n <= 20) and spherical for very long ones (n = 100), we find the emergence of ordered helical structures at intermediate lengths (n similar to 40). We find different types of disordered helices and toroid-like structures at n = 60. We also report a sharp transition in the stability of the collapsed state as a function of the chain length through relevant free energy calculations. While the collapsed state is only marginally metastable for C20H42, a clear bistable free energy surface emerges only when the chain is about 30 monomers long. For n = 30, the polymer exhibits an intermittent oscillation between the collapsed and the coil structures, characteristic of two stable states separated by a small barrier.

Jump Reorientation of Water Molecules Confined in Narrow Carbon Nanotubes

We used molecular dynamics (MD) simulations to study the reorientational dynamics of water molecules confined inside narrow carbon nanotubes immersed in a bath of water. Our simulations show that the confined water molecules exhibit bistability in their reorientational relaxation, which proceeds by angular jumps between the two stable states. The angular jump of a water molecule in the bulk involves the breaking of a hydrogen bond with one of its neighbors and the formation of a hydrogen bond with a different neighbor. In contrast, the angular jump of a confined water molecule corresponds to an interchange of the two hydrogen atoms that can form a hydrogen bond with the same neighbor. The free energy barrier between these two states is a few k(B)T. The analytic solution of a simplified two-state jump model that qualitatively explains the reorientational behavior observed in simulations is also presented.

Role of conformational dynamics in kinetics of an enzymatic cycle in a nonequilibrium steady state

Enzyme is a dynamic entity with diverse time scales, ranging from picoseconds to seconds or even longer. Here we develop a rate theory for enzyme catalysis that includes conformational dynamics as cycling on a two-dimensional (2D) reaction free energy surface involving an intrinsic reaction coordinate (X) and an enzyme conformational coordinate (Q). The validity of Michaelis-Menten (MM) equation, i.e., substrate concentration dependence of enzymatic velocity, is examined under a nonequilibrium steady state. Under certain conditions, the classic MM equation holds but with generalized microscopic interpretations of kinetic parameters. However, under other conditions, our rate theory predicts either positive (sigmoidal-like) or negative (biphasic-like) kinetic cooperativity due to the modified effective 2D reaction pathway on X-Q surface, which can explain non-MM dependence previously observed on many monomeric enzymes that involve slow or hysteretic conformational transitions. Furthermore, we find that a slow conformational relaxation during product release could retain the enzyme in a favorable configuration, such that enzymatic turnover is dynamically accelerated at high substrate concentrations. The effect of such conformation retainment in a nonequilibrium steady state is evaluated.

A semi-empirical theory of molten salts

A semi-empirical model is presented for describing the interionic interactions in molten salts using the experimentally available structure data. An extension of Bertaut's method of non-overlapping charges is used to estimate the electrostatic interaction energy in ionic melts. It is shown, in agreement with earlier computer simulation studies, that this energy increases when an ionic salt melts. The repulsion between ions is described using a compressible ion theory which uses structure-independent parameters. The van der Waals interactions and the thermal free energy are also included in the total energy, which is minimised with respect to isostructural volume variations to calculate the equilibrium density. Detailed results are presented for three molten systems, NaCl, CaCl2 and ZnCl2, and are shown to be in satisfactory agreement with experiments. With reliable structural data now being reported for several other molten salts, the present study gains relevance.

Adsopption isotherms-A critique

The possibility or the impossibility of separating the particle and the electrode interactions is discussed in a wider context of the effects due to any two interaction potentials on the equation of state. The involved nature of the pressure dependence on two individually definable forces is illustrated through the Percus Yevick results for the adhesive hard spheres. An alternative form of the adsorption isotherm is given to bring home the intimate relationship between the actual equation of state and the free energy of adsorption. Thermodynamic consequences of congruence with respect to E (or q) as reflected through the linear plots of q (or E) vs. θ are well known. Mathematical consequences of simultaneous congruence have been pointed out recently. In this paper, the physical nature of congruence hypothesis is revealed. In passing "the pseudo-congruence" is also discussed. It is emphasised that the problem is no less ambiguous with regard to modelling the particle/particle interaction. The ad hoc nature of our dependence of the available equations of state is emphasised through a discussion on the HFL theory. Finally, a heuristic method for modelling ΔG mathematically-incorporating its behaviour at saturation coverages-is advanced. The more interesting aspects of this approach, which generalises almost all isotherms hitherto known, are sketched.

Phase diagram of vortex matter in layered superconductors with tilted columnar pinning centers

We study the vortex matter phase diagram of a layered superconductor in the presence of columnar pinning defects, tilted with respect to the normal to the layers. We use numerical minimization of the free energy written as a functional of the time-averaged vortex density of the Ramakrishnan-Yussouff form, supplemented by the appropriate pinning potential. We study the case where the pin density is smaller than the areal vortex density. At lower pin concentrations, we find, for temperatures of the order of the melting temperature of the unpinned lattice, a Bose glass type phase which at lower temperatures converts, via a first-order transition, to a Bragg glass, while, at higher temperatures, it crosses over to an interstitial liquid. At somewhat higher concentrations, no transition to a Bragg glass is found even at the lowest temperatures studied. While qualitatively the behavior we find is similar to that obtained using the same procedures for columnar pins normal to the layers, there are important and observable quantitative differences, which we discuss.