PLUMED 2: New feathers for an old bird

Enhancing sampling and analyzing simulations are central issues in molecular simulation. Recently, we introduced PLUMED, an open-source plug-in that provides some of the most popular molecular dynamics (MD) codes with implementations of a variety of different enhanced sampling algorithms and collective variables (CVs). The rapid changes in this field, in particular new directions in enhanced sampling and dimensionality reduction together with new hardware, require a code that is more flexible and more efficient. We therefore present PLUMED 2 here a,complete rewrite of the code in an object-oriented programming language (C++). This new version introduces greater flexibility and greater modularity, which both extends its core capabilities and makes it far easier to add new methods and CVs. It also has a simpler interface with the MD engines and provides a single software library containing both tools and core facilities. Ultimately, the new code better serves the ever-growing community of users and contributors in coping with the new challenges arising in the field.

Program summary

Program title: PLUMED 2

Catalogue identifier: AEEE_v2_0

Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEEE_v2_0.html

Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland

Licensing provisions: Yes

No. of lines in distributed program, including test data, etc.: 700646

No. of bytes in distributed program, including test data, etc.: 6618136

Distribution format: tar.gz

Programming language: ANSI-C++.

Computer: Any computer capable of running an executable produced by a C++ compiler.

Operating system: Linux operating system, Unix OSs.

Has the code been vectorized or parallelized?: Yes, parallelized using MPI.

RAM: Depends on the number of atoms, the method chosen and the collective variables used.

Classification: 3, 7.7, 23. Catalogue identifier of previous version: AEEE_v1_0.

Journal reference of previous version: Comput. Phys. Comm. 180 (2009) 1961.

External routines: GNU libmatheval, Lapack, Bias, MPI. (C) 2013 Elsevier B.V. All rights reserved.

Second-order NMR spectra at high field of common organic functional groups

The proton NMR spectra of aryl n-propyl sulfides gave rise to what may appear to be first-order proton NMR spectra. Upon oxidation to the corresponding sulfone, the spectra changed appearance dramatically and were clearly second-order. A detailed analysis of these second-order spectra, in the sulfone series, provided vicinal coupling constants which indicated that these compounds had a moderate preference for the anti-conformer, reflecting the much greater size of the sulfone over the sulfide. It also emerged, from this study, that the criterion for observing large second-order effects in the proton NMR spectra of 1,2-disubstituted ethanes was that the difference in vicinal coupling constants must be large and the difference in geminal coupling constants must be small. n-Propyl triphenylphosphonium bromide and 2-trimethylsilylethanesulfonyl chloride, and derivatives thereof, also exhibited second-order spectra, again due to the bulky substituents. Since these spectra are second-order due to magnetic nonequivalence of the nuclei in question, not chemical shifts, the proton spectra are perpetually second-order and can never be rendered first-order by using higher field NMR spectrometers.

Chlorostannate(II) Ionic Liquids: Speciation, Lewis Acidity, and Oxidative Stability

The anionic speciation of chlorostannate(II) ionic liquids, prepared by mixing 1-alkyl-3-methylimidazolium chloride and tin(II) chloride in various molar ratios, chi(SnCl2), was investigated in both solid and liquid states. The room temperature ionic liquids were investigated by Sn-119 NMR spectroscopy, X-ray photoelectron spectroscopy, and viscometry. Crystalline samples were studied using Raman spectroscopy, single-crystal X-ray crystallography, and differential scanning calorimetry. Both liquid and solid systems (crystallized from the melt) contained [SnCl3](-) in equilibrium with Cl- when chi(SnCl2) < 0.50, [SnCl3](-) in equilibrium with [Sn2Cl5](-) when chi(SnCl2) > 0.50, and only [SnCl3](-) when chi(SnCl2) = 0.50. Tin(II) chloride was found to precipitate when chi(SnCl2) > 0.63. No evidence was detected for the existence of [SnCl4](-) across the entire range of chi(SnCl2) although such anions have been reported in the literature for chlorostannate(II) organic salts crystallized from organic solvents. Furthermore, the Lewis acidity of the chlorostannate(II)-based systems, expressed by their Gutmann acceptor number, has been determined as a function of the composition, chi(SnCl2), to reveal Lewis acidity for chi(SnCl2) > 0.50 samples comparable to the analogous systems based on zinc(II). A change of the Lewis basicity of the anion was estimated using H-1 NMR spectroscopy, by comparison of the measured chemical shifts of the C-2 hydrogen in the imidazolium ring. Finally, compositions containing free chloride anions (chi(SnCl2) < 0.50) were found to oxidize slowly in air to form a chlorostannate(IV) ionic liquid containing the [SnCl6](2-) anion.

