Flexibility and inhibitor binding in Cdc25 phosphatases

Cdc25 phosphatases involved in cell cycle checkpoints are now active targets for the development of anti-cancer therapies. Rational drug design would certainly benefit from detailed structural information for Cdc25s. However, only apo- or sulfate-bound crystal structures of the Cdc25 catalytic domain have been described so far. Together with previously available crystalographic data, results from molecular dynamics simulations, bioinformatic analysis, and computer-generated conformational ensembles shown here indicate that the last 30-40 residues in the C-terminus of Cdc25B are partially unfolded or disordered in solution. The effect of C-terminal flexibility upon binding of two potent small molecule inhibitors to Cdc25B is then analyzed by using three structural models with variable levels of flexibility, including an equilibrium distributed ensemble of Cdc25B backbone conformations. The three Cdc25B structural models are used in combination with flexible docking, clustering, and calculation of binding free energies by the linear interaction energy approximation to construct and validate Cdc25B-inhibitor complexes. Two binding sites are identified on top and beside the Cdc25B active site. The diversity of interaction modes found increases with receptor flexibility. Backbone flexibility allows the formation of transient cavities or compact hydrophobic units on the surface of the stable, folded protein core that are unexposed or unavailable for ligand binding in rigid and densely packed crystal structures. The present results may help to speculate on the mechanisms of small molecule complexation to partially unfolded or locally disordered proteins.

Effect of SO(2) on the Transport Properties of an Imidazolium Ionic Liquid and Its Lithium Solution

Transport coefficients have been measured as a function of the concentration of sulfur dioxide, SO(2), dissolved in 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)-imide, [BMMI][Tf(2)N], as well as in its lithium salt solution, Li[Tf(2)N]. The SO(2) reduces viscosity and density and increases conductivity and diffusion coefficients in both the neat [BMMI] [Tf(2)N] and the [BMMI][Tf(2)N]-Li[Tf(2)N] solution. The conductivity enhancement is not assigned to a simple viscosity effect; the weakening of ionic interactions upon SO(2) addition also plays a role. Microscopic details of the SO(2) effect were unraveled using Raman spectroscopy and molecular dynamics (MD) simulations. The Raman spectra suggest that the Li(+)-[Tf(2)N] interaction is barely affected by SO(2), and the SO(2)-[Tf(2)N] interaction is weaker than previously observed in an investigation of an ionic liquid containing the bromide anion. Transport coefficients calculated by MD simulations show the same trend as the experimental data with respect to SO(2) content. The MD simulations provide structural information on SO(2) molecules around [Tf(2)N], in particular the interaction of the sulfur atom of SO(2) with oxygen and fluorine atoms of the anion. The SO(2)-[BMMI] interaction is also important because the [BMMI] cations with above-average mobility have a larger number of nearest-neighbor SO(2) molecules.

Unraveling Dynamical Heterogeneity in the Ionic Liquid 1-Butyl-3-methylimidazolium Chloride

Heterogeneous dynamics within a time range of nanoseconds was investigated by molecular dynamics (MD) simulations of 1 -butyl-3-methylimidazolium chloride ([bmim]Cl). After identifying groups of fast and slow ions, it was shown that the separation between the location of the center of mass and the center of charge of cations, d(CMCC), is a signature of such difference in ionic mobility. The distance d(CMCC) can be used as a signature because it relaxes in the time window of the dynamical heterogeneity. The relationship between the parameter dcmcc and conformations of the side alkyl chain in [bmim] is discussed. Since the relatively slow relaxation of dcmcc is a consequence of coexisting polar and nonpolar domains in the bulk, the MD simulations reveal a subtle interplay between structural and dynamical heterogeneity in ionic liquids.

Host-guest system of 4-nerolidylcatechol in 2-hydroxypropyl-beta-cyclodextrin: preparation, characterization and molecular modeling

