Effect of coordinated ions on structure and flexibility of parallel G-quadruplexes: A molecular dynamics study

ingle 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(G(7)) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the O6 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. There 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 Angstrom 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 Nai 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, O6 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.

Molecular analysis of Friedreich's ataxia locus in the Indian population

Friedreich's ataxia (FRDA) is an autosomal recessive neurodegenerative disorder caused by expansion of GAA repeats in the frataxin gene. We have carried out the first molecular analysis at the Friedreich's ataxia locus in the Indian population. Materials and methods - Three families clinically diagnosed for Friedreich's ataxia were analyzed for GAA expansion at the FRDA locus. The distribution of GAA repeats was also estimated in normal individuals of Indian origin. Results - All patients clinically diagnosed for Friedreich's ataxia were found to be homozygous for GAA repeat expansion. The GAA repeat in the normal population show a bimodal distribution with 94% of alleles ranging from 7-16 repeats. Conclusion - Indian patients with expansion at the FRDA locus showed typical clinical features of Friedreich's ataxia. The low frequency of large normal alleles (6%) could indicate that the prevalence of this disease in the Indian population is likely to be low.

Asymmetry in structural and thermo-mechanical behavior of intermetallic NiAl nanowire under tensile/compressive loading: A molecular dynamics study

The asymmetric stress strain behavior under tension/compression in an initial < 100 > B-2-NiAl nanowire is investigated considering two different surface configurations i.e., < 100 >/(0 1 0) (0 0 1) and < 100 >/(0 1 1) (0 - 1 1). This behavior is attributed to two different deformation mechanisms namely a slip dominated deformation under compression and a known twinning dominated deformation under tension. It is also shown that B2 -> BCT (body-centered-tetragonal) phase transformation under tensile loading is independent of the surface configurations for an initial < 100 > oriented NiAl nanowire. Under tensile loading, the nanowire undergoes a stress-induced martensiticphase transformation from an initial B2 phase to BCT phase via twinning along {110} plane with failure strain of similar to 0.30. On the other hand, a compressive loading causes failure of these nanowires via brittle fracture after compressive yielding, with a maximum failure strain of similar to-0.12. Such brittle fracture under compressive loading occurs via slip along {110} plane without any phase transformations. Softening/hardening behavior is also reported for the first time in these nanowires under tensile/compressive loadings, which cause asymmetry in their yield strength behavior in the stress strain space. Result shows that a sharp increase in energy with increasing strain under compressive loading causes hardening of the nanowire, and hence, gives improved yield strength as compared to tensile loading. (C) 2010 Elsevier Ltd. All rights reserved.

Structure and dynamics of benzene in one-dimensional channels

Molecular dynamics investigation of benzene in one-dimensional channel systems A1PO(4)-5, VPI-5, and carbon nanotube is reported. The results suggest that, in all the three host systems, the plane of benzene is almost perpendicular to the channel axis when the molecule is near the center of the channel and the plane of benzene is parallel to the channel axis when the molecule is near the wall of the channel. The density distribution of benzene as a function of channel length, z and the radial distance, r, from the channel axis is also different in the three host structures. Anisotropy in translational diffusion coefficient, calculated in body-fixed frame of benzene, suggests that benzene prefers to move with its plane parallel to the direction of motion in A1PO(4)-5 and VPI-5 whereas in carbon nanotube the motion occurs predominantly with the plane of the benzene perpendicular to the direction of motion.;Anisotropy associated with the rotational motion is seen to alter significantly in confinement as compared to liquid benzene. In A1PO(4)-5, the rotational anisotropy is reversed as compared to liquid benzene thereby suggesting that anisotropy arising out of molecular geometry can be reduced. Reorientational correlation times for C-6 and C-2 axes Of benzene are reported. Apart from the inertial decay of reorientational correlation function due to free, rotation, two other distinct regimes of decay are observed in narrower channels (AIPO(4)-5 and carbon nanotube): (i) an initial fast decay (0.5-2 ps) and (ii) a slower decay (>2 ps) of reorientational correlation function where C-6 decays slower than C-2 Similar to what is observed in liquid benzene. In the initial fast decay, it is seen that the decay for C-6 is faster than C-2 which is in contrast to what is observed in liquid benzene or for benzene confined in VPI-5.

Molecular theory of solvation and solvation dynamics of a classical ion in a dipolar liquid

