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

Serotonin and corticosterone rhythms in mice exposed to cigarette smoke and in patients with COPD:implication for COPD-associated neuropathogenesis

The circadian timing system controls daily rhythms of physiology and behavior, and disruption of clock function can trigger stressful life events. Daily exposure to cigarette smoke (CS) can lead to alteration in diverse biological and physiological processes. Smoking is associated with mood disorders, including depression and anxiety. Patients with chronic obstructive pulmonary disease (COPD) have abnormal circadian rhythms, reflected by daily changes in respiratory symptoms and lung function. Corticosterone (CORT) is an adrenal steroid that plays a considerable role in stress and anti-inflammatory responses. Serotonin (5-hydroxytryptamine; 5HT) is a neurohormone, which plays a role in sleep/wake regulation and affective disorders. Secretion of stress hormones (CORT and 5HT) is under the control of the circadian clock in the suprachiasmatic nucleus. Since smoking is a contributing factor in the development of COPD, we hypothesize that CS can affect circadian rhythms of CORT and 5HT secretion leading to sleep and mood disorders in smokers and patients with COPD. We measured the daily rhythms of plasma CORT and 5HT in mice following acute (3 d), sub-chronic (10 d) or chronic (6 mo) CS exposure and in plasma from non-smokers, smokers and patients with COPD. Acute and chronic CS exposure affected both the timing (peak phase) and amplitude of the daily rhythm of plasma CORT and 5HT in mice. Acute CS appeared to have subtle time-dependent effects on CORT levels but more pronounced effects on 5HT. As compared with CORT, plasma 5HT was slightly elevated in smokers but was reduced in patients with COPD. Thus, the effects of CS on plasma 5HT were consistent between mice and patients with COPD. Together, these data reveal a significant impact of CS exposure on rhythms of stress hormone secretion and subsequent detrimental effects on cognitive function, depression-like behavior, mood/anxiety and sleep quality in smokers and patients with COPD

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

The Dynamics of Ca2+ Ions within the Solvation Shell of Calbindin D9k

The encounter of a Ca2+ ion with a protein and its subsequent binding to specific binding sites is an intricate process that cannot be fully elucidated from experimental observations. We have applied Molecular Dynamics to study this process with atomistic details, using Calbindin D9k (CaB) as a model protein. The simulations show that in most of the time the Ca2+ ion spends within the Debye radius of CaB, it is being detained at the 1st and 2nd solvation shells. While being detained near the protein, the diffusion coefficient of the ion is significantly reduced. However, due to the relatively long period of detainment, the ion can scan an appreciable surface of the protein. The enhanced propagation of the ion on the surface has a functional role: significantly increasing the ability of the ion to scan the protein's surface before being dispersed to the bulk. The contribution of this mechanism to Ca2+ binding becomes significant at low ion concentrations, where the intervals between successive encounters with the protein are getting longer. The efficiency of the surface diffusion is affected by the distribution of charges on the protein's surface. Comparison of the Ca2+ binding dynamics in CaB and its E60D mutant reveals that in the wild type (WT) protein the carboxylate of E60 function as a preferred landing-site for the Ca2+ arriving from the bulk, followed by delivering it to the final binding site. Replacement of the glutamate by aspartate significantly reduced the ability to transfer Ca2+ ions from D60 to the final binding site, explaining the observed decrement in the affinity of the mutated protein to Ca2+

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

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

General Anesthetics Predicted to Block the GLIC Pore with Micromolar Affinity

Although general anesthetics are known to modulate the activity of ligand-gated ion channels in the Cys-loop superfamily, there is at present neither consensus on the underlying mechanisms, nor predictive models of this modulation. Viable models need to offer quantitative assessment of the relative importance of several identified anesthetic binding sites. However, to date, precise affinity data for individual sites has been challenging to obtain by biophysical means. Here, the likely role of pore block inhibition by the general anesthetics isoflurane and propofol of the prokaryotic pentameric channel GLIC is investigated by molecular simulations. Microscopic affinities are calculated for both single and double occupancy binding of isoflurane and propofol to the GLIC pore. Computations are carried out for an open-pore conformation in which the pore is restrained to crystallographic radius, and a closed-pore conformation that results from unrestrained molecular dynamics equilibration of the structure. The GLIC pore is predicted to be blocked at the micromolar concentrations for which inhibition by isofluorane and propofol is observed experimentally. Calculated affinities suggest that pore block by propofol occurs at signifcantly lower concentrations than those for which inhibition is observed: we argue that this discrepancy may result from binding of propofol to an allosteric site recently identified by X-ray crystallography, which may cause a competing gain-of-function effect. Affinities of isoflurane and propofol to the allosteric site are also calculated, and shown to be 3 mM for isoflurane and for propofol; both anesthetics have a lower affinity for the allosteric site than for the unoccupied pore