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.

Effect of coordinated ions on structure and flexiblity of parallel G-quandruplexes: a molecular dynamics study

Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G7) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the 06 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 A from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na+ counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, 06 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures

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 new approach to the kinetics of size reduction in ball milling

The rate of breakage of feed in ball milling is usually represented in the form of a first-order rate equation. The equation was developed by treating a simple batch test mill as a well mixed reactor. Several case of deviation from the rule have been reported in the literature. This is attributed to the fact that accumulated fines interfere with the feed material and breaking events are masked by these fines. In the present paper, a new rate equation is proposed which takes into account the retarding effect of fines during milling. For this purpose the analogy of diffusion of ions through permeable membranes is adopted, with suitable modifications. The validity of the model is cross checked with the data obtained in batch grinding of ?850/+600 ?m size quartz. The proposed equation enables calculation of the rate of breakage of the feed at any instant of time.

Homonuclear NOE in spatially periodic spin systems

The nuclear Overhauser effect equations are solved analytically for a homonuclear group of spins whose sites are periodically arranged, including the special cases where the spins lie at the vertices of a regular polygon and on a one-dimensional lattice. t is shown that, for long correlation times, the equations governing magnetization transfer resemble a diffusion equation. Furthermore the deviation from exact diffusion is quantitatively related to the molecular tumbling correlation time. Equations are derived for the range of magnetization travel subsequent to the perturbation of a single spin in a lattice for both the case of strictly dipolar relaxation and the more general situation where additional T1 mechanisms may be active. The theory given places no restrictions on the delay (or mixing) times, and it includes all the spins in the system. Simulations are presented to confirm the theory.

Molecular dynamics simulations on parallel and antiparallel C.G*G triplexes

Molecular dynamics (MD) studies have been carried out on the Hoogsteen hydrogen bonded parallel and the reverse Hoogsteen hydrogen banded antiparallel C.G*G triplexes. Earlier, the molecular mechanics studies had shown that the parallel structure was energetically more favourable than the antiparallel structure. To characterize the structural stability of the two triplexes and to investigate whether the antiparallel structure can transit to an energetically more favourable structure, due to the local fluctuations in the structure during the MD simulation, the two structures were subjected to 200ps of constant temperature vacuum MD simulations at 300K. Initially no constraints were applied to the structures and it was observed that for the antiparallel tripler, the structure showed a large root mean square deviation from the starting structure within the first 12ps and the N4-H41-O6 hydrogen bond in the WC duplex got distorted due to a high propeller twist and a moderate increase in the opening angle in the basepairs. Starting from an initial value of 30 degrees, helical twist of the average structure from this simulation had a value of 36 degrees, while the parallel structure stabilized at a twist of 33 degrees. In spite of the hydrogen bond distortions in the antiparallel tripler, it was energetically comparable to the parallel tripler. To examine the structural characteristics of an undistorted structure, another MD simulation was performed on the antiparallel tripler by constraining all the hydrogen bonds. This structure stabilized at an average twist of 33 degrees. In the course of the dynamics though the energy of the molecule - compared to the initial structure - improved, it did not become comparable to the parallel structure. Energy minimization studies performed in the presence of explicit water and counterions also showed the two structures to be equally favourable energetically Together these results indicate that the parallel C.G*G tripler with Hoogsteen hydrogen bonds also represents a stereochemically and energetically favourable structure for this class of triplexes.

(2E)-3-(4-Ethoxyphenyl)-1-(2-methyl-4-phenylquinolin-3-yl)prop-2-en-1-on e monohydrate

The title hydrate, C27H23NO2 center dot H2O, features an almost planar quinoline residue (r.m.s. deviation = 0.015 angstrom) with the benzene dihedral angle = 63.80 (7)degrees] and chalcone C-C-C-O torsion angle = -103.38 (18)degrees] substituents twisted significantly out of its plane. The configuration about the C=C bond 1.340 (2) angstrom] is E. In the crystal, molecules related by the 21 symmetry operation are linked along the b axis via water molecules that form O-H center dot center dot center dot O-c and O-H center dot center dot center dot N-q hydrogen bonds (c = carbonyl and q = quinoline). A C-H center dot center dot center dot O interaction also occurs.

Phase diagram of the two-dimensional disordered Hubbard model in the Hartree-Fock approximation

Ground-state properties of the two-dimensional Hubbard model with point-defect disorder are investigated numerically in the Hartree-Fock approximation. The phase diagram in the p(point defect concentration)-delta(deviation from half filling) plane exhibits antiferromagnetic, spin-density-wave, paramagnetic, and spin-glass-like phases. The disorder stabilizes the antiferromagnetic phase relative to the spin-density-wave phase. The presence of U strongly enhances the localization in the antiferromagnetic phase. The spin-density-wave and spin-glass-like phases are weakly localized.

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]

Surface tension of some binary organic solutions of alkylbenzenes and 1-alkanols