Chlorometallate(III) ionic liquids as Lewis acidic catalysts - a quantitative study of acceptor properties

The Gutmann Acceptor Number (AN), which is a quantitative measure of Lewis acidity, has been estimated using the P-31 NMR chemical shift of a probe molecule, triethylphosphine oxide, for a range of chlorometallate(III) ionic liquids, based on Group 13 metals (aluminium(III), gallium(III) and indium(III)) and the 1-octyl-3-methylimidazolium cation, at different compositions. The results were interpreted in terms of extant speciation studies of chlorometallate(III) ionic liquids, and compared with a range of standard molecular solvents and acids. The value of these data were illustrated in terms of the selection of appropriate ionic liquids for specific applications.

Rational design of a dual-mode optical and chemical prodrug

Purpose. The purpose of this study is to demonstrate the rational design and behaviour of the first dual mode optical and chemical prodrug, exemplified by an acetyl salicylic acid-based system. Methods. A cyclic 1,4-benzodioxinone prodrug was synthesised by reaction of 3,5-dimethoxybenzoin and acetyl salicoyl chloride with pyridine. After purification by column chromatography and recrystallization, characterization was achieved using infrared and NMR spectroscopies, mass spectrometry, elemental analysis and single crystal X-ray diffraction. Light-triggered drug liberation was characterised via UV-visible spectroscopy following low-power 365 nm irradiation for controlled times. Chemical drug liberation was characterised via UV-visible spectroscopy in pH 5.5 solution. Results. The synthetic method yielded pure prodrug, with full supporting characterisation. Light-triggered drug liberation proceeded at a rate of 8.30 10j2 sj1, while chemical, hydrolytic liberation proceeded independently at 1.89 10j3 sj1. The photochemical and hydrolytic reactions were both quantitative. Conclusions. This study demonstrates the first rational dual-mode optical and chemical prodrug, using acetyl salicylic acid as a model, acting as a paradigm for future dual-mode systems. Photochemical drug liberation proceeds 44 times faster than chemical liberation, suggesting potential use in drug-eluting medical devices where an additional burst of drug is required at the onset of infection.

Neutron diffraction, NMR and molecular dynamics study of glucose dissolved in the ionic liquid 1-ethyl-3-methylimidazolium acetate

beta-D-glucose dissolved in the ionic liquid 1-ethyl-3-methylimidazolium acetate in a 6 : 1 molar ratio (ionic liquid : glucose) has been studied by neutron scattering, NMR and molecular dynamics simulations. Good agreement was found between simulated neutron scattering profiles generated for isotopically substituted liquid systems and those experimentally determined as well as between simulated and experimental diffusion coefficients obtained by Pulsed Field Gradient NMR spectroscopy. The overriding glucose-ionic liquid interactions in the liquid are hydrogen-bonding between acetate oxygens and sugar hydroxyl groups. The ionic liquid cation was found to play only a minor role in the solvation of the sugar and does not participate in hydrogen-bonding with the sugar to any significant degree. NOESY experiments lend further evidence that there is no direct interaction between sugar hydroxyl groups and acidic hydrogens on the ionic liquid cation.

OXOIRON(IV) PORPHYRINS DERIVED FROM CHARGED IRON(III) TETRAARYLPORPHYRINS AND CHEMICAL OXIDANTS IN AQUEOUS AND METHANOLIC SOLUTIONS

