The Structure of C2H4 Clusters from Theoretical Interaction Potentials and Vibrational Predissociation Data

Optimized geometries and binding energies are calculated for ethene (ethylene) dimers, trimers, and tetramers based on a pairwise additive dimer potential. From these results intermolecular frequencies and relative abundancies (catchment areas) of the different isomers are obtained and compared with the results of accurate measurements of the photodissociation upon absorption of one photon of a CO2 laser in the region of the ν7 monomer absorption band at 949 cm-1. The clusters are size selected in a scattering experiment and show for a cluster size from n=2 to n=6 a frequency maximum shifted by 3 cm-1 to the blue compared with the monomer. The result is explained by the predominance of chains and chain-like structures of the clusters in the photodissociation process. The chains consist of cross-like dimer sub-units

Electron correlations for ground state properties of group IV semiconductors

Valence energies for crystalline C, Si, Ge, and Sn with diamond structure have been determined using an ab-initio approach based on information from cluster calculations. Correlation contributions, in particular, have been evaluated in the coupled electron pair approximation (CEPA), by means of increments obtained for localized bond orbitals and for pairs and triples of such bonds. Combining these results with corresponding Hartree-Fock (HF) data, we recover about 95 % of the experimental cohesive energies. Lattice constants are overestimated at the HF level by about 1.5 %; correlation effects reduce these deviations to values which are within the error bounds of this method. A similar behavior is found for the bulk modulus: the HF values which are significantly too high are reduced by correlation effects to about 97 % of the experimental values.Comment: 22 pages, latex, 2 figure

Optimizing end-labeled free-solution electrophoresis by increasing the hydrodynamic friction of the drag-tag

We study the electrophoretic separation of polyelectrolytes of varying lengths by means of end-labeled free-solution electrophoresis (ELFSE). A coarse-grained molecular dynamics simulation model, using full electrostatic interactions and a mesoscopic Lattice Boltzmann fluid to account for hydrodynamic interactions, is used to characterize the drag coefficients of different label types: linear and branched polymeric labels, as well as transiently bound micelles. It is specifically shown that the label's drag coefficient is determined by its hydrodynamic size, and that the drag per label monomer is largest for linear labels. However, the addition of side chains to a linear label offers the possibility to increase the hydrodynamic size, and therefore the label efficiency, without having to increase the linear length of the label, thereby simplifying synthesis. The third class of labels investigated, transiently bound micelles, seems very promising for the usage in ELFSE, as they provide a significant higher hydrodynamic drag than the other label types. The results are compared to theoretical predictions, and we investigate how the efficiency of the ELFSE method can be improved by using smartly designed drag-tags.Comment: 32 pages, 11 figures, submitted to Macromolecule

Multi-Particle Collision Dynamics -- a Particle-Based Mesoscale Simulation Approach to the Hydrodynamics of Complex Fluids

In this review, we describe and analyze a mesoscale simulation method for fluid flow, which was introduced by Malevanets and Kapral in 1999, and is now called multi-particle collision dynamics (MPC) or stochastic rotation dynamics (SRD). The method consists of alternating streaming and collision steps in an ensemble of point particles. The multi-particle collisions are performed by grouping particles in collision cells, and mass, momentum, and energy are locally conserved. This simulation technique captures both full hydrodynamic interactions and thermal fluctuations. The first part of the review begins with a description of several widely used MPC algorithms and then discusses important features of the original SRD algorithm and frequently used variations. Two complementary approaches for deriving the hydrodynamic equations and evaluating the transport coefficients are reviewed. It is then shown how MPC algorithms can be generalized to model non-ideal fluids, and binary mixtures with a consolute point. The importance of angular-momentum conservation for systems like phase-separated liquids with different viscosities is discussed. The second part of the review describes a number of recent applications of MPC algorithms to study colloid and polymer dynamics, the behavior of vesicles and cells in hydrodynamic flows, and the dynamics of viscoelastic fluids

A combined computational and experimental investigation of the [2Fe–2S] cluster in biotin synthase