The surface tensions of binary mixtures of 1-alkanols (Cl-Cd with benzene, toluene, or xylene were measured. The results were correlated with the activity coefficients calculated through the group contribution method such as UNIFAC, with the maximum deviation from the experimental results less that 5%. The coefficients of the correlation are correlated with the chain length.

Dependence of magnetism on GdFeO3 distortion in the t(2g) system ARuO(3) (A = Sr, Ca)

We have examined the stability of the ferromagnetic (FM) state in CaRuO3 and SrRuO3 as a function of the GdFeO3 distortion. Model calculations predict the dependence of the FM transition temperature (T-c) on the rotation angle theta to vary as cos(2)(2 theta) for e(g)-electron systems. However, here, we find an initial increase and then the expected decrease. Furthermore, a much faster decrease is found than predicted for e(g)-electron systems. Considering the specific case of CaRuO3, a larger deviation of the Ru-O-Ru angle from 180 degrees in CaRuO3 as compared to SrRuO3 should result in a more reduced bandwidth, thereby making the former more correlated. The absence of long-range magnetic order in the more correlated CaRuO3 is traced to the strong collapse of various exchange interaction strengths that arises primarily from the volume reduction and increased distortion of the RuO6 octahedra network that accompanies the presence of a smaller ion at the A site.

Thermodynamic Study of Mixed Anionic Solid Solutions Using Gradient Solid Electrolytes System K2CO3 - K2SO4

The thermodynamic properties of K2CO3 -KSO, solid solutions with hexagonal structure have been measured using a solid-state cell, incorporating a composite solid electrolyte with step-changes in composition. The cell with the configuration Pt, CO2' + O2' || K2CO3 | K2(CO3)x(SO4)1-x || CO2'' + O2'' + Pt X =1 X=X was investigated in the temperature range of 925 to 1165 K. The composite gradient solid electrolyte consisted of pure K2CO3 at one extremity and the solid solution under study at the other. The Nernstian response of the cell to changes in partial pressures of CO2 and O2 at the electrodes and temperature was demonstrated. The activity of K2CO3 in the solid solution was measured by three techniques. All three methods gave identical results, indicating unit transport number for K+ ions and negligible diffusion potential due to concentration gradients of carbonate and sulfate ions. The activity of K2CO3 exhibits positive deviation from Raoult's law. The excess Gibbs energy of mixing of the solid solution can be represented using a subregular solution model DELTAG(E) = X(1 - X)[5030X + 4715(1 - X)] J mol-1 By combining this information with the phase diagram, mixing properties of the liquid phase were obtained.

Guessing Revisited: A Large Deviations Approach

The problem of guessing a random string is revisited. A close relation between guessing and compression is first established. Then it is shown that if the sequence of distributions of the information spectrum satisfies the large deviation property with a certain rate function, then the limiting guessing exponent exists and is a scalar multiple of the Legendre-Fenchel dual of the rate function. Other sufficient conditions related to certain continuity properties of the information spectrum are briefly discussed. This approach highlights the importance of the information spectrum in determining the limiting guessing exponent. All known prior results are then re-derived as example applications of our unifying approach.

2-(4-Methylsulfanylphenyl)-1H-benzimidazol-3-ium bromide

In the cation of the title compound, C14H13N2S+center dot Br-, the essentially planar benzimidazole system (r.m.s. deviation = 0.0082 angstrom) is substituted with a 4-methylsulfanylphenyl ring. The dihedral angle between the benzimidazole system and the 4-methylsulfanylphenyl ring is 2.133 (2)degrees. The crystal structure is characterized by strong and highly directional intermolecular N-H center dot center dot center dot Br hydrogen bonds involving the bromide ion. Moreover, C-H center dot center dot center dot S interactions result in chains of molecules along the c axis. The supramolecular assembly is further stabilized by pi-pi stacking interactions between the benzimidazole system and 4-methylsulfanylphenyl rings centroid-centroid distance = 3.477 (4) angstrom].

Signals of additional Z ` boson in e(+)e(-) -> W+W- at the ILC with polarized beams

We consider the possibility of fingerprinting the presence of heavy additional Z' bosons that arise naturally in extensions of the standard model such as E-6 models and left-right symmetric models, through their mixing with the standard model Z boson. By considering a class of observables including total cross sections, energy distributions and angular distributions of decay leptons we find significant deviation from the standard model predictions for these quantities with right-handed electrons and left-handed positrons at root s= 800GeV. The deviations being less pronounced at smaller centre of mass energies as the models are already tightly constrained. Our work suggests that the ILC should have a strong beam polarization physics program particularly with these configurations. On the other hand, a forward backward asymmetry and lepton fraction in the backward direction are more sensitive to new physics with realistic polarization due to interesting interplay with the neutrino t-channel diagram. This process complements the study of fermion pair production processes that have been considered for discrimination between these models.