PACKMOL: A Package for Building Initial Configurations for Molecular Dynamics Simulations

Adequate initial configurations for molecular dynamics simulations consist of arrangements of molecules distributed in space in such a way to approximately represent the system`s overall structure. In order that the simulations are not disrupted by large van der Waals repulsive interactions, atoms from different molecules Must keep safe pairwise distances. Obtaining Such a molecular arrangement can be considered it packing problem: Each type molecule must satisfy spatial constraints related to the geometry of the system, and the distance between atoms of different molecules Must be greater than some specified tolerance. We have developed a code able to pack millions of atoms. grouped in arbitrarily complex molecules, inside a variety of three-dimensional regions. The regions may be intersections of spheres, ellipses, cylinders, planes, or boxes. The user must provide only the structure of one molecule of each type and the geometrical constraints that each type of molecule must satisfy. Building complex mixtures, interfaces, solvating biomolecules in water, other solvents, or mixtures of solvents, is straight forward. In addition. different atoms belonging to the same molecule may also be restricted to different spatial regions, in Such a way that more ordered molecular arrangements call be built, as micelles. lipid double-layers, etc. The packing time for state-of-the-art molecular dynamics systems varies front a few seconds to a few Minutes in a personal Computer. The input files are simple and Currently compatible with PDB, Tinker, Molden, or Moldy coordinate files. The package is distributed as free software and call be downloaded front http://www.ime.unicamp.br/similar to martinez/packmol/. (C) 2009 Wiley Periodicals. Inc. J Comput Chem 30: 2157-2164, 2009

SUPPORT OF MAXIMIZING MEASURES FOR TYPICAL C(0) DYNAMICS ON COMPACT MANIFOLDS

Given a compact manifold X, a continuous function g : X -> IR, and a map T : X -> X, we study properties of the T-invariant Borel probability measures that maximize the integral of g. We show that if X is a n-dimensional connected Riemaniann manifold, with n >= 2, then the set of homeomorphisms for which there is a maximizing measure supported on a periodic orbit is meager. We also show that, if X is the circle, then the ""topological size"" of the set of endomorphisms for which there are g maximizing measures with support on a periodic orbit depends on properties of the function g. In particular, if g is C(1), it has interior points.

A DISCRETIZATION SCHEME FOR AN ONE-DIMENSIONAL REACTION-DIFFUSION EQUATION WITH DELAY AND ITS DYNAMICS

The goal of this paper is to present an approximation scheme for a reaction-diffusion equation with finite delay, which has been used as a model to study the evolution of a population with density distribution u, in such a way that the resulting finite dimensional ordinary differential system contains the same asymptotic dynamics as the reaction-diffusion equation.

Structure, Dynamics, and RNA Interaction Analysis of the Human SBDS Protein

Shwachman-Bodian-Diamond syndrome is an autosomal recessive genetic syndrome with pleiotropic phenotypes, including pancreatic deficiencies, bone marrow dysfunctions with increased risk of myelodysplasia or leukemia, and skeletal abnormalities. This syndrome has been associated with mutations in the SBDS gene, which encodes a conserved protein showing orthologs in Archaea and eukaryotes. The Shwachman-Bodian-Diamond syndrome pleiotropic phenotypes may be an indication of different cell type requirements for a fully functional SBDS protein. RNA-binding activity has been predicted for archaeal and yeast SBDS orthologs, with the latter also being implicated in ribosome biogenesis. However, full-length SBDS orthologs function in a species-specific manner, indicating that the knowledge obtained from model systems may be of limited use in understanding major unresolved issues regarding SBDS function, namely, the effect of mutations in human SBDS on its biochemical function and the specificity of RNA interaction. We determined the solution structure and backbone dynamics of the human SBDS protein and describe its RNA binding site using NMR spectroscopy. Similarly to the crystal structures of Archaea, the overall structure of human SBDS comprises three well-folded domains. However, significant conformational exchange was observed in NMR dynamics experiments for the flexible linker between the N-terminal domain and the central domain, and these experiments also reflect the relative motions of the domains. RNA titrations monitored by heteronuclear correlation experiments and chemical shift mapping analysis identified a classic RNA binding site at the N-terminal FYSH (fungal, Yhr087wp, Shwachman) domain that concentrates most of the mutations described for the human SBDS. (C) 2010 Elsevier Ltd. All rights reserved.

Picosecond Dynamics of the Prototropic Reactions of 7-Hydroxyflavylium Photoacids Anchored at an Anionic Micellar Surface

Three water-insoluble, micelle-anchored flavylium salts, 7-hydroxy-3-octyl-flavylium chloride, 4`-hexyl-7-hydroxyflavylium chloride, and 6-hexyl-7-hydroxy-4-methyl-flavylium chloride, have been employed to probe excited-state prototropic reactions in micellar sodium dodecyl sulfate (SDS). In SDS micelles, the fluorescence decays of these three flavylium salts are tetraexponential functions in the pH range from 1.0 to 4.6 at temperatures from 293 to 318 K. The four components of the decays are assigned to Four kinetically coupled excited species in the micelle: specifically, promptly deprotonable (AH(+)*) and nonpromptly deprotonable (AH(h)(+)*) orientations of the acid in the micelle. the base-proton geminate pair (A*center dot center dot center dot H(+)), and the free conjugate base (A*). The initial prompt deprotonation to form the germinate pair occurs at essentially the same rate (k(d) similar to 6-7 x 10(10) s(-1)) for all three photoacids. Recombination of the germinate pair is similar to 3-fold faster than the rate of proton escape from the pair (k(rec) similar to 3 x 10(10) s(-1) and k(diss) similar to 1 x 10(10) s(-1)), corresponding to an intrinsic recombination efficiency of the pair of similar to 75%. Finally, the reprotonation of the short-lived free A* (200-350 ps, depending oil the photoacid) has two components, only one of which depends oil the proton concentration in the intermicellar aqueous phase. Ultrafast transfer of the proton to water and substantial compartmentalization of the photogenerated proton at the micelle surface Oil the picosecond time scale strongly suggest preferential transfer of the proton to preformed hydrogen-bonded water bridges between the photoacid and the anionic headgroups. This localizes the proton in the vicinity of the excited base much more efficiently than ill bulk water, resulting ill the predominance of geminate re reprotonation at the micelle surface.

