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

Treatment of two postoperative endophthalmitis cases due to Aspergillus flavus and Scopulariopsis spp. with local and systemic antifungal therapy

<p>Abstract</p> <p>Background</p> <p>Endophthalmitis is the inflammatory response to invasion of the eye with bacteria or fungi. The incidence of endophthalmitis after cataract surgery varies between 0.072–0.13 percent. Treatment of endophthalmitis with fungal etiology is difficult.</p> <p>Case Presentation</p> <p><b>Case 1: </b>A 71-year old male diabetic patient developed postoperative endophthalmitis due to <it>Aspergillus flavus</it>. The patient was treated with topical amphotericin B ophthalmic solution, intravenous (IV) liposomal amphotericin-B and caspofungin following vitrectomy.</p> <p><b>Case 2: </b>A 72-year old male cachectic patient developed postoperative endophthalmitis due to <it>Scopulariopsis </it>spp. The patient was treated with topical and IV voriconazole and caspofungin.</p> <p>Conclusion</p> <p><it>Aspergillus </it>spp. are responsible of postoperative fungal endophthalmitis. Endophthalmitis caused by <it>Scopulariopsis </it>spp. is a very rare condition. The two cases were successfully treated with local and systemic antifungal therapy.</p

Insights into the Complex Formed by Matrix Metalloproteinase-2 and Alloxan Inhibitors: Molecular Dynamics Simulations and Free Energy Calculations

Matrix metalloproteinases (MMP) are well-known biological targets implicated in tumour progression, homeostatic regulation, innate immunity, impaired delivery of pro-apoptotic ligands, and the release and cleavage of cell-surface receptors. Hence, the development of potent and selective inhibitors targeting these enzymes continues to be eagerly sought. In this paper, a number of alloxan-based compounds, initially conceived to bias other therapeutically relevant enzymes, were rationally modified and successfully repurposed to inhibit MMP-2 (also named gelatinase A) in the nanomolar range. Importantly, the alloxan core makes its debut as zinc binding group since it ensures a stable tetrahedral coordination of the catalytic zinc ion in concert with the three histidines of the HExxHxxGxxH metzincin signature motif, further stabilized by a hydrogen bond with the glutamate residue belonging to the same motif. The molecular decoration of the alloxan core with a biphenyl privileged structure allowed to sample the deep S1′ specificity pocket of MMP-2 and to relate the high affinity towards this enzyme with the chance of forming a hydrogen bond network with the backbone of Leu116 and Asn147 and the side chains of Tyr144, Thr145 and Arg149 at the bottom of the pocket. The effect of even slight structural changes in determining the interaction at the S1′ subsite of MMP-2 as well as the nature and strength of the binding is elucidated via molecular dynamics simulations and free energy calculations. Among the herein presented compounds, the highest affinity (pIC50 = 7.06) is found for BAM, a compound exhibiting also selectivity (>20) towards MMP-2, as compared to MMP-9, the other member of the gelatinases

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

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

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

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

A Kinetic Model of Trp-Cage Folding from Multiple Biased Molecular Dynamics Simulations

Trp-cage is a designed 20-residue polypeptide that, in spite of its size, shares several features with larger globular proteins. Although the system has been intensively investigated experimentally and theoretically, its folding mechanism is not yet fully understood. Indeed, some experiments suggest a two-state behavior, while others point to the presence of intermediates. In this work we show that the results of a bias-exchange metadynamics simulation can be used for constructing a detailed thermodynamic and kinetic model of the system. The model, although constructed from a biased simulation, has a quality similar to those extracted from the analysis of long unbiased molecular dynamics trajectories. This is demonstrated by a careful benchmark of the approach on a smaller system, the solvated Ace-Ala3-Nme peptide. For the Trp-cage folding, the model predicts that the relaxation time of 3100 ns observed experimentally is due to the presence of a compact molten globule-like conformation. This state has an occupancy of only 3% at 300 K, but acts as a kinetic trap. Instead, non-compact structures relax to the folded state on the sub-microsecond timescale. The model also predicts the presence of a state at of 4.4 Å from the NMR structure in which the Trp strongly interacts with Pro12. This state can explain the abnormal temperature dependence of the and chemical shifts. The structures of the two most stable misfolded intermediates are in agreement with NMR experiments on the unfolded protein. Our work shows that, using biased molecular dynamics trajectories, it is possible to construct a model describing in detail the Trp-cage folding kinetics and thermodynamics in agreement with experimental data