High-quality annotation of promoter regions for 913 bacterial genomes

Motivation: The number of bacterial genomes being sequenced is increasing very rapidly and hence, it is crucial to have procedures for rapid and reliable annotation of their functional elements such as promoter regions, which control the expression of each gene or each transcription unit of the genome. The present work addresses this requirement and presents a generic method applicable across organisms. Results: Relative stability of the DNA double helical sequences has been used to discriminate promoter regions from non-promoter regions. Based on the difference in stability between neighboring regions, an algorithm has been implemented to predict promoter regions on a large scale over 913 microbial genome sequences. The average free energy values for the promoter regions as well as their downstream regions are found to differ, depending on their GC content. Threshold values to identify promoter regions have been derived using sequences flanking a subset of translation start sites from all microbial genomes and then used to predict promoters over the complete genome sequences. An average recall value of 72% (which indicates the percentage of protein and RNA coding genes with predicted promoter regions assigned to them) and precision of 56% is achieved over the 913 microbial genome dataset.

Phase relationships in the system Ni-W-O and thermodynamic properties of NiWO,

The phase diagram of the Ni-W-O system at 1200 K was established by metallographic and X-ray identification of the phases present after equilibration at controlled oxygen potentials. The oxygen partial pressures over the samples were fixed by metered streams of CO+CO2 gas mixtures. There was only one ternary oxide, nickel tungstate (NiWO4), in the Ni-W-O system at a total pressure of 1 atm, and this compound decomposed to a mixture of Ni+WO2.72 on lowering the oxygen potential. The Gibbs' free energy of formation of NiWO4 was determined from the measurement of the e.m.f. of the solid oxide galvanic cell, Pt, Ni+NiWO4+WO2.72/CaO-ZrO2/Ni+NiO, Pt and thermodynamic properties of tungsten and nickel oxides available in the literature. For the reaction, NiO(s)+WO3(s)rarrNiWO4(s) DeltaG°=–10500–0.708 T (±250) cal mol–1.

A nanosecond molecular dynamics study of antiparallel d(G)(7) quadruplex structures: Effect of the coordinated cations

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G(7)) quadruplex structures with different coordinated ions, namely Na+ and K+ ion, water and Na+ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of similar to1.5 Angstrom from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K+ ions are quite immobile. MD studies at 400 K indicate that K+ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K+ ion coordinated structure is more favourable than that for the Na+ ion coordinated antiparallel quadruplex structure.

Padé approach to potential transients: Part 1. Electron transfer without and with coupling to first-order homogeneous reactions at planar electrodes

Potential transients are obtained by using “Padé approximants” (an accurate approximation procedure valid globally — not just perturbatively) for all amplitudes of concentration polarization and current densities. This is done for several mechanistic schemes under constant current conditions. We invert the non-linear current-potential relationship in the form (using the Lagrange or the Ramanujan method) of power series appropriate to the two extremes, namely near reversible and near irreversible. Transforming both into the Pad́e expressions, we construct the potential-time profile by retaining whichever is the more accurate of the two. The effectiveness of this method is demonstrated through illustrations which include couplings of homogeneous chemical reactions to the electron-transfer step.

On the evaluation of stability of rare earth oxides as face coats for investment casting of titanium

Attempts have been made to evaluate the thermal stability of rare earth oxide face coats against liquid titanium. Determination of microhardness profiles and concentration profiles of oxygen and metallic constituents of oxide in investment cast titanium rods has allowed grActation of relative stability of rare earth oxides. The relative stability of evaluated oxides in the order of increasing stability follows the sequence CeO2 — ZrO2 — Gd2O3 — didymium oxide — Sm2O3 —Nd2O3 — Y2O3. The grading does not follow the free energy data of the formation of these oxides. A better correlation with the experimental observations is obtained when the solubility of the metallic species in titanium is also taken into consideration.

Improved Inverted V-Notch or Chimney Weir

This paper is devoted to the improvement of the measuring range of inverted V-notch (IVN) weir, a practical linear sharp-crested weir, designed earlier by the writers. The range of linearity of IVN can be considerably enhanced (by more than 200%) by the addition of a retangular weir of width 0.265W (W = half crest width) at a depth of 0.735d (d = altitude of IVN), above the crest of the weir, which is equivalent to providing at this depth two vertical straight lines to the IVN, resulting in a chimney-shaped profile; hence, the modified weir is named chimney weir. The design parameters of the weir, that is, the linearity range, base flow depth, and datum constant, which fixes the reference plane of the weir, are estimated by solving the nonlinear programming problem using a numerical optimization procedure. For flows through this weir above a depth of 0.22d, the discharges are proportional to the depth of flow measured above a reference plane situated at 0.08d above the weir crest for all heads in the range 0.22d <= h <= 2.43d, within a maximum percentage deviation of ±1.5 from the theoretical discharge. A significant result of the analysis is that the same linear head-discharge relationship governing the flow through the IVN is also valid for the extended chimney weir. Experiments with three different chimney weirs show excellent agreement with the theory by giving a constant average coefficient of discharge for each weir.

