Conformational and Structural Relaxations of Poly(ethylene oxide) and Poly(propylene oxide) Melts: Molecular Dynamics Study of Spatial Heterogeneity, Cooperativity, and Correlated Forward-Backward Motion

Performing molecular dynamics simulations for all-atom models, we characterize the conformational and structural relaxations of poly(ethylene oxide) and poly(propylene oxide) melts. The temperature dependence of these relaxation processes deviates from an Arrhenius law for both polymers. We demonstrate that mode-coupling theory captures some aspects of the glassy slowdown, but it does not enable a complete explanation of the dynamical behavior. When the temperature is decreased, spatially heterogeneous and cooperative translational dynamics are found to become more important for the structural relaxation. Moreover, the transitions between the conformational states cease to obey Poisson statistics. In particular, we show that, at sufficiently low temperatures, correlated forward-backward motion is an important aspect of the conformational relaxation, leading to strongly nonexponential distributions for the waiting times of the dihedrals in the various conformational statesComment: 13 pages, 13 figure

A critical assessment of methods for the intrinsic analysis of liquid interfaces: 2. density profiles

Substantial improvements in the molecular level understanding of fluid interfaces have recently been achieved by recognizing the importance of detecting the intrinsic surface of the coexisting condensed phases in computer simulations (i.e., after the removal of corrugations caused by capillary waves) and by developing several methods for identifying the molecules that are indeed located at the boundary of the two phases. In our previous paper [J. Phys. Chem. C 2010, 114, 11169], we critically compared those methods in terms of reliability, robustness, and computation speed. Once the intrinsic surface of a given phase is detected, various profiles, such as the density profiles of the components, can be calculated relative to this intrinsic surface rather than to the macroscopically planar Gibbs dividing surface. As a continuation of our previous study, here we present a detailed and critical comparison of various methods that can be used to calculate intrinsic density profiles once the full set of truly interfacial molecules has been identified. Two of the methods, the Fourier function and the Voronoi tessellation, are already described in the literature; two other methods, the covering surface and the triangular interpolation, are newly proposed algorithms; one method, the modified grid-based intrinsic profile (GIP) method, is an improvement over an existing procedure. The different methods are again compared in terms of accuracy and computational cost. On the basis of this comparison, we propose a fast and accurate protocol to be routinely used for intrinsic surface analyses in computer simulations

A critical assessment of methods for the intrinsic analysis of liquid interfaces. 1. surface site distributions

Substantial progress in our understanding of interfacial structure and dynamics has stemmed from the recent development of algorithms that allow for an intrinsic analysis of fluid interfaces. These work by identifying the instantaneous location of the interface, at the atomic level, for each molecular configuration and then computing properties relative to this location. Such a procedure eliminates the broadening of the interface caused by capillary waves and reveals the underlying features of the system. However, a precise definition of which molecules actually belong to the interfacial layer is difficult to achieve in practice. Furthermore, it is not known if the different intrinsic analysis methods are consistent with each other and yield similar results for the interfacial properties. In this paper, we carry out a systematic and detailed comparison of the available methods for intrinsic analysis of fluid interfaces, based on a molecular dynamics simulation of the interface between liquid water and carbon tetrachloride. We critically assess the advantages and shortcomings of each method, based on reliability, robustness, and speed of computation, and establish consistent criteria for determining which molecules belong to the surface layer. We believe this will significantly contribute to make intrinsic analysis methods widely and routinely applicable to interfacial systems

Solvation free energy profile of the SCN- ion across the water-1,2-dichloroethane liquid/liquid interface. A computer simulation study

The solvation free energy profile of a single SCN- ion is calculated across the water-1,2-dichloroethane liquid/liquid interface at 298 K by the constraint force method. The obtained results show that the free energy cost of transferring the ion from the aqueous to the organic phase is about 70 kJ/mol, The free energy profile shows a small but clear well at the aqueous side of the interface, in the subsurface region of the water phase, indicating the ability of the SCN- ion to be adsorbed in the close vicinity of the interface. Upon entrance of the SCN- ion to the organic phase a coextraction of the water molecules of its first hydration shell occurs. Accordingly, when it is located at the boundary of the two phases the SCN- ion prefers orientations in which its bulky S atom is located at the aqueous side, and the small N atom, together with its first hydration shell, at the organic side of the interface

Properties of the Liquid-Vapor Interface of Acetone-Water Mixtures. A Computer Simulation and ITIM Analysis Study

