Computer simulations of structures, energetics and dynamics of myoglobin center dot center dot center dot ligand complexes

Myoglobin has been studied in considerable detail using different experimental and computational techniques over the past decades. Recent developments in time-resolved spectroscopy have provided experimental data amenable to detailed atomistic simulations. The main theme of the present review are results on the structures, energetics and dynamics of ligands ( CO, NO) interacting with myoglobin from computer simulations. Modern computational methods including free energy simulations, mixed quantum mechanics/molecular mechanics simulations, and reactive molecular dynamics simulations provide insight into the dynamics of ligand dynamics in confined spaces complementary to experiment. Application of these methods to calculate and understand experimental observations for myoglobin interacting with CO and NO are presented and discussed.

A complete atomistic model of molten polyethylene from neutron scattering data: a new methodology for polymer structure

A new approach to the study of the local organization in amorphous polymer materials is presented. The method couples neutron diffraction experiments that explore the structure on the spatial scale 1–20 Å with the reverse Monte Carlo fitting procedure to predict structures that accurately represent the experimental scattering results over the whole momentum transfer range explored. Molecular mechanics and molecular dynamics techniques are also used to produce atomistic models independently from any experimental input, thereby providing a test of the viability of the reverse Monte Carlo method in generating realistic models for amorphous polymeric systems. An analysis of the obtained models in terms of single chain properties and of orientational correlations between chain segments is presented. We show the viability of the method with data from molten polyethylene. The analysis derives a model with average C-C and C-H bond lengths of 1.55 Å and 1.1 Å respectively, average backbone valence angle of 112, a torsional angle distribution characterized by a fraction of trans conformers of 0.67 and, finally, a weak interchain orientational correlation at around 4 Å.

Structure and magnetic properties of mono- and bi-dispersed ferrofluids as revealed by simulations

Using a Langevin molecular dynamics simulation, we show that the magnetic properties of a mono- and bi-dispersed ferrofluid system depend on the volume fraction and the dipolar coupling parameter. For the bi-dispersed system, most of the chains are formed by the large particles, but the aggregation behavior of the large particles is hindered by the presence of the small particles, which are predominantly attached to the end of the particle chain. To further elucidate the microscopic fluid structure, anisotropic structure factors are calculated.

Influence of the size distribution of particles on the viscous property of an electrorheological fluid

The influence of the size distribution of particles on the viscous property of an electrorheological fluid has been investigated by the molecular dynamic simulation method. The shear stress of the fluid is found to decrease with the increase of the variance sigma(2) of the Gaussian distribution of the particle size, and then reach a steady value when sigma is larger than 0.5. This phenomenon is attributed to the influence of the particle size distribution on the dynamic structural evolution in the fluid as well as the strength of the different chain-like structures formed by the particles.

Many-body effect in electrorheological responses

The many-body effect in the kinetic responses of ER fluids is studied by a molecular-dynamic simulation method. The mutual polarization effects of the particles are considered by self-consistently calculating the dipole strength on each particle according to the external field and the dipole field due to all the other particles in the fluids. The many-body effect is found to increase with the enhancement of the particle concentration and the permittivity ratio between the solvent and the particles. The calculated response times are shorter than that predicted with the 'point-dipole' model and agree very well with experimental results. The many-body effect enhances the shear stresses of the fluids by several times. But they are not proportional to the many-body correction factor lambda as expected. This is due to the fact that larger interaction forces between the particles lead to coarsening of the fibers formed in the suspensions. The results show that the many-body and multipolar interaction between the particles must be treated comprehensively in the simulations in order to get more reliable results.

The effect of electric field direction on the-shear stress of ER fluids

The effect of the direction of external electric field on the shear stress of an ER fluid has been studied by molecular-dynamics simulation. Due to the formation of inclined chains, the shear stress strongly depends on the direction of the field, and it may be very large under some special field direction. And theoretical model of ideal microstructure of ER fluids has proved this result. Thus the ER effect may be greatly enhanced just by choosing an optimum direction for the field without any additional requirement, suggesting a promising way to the practical application of ER fluids.

First structural analysis of a naphthalene-based poly(ether ketone): Crystal and molecular simulation from X-ray powder data and diffraction modeling

Polycondensation of 2,6-dihydroxynaphthalene with 4,4'-bis(4"-fluorobenzoyl)biphenyl affords a novel, semicrystalline poly(ether ketone) with a melting point of 406 degreesC and glass transition temperature (onset) of 168 degreesC. Molecular modeling and diffraction-simulation studies of this polymer, coupled with data from the single-crystal structure of an oligomer model, have enabled the crystal and molecular structure of the polymer to be determined from X-ray powder data. This structure-the first for any naphthalene-containing poly(ether ketone)-is fully ordered, in monoclinic space group P2(1)/b, with two chains per unit cell. Rietveld refinement against the experimental powder data gave a final agreement factor (R-wp) of 6.7%.

Simple model of adsorption on external surface of carbon nanotubes: a new analytical approach basing on molecular simulation data

Nitrogen adsorption on carbon nanotubes is wide- ly studied because nitrogen adsorption isotherm measurement is a standard method applied for porosity characterization. A further reason is that carbon nanotubes are potential adsorbents for separation of nitrogen from oxygen in air. The study presented here describes the results of GCMC simulations of nitrogen (three site model) adsorption on single and multi walled closed nanotubes. The results obtained are described by a new adsorption isotherm model proposed in this study. The model can be treated as the tube analogue of the GAB isotherm taking into account the lateral adsorbate-adsorbate interactions. We show that the model describes the simulated data satisfactorily. Next this new approach is applied for a description of experimental data measured on different commercially available (and characterized using HRTEM) carbon nanotubes. We show that generally a quite good fit is observed and therefore it is suggested that the observed mechanism of adsorption in the studied materials is mainly determined by adsorption on tubes separated at large distances, so the tubes behave almost independently.