A Numerical-Simulation Of The Gneiss-Charnockite Transformations In Southern India

This paper presents a numerical simulation of the well-documented, fluid-controlled Kabbal and Ponmudi type gneiss-chamockite transformations in southern India using a free energy minimization method. The computations have considered all the major solid phases and important fluid species in the rock - C-O-H and rock - C-O-H-N systems. Appropriate activity-composition relations for the solid solutions and equations of state for the fluids have been included in order to evaluate the mineral-fluid equilibria attending the incipient chamockite development in the gneisses. The C-O-H fluid speciation pattern in both the Kabbal and Ponmudi type systems indicates that CO2 and H2O make up the bulk of the fluid phase with CO, CH4, H-2 and O2 as minor constituents. In the graphite-buffered Ponmudi-system, the abundance of CO, CH4 and H-2 is orders of magnitude higher than that in the graphite-free Kabbal system. Simulation with C-O-H-N fluids of varying composition demonstrates the complementary role of CO2 and N2 as rather inert dilutants of H2O in the fluid phase. The simulation, carried out on available whole-rock data, has demonstrated the dependence of the transformation X(H2O) on P,T, and phase and chemical composition of the precursor gneiss.

Direct tests of the entropic model of quasicrystals

Assuming an entropic origin for phason elasticity in quasicrystals, we derive predictions for the temperature dependence of grain-boundary structure and free energy, the nature of the elastic instability in these systems, and the behavior of sound damping near the instability. We believe that these will provide decisive tests of the entropic model for quasicrystals.

Mechanical alloying—a novel synthesis route for amorphous phases

Mechanical alloying (MA) pioneered by Benjamin is a technique for the extension of solid solubility in systems where the equilibrium solid solubility is limited. This technique has, in recent years, emerged as a novel alternate route for rapid solidification processing (RSP) for the production of metastable crystalline, quasicrystalline, amorphous phases and nanocrystalline materials. The glass-forming composition range (GFR), in general, is found to be much wider in case of MA in comparison with RSP. The amorphous powders produced by MA can be compacted to bulk shapes and sizes and can be used as precursors to obtain high strength materials. This paper reports the work done on solid state amorphization by MA in Ti-Ni-Cu and Al-Ti systems where a wide GFR has been obtained. Al-Ti is a classic case where no glass formation has been observed by RSP, while a GFR of 25–90 at.% Ti has been obtained in this system, thus demonstrating the superiority of MA over RSP. The free energy calculations made to explain GFR are also presented.

Nonlinear hydrodynamics of a hard-sphere fluid near the glass transition

We conduct a numerical study of the dynamic behavior of a dense hard-sphere fluid by deriving and integrating a set of Langevin equations. The statics of the system is described by a free-energy functional of the Ramakrishnan-Yussouff form. We find that the system exhibits glassy behavior as evidenced through a stretched exponential decay and a two-stage relaxation of the density correlation function. The characteristic times grow with increasing density according to the Vogel-Fulcher law. The wave-number dependence of the kinetics is extensively explored. The connection of our results with experiment, mode-coupling theory, and molecular-dynamics results is discussed.

Dynamic coupling between the LID and NMP domain motions in the catalytic conversion of ATP and AMP to ADP by adenylate kinase

The catalytic conversion of adenosine triphosphate (ATP) and adenosine monophosphate (AMP) to adenosine diphosphate (ADP) by adenylate kinase (ADK) involves large amplitude, ligand induced domain motions, involving the opening and the closing of ATP binding domain (LID) and AMP binding domain (NMP) domains, during the repeated catalytic cycle. We discover and analyze an interesting dynamical coupling between the motion of the two domains during the opening, using large scale atomistic molecular dynamics trajectory analysis, covariance analysis, and multidimensional free energy calculations with explicit water. Initially, the LID domain must open by a certain amount before the NMP domain can begin to open. Dynamical correlation map shows interesting cross-peak between LID and NMP domain which suggests the presence of correlated motion between them. This is also reflected in our calculated two-dimensional free energy surface contour diagram which has an interesting elliptic shape, revealing a strong correlation between the opening of the LID domain and that of the NMP domain. Our free energy surface of the LID domain motion is rugged due to interaction with water and the signature of ruggedness is evident in the observed root mean square deviation variation and its fluctuation time correlation functions. We develop a correlated dynamical disorder-type theoretical model to explain the observed dynamic coupling between the motion of the two domains in ADK. Our model correctly reproduces several features of the cross-correlation observed in simulations. (C) 2011 American Institute of Physics. doi:10.1063/1.3516588]