Picosecond Dynamics of Proton Transfer of a 7-Hydroxyflavylium Salt in Aqueous-Organic Solvent Mixtures

The intermediacy of the geminate base proton pair (A*center dot center dot center dot H(+)) in excited-state proton-transfer (ESPT) reactions (two-step mechanism) has been investigated employing the synthetic flavylium salt 7-hydroxy-4-methyl-flavylium chloride (HMF). In aqueous solution, the ESPT mechanism involves solely the excited acid AH* and base A* forms of HMF as indicated by the fluorescence spectra and double-exponential fluorescence decays (two species, two decay times). However, upon addition of either 1,4-dioxane or 1,2-propylene glycol, the decays become triple-exponential with a term consistent with the presence of the geminate base proton pair A*center dot center dot center dot H(+). The geminate pair becomes detectable because of the increase in the recombination rate constant, k(rec), of (A*center dot center dot center dot H(+)) with increasing the mole fraction of added organic cosolvent. Because the two-step ESPT mechanism splits the intrinsic prototropic reaction rates (deprotonation of AH(+)*, k(d), and recombination, k(rec) of A*center dot center dot center dot H(+)) from the diffusion controlled rates (dissociation, k(diss) and formation, k(diff)[H(+)], of A*center dot center dot center dot H+), the experimental detection of the geminate pair provides a wealth of information on the proton-transfer reaction (k(d) and k(rec)) as well as on proton diffusion/migration (k(diss) and k(diff)).

The shoving model for the glass-former LiCl center dot 6H(2)O: A molecular dynamics simulation study

Molecular dynamics (MD) simulations of LiCl center dot 6H(2)O Showed that the diffusion coefficient D, and also I lie structural relaxation time , follow a power law at high temperatures, D(-1) proportional to (T - T(0))(-mu), with the same experimental parameters for viscosity (T(0) = 207 K, mu = 2.08). Decoupling between D and occurs at T(x) similar to 1.1 T(0). High frequency acoustic excitations for the LiCl center dot 6H(2)O model were obtained by the calculation of time correlation functions of mass current fluctuations. The temperature dependence of the instantaneous shear modulus, G,(T), was considered in the shoving model for supercooled liquids [J.C. Dyre, T. Christensen, N.B. Olsen, J. Non-Cryst. Solids 352 (2006) 4635] resulting in a linear relationship log (D(-1)) vs. G root T. (C) 2009 Elsevier B.V. All rights reserved.

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.