GSHMC: an efficient method for molecular simulation

The hybrid Monte Carlo (HMC) method is a popular and rigorous method for sampling from a canonical ensemble. The HMC method is based on classical molecular dynamics simulations combined with a Metropolis acceptance criterion and a momentum resampling step. While the HMC method completely resamples the momentum after each Monte Carlo step, the generalized hybrid Monte Carlo (GHMC) method can be implemented with a partial momentum refreshment step. This property seems desirable for keeping some of the dynamic information throughout the sampling process similar to stochastic Langevin and Brownian dynamics simulations. It is, however, ultimate to the success of the GHMC method that the rejection rate in the molecular dynamics part is kept at a minimum. Otherwise an undesirable Zitterbewegung in the Monte Carlo samples is observed. In this paper, we describe a method to achieve very low rejection rates by using a modified energy, which is preserved to high-order along molecular dynamics trajectories. The modified energy is based on backward error results for symplectic time-stepping methods. The proposed generalized shadow hybrid Monte Carlo (GSHMC) method is applicable to NVT as well as NPT ensemble simulations.

New insights on growth mechanisms of protein clusters at surfaces: an AFM and simulation study

Despite its relevance to a wide range of technological and fundamental areas, a quantitative understanding of protein surface clustering dynamics is often lacking. In inorganic crystal growth, surface clustering of adatoms is well described by diffusion-aggregation models. In such models, the statistical properties of the aggregate arrays often reveal the molecular scale aggregation processes. We investigate the potential of these theories to reveal hitherto hidden facets of protein clustering by carrying out concomitant observations of lysozyme adsorption onto mica surfaces, using atomic force microscopy. and Monte Carlo simulations of cluster nucleation and growth. We find that lysozyme clusters diffuse across the substrate at a rate that varies inversely with size. This result suggests which molecular scale mechanisms are responsible for the mobility of the proteins on the substrate. In addition the surface diffusion coefficient of the monomer can also be extracted from the comparison between experiments and simulations. While concentrating on a model system of lysozyme-on-mica, this 'proof of concept' study successfully demonstrates the potential of our approach to understand and influence more biomedically applicable protein-substrate couples.

Notations and conventions in molecular spectroscopy: part 1. General spectroscopic notation

The field of Molecular Spectroscopy was surveyed in order to determine a set of conventions and symbols which are in common use in the spectroscopic literature. This document, which is Part I in a series, establishes the notations and conventions used for general spectroscopic notations and deals with quantum mechanics, quantum numbers (vibrational states, angular momentum and energy levels), spectroscopic transitions, and miscellaneous notations (e.g. spectroscopic terms). Further parts will follow, dealing inter alia with symmetry notation, permutation and permutation-inversion symmetry notation, vibration-rotation spectroscopy and electronic spectroscopy.

Crystal and molecular structure of the thioether-analogue of PEEK from X-ray powder data and diffraction-modelling

Analysis of X-ray powder data for the melt-crystallisable aromatic poly(thioether thioether ketone) [-S-Ar-S-Ar-CO-Ar](n), ('PTTK', Ar= 1,4-phenylene), reveals that it adopts a crystal structure very different from that established for its ether-analogue PEEK. Molecular modelling and diffraction-simulation studies of PTTK show that the structure of this polymer is analogous to that of melt-crystallised poly(thioetherketone) [-SAr-CO-Ar](n) in which the carbonyl linkages in symmetry-related chains are aligned anti-parallel to one another. and that these bridging units are crystallographically interchangeable. The final model for the crystal structure of PTTK is thus disordered, in the monoclinic space group 121a (two chains per unit cell), with cell dimensions a = 7.83, b = 6.06, c = 10.35 angstrom, beta = 93.47 degrees. (c) 2005 Elsevier Ltd. All rights reserved.

Quantum mechanics studies of the tautomers of dlhydrodiacetylformoin, an important Maillard intermediate and the synthesis of furaneol

Conformational analyses have been carried out on the acyclic and cyclic forms of dihydrodiacetylformoin, an important Maillard intermediate and precursor for furaneol. For the acyclic forms, the 2,5-dicarbonyl isomers have the lowest energy, while for the cyclic forms, the 3-carbonyl are favoured over the 4-carbonyl isomers. The likely path for cyclisation is investigated and it is shown that the favoured path is dependent upon the relative chiralities of the carbon atoms and in particular that the reaction proceeds more readily if C2 and C3 have different chiralities. After cyclisation, the reaction path to produce furaneol proceeds via the loss of a water molecule. This reaction has been studied with a model including two water molecules and a hydroxide anion and shows relatively low-energy barriers. (C) 2008 Published by Elsevier B.V.

Solid phases of cyclopentane: combined experimental and simulation study

The phase diagram of cyclopentane has been studied by powder neutron diffraction, providing diffraction patterns for phases I, II, and III, over a range of temperatures and pressures. The putative phase IV was not observed. The structure of the ordered phase III has been solved by single-crystal diffraction. Computational modeling reveals that there are many equienergetic ordered structures for cyclopentane within a small energy range. Molecular dynamics simulations reproduce the structures and diffraction patterns for phases I and III and also show an intermediate disordered phase, which is used to interpret phase II.