Phenomenological Ginzburg-Landau-like theory for superconductivity in the cuprates

We propose and develop here a phenomenological Ginzburg-Landau-like theory of cuprate high-temperature superconductivity. The free energy of a cuprate superconductor is expressed as a functional F of the complex spin-singlet pair amplitude psi(ij) equivalent to psi(m) = Delta(m) exp(i phi(m)), where i and j are nearest-neighbor sites of the square planar Cu lattice in which the superconductivity is believed to primarily reside, and m labels the site located at the center of the bond between i and j. The system is modeled as a weakly coupled stack of such planes. We hypothesize a simple form FDelta, phi] = Sigma(m)A Delta(2)(m) + (B/2)Delta(4)(m)] + C Sigma(< mn >) Delta(m) Delta(n) cos(phi(m) - phi(n)) for the functional, where m and n are nearest-neighbor sites on the bond-center lattice. This form is analogous to the original continuum Ginzburg-Landau free-energy functional; the coefficients A, B, and C are determined from comparison with experiments. A combination of analytic approximations, numerical minimization, and Monte Carlo simulations is used to work out a number of consequences of the proposed functional for specific choices of A, B, and C as functions of hole density x and temperature T. There can be a rapid crossover of from small to large values as A changes sign from positive to negative on lowering T; this crossover temperature T-ms(x) is identified with the observed pseudogap temperature T*(x). The thermodynamic superconducting phase-coherence transition occurs at a lower temperature T-c(x), and describes superconductivity with d-wave symmetry for positive C. The calculated T-c(x) curve has the observed parabolic shape. The results for the superfluid density rho(s)(x, T), the local gap magnitude , the specific heat C-v(x, T) (with and without a magnetic field), as well as vortex properties, all obtained using the proposed functional, are compared successfully with experiments. We also obtain the electron spectral density as influenced by the coupling between the electrons and the correlation function of the pair amplitude calculated from the functional, and compare the results successfully with the electronic spectrum measured through angle resolved photoemission spectroscopy (ARPES). For the specific heat, vortex structure, and electron spectral density, only some of the final results are reported here; the details are presented in subsequent papers.

Thermodynamics and phase equilibria in the Ga---Sb system

The optimum values of the solution parameters of a multiparameter integral free-energy function have been determined using experimental data from the Ga-Sb system. The equation is represented as DELTAG(xs) = x(1 - x)[(1 - x)(a1 + a2T + a3T ln T) + x(a4 + a5T + a6T ln T) + x(1 - x)(a7 + a8T + a9xT)].The integral and the corresponding partial form of the free energy function have been found to be of use when interpreting the high temperature thermodynamic data, atomic interactions and phase equilibria in the Ga-Sb system.

Glass Transition in the Density Functional Theory of Freezing

A mean-field description of the glass transition in the hard-sphere system is obtained by numerically locating "glassy" minima of a model free-energy functional. These minima, characterized by inhomogeneous but aperiodic density distributions, appear as the average density is increased above the value at which equilibrium crystallization takes place. Investigations of the density distribution and local bond-orientational order at these minima yield results similar to those obtained from simulations.

Barrier crossing in one and three dimensions by a long chain

We consider the Kramers problem for a long chain polymer trapped in a biased double-well potential. Initially the polymer is in the less stable well and it can escape from this well to the other well by the motion of its N beads across the barrier to attain the configuration having lower free energy. In one dimension we simulate the crossing and show that the results are in agreement with the kink mechanism suggested earlier. In three dimensions, it has not been possible to get an analytical `kink solution' for an arbitrary potential; however, one can assume the form of the solution of the nonlinear equation as a kink solution and then find a double-well potential in three dimensions. To verify the kink mechanism, simulations of the dynamics of a discrete Rouse polymer model in a double well in three dimensions are carried out. We find that the time of crossing is proportional to the chain length, which is in agreement with the results for the kink mechanism. The shape of the kink solution is also in agreement with the analytical solution in both one and three dimensions.