A microscopic theory of equilibrium solvation and solvation dynamics of a classical, polar, solute molecule in dipolar solvent is presented. Density functional theory is used to explicitly calculate the polarization structure around a solvated ion. The calculated solvent polarization structure is different from the continuum model prediction in several respects. The value of the polarization at the surface of the ion is less than the continuum value. The solvent polarization also exhibits small oscillations in space near the ion. We show that, under certain approximations, our linear equilibrium theory reduces to the nonlocal electrostatic theory, with the dielectric function (c(k)) of the liquid now wave vector (k) dependent. It is further shown that the nonlocal electrostatic estimate of solvation energy, with a microscopic c(k), is close to the estimate of linearized equilibrium theories of polar liquids. The study of solvation dynamics is based on a generalized Smoluchowski equation with a mean-field force term to take into account the effects of intermolecular interactions. This study incorporates the local distortion of the solvent structure near the ion and also the effects of the translational modes of the solvent molecules.The latter contribution, if significant, can considerably accelerate the relaxation of solvent polarization and can even give rise to a long time decay that agrees with the continuum model prediction. The significance of these results is discussed.

Molecular and crystal structure of 4-ethynylcyanobenzene. An example of Csingle bondHcdots, three dots, centeredN hydrogen bonding

The molecular and crystal structures of 4-ethynylcyanobenzene arereported. The packing of molecules in the crystal is found to be homologous with the crystal structures of HCN, cyanoacetylene and 4-cyano-4'-ethynylbiphenyl. Alternatively, these four crystals could be said to constitute a structural homologous series. The influence of C-H center dot center dot center dot N hydrogen bonding in directing a linear supramolecular arrangement of molecules with ethynyl and cyano groups at opposite ends, is illustrated. (C) 2010 Elsevier B.V. All rights reserved.

A molecular theory of collective orientational relaxation in pure and binary dipolar liquids

A molecular theory of collective orientational relaxation of dipolar molecules in a dense liquid is presented. Our work is based on a generalized, nonlinear, Smoluchowski equation (GSE) that includes the effects of intermolecular interactions through a mean‐field force term. The effects of translational motion of the liquid molecules on the orientational relaxation is also included self‐consistently in the GSE. Analytic expressions for the wave‐vector‐dependent orientational correlation functions are obtained for one component, pure liquid and also for binary mixtures. We find that for a dipolar liquid of spherical molecules, the correlation function ϕ(k,t) for l=1, where l is the rank of the spherical harmonics, is biexponential. At zero wave‐vector, one time constant becomes identical with the dielectric relaxation time of the polar liquid. The second time constant is the longitudinal relaxation time, but the contribution of this second component is small. We find that polar forces do not affect the higher order correlation functions (l>1) of spherical dipolar molecules in a linearized theory. The expression of ϕ(k,t) for a binary liquid is a sum of four exponential terms. We also find that the wave‐vector‐dependent relaxation times depend strongly on the microscopic structure of the dense liquid. At intermediate wave vectors, the translational diffusion greatly accelerates the rate of orientational relaxation. The present study indicates that one must pay proper attention to the microscopic structure of the liquid while treating the translational effects. An analysis of the nonlinear terms of the GSE is also presented. An interesting coupling between the number density fluctuation and the orientational fluctuation is uncovered.

Ab-initio molecular orbital studies on C---H…X hydrogen bonded systemsstar

C---H…X hydrogen bonded systems are studied by the STO-3G method. The proton donor ability of carbon is analysed in terms of its hybridization states and the substituents.

Chlorine-35 nuclear quadrupole resonance study of molecular torsional motion in some chloropyridine-N-oxides and chloroquinoline-N-oxides

Chlorine has been substituted at the 2- and 4-positions in the pyridine and quinoline rings of the corresponding N-oxides and 35Cl n.q.r. spectra have been studied in the temperature range 77–300 K. The change in the n.q.r. frequencies in N-oxides as compared to their parent compounds are interpreted in terms of the conjugative effect and the inductive effect of the N+—O– group. The negative temperature coefficients of the resonance frequencies in chloropyridine-N-oxides have been analysed using the Bayer, Kushida and Brown equations. The calculated torsional frequencies, which are in the range 52–78 cm–1, are found to be only slightly temperature dependent.

A self-consistent formulation for analysis and generation of non-uniform DNA structure

A fully self-consistent formulation is described here for the analysis and generation of base-pairs in non-uniform DNA structures, in terms of various local parameters. It is shown that the internal "wedge parameters" are mathematically related to the parameters describing the base-pair orientation with respect to an external helix axis. Hence any one set of three translation and three rotation parameters are necessary and sufficient to completely describe the relative orientation of the base-pairs comprising a step (or doublet). A general procedure is outlined for obtaining an average or global helix axis from the local helix axes for each step. A graphical representation of the local helix axes in the form of a polar plot is also shown and its application for estimating the curvature of oligonucleotide structures is illustrated, with examples of both A and B type structures.

Conformations of three heterocyclicperhydropyrrolobenzofurans and polymeric assembly via co-operative inter-molecular C-H center dot center dot center dot O hydrogen bonds

