Temperature dependence of molecular dynamics and supramolecular aggregation in MEH-PPV films: A solid-state NMR, X-ray and fluorescence spectroscopy study

This article presents an investigation of the temperature induced modification in the microstructure and dynamics of poly[2-methoxy-5-(2`-ethylhexyloxy)-1,4-phenylenevinylene] (MEH-PPV) cast films using Wide-Angle X-ray Scattering (WAXS), solid-state Nuclear Magnetic Resonance (NMR), and Fluorescence Spectroscopy (PL). MEH-PPV chain motions were characterized as a function of temperature by NMR. The results indicated that the solvent used to cast the films influences the activation energy of the side-chain motions. This was concluded from the comparison of the activation energy of the toluene cast film, E(a) = (54 +/- 8) kJ/mol, and chloroform cast film, E(a) = (69 +/- 5) kJ/mol, and could be attributed to the higher side-chain packing provided by chloroform, that preferentially solvates the side chain in contrast to toluene that solvates mainly the backbone. Concerning the backbone mobility, it was observed that the torsional motions in the MEH-PPV have average amplitude of similar to 10 degrees at 300 K, which was found to be independent of the solvent used to cast the films. In order to correlate the molecular dynamics processes with the changes in the microstructure of the polymer, in situ WAXS experiments as a function of temperature were performed and revealed that the interchain spacing in the MEH-PPV molecular aggregates increases as a function of temperature, particularly at temperatures where molecular relaxations occur. It was also observed that the WAXS peak associated with the bilayer spacing becomes narrower and its intensity increases whereas the peak associated with the inter-backbone planes reduces its intensity for higher temperatures. This last result Could be interpreted as a decrease in the number of aggregates and the reduction of the interchain species during the MEH-PPV relaxation processes. These WAXS results were correlated with PL spectra modifications observed upon temperature treatments. (C) 2008 Elsevier Ltd. All rights reserved.

Combined Crystallographic and Solution Molecular Dynamics Study of Allosteric Effects in Ester and Ketone p-tert-Butylcalix[4]arene Derivatives and Their Complexes with Acetonitrile, Cd(II), and Pb(II)

We describe here a procedure to bridge the gap in the field of calixarene physicochemistry between solid-state atomic-resolution structural information and the liquid-state low-resolution thermodynamics and spectroscopic data. We use MD simulations to study the kinetics and energetics involved in the complexation of lower rim calix[4]arene derivatives (L), containing bidentate ester (1) and ketone (2) pendant groups, with acetonitrile molecule (MeCN) and Cd2+ and Pb2+ ions (M2+) in acetonitrile solution. On one hand, we found that the prior inclusion of MeCN into the calix to form a L(MeCN) adduct has only a weak effect in preorganizing the hydrophilic cavity toward metal ion binding. On the other hand, the strong ion-hydrophilic cavity interaction produces a wide open calix which enhances the binding of one MeCN molecule (allosteric effect) to stabilize the whole (M2+)1(MeCN) bifunctional complex. We reach two major conclusions: (i) the MD results for the (M2+)1(MeCN) binding are in close agreement with the ""endo"", fully encapsulated, metal complex found by X-ray diffraction and in vacuo MD calculations, and (ii) the MD structure for the more flexible 2 ligand, however, differs from the also endo solid-state molecule. In fact, it shows strong solvation effects at the calixarene lower bore by competing MeCN molecules that share the metal coordination sphere with the four C=O oxygens of an ""exo"" (M2+)2(MeCN) complex.

Only subtle protein conformational adaptations are required for ligand binding to thyroid hormone receptors: Simulations using a novel multipoint steered molecular dynamics approach

Thyroid hormone receptors (TR) are hormone-dependent transcription regulators that play a major role in human health, development, and metabolic functions. The thyroid hormone resistance syndrome, diabetes, obesity, and some types of cancer are just a few examples of important diseases that are related to TR malfunctioning, particularly impaired hormone binding. Ligand binding to and dissociation from the receptor ultimately control gene transcription and, thus, detailed knowledge of binding and release mechanisms are fundamental for the comprehension of the receptor`s biological function and development of pharmaceuticals. In this work, we present the first computational study of ligand entry into the ligand binding domain (LBD) of a nuclear receptor. We report molecular dynamics simulations of ligand binding to TRs using a generalization of the steered molecular dynamics technique designed to perform single-molecule pulling simulations along arbitrarily nonlinear driving pathways. We show that only gentle protein movements and conformational adaptations are required for ligand entry into the LBDs and that the magnitude of the forces applied to assist ligand binding are of the order of the forces involved in ligand dissociation. Our simulations suggest an alternative view for the mechanisms ligand binding and dissociation of ligands from nuclear receptors in which ligands can simply diffuse through the protein surface to reach proper positioning within the binding pocket. The proposed picture indicates that the large-amplitude protein motions suggested by the apo- and holo-RXR alpha crystallographic structures are not required, reconciling conformational changes of LBDs required for ligand entry with other nuclear receptors apo-structures that resemble the ligand-bound LBDs.

Ligand dissociation from estrogen receptor is mediated by receptor dimerization: Evidence from molecular dynamics simulations

Estrogen Receptor (ER) is an important target for pharmaceutical design. Like other ligand-dependent transcription factors, hormone binding regulates ER transcriptional activity. Nevertheless, the mechanisms by which ligands enter and leave ERs and other nuclear receptors remain poorly understood. Here, we report results of locally enhanced sampling molecular dynamics simulations to identify dissociation pathways of two ER ligands [the natural hormone 17 beta-estradiol (E-2) and the selective ER modulator raloxifene (RAL)] from the human ER alpha ligand-binding domain in monomeric and dimeric forms. E-2 dissociation occurs via three different pathways in ER monomers. One resembles the mousetrap mechanism (Path I), involving repositioning of helix 12 (H12), others involve the separation of H8 and H11 (Path II), and a variant of this pathway at the bottom of the ligand-binding domain (Path II`). RAL leaves the receptor through Path I and a Path I variant in which the ligand leaves the receptor through the loop region between H11 and H12 (Path I`). Remarkably, ER dimerization strongly suppresses Paths II and II` for E-2 dissociation and modifies RAL escape routes. We propose that differences in ligand release pathways detected in the simulations for ER monomers and dimers provide an explanation for previously observed effects of ER quaternary state on ligand dissociation rates and suggest that dimerization may play an important, and hitherto unexpected, role in regulation of ligand dissociation rates throughout the nuclear receptor family.

The effects of anticalcification treatments and hydration on the molecular dynamics of bovine pericardium collagen as revealed by (13)C solid-state NMR

This article describes a solid-state NMR (SSNMR) investigation of the influence of hydration and chemical cross-linking on the molecular dynamics of the constituents of the bovine pericardium (BP) tissues and its relation to the mechanical properties of the tissue. Samples of natural phenetylamine-diepoxide (DE)- and glutaraldehyde (GL)-fixed BP were investigated by (13)C cross-polarization SSNMR to probe the dynamics of the collagen, and the results were correlated to the mechanical properties of the tissues, probed by dynamical mechanical analysis. For samples of natural BP, the NMR results show that the higher the hydration level the more pronounced the molecular dynamics of the collagen backbone and sidechains, decreasing the tissue`s elastic modulus. In contrast, in DE- and GL-treated samples, the collagen molecules are more rigid, and the hydration seems to be less effective in increasing the collagen molecular dynamics and reducing the mechanical strength of the samples. This is mostly attributed to the presence of cross-links between the collagen plates, which renders the collagen mobility less dependent on the water absorption in chemically treated samples. Copyright (C) 2010 John Wiley & Sons, Ltd.

Prigogine-Defay Ratio for an Ionic Glass-Former: Molecular Dynamics Simulations

The pressure dependence of the glass-transition temperature, T(g)(P), of the ionic glass-former 2Ca(NO(3))(2) center dot 3KNO(3), CKN, has been obtained by molecular dynamics (MD) simulations The liquid-glass difference of thermal expansivity, Delta alpha, heat capacity, Delta C(p), and isothermal compressibility, Delta kappa, have been calculated as a function of pressure. It has been found that the Ehrenfest relation dT(g)/dP = TV Delta alpha/Delta C(p) predicts the pressure dependence of T, but the other Ehrenfest relation, dT(g)/dP = Delta kappa/Delta alpha, does not. Consequently, the Prigogine-Defay ratio, Pi = Delta C(p)Delta kappa/TV Delta alpha(2), is Pi similar to 1.2 at low pressures, but increases 1 order of magnitude at high pressures. The pressure dependence of the Prigogine-Defay ratio is interpreted in light of recent explanations for the finding Pi > 1.

Docking and molecular dynamics simulation of quinone compounds with trypanocidal activity

In this work, two different docking programs were used, AutoDock and FlexX, which use different types of scoring functions and searching methods. The docking poses of all quinone compounds studied stayed in the same region in the trypanothione reductase. This region is a hydrophobic pocket near to Phe396, Pro398 and Leu399 amino acid residues. The compounds studied displays a higher affinity in trypanothione reductase (TR) than glutathione reductase (GR), since only two out of 28 quinone compounds presented more favorable docking energy in the site of human enzyme. The interaction of quinone compounds with the TR enzyme is in agreement with other studies, which showed different binding sites from the ones formed by cysteines 52 and 58. To verify the results obtained by docking, we carried out a molecular dynamics simulation with the compounds that presented the highest and lowest docking energies. The results showed that the root mean square deviation (RMSD) between the initial and final pose were very small. In addition, the hydrogen bond pattern was conserved along the simulation. In the parasite enzyme, the amino acid residues Leu399, Met400 and Lys402 are replaced in the human enzyme by Met406, Tyr407 and Ala409, respectively. In view of the fact that Leu399 is an amino acid of the Z site, this difference could be explored to design selective inhibitors of TR.

Tayloring Brazil: a system dynamics model for monetary policy feedback

The thesis introduces a system dynamics Taylor rule model of new Keynesian nature for monetary policy feedback in Brazil. The nonlinear Taylor rule for interest rate changes con-siders gaps and dynamics of GDP growth and inflation. The model closely tracks the 2004 to 2011 business cycle and outlines the endogenous feedback between the real interest rate, GDP growth and inflation. The model identifies a high degree of endogenous feedback for monetary policy and inflation, while GDP growth remains highly exposed to exogenous eco-nomic conditions. The results also show that the majority of the monetary policy moves during the sample period was related to GDP growth, despite higher coefficients of inflation parameters in the Taylor rule. This observation challenges the intuition that inflation target-ing leads to a dominance of monetary policy moves with respect to inflation. Furthermore, the results suggest that backward looking price-setting with respect to GDP growth has been the dominant driver of inflation. Moreover, simulation exercises highlight the effects of the new BCB strategy initiated in August 2011 and also consider recession and inflation avoid-ance versions of the Taylor rule. In methodological terms, the Taylor rule model highlights the advantages of system dynamics with respect to nonlinear policies and to the stock-and-flow approach. In total, the strong historical fit and some counterintuitive observations of the Taylor rule model call for an application of the model to other economies.

Preferential location of prilocaine and etidocaine in phospholipid bilayers: A molecular dynamics study

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Molecular dynamics simulation on devitrification: Isothermal devitrification and thermodynamics of PbF2 glasses

The vitrification and devitrification features of lead fluoride are investigated by means of molecular dynamic simulations. The influence of heating rate on the devitrification temperature as well as the dependence of the glass properties on its thermal history, i.e., the cooling rate employed, is identified. As expected, different glasses are obtained when the cooling rates differ. Diffusion coefficient analysis during heating of glass and crystal, indicates that the presence of defects on the glassy matrix favors the transition processes from the ionic to a superionic state, with high mobility of fluorine atoms, responsible for the high anionic conduction of lead fluoride. Nonisothermal and isothermal devitrification processes are simulated in glasses obtained at different cooling rates and structural organizations occurring during the heat treatments are clearly observed. When a fast cooling rate is employed during the glass formation, the devitrification of a single crystal (limited by the cell dimensions) is observed, while the glass obtained with slower cooling rate, allowing relaxations and organization of various regions on the glass bulk during the cooling process, devitrifies in more than one crystalline plane. (C) 2004 American Institute of Physics.

Molecular dynamics simulations on devitrification: the PbF2 case

In this work molecular dynamics simulations were performed to reproduce the kinetic and thermodynamic transformations occurring during melt crystallization, vitrification, and glass crystallization (devitrification) of PbF2. Two potential parameters were analyzed in order to access the possibility of modeling these properties. These interionic potentials are models developed to describe specific characteristic of PbF2, and thermodynamic properties were well reproduced by one of them, while the other proved well adapted to simulate the crystalline structure of this fluoride. By a modeled nonisothermal heat treatment of the glass, it was shown that the devitrification of a cubic structure in which the Pb-Pb distances are in good agreement with theory and experiment. (C) 2002 American Institute of Physics.

On the unzipping mechanisms of carbon nanotubes: Insights from reactive molecular dynamics simulations

Unzipping carbon nanotubes (CNTs) is considered one of the most promising approaches for the controlled and large-scale production of graphene nanoribbons (GNR). These structures are considered of great importance for the development of nanoelectronics because of its dimensions and intrinsic nonzero band gap value. Despite many years of investigations some details on the dynamics of the CNT fracture/unzipping processes remain unclear. In this work we have investigated some of these process through molecular dynamics simulations using reactive force fields (ReaxFF), as implemented in the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code. We considered multi-walled CNTs of different dimensions and chiralities and under induced mechanical stretching. Our preliminary results show that the unzipping mechanisms are highly dependent on CNT chirality. Well-defined and distinct fracture patterns were observed for the different chiralities. Armchair CNTs favor the creation of GNRs with well-defined armchair edges, while zigzag and chiral ones produce GNRs with less defined and defective edges. © 2012 Materials Research Society.

Controlled route to the fabrication of carbon and boron nitride nanoscrolls: A molecular dynamics investigation

Carbon nanoscrolls (graphene layers rolled up into papyrus-like tubular structures) are nanostructures with unique and interesting characteristics that could be exploited to build several new nanodevices. However, an efficient and controlled synthesis of these structures was not achieved yet, making its large scale production a challenge to materials scientists. Also, the formation process and detailed mechanisms that occur during its synthesis are not completely known. In this work, using fully atomistic molecular dynamics simulations, we discuss a possible route to nanoscrolls made from graphene layers deposited over silicon oxide substrates containing chambers/pits. The scrolling mechanism is triggered by carbon nanotubes deposited on the layers. The process is completely general and can be used to produce scrolls from other lamellar materials, like boron nitride, for instance. © 2013 American Institute of Physics.