The interaction of 4-nerolidylcatechol (4-NRC), a potent antioxidant agent, and 2-hydroxypropyl-beta-cyclodextrin (HP-beta-CD) was investigated by the solubility method using Fourier transform infrared (FTIR) methods in addition to UV-Vis, (1)H-nuclear magnetic resonance (NMR) spectroscopy and molecular modeling. The inclusion complexes were prepared using grinding, kneading and freeze-drying methods. According to phase solubility studies in water a B(S)-type diagram was found, displaying a stoichiometry complexation of 2:1 (drug:host) and stability constant of 6494 +/- A 837 M(-1). Stoichiometry was established by the UV spectrophotometer using Job`s plot method and, also confirmed by molecular modeling. Data from (1)H-NMR, and FTIR, experiments also provided formation evidence of an inclusion complex between 4-NRC and HP-beta-CD. 4-NRC complexation indeed led to higher drug solubility and stability which could probably be useful to improve its biological properties and make it available to oral administration and topical formulations.

Shielding of ionic interactions by sulfur dioxide in an ionic liquid

The effect of adding SO(2) on the structure and dynamics of 1-butyl-3-methylimidazolium bromide (BMIBr) was investigated by low-frequency Raman spectroscopy and molecular dynamics (MD) simulations. The MD simulations indicate that the long-range structure of neat BMIBr is disrupted resulting in a liquid with relatively low viscosity and high conductivity, but strong correlation of ionic motion persists in the BMIBr-SO(2) mixture due to ionic pairing. Raman spectra within the 5 < omega < 200 cm(-1) range at low temperature reveal the short-time dynamics, which is consistent with the vibrational density of states calculated by MD simulations. Several time correlation functions calculated by MD simulations give further insights on the structural relaxation of BMIBr-SO(2).

Transport coefficients, Raman spectroscopy, and computer simulation of lithium salt solutions in an ionic liquid

Lithium salt solutions of Li(CF3SO2)(2)N, LiTFSI, in a room-temperature ionic liquid (RTIL), 1-butyl-2,3-dimethyl-imidazolium cation, BMMI, and the (CF3SO2)(2)N-, bis(trifluoromethanesulfonyl)imide anion, [BMMI][TFSI], were prepared in different concentrations. Thermal properties, density, viscosity, ionic conductivity, and self-diffusion coefficients were determined at different temperatures for pure [BMMI][TFSI] and the lithium solutions. Raman spectroscopy measurements and computer simulations were also carried out in order to understand the microscopic origin of the observed changes in transport coefficients. Slopes of Walden plots for conductivity and fluidity, and the ratio between the actual conductivity and the Nernst-Einstein estimate for conductivity, decrease with increasing LiTFSI content. All of these studies indicated the formation of aggregates of different chemical nature, as it is corroborated by the Raman spectra. In addition, molecular dynamics (MD) simulations showed that the coordination of Li+ by oxygen atoms of TFSI anions changes with Li+ concentration producing a remarkable change of the RTIL structure with a concomitant reduction of diffusion coefficients of all species in the solutions.

Raman spectra of polymer electrolytes based on poly(ethylene glycol) dimethyl ether, lithium perchlorate, and the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate

Raman spectra of polymer electrolytes based on poly(ethylene glycol) dimethyl ether (PEGdME) with LiClO(4), PEGdME/LiClO(4), and the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate, PEGdME/[bmim]PF(6), are compared. Raman spectroscopy suggests stronger interactions in PEGdME/LiClO(4) than PEGdmE/[bmim]PF(6), thus corroborating previous results obtained by molecular dynamics simulations. Quantum Chemistry methods have been used to calculate vibrational frequencies and the equilibrium structure of segments of the polymer chain around the cation. A consistent picture has been obtained from Raman spectroscopy, density functional theory (DFT) calculations, and molecular dynamics simulations for these polymer electrolytes. (C) 2010 Elsevier B.V. All rights reserved.

Effect of the Y(2)O(3)-ZrO(2) support composition on nickel catalyst evaluated in dry reforming of methane

Nickel catalysts with a load of 5 wt.% Ni, supported on pure ZrO(2) and ZrO(2) stabilized with 4 mol%, 8 mol% and 12 mol% of Y(2)O(3), were prepared by the polymerization method. The samples were characterized by X-ray diffraction (XRD), temperature-programmed reduction with hydrogen (TPR-H(2)), specific surface area (BET) and electronic paramagnetic resonance (EPR) and tested as catalysts for carbon dioxide reforming of methane. The XRD patterns showed the presence of the oxide precursor (NiO) and the tetragonal phase of a Y(2)O(3)-ZrO(2) solid solution. According to the TPR-H(2) analysis, the reduction of various NiO species was influenced by the composition of the support. Catalytic tests were conducted at 800 degrees C for 6 h, and the composition of the gaseous products and the catalytic conversion rate depended on the composition of the Y(2)O(3)-ZrO(2) solid solution and its influence on the supported NiO species. A direct relation was observed between the variation in the support, the nickel species supported on it and the performance in the catalytic tests. (C) 2008 Elsevier B.V. All rights reserved.

Modelling water adsorption on Au(210) surfaces: II. Monte Carlo simulations

Canonical Monte Carlo simulations for the Au(210)/H(2)O interface, using a force field recently proposed by us, are reported. The results exhibit the main features normally observed in simulations of water molecules in contact with different noble metal surfaces. The calculations also assess the influence of the surface topography on the structural aspects of the adsorbed water and on the distribution of the water molecules in the direction normal to the metal surface plane. The adsorption process is preferential at sites in the first layer of the metal. The analysis of the density profiles and dipole moment distributions points to two predominant orientations. Most of the molecules are adsorbed with the molecular plane parallel to surface, while others adsorb with one of the O-H bonds parallel to the surface and the other bond pointing towards the bulk liquid phase. There is also evidence of hydrogen bond formation between the first and second solvent layers at the interface. (c) 2007 Elsevier B.V. All rights reserved.

Quasiharmonic elastic constants corrected for deviatoric thermal stresses

The quasiharmonic approximation (QHA), in its simplest form also called the statically constrained (SC) QHA, has been shown to be a straightforward method to compute thermoelastic properties of crystals. Recently we showed that for noncubic solids SC-QHA calculations develop deviatoric thermal stresses at high temperatures. Relaxation of these stresses leads to a series of corrections to the free energy that may be taken to any desired order, up to self-consistency. Here we show how to correct the elastic constants obtained using the SC-QHA. We exemplify the procedure by correcting to first order the elastic constants of MgSiO(3) perovskite and MgSiO(3) postperovskite, the major phases of the Earth's lower mantle. We show that this first-order correction is quite satisfactory for obtaining the aggregated elastic averages of these minerals and their velocities in the lower mantle. This type of correction is also shown to be applicable to experimental measurements of elastic constants in situations where deviatoric stresses can develop, such as in diamond-anvil cells.

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: A statistical model

We propose a statistical model to account for the gel-fluid anomalous phase transitions in charged bilayer- or lamellae-forming ionic lipids. The model Hamiltonian comprises effective attractive interactions to describe neutral-lipid membranes as well as the effect of electrostatic repulsions of the discrete ionic charges on the lipid headgroups. The latter can be counterion dissociated (charged) or counterion associated (neutral), while the lipid acyl chains may be in gel (low-temperature or high-lateral-pressure) or fluid (high-temperature or low-lateral-pressure) states. The system is modeled as a lattice gas with two distinct particle types-each one associated, respectively, with the polar-headgroup and the acyl-chain states-which can be mapped onto an Ashkin-Teller model with the inclusion of cubic terms. The model displays a rich thermodynamic behavior in terms of the chemical potential of counterions (related to added salt concentration) and lateral pressure. In particular, we show the existence of semidissociated thermodynamic phases related to the onset of charge order in the system. This type of order stems from spatially ordered counterion association to the lipid headgroups, in which charged and neutral lipids alternate in a checkerboard-like order. Within the mean-field approximation, we predict that the acyl-chain order-disorder transition is discontinuous, with the first-order line ending at a critical point, as in the neutral case. Moreover, the charge order gives rise to continuous transitions, with the associated second-order lines joining the aforementioned first-order line at critical end points. We explore the thermodynamic behavior of some physical quantities, like the specific heat at constant lateral pressure and the degree of ionization, associated with the fraction of charged lipid headgroups.

Study of the optical and magnetic properties of pyrimidine in water combining PCM and QM/MM methodologies

The solvent effects on the low-lying absorption spectrum and on the (15)N chemical shielding of pyrimidine in water are calculated using the combined and sequential Monte Carlo simulation and quantum mechanical calculations. Special attention is devoted to the solute polarization. This is included by an iterative procedure previously developed where the solute is electrostatically equilibrated with the solvent. In addition, we verify the simple yet unexplored alternative of combining the polarizable continuum model (PCM) and the hybrid QM/MM method. We use PCM to obtain the average solute polarization and include this in the MM part of the sequential QM/MM methodology, PCM-MM/QM. These procedures are compared and further used in the discrete and the explicit solvent models. The use of the PCM polarization implemented in the MM part seems to generate a very good description of the average solute polarization leading to very good results for the n-pi* excitation energy and the (15)N nuclear chemical shield of pyrimidine in aqueous environment. The best results obtained here using the solute pyrimidine surrounded by 28 explicit water molecules embedded in the electrostatic field of the remaining 472 molecules give the statistically converged values for the low lying n-pi* absorption transition in water of 36 900 +/- 100 (PCM polarization) and 36 950 +/- 100 cm(-1) (iterative polarization), in excellent agreement among one another and with the experimental value observed with a band maximum at 36 900 cm(-1). For the nuclear shielding (15)N the corresponding gas-water chemical shift obtained using the solute pyrimidine surrounded by 9 explicit water molecules embedded in the electrostatic field of the remaining 491 molecules give the statistically converged values of 24.4 +/- 0.8 and 28.5 +/- 0.8 ppm, compared with the inferred experimental value of 19 +/- 2 ppm. Considering the simplicity of the PCM over the iterative polarization this is an important aspect and the computational savings point to the possibility of dealing with larger solute molecules. This PCM-MM/QM approach reconciles the simplicity of the PCM model with the reliability of the combined QM/MM approaches.

Hydrogen bond interactions between acetone and supercritical water

Hydrogen bond interactions between acetone and supercritical water are investigated using a combined and sequential Monte Carlo/quantum mechanics (S-MC/QM) approach. Simulation results show a dominant presence of con. gurations with one hydrogen bond for different supercritical states, indicating that this specific interaction plays an important role on the solvation properties of acetone in supercritical water. Using QM MP2/aug-cc-pVDZ the calculated average interaction energy reveals that the hydrogen-bonded acetone-water complex is energetically more stable under supercritical conditions than ambient conditions and its stability is little affected by variations of temperature and/or pressure. All average results reported here are statistically converged.

Formation of atomic carbon chains from graphene nanoribbons

The formation of one-dimensional carbon chains from graphene nanoribbons is investigated using ab initio molecular dynamics. We show under what conditions it is possible to obtain a linear atomic chain via pulling of the graphene nanoribbons. The presence of dimers composed of two-coordinated carbon atoms at the edge of the ribbons is necessary for the formation of the linear chains, otherwise there is simply the full rupture of the structure. The presence of Stone-Wales defects close to these dimers may lead to the formation of longer chains. The local atomic configuration of the suspended atoms indicates the formation of single and triple bonds, which is a characteristic of polyynes.

Contrasting LH-HH subband splitting of strained quantum wells grown along [001] and [113] directions

Contrasting responses for the temperature tuning of the electronic structure in semiconductor quantum wells are discussed for heterolayered structures grown along (001) and (113) directions. The temperature affects the strain modulation of the deformation potentials and the effective optical gap is tuned along with the intersub-band splitting in the valence band. A multiband theoretical model accounts for the characterization of the electronic structure, highlighting the main qualitative and quantitative differences between the two systems under study. The microscopic source of strain fields and the detailed mapping of their distribution are provided by a simulation using classical molecular-dynamics technics.