Molecular dynamics simulations of the liquid-vapor interface of acetone-water mixtures of different compositions, covering the entire composition range have been performed on the canonical (N, V, T) ensemble at 298 K, using a model combination that excellently describes the mixing properties of these compounds. The properties of the intrinsic liquid surfaces have been analyzed in terms of the Identification of the Truly Interfacial Molecules (ITIM) method. Thus, the composition, width, roughness, and separation of the subsurface molecular layers, as well as self-association, orientation, and dynamics of exchange with the bulk phase of the surface molecules have been analyzed in detail. Our results show that acetone molecules are strongly adsorbed at the liquid surface, and this adsorption extends to several molecular layers. Like molecules in the surface layer are found to form relatively large lateral self-associates. The effect of the vicinity of the vapor phase on a number of properties of the liquid phase vanishes beyond the first molecular layer, with the second subsurface layer already part of the bulk liquid phase in these respects. The orientational preferences of the surface molecules are governed primarily by the dipole-dipole interaction of the neighboring acetone molecules, and hydrogen bonding interaction of the neighboring acetone-water pairs. (Figure Presented). © 2015 American Chemical Society

Active Nuclear Receptors Exhibit Highly Correlated AF-2 Domain Motions

Nuclear receptor ligand binding domains (LBDs) convert ligand binding events into changes in gene expression by recruiting transcriptional coregulators to a conserved activation function-2 (AF-2) surface. While most nuclear receptor LBDs form homo- or heterodimers, the human nuclear receptor pregnane X receptor (PXR) forms a unique and essential homodimer and is proposed to assemble into a functional heterotetramer with the retinoid X receptor (RXR). How the homodimer interface, which is located 30 Å from the AF-2, would affect function at this critical surface has remained unclear. By using 20- to 30-ns molecular dynamics simulations on PXR in various oligomerization states, we observed a remarkably high degree of correlated motion in the PXR–RXR heterotetramer, most notably in the four helices that create the AF-2 domain. The function of such correlation may be to create “active-capable” receptor complexes that are ready to bind to transcriptional coactivators. Indeed, we found in additional simulations that active-capable receptor complexes involving other orphan or steroid nuclear receptors also exhibit highly correlated AF-2 domain motions. We further propose a mechanism for the transmission of long-range motions through the nuclear receptor LBD to the AF-2 surface. Taken together, our findings indicate that long-range motions within the LBD scaffold are critical to nuclear receptor function by promoting a mobile AF-2 state ready to bind coactivators

Floating Patches of HCN at the Surface of Their Aqueous Solutions - Can They Make "HCN World" Plausible?

The liquid/vapor interface of the aqueous solutions of HCN of different concentrations has been investigated using molecular dynamics simulation and intrinsic surface analysis. Although HCN is fully miscible with water, strong interfacial adsorption of HCN is observed at the surface of its aqueous solutions, and, at the liquid surface, the HCN molecules tend to be located even at the outer edge of the surface layer. It turns out that in dilute systems the HCN concentration can be about an order of magnitude larger in the surface layer than in the bulk liquid phase. Furthermore, HCN molecules show a strong lateral self-association behavior at the liquid surface, forming thus floating HCN patches at the surface of their aqueous solutions. Moreover, HCN molecules are staying, on average, an order of magnitude longer at the liquid surface than water molecules, and this behavior is more pronounced at smaller HCN concentrations. Because of this enhanced dynamical stability, the floating HCN patches can provide excellent spots for polymerization of HCN, which can be the key step in the prebiotic synthesis of partially water-soluble adenine. All of these findings make the hypothesis of "HCN world" more plausible

Intrinsic structure and dynamics of the water/nitrobenzene interface

In this paper we present results of a detailed and systematic molecular dynamics study of the water/nitrobenzene interface. Using a simple procedure to eliminate fluctuations of the interface position, we are able to obtain true intrinsic profiles for several properties (density, hydrogen bonds, molecular orientation, etc.) in the direction perpendicular to the interfacial plane. Our results show that both water and organic inter-facial molecules form a tightly packed layer oriented parallel to the interface, with reduced mobility in the perpendicular direction. Beyond this layer, water quickly restores its bulk structure, while nitrobenzene exhibits structural anisotropies that extend further into the bulk region: Water molecules that protrude farthest into the organic phase point one hydrogen atom in the direction perpendicular to the interface, forming a hydrogen bond with a nitrobenzene oxygen. By fitting both the global and the intrinsic density profiles, we obtain estimates for the total and intrinsic interface widths, respectively. These are combined with capillary wave theory to produce a self-consistent method for the calculation of the inter-facial tension. Values calculated using this method are in very good agreement with direct calculations from the components of the pressure tensor