In 1-cyclo-hexyl-6,6,8a-trimethyl-3a,6,7,8a-tetra-hydro-1H-1-benzofuro[2,3-b]pyrrole-2,4(3H,5H)-dione, C19H27NO3, (I), and the isomorphous compounds 6,6,8a-trimethyl-1-phenyl-3a,6,7,8a-tetra-hydro-1H-1-benzofuro[2,3-b]pyrrole-2,4(3H,5H)-dione, C19H21NO3, (II), and 6,6,8a-trimethyl-1-(3-pyridyl)-3a,6,7,8a-tetra-hydro-1H-1-benzofuro[2,3-b]pyrrole-2,4(3H,5H)-dione, C18H20N2O3, (III), the tetra-hydro-benzo-dihydro-furo-pyrrolidine ring systems are folded at the cis junction of the five-membered rings, giving rise to a non-planar shape of the tricyclic cores. The dihydro-furan and pyrrolidine rings in (I) are puckered and adopt an envelope conformation. The cyclo-hexene rings adopt a half-chair conformation in all the mol-ecules, while the substituent N-cyclo-hexyl ring in (I) assumes a chair form. Short intra-molecular C-HcO contacts form S(5) and S(6) motifs. The isomorphous compounds (II) and (III) are effectively isostructural, and aggregate into chains via inter-molecular C-HcO hydrogen bonds.

Crystal and molecular structure of the ammonium salt of the dinucleoside monophosphate d(CpG)

The crystal and molecular structure of the ammonium salt of deoxycytidylyl-(3'-5')-deoxyguanosine has been determined from 0.85 A resolution single crystal X-ray diffraction data. The crystals obtained by acetone diffusion technique at -20 degrees C, are orthorhombic, P212121, a = 12.880(2), b = 17444(2) and c = 27.642(2) A. The structure was solved by high resolution Patterson and Fourier methods and refined to R = 0.136. There are two d(CpG) molecules in the asymmetric unit forming a mini left handed Z-DNA helix. This is in contrast to the earlier reported forms of d(CpG) where the molecules form self base paired duplexes. There are two ammonium ions in the asymmetric unit. The major groove NH+4 ion interacts with N7 of guanines through water bridges besides making H-bonded interactions directly with the phosphate oxygen atoms. A second NH+4 ion is found in the minor groove interacting directly with the phosphate oxygen atoms. Symmetry related molecules pack in such a way that the cytosine base stacks on cytosine and guanine base on guanine. Our structure demonstrates that alternating d(CpG) sequences have the ability to adopt the left handed Z-DNA structure even at the dimer level i.e., in a sequence which is only two base pairs long.

Molecular structure-activity relationship study of some non-steroidal antiinflammatory agents using electrostatic potential mapping

A series of 6,11-dihydro-11-oxodibenz[b,e]oxepin-2-acetic acids (DOAA) which are known to be anti-inflammatory agents were studied. The geometries of some of the molecules obtained from X-ray crystallography were used in the calculations as such while the geometries of their derivatives were obtained by local, partial geometry optimization around the Sites of substitution employing the AMI method, keeping the remaining parts of the geometries the same as those in the parent molecules. Molecular electrostatic potential (MEP) mapping was performed for the molecules using optimized hybridization displacement charges (HDC) combined with Lowdin charges, as this charge distribution has been shown earlier to yield near ab initio quality results. A good correlation has been found between the MEP values near the oxygen atoms of the hydroxyl groups of the carboxy groups of the molecules and their anti-inflammatory activities. The result is broadly in agreement with the model proposed earlier by other authors regarding the structure-activity relationship for other similar molecules.

Molecular Complexes of Alprazolam with Carboxylic Acids, Boric Acid, Boronic Acids, and Phenols. Evaluation of Supramolecular Heterosynthons Mediated by a Triazole Ring

A series of molecular complexes, both co-crystals and salts, of a triazole drug-alprazolam-with carboxylic acids, boric acid, boronic acids, and phenols have been analyzed with respect to heterosynthons present in the crystal structures. In all cases, the triazole ring behaves as an efficient hydrogen bond acceptor with the acidic coformers. The hydrogen bond patterns exhibited with aromatic carboxylic acids were found to depend on the nature and position of the substituents. Being a strong acid, 2,6-dihydroxybenzoic acid forms a salt with alprazolam. With aliphatic dicarboxylic acids alprazolam forms hydrates and the water molecules play a central role in synthon formation and crystal packing. The triazole ring makes two distinct heterosynthons in the molecular complex with boric acid. Boronic acids and phenols form consistent hydrogen bond patterns, and these are seemingly independent of the substitutional effects. Boronic acids form noncentrosymmetric cyclic synthons, while phenols form O-H center dot center dot center dot N hydrogen bonds with the triazole ring. (C) 2010 Wiley-Liss, Inc. and the American Pharmacists Association J Pharm Sci 99:3743-3753, 2010.