Alkoxy Chain Effect on the Viscosity of a Quaternary Ammonium Ionic Liquid: Molecular Dynamics Simulations

The viscosity of ionic liquids based on quaternary ammonium cations is reduced when one of the alkyl chains is replaced by an alkoxy chain (Zhou et al. Chem. Eur. J. 2005, 11, 752.). A microscopic picture of the role played by the ether function in decreasing the viscosity of quaternary ammonium ionic liquids is provided here by molecular dynamics (MD) simulations. A model for the ionic liquid N-ethyl-N,N-dimethyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide, MOENM(2)E TFSI, is compared to the tetraalky-lammonium counterpart. The alkoxy derivative has lower viscosity, higher ionic diffusion coefficients, and higher conductivity than the tetraalkyl system at the same density and temperature. A clear signature of the ether function on the liquid structure is observed in cation-cation correlations, but not in anion-anion or anion-cation correlations. In both the alkyl and the alkoxy ionic liquids, there is aggregation of long chains of neighboring cations within micelle-like structures. The MD simulations indicate that the less effective assembly between the more flexible alkoxy chains, in comparison to alkyl chains, is the structural reason for higher ionic mobility in MOENM(2)E TFSI.

Interfacial Electron Transfer Dynamics Following Laser Flash Photolysis of [Ru(bpy)(2)((4,4 `-PO(3)H(2))(2)bpy)](2+) in TiO(2) Nanoparticle Films in Aqueous Environments

Nanosecond laser flash photolysis has been used to investigate injection and back electron transfer from the complex [(Ru-(bpy)(2)(4,4`-(PO(3)H(2))(2)bpy)](2+) surface-bound to TiO(2) (TiO(2)-Ru(II)). The measurements were conducted under conditions appropriate for water oxidation catalysis by known single-site water oxidation catalysts. Systematic variations in average lifetimes for back electron transfer, - were observed with changes in pH, surface coverage, incident excitation intensity, and applied bias. The results were qualitatively consistent with a model involving rate-limiting thermal activation of injected electrons from trap sites to the conduction band or shallow trap sites followed by site-to-site hopping and interfacial electron transfer, TiO(2)(e(-))-Ru(3+) -> TiO(2)-Ru(2+). The appearance of pH-dependent decreases in the efficiency of formation of TiO(2)-Ru(3+) and in incident-photon-to-current efficiencies with the added reductive scavenger hydroquinone point to pH-dependent back electron transfer processes on both the sub-nanosecond and millisecond-microsecond time scales, which could be significant in limiting long-term storage of multiple redox equivalents.

Excited-State Dynamics in fac-[Re(CO)(3)(Me(4)phen)(L)](+)

Excited-state dynamics in fac-[Re(CO)(3)(Me(4)phen)(cis-L)](+) (Me(4)phen = 3,4,7,8-tetramethyl-1,10-phenanthroline, L = 4-styrylpyridine (stpy) or 1,2-bis(4-pyridyl)ethylene (bpe)) were investigated by steady-state and time-resolved techniques. A complex equilibrium among three closely lying excited states, 3IL(cis-L), (3)MLCT(Re -> me4phen), and (3)IL(Me4phen), has been established. Under UV irradiation, cis-to-trans isomerization of coordinated cis-L is observed with a quantum yield of 0.15 in acetonitrile solutions. This photoreaction competes with radiative decay from (3)MLCT(Re -> Me4phen) and (3)IL(Me4phen) excited states, leading to a decrease in the emission quantum yield relative to the nonisomerizable complex fac-[Re(CO)(3)(Me(4)phen)(bpa)](+) (bpa = 1,2-bis(4-pyridyl)ethane). From temperature-dependent time-resolved emission measurements in solution and in poly(methyl methacrylate) (PMMA) films, energy barriers (Delta E(a)) for interconversion between (3)MLCT(Re -> me4Phen) and (3)IL(Me4phen) emitting states were determined. For L = cis-stpy, Delta E(a) = 11 (920 cm(-1)) and 15 kJ mol(-1) (1254 cm(-1)) in 5:4 propionitrile/butyronitrile and PMMA, respectively. For L = cis-bpe, Delta E(a) = 13 kJ mol(-1) (1087 cm(-1)) in 5:4 propionitrile/butyronitrile. These energy barriers are sufficient to decrease the rate constant for internal conversion from higher-lying (3)IL(me4phen) state to (3)MLCT(Re -> Me4phen), k(i) congruent to 10(6) s(-1). The decrease in rate allows for the observation of intraligand phosphorescence, even in fluid medium at room temperature. Our results provide additional insight into the role of energy gap and excited-state dynamics on the photochemical and photophysical properties of Re(I) polypyridyl complexes.

Excited-state dynamics of polyfluorene derivatives in solution

The excited-state dynamics of two polyfluorene copolymers, one fully conjugated containing phenylene vinylene units alternated with 9,9`-dihexylfluorenyl groups and the other segmented by -(CH2)(8)- spacer, were studied in dilute solution of different solvents using a picosecond single-photon timing technique. The excited-state dynamics of the segmented copolymer follows the Forster resonant energy-transfer model which describes intrachain energy-transfer kinetics among random oriented chromophores. Energy transfer is confirmed by analysis of fluorescence anisotropy relaxation with the measurement of a short decay component of about 60 ps. The fluorescence decay surface of the fully conjugated copolymer is biexponential with decay times of about 470 and 900 ps, ascribed to deactivation of chain moieties containing trans and cis isomers already in a photostationary condition. Thus, energy transfer is very fast due to the conjugated nature and rigid-rod-like structure of this copolymer chain.

On the interaction of bovine serum albumin with ionic surfactants: Temperature induced EPR changes of a maleimide nitroxide reflect local protein dynamics and probe solvent accessibility

The interaction of bovine serum albumin (BSA) with the ionic surfactants sodium dodecylsulfate (SDS, anionic), cetyltrimethylammonium chloride (CTAC, cationic) and N-hexadecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (HPS, zwitterionic) was studied by electron paramagnetic resonance (EPR) spectroscopy of spin label covalently bound to the single free thiol group of the protein. EPR spectra simulation allows to monitor the protein dynamics at the labeling site and to estimate the changes in standard Gibbs free energy, enthalpy and entropy for transferring the nitroxide side chain from the more motionally restricted to the less restricted component. Whereas SDS and CTAC showed similar increases in the dynamics of the protein backbone for all measured concentrations. HPS presented a smaller effect at concentrations above 1.5 mM. At 10 mM of surfactants and 0.15 mM BSA, the standard Gibbs free energy change was consistent with protein backbone conformations more expanded and exposed to the solvent as compared to the native protein, but with a less pronounced effect for HPS. In the presence of the surfactants, the enthalpy change, related to the energy required to dissociate the nitroxide side chain from the protein, was greater, suggesting a lower water activity. The nitroxide side chain also detected a higher viscosity environment in the vicinity of the paramagnetic probe induced by the addition of the surfactants. The results suggest that the surfactant-BSA interaction, at higher surfactant concentration, is affected by the affinities of the surfactant to its own micelles and micelle-like aggregates. Complementary DLS data suggests that the temperature induced changes monitored by the nitroxide probe reflects local changes in the vicinity of the single thiol group of Cys-34 BSA residue. (C) 2011 Elsevier B.V. All rights reserved.

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.

Towards Presidential Rule in Ukraine : Hybrid Regime Dynamics under Semi-Presidentialism

This article sets out to analyse recent regime developments in Ukraine in relation to semi-presidentialism. The article asks: to what extent and in what ways theoretical arguments against semi-presidentialism (premier-presidential and president-parliamentary systems) are relevant for understanding the changing directions of the Ukrainian regime since the 1990s? The article also reviews the by now overwhelming evidence suggesting that President Yanukovych is turning Ukraine into a more authoritarian hybrid regime and raises the question to what extent the president-parliamentary system might serve this end. The article argues that both kinds of semi-presidentialism have, in different ways, exacerbated rather than mitigated institutional conflict and political stalemate. The return to the president-parliamentary system in 2010 – the constitutional arrangement with the most dismal record of democratisation – was a step in the wrong direction. The premier-presidential regime was by no means ideal, but it had at least two advantages. It weakened the presidential dominance and it explicitly anchored the survival of the government in parliament. The return to the 1996 constitution ties in well with the notion that President Viktor Yanukovych has embarked on an outright authoritarian path.