The oxidation of six charged iron(III) tetraarylporphyrins with chemical oxidants has been investigated. In aqueous solution each can be converted by tert-butyl hydroperoxide or monopersulphate into its corresponding oxoiron(IV) porphyrin, whereas in methanol only the iron(III) tetra(N-methylpyridyl)porphyrins form detectable ferryl porphyrins at ambient temperatures. On standing, the iron species revert to the parent porphyrin with a small loss due to non-reversible oxidative destruction. That the oxidised porphyrin intermediates are oxoiron(IV) species has been determined using UV-VIS, resonance Raman, H1 NMR and EPR spectroscopy.

NMR analysis of Nile Blue (C I Basic Blue 12) and Thionine (C I 52000) in solution

The dyes Nile Blue (C I Basic Blue 12) and Thionine (C I 52000) were examined in both ionic and neutral forms in different solvents using NMR and UV-visible spectroscopy to firmly establish the structures of the molecules and to assess the nature and extent of their aggregation H-1 and C-13 NMR assignments and chemical shift data were used together with nuclear Overhauser effect information to propose a self-assembly structure These data were supplemented with variable temperature dilution and diffusion-based experimental results using H-1 NMR spectroscopy thereby enabling extended aggregate structures to be assessed in terms of the relative strength of self-association and the extent to which extended aggregates could form (C) 2010 Elsevier Ltd All rights reserved

Chemical structural effects on γ-ray spectra of positron annihilation in fluorobenzenes

Spectra of ?-ray Doppler shifts for positron annihilation in benzene and its fluoro-derivatives are simulated using low energy plane wave positron (LEPWP) approximation. The results are compared with available measurements. It is found that the Doppler shifts in these larger aromatic compounds are dominated by the contributions of the valence electrons and that the LEPWP model overestimates the measurements by approximately 30%, in agreement with previous findings in noble gases and small molecules. It is further revealed that the halogen atoms not only switch the sign of the charges on carbon atoms that they bond to, but that they also polarize other C-H bonds in the molecule leading to a redistribution of the molecular electrostatic potentials. As a result, it is likely that the halogen atoms contribute more significantly to the annihilation process. The present study also suggests that, while the Doppler shifts are sensitive to the number of valence electrons in the molecules, they are less sensitive to the chemical structures of isomers that have the same numbers and type of atoms and, hence, the same numbers of electrons. Further investigation of this effect is warranted. © EDP Sciences, Società Italiana di Fisica, Springer-Verlag 2012.

Structure and dynamics of aqueous 2-propanol: A THz-TDS, NMR and neutron diffraction study

Aqueous liquid mixtures, in particular, those involving amphiphilic species, play an important role in many physical, chemical and biological processes. Of particular interest are alcohol/water mixtures; however, the structural dynamics of such systems are still not fully understood. Herein, a combination of terahertz time-domain spectroscopy (THz-TDS) and NMR relaxation time analysis has been applied to investigate 2-propanol/water mixtures across the entire composition range; while neutron diffraction studies have been carried out at two specific concentrations. Excellent agreement is seen between the techniques with a maximum in both the relative absorption coefficient and the activation energy to molecular motion occurring at ∼90 mol% H₂O. Furthermore, this is the same value at which well-established excess thermodynamic functions exhibit a maximum/minimum. Additionally, both neutron diffraction and THz-TDS have been used to provide estimates of the size of the hydration shell around 2-propanol in solution. Both methods determine that between 4 and 5 H₂O molecules per 2-propanol are found in the 2-propanol/water clusters at 90 mol% H₂O. Based on the acquired data, a description of the structure of 2-propanol/water across the composition range is presented.

Prevotella denticola Lipopolysaccharide from a Cystic Fibrosis Isolate Possesses a Unique Chemical Structure

We report the first complete structural characterization of the lipopolysaccharide (LPS) from a cystic fibrosis (CF) clinical isolate of Prevotella denticola (B003V1S1X). Chemical, spectroscopic, and spectrometric analyses revealed a unique rough-type LPS (LOS) structure. The structure has a highly negatively charged heptasaccharide core region containing hexoses, with the first two sugars, 3-deoxy-D-manno-oct-2-ulosonic acid (Kdo) and mannose, highly phosphorylated. Furthermore, the lipid A moiety has the typical structure for the genus Prevotella, and was also highly phosphorylated.