Biotin synthase was the first example of what is now regarded as a distinctive enzyme class within the radical S-adenosylmethionine superfamily, the members of which use Fe/S clusters as the sulphur source in radical sulphur insertion reactions. The crystal structure showed that this enzyme contains a [2Fe–2S] cluster with a highly unusual arginine ligand, besides three normal cysteine ligands. However, the crystal structure is at such a low resolution that neither the exact coordination mode nor the role of this exceptional ligand has been elucidated yet, although it has been shown that it is not essential for enzyme activity. We have used quantum refinement of the crystal structure and combined quantum mechanical and molecular mechanical calculations to explore possible coordination modes and their influences on cluster properties. The investigations show that the protonation state of the arginine ligand has little influence on cluster geometry, so even a positively charged guanidinium moiety would be in close proximity to the iron atom. Nevertheless, the crystallised enzyme most probably contains a deprotonated (neutral) arginine coordinating via the NH group. Furthermore, the Fe···Fe distance seems to be independent of the coordination mode and is in perfect agreement with distances in other structurally characterised [2Fe–2S] clusters. The exceptionally large Fe···Fe distance found in the crystal structure could not be reproduced

Structure of the Nucleotide Radical Formed during Reaction of CDP/TTP with the E441Q-α2β2 of E. coli Ribonucleotide Reductase

The Escherichia coli ribonucleotide reductase (RNR) catalyzes the conversion of nucleoside diphosphates to deoxynucleotides and requires a diferric-tyrosyl radical cofactor for catalysis. RNR is composed of a 1:1 complex of two homodimeric subunits: α and β. Incubation of the E441Q-α mutant RNR with substrate CDP and allosteric effector TTP results in loss of the tyrosyl radical and formation of two new radicals on the 200 ms to min time scale. The first radical was previously established by stopped flow UV/vis spectroscopy and pulsed high field EPR spectroscopy to be a disulfide radical anion. The second radical was proposed to be a 4′-radical of a 3′-keto-2′-deoxycytidine 5′-diphosphate. To identify the structure of the nucleotide radical [1′-[superscript 2]H], [2′-[superscript 2]H], [4′-[superscript 2]H], [5′-[superscript 2]H], [U−[superscript 13]C, [superscript 15]N], [U−[superscript 15]N], and [5,6 -[superscript 2]H] CDP and [β-[superscript 2]H] cysteine-α were synthesized and incubated with E441Q-α2β2 and TTP. The nucleotide radical was examined by 9 GHz and 140 GHz pulsed EPR spectroscopy and 35 GHz ENDOR spectroscopy. Substitution of [superscript 2]H at C4′ and C1′ altered the observed hyperfine interactions of the nucleotide radical and established that the observed structure was not that predicted. DFT calculations (B3LYP/IGLO-III/B3LYP/TZVP) were carried out in an effort to recapitulate the spectroscopic observations and lead to a new structure consistent with all of the experimental data. The results indicate, unexpectedly, that the radical is a semidione nucleotide radical of cytidine 5′-diphosphate. The relationship of this radical to the disulfide radical anion is discussed.National Institutes of Health (U.S.) (GM29595)(EB002804)(EB002026

Many-particle Brownian and Langevin Dynamics Simulations with the Brownmove package

<p>Abstract</p> <p>Background</p> <p>Brownian Dynamics (BD) is a coarse-grained implicit-solvent simulation method that is routinely used to investigate binary protein association dynamics, but due to its efficiency in handling large simulation volumes and particle numbers it is well suited to also describe many-protein scenarios as they often occur in biological cells.</p> <p>Results</p> <p>Here we introduce our "brownmove" simulation package which was designed to handle many-particle problems with varying particle numbers and allows for a very flexible definition of rigid and flexible protein and polymer models. Both a Brownian and a Langevin dynamics (LD) propagation scheme can be used and hydrodynamic interactions are treated efficiently with our recently introduced TEA-HI ansatz [Geyer, Winter, JCP 130 (2009) 114905]. With simulations of constrained polymers and flexible models of spherical proteins we demonstrate that it is crucial to include hydrodynamics when multi-bead models are used in BD or LD simulations. Only then both the translational and the rotational diffusion coefficients and the timescales of the internal dynamics can be reproduced correctly. In the third example project we show how constant density boundary conditions [Geyer et al, JCP 120 (2004) 4573] can be used to set up a non-equilibrium simulation of diffusional transport across an array of fixed obstacles. Finally, we demonstrate how the agglomeration dynamics of multiple particles with attractive patches can be analysed conveniently with the help of a dynamic interaction network.</p> <p>Conclusions</p> <p>Combining BD and LD propagation, fast hydrodynamics, a flexible protein model, and interfaces for "open" simulation settings, our freely available "brownmove" simulation package constitutes a new platform for coarse-grained many-particle simulations of biologically relevant diffusion and transport processes.</p