Structure of the dimeric N-glycosylated form of fungal β-N-acetylhexosaminidase revealed by computer modeling, vibrational spectroscopy, and biochemical studies

<p>Abstract</p> <p>Background</p> <p>Fungal β-<it>N</it>-acetylhexosaminidases catalyze the hydrolysis of chitobiose into its constituent monosaccharides. These enzymes are physiologically important during the life cycle of the fungus for the formation of septa, germ tubes and fruit-bodies. Crystal structures are known for two monomeric bacterial enzymes and the dimeric human lysosomal β-<it>N</it>-acetylhexosaminidase. The fungal β-<it>N</it>-acetylhexosaminidases are robust enzymes commonly used in chemoenzymatic syntheses of oligosaccharides. The enzyme from <it>Aspergillus oryzae </it>was purified and its sequence was determined.</p> <p>Results</p> <p>The complete primary structure of the fungal β-<it>N</it>-acetylhexosaminidase from <it>Aspergillus oryzae </it>CCF1066 was used to construct molecular models of the catalytic subunit of the enzyme, the enzyme dimer, and the <it>N</it>-glycosylated dimer. Experimental data were obtained from infrared and Raman spectroscopy, and biochemical studies of the native and deglycosylated enzyme, and are in good agreement with the models. Enzyme deglycosylated under native conditions displays identical kinetic parameters but is significantly less stable in acidic conditions, consistent with model predictions. The molecular model of the deglycosylated enzyme was solvated and a molecular dynamics simulation was run over 20 ns. The molecular model is able to bind the natural substrate – chitobiose with a stable value of binding energy during the molecular dynamics simulation.</p> <p>Conclusion</p> <p>Whereas the intracellular bacterial β-<it>N</it>-acetylhexosaminidases are monomeric, the extracellular secreted enzymes of fungi and humans occur as dimers. Dimerization of the fungal β-<it>N</it>-acetylhexosaminidase appears to be a reversible process that is strictly pH dependent. Oligosaccharide moieties may also participate in the dimerization process that might represent a unique feature of the exclusively extracellular enzymes. Deglycosylation had only limited effect on enzyme activity, but it significantly affected enzyme stability in acidic conditions. Dimerization and <it>N</it>-glycosylation are the enzyme's strategy for catalytic subunit stabilization. The disulfide bridge that connects Cys⁴⁴⁸with Cys⁴⁸³stabilizes a hinge region in a flexible loop close to the active site, which is an exclusive feature of the fungal enzymes, neither present in bacterial nor mammalian structures. This loop may play the role of a substrate binding site lid, anchored by a disulphide bridge that prevents the substrate binding site from being influenced by the flexible motion of the loop.</p

Roles of Electrostatics and Conformation in Protein-Crystal Interactions

In vitro studies have shown that the phosphoprotein osteopontin (OPN) inhibits the nucleation and growth of hydroxyapatite (HA) and other biominerals. In vivo, OPN is believed to prevent the calcification of soft tissues. However, the nature of the interaction between OPN and HA is not understood. In the computational part of the present study, we used molecular dynamics simulations to predict the adsorption of 19 peptides, each 16 amino acids long and collectively covering the entire sequence of OPN, to the {100} face of HA. This analysis showed that there is an inverse relationship between predicted strength of adsorption and peptide isoelectric point (P<0.0001). Analysis of the OPN sequence by PONDR (Predictor of Naturally Disordered Regions) indicated that OPN sequences predicted to adsorb well to HA are highly disordered. In the experimental part of the study, we synthesized phosphorylated and non-phosphorylated peptides corresponding to OPN sequences 65–80 (pSHDHMDDDDDDDDDGD) and 220–235 (pSHEpSTEQSDAIDpSAEK). In agreement with the PONDR analysis, these were shown by circular dichroism spectroscopy to be largely disordered. A constant-composition/seeded growth assay was used to assess the HA-inhibiting potencies of the synthetic peptides. The phosphorylated versions of OPN65-80 (IC50 = 1.93 µg/ml) and OPN220-235 (IC50 = 1.48 µg/ml) are potent inhibitors of HA growth, as is the nonphosphorylated version of OPN65-80 (IC50 = 2.97 µg/ml); the nonphosphorylated version of OPN220-235 has no measurable inhibitory activity. These findings suggest that the adsorption of acidic proteins to Ca2+-rich crystal faces of biominerals is governed by electrostatics and is facilitated by conformational flexibility of the polypeptide chain