Complexation of siRNA with Dendrimer: A Molecular Modeling Approach

This paper reports the structural behavior and thermodynamics of the complexation of siRNA with poly(amidoamine) (PAMAM) dendrimers of generation 3 (G3) and 4 (G4) through fully atomistic molecular dynamics (MD) simulations accompanied by free energy calculations and inherent structure determination. We have also done simulation with one siRNA and two dendrimers (2 x G3 or 2xG4) to get the microscopic picture of various binding modes. Our simulation results reveal the formation of stable siRNA-dendrimer complex over nanosecond time scale. With the increase in dendrimcr generation, the charge ratio increases and hence the binding energy between siRNA and dendrimer also increases in accordance with available experimental measurements. Calculated radial distribution functions of amines groups of various subgenerations in a given generation of dendrimer and phosphate in backbone of siRNA reveals that one dendrimer of generation 4 shows better binding with siRNA almost wrapping the dendrimer when compared to the binding with lower generation dendrimer like G3. In contrast, two dendrimers of generation 4 show binding without siRNA wrapping the den-rimer because of repulsion between two dendrimers. The counterion distribution around the complex and the water molecules in the hydration shell of siRNA give microscopic picture of the binding dynamics. We see a clear correlation between water. counterions motions and the complexation i.e. the water molecules and counterions which condensed around siRNA are moved away from the siRNA backbone when dendrimer start binding to the siRNA back hone. As siRNA wraps/bind to the dendrimer counterions originally condensed onto siRNA (Na-1) and dendrimer (Cl-) get released. We give a quantitative estimate of the entropy of counterions and show that there is gain in entropy due to counterions release during the complexation. Furthermore, the free energy of complexation of IG3 and IG4 at two different salt concentrations shows that increase in salt concentration leads to the weakening of the binding affinity of siRNA and dendrimer.

Steady incompressible flow of cohesionless granular materials through a wedge-shaped hopper: Frictional 14kinetic solution to the smooth wall, radial gravity problem

Hybrid frictional-kinetic equations are used to predict the velocity, grain temperature, and stress fields in hoppers. A suitable choice of dimensionless variables permits the pseudo-thermal energy balance to be decoupled from the momentum balance. These balances contain a small parameter, which is analogous to a reciprocal Reynolds number. Hence an approximate semi-analytical solution is constructed using perturbation methods. The energy balance is solved using the method of matched asymptotic expansions. The effect of heat conduction is confined to a very thin boundary layer near the exit, where it causes a marginal change in the temperature. Outside this layer, the temperature T increases rapidly as the radial coordinate r decreases. In particular, the conduction-free energy balance yields an asymptotic solution, valid for small values of r, of the form T proportional r-4. There is a corresponding increase in the kinetic stresses, which attain their maximum values at the hopper exit. The momentum balance is solved by a regular perturbation method. The contribution of the kinetic stresses is important only in a small region near the exit, where the frictional stresses tend to zero. Therefore, the discharge rate is only about 2.3% lower than the frictional value, for typical parameter values. As in the frictional case, the discharge rate for deep hoppers is found to be independent of the head of material.

Entropy and dynamics of water in hydration layers of a bilayer

We compute the entropy and transport properties of water in the hydration layer of dipalmitoylphosphatidylcholine bilayer by using a recently developed theoretical scheme two-phase thermodynamic model, termed as 2PT method; S.-T. Lin et al., J. Chem. Phys. 119, 11792 (2003)] based on the translational and rotational velocity autocorrelation functions and their power spectra. The weights of translational and rotational power spectra shift from higher to lower frequency as one goes from the bilayer interface to the bulk. Water molecules near the bilayer head groups have substantially lower entropy (48.36 J/mol/K) than water molecules in the intermediate region (51.36 J/mol/K), which have again lower entropy than the molecules (60.52 J/mol/K) in bulk. Thus, the entropic contribution to the free energy change (T Delta S) of transferring an interface water molecule to the bulk is 3.65 kJ/mol and of transferring intermediate water to the bulk is 2.75 kJ/mol at 300 K, which is to be compared with 6.03 kJ/mol for melting of ice at 273 K. The translational diffusion of water in the vicinity of the head groups is found to be in a subdiffusive regime and the rotational diffusion constant increases going away from the interface. This behavior is supported by the slower reorientational relaxation of the dipole vector and OH bond vector of interfacial water. The ratio of reorientational relaxation time for Legendre polynomials of order 1 and 2 is approximately 2 for interface, intermediate, and bulk water, indicating the presence of jump dynamics in these water molecules. (C) 2010 American Institute of Physics. doi:10.1063/1.3494115]

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.

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.

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.

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.