Crystalline ethane-1,2-diol does not have intra-molecular hydrogen bonding: Experimental and theoretical charge density studies

The X-ray structure and electron density distribution of ethane-1,2-diol (ethylene glycol), obtained at a resolution extending to 1.00 Å−1 in sin θ/λ (data completion = 100% at 100 K) by in situ cryocrystallization technique is reported. The diol is in the gauche (g′Gt) conformation with the crystal structure stabilised by a network of inter-molecular hydrogen bonds. In addition to the well-recognized O–H···O hydrogen bonds there is topological evidence for C–H···O inter-molecular interactions. There is no experimental electron density based topological evidence for the occurrence of an intra-molecular hydrogen bond. The O···H spacing is not, vert, similar0.45 Å greater than in the gas-phase with an O–H···O angle close to 90°, calling into question the general assumption that the gauche conformation of ethane-1,2-diol is stabilised by the intra-molecular oxygen–hydrogen interaction.

Molecular Dynamics Simulations of Homogeneous Nucleation

Nucleation is the first step in a phase transition where small nuclei of the new phase start appearing in the metastable old phase, such as the appearance of small liquid clusters in a supersaturated vapor. Nucleation is important in various industrial and natural processes, including atmospheric new particle formation: between 20 % to 80 % of atmospheric particle concentration is due to nucleation. These atmospheric aerosol particles have a significant effect both on climate and human health. Different simulation methods are often applied when studying things that are difficult or even impossible to measure, or when trying to distinguish between the merits of various theoretical approaches. Such simulation methods include, among others, molecular dynamics and Monte Carlo simulations. In this work molecular dynamics simulations of the homogeneous nucleation of Lennard-Jones argon have been performed. Homogeneous means that the nucleation does not occur on a pre-existing surface. The simulations include runs where the starting configuration is a supersaturated vapor and the nucleation event is observed during the simulation (direct simulations), as well as simulations of a cluster in equilibrium with a surrounding vapor (indirect simulations). The latter type are a necessity when the conditions prevent the occurrence of a nucleation event in a reasonable timeframe in the direct simulations. The effect of various temperature control schemes on the nucleation rate (the rate of appearance of clusters that are equally able to grow to macroscopic sizes and to evaporate) was studied and found to be relatively small. The method to extract the nucleation rate was also found to be of minor importance. The cluster sizes from direct and indirect simulations were used in conjunction with the nucleation theorem to calculate formation free energies for the clusters in the indirect simulations. The results agreed with density functional theory, but were higher than values from Monte Carlo simulations. The formation energies were also used to calculate surface tension for the clusters. The sizes of the clusters in the direct and indirect simulations were compared, showing that the direct simulation clusters have more atoms between the liquid-like core of the cluster and the surrounding vapor. Finally, the performance of various nucleation theories in predicting simulated nucleation rates was investigated, and the results among other things highlighted once again the inadequacy of the classical nucleation theory that is commonly employed in nucleation studies.

Molecular Basis for the Role of Staphylococcus aureus Penicillin Binding Protein 4 in Antimicrobial Resistance

Penicillin binding proteins (PBPs) are membrane-associated proteins that catalyze the final step of murein biosynthesis. These proteins function as either transpeptidases or carboxypeptidases and in a few cases demonstrate transglycosylase activity. Both transpeptidase and carboxypeptidase activities of PBPs occur at the D-Ala-D-Ala terminus of a murein precursor containing a disaccharide pentapeptide comprising N-acetyl-glucosamine and N-acetyl-muramic acid-L-Ala-D-Glu-L-Lys-D-Ala-D-Ala. beta-Lactam antibiotics inhibit these enzymes by competing with the pentapeptide precursor for binding to the active site of the enzyme. Here we describe the crystal structure, biochemical characteristics, and expression profile of PBP4, a low-molecular-mass PBP from Staphylococcus aureus strain COL. The crystal structures of PBP4-antibiotic complexes reported here were determined by molecular replacement, using the atomic coordinates deposited by the New York Structural Genomics Consortium. While the pbp4 gene is not essential for the viability of S. aureus, the knockout phenotype of this gene is characterized by a marked reduction in cross-linked muropeptide and increased vancomycin resistance. Unlike other PBPs, we note that expression of PBP4 was not substantially altered under different experimental conditions, nor did it change across representative hospital- or community-associated strains of S. aureus that were examined. In vitro data on purified recombinant S. aureus PBP4 suggest that it is a beta-lactamase and is not trapped as an acyl intermediate with beta-lactam antibiotics. Put together, the expression analysis and biochemical features of PBP4 provide a framework for understanding the function of this protein in S. aureus and its role in antimicrobial resistance.

Two-Phase Thermodynamic Model for Efficient and Accurate Absolute Entropy of Water from Molecular Dynamics Simulations

Presented here is the two-phase thermodynamic (2PT) model for the calculation of energy and entropy of molecular fluids from the trajectory of molecular dynamics (MD) simulations. In this method, the density of state (DoS) functions (including the normal modes of translation, rotation, and intramolecular vibration motions) are determined from the Fourier transform of the corresponding velocity autocorrelation functions. A fluidicity parameter (f), extracted from the thermodynamic state of the system derived from the same MD, is used to partition the translation and rotation modes into a diffusive, gas-like component (with 3Nf degrees of freedom) and a nondiffusive, solid-like component. The thermodynamic properties, including the absolute value of entropy, are then obtained by applying quantum statistics to the solid component and applying hard sphere/rigid rotor thermodynamics to the gas component. The 2PT method produces exact thermodynamic properties of the system in two limiting states: the nondiffusive solid state (where the fluidicity is zero) and the ideal gas state (where the fluidicity becomes unity). We examine the 2PT entropy for various water models (F3C, SPC, SPC/E, TIP3P, and TIP4P-Ew) at ambient conditions and find good agreement with literature results obtained based on other simulation techniques. We also validate the entropy of water in the liquid and vapor phases along the vapor-liquid equilibrium curve from the triple point to the critical point. We show that this method produces converged liquid phase entropy in tens of picoseconds, making it an efficient means for extracting thermodynamic properties from MD simulations.

Electronic structures of electron donor-acceptor complexes: results from ultraviolet photoelectron spectroscopy and molecular orbital calculations

HeI photoelectron spectra of 1:1 electron donor-acceptor complexes are discussed in the light of molecular orbital calculations. The complexes discussed include those formed by BH3, BF3 and SO2. Some systematics have been found in the ionization energy shifts of the complexes compared to the free components and these are related to the strength of the donor-acceptor bond. Hel spectra of hydrogen bonded complexes are discussed in comparison with results from MO calculations. Limitations of such studies as well as scope for further investigations are indicated.

Molecular Caches: A caching structure for dynamic creation of application-specific Heterogeneous cache regions

CMPs enable simultaneous execution of multiple applications on the same platforms that share cache resources. Diversity in the cache access patterns of these simultaneously executing applications can potentially trigger inter-application interference, leading to cache pollution. Whereas a large cache can ameliorate this problem, the issues of larger power consumption with increasing cache size, amplified at sub-100nm technologies, makes this solution prohibitive. In this paper in order to address the issues relating to power-aware performance of caches, we propose a caching structure that addresses the following: 1. Definition of application-specific cache partitions as an aggregation of caching units (molecules). The parameters of each molecule namely size, associativity and line size are chosen so that the power consumed by it and access time are optimal for the given technology. 2. Application-Specific resizing of cache partitions with variable and adaptive associativity per cache line, way size and variable line size. 3. A replacement policy that is transparent to the partition in terms of size, heterogeneity in associativity and line size. Through simulation studies we establish the superiority of molecular cache (caches built as aggregations of molecules) that offers a 29% power advantage over that of an equivalently performing traditional cache.

Am unexpected entry into the silver(I) chemistry of diphosphazane ligands: Crystal and molecular structures of [Ag{Ph2PN(Pr-i)PPh2}(CF3SO3)](2), [Ag{Ph2PN(Pr-i)PPh(OC6H3Me2-2,6)}(2)]PF6, and [Ag-3(mu-Cl)(2){Ph2PN(R)PPh2}(3)]PF6 (R = H, Pr-i)

Mononuclear, binuclear and trinuclear silver(l) complexes were obtained unexpectedly while probing the reactivity of diphosphazane ligands of the type X2PN(Pr-i)PXY towards the ruthenium-based precursor Ru(bipy)(2)Cl-2 center dot 2H(2)O, in the presence of a silver salt as a chloride scavenger. Subsequently, the reactions of AgX [X = Cl, NO3 or CF3SO3] with Ph2PN(R)PPh(Y) [R = H, Y = Ph; R = Pr-i, Y = Ph or OC6H3Me2-2,6] in a 1: 1 or 1:2 molar ratio have been investigated. Mononuclear or binuclear Ag(I) complexes containing either chelating or bridging diphosphazane ligands are obtained. Trinuclear silver(l) complexes are accessible by the treatment of diphosphazane ligands, Ph2PN(R)PPh2 [R = H, Pr-i] with AgCl using piperidine as the solvent. In the presence of a suitable chloride donor species, the mononuclear and binuclear complexes of Ph2PN(Pr-i)PPh2 are transformed slowly to the trinuclear complex [Ag-3(mu-Cl)(2){Ph2PN(Pr-i)PPh2}(3)]X, over a period 20 h. The structures of representative complexes have been confirmed by X-ray crystallography and the salient structural features are discussed

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.

Phylogenetics and biogeography of a spectacular Old World radiation of butterflies: the subtribe Mycalesina (Lepidoptera: Nymphalidae: Satyrini)

Background: Butterflies of the subtribe Mycalesina (Nymphalidae: Satyrinae) are important model organisms in ecology and evolution. This group has radiated spectacularly in the Old World tropics and presents an exciting opportunity to better understand processes of invertebrate rapid radiations. However, the generic-level taxonomy of the subtribe has been in a constant state of flux, and relationships among genera are unknown. There are six currently recognized genera in the group. Mycalesis, Lohora and Nirvanopsis are found in the Oriental region, the first of which is the most speciose genus among mycalesines, and extends into the Australasian region. Hallelesis and Bicyclus are found in mainland Africa, while Heteropsis is primarily Madagascan, with a few species in Africa. We infer the phylogeny of the group with data from three genes (total of 3139 bp) and use these data to reconstruct events in the biogeographic history of the group.,Results: The results indicate that the group Mycalesina radiated rapidly around the Oligocene-Miocene boundary. Basal relationships are unresolved, but we recover six well-supported clades. Some species of Mycalesis are nested within a primarily Madagascan clade of Heteropsis,while Nirvanopsis is nested within Lohora. The phylogeny suggests that the group had its origin either in Asia or Africa, and diversified through dispersals between the two regions, during the late Oligocene and early Miocene. The current dataset tentatively suggests that the Madagascan fauna comprises two independent radiations. The Australasian radiation shares a common ancestor derived from Asia. We discuss factors that are likely to have played a key role in the diversification of the group. Conclusions: We propose a significantly revised classification scheme for Mycalesina. We conclude that the group originated and radiated from an ancestor that was found either in Asia or Africa, with dispersals between the two regions and to Australasia. Our phylogeny paves the way for further comparative studies on this group that will help us understand the processes underlying diversification in rapid radiations of invertebrates.

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.