Energetic basis of molecular recognition in a DNA aptamer

The thermal stability and ligand binding properties of the L-argininamide-binding DNA aptamer (5'-GATCGAAACGTAGCGCCTTCGATC3') were studied by spectroscopic and calorimetric methods. Differential calorimetric studies showed that the uncomplexed aptamer melted in a two-state reaction with a melting temperature T-m = 50.2 +/- 0.2 degrees C and a folding enthalpy Delta H degrees(fold) = -49.0 +/- 2.1 kcal mol(-1). These values agree with values of T-m = 49.6 degrees C and Delta H degrees(fold) = -51.2 kcal mol(-1) predicted for a simple hairpin structure. Melting of the uncomplexed aptamer was dependent upon salt concentration, but independent of strand concentration. The T of aptamer melting was found to increase as L-argininamide concentrations increased. Analysis of circular dichroism titration data using a single-site binding model resulted in the determination of a binding free energy Delta G degrees(bind) = -5.1 kcal mol(-1). Isothermal titration calorimetry studies revealed an exothermic binding reaction with Delta H degrees(bind) = -8.7 kcal mol(-1). Combination of enthalpy and free energy produce ail unfavorable entropy of -T Delta S degrees = +3.6 kcal mol(-1). A molar heat capacity change of -116 cal mol(-1) K-1 was determined from calorimetric measurements at four temperatures over the range of 15-40 degrees C. Molecular dynamics simulations were used to explore the structures of the unligated and ligated aptamer structures.

X-ray crystal structure and molecular mechanics calculations of (2,6-iso-dipropyl-phenylamide) dimethyl (tetra-methylcyclopenta-dienyl) silane titanium dichloride

Crystal and molecular structure of (2.6-dipropylphenylamide) dimethyl (tetra-methyl cyclopentadienyl) silane titanium dichloride (I) was fully characterized by X-ray diffraction. The crystal is obtained from a mixture of ether/hexane as orthorhombic. with a = 12.658 (3) Angstrom. b = 16.62 (3) Angstrom. c = 11.760 (2) Angstrom. V = 2474.2 (9) Angstrom(3). Z = 4, space group Pnma. R = 0.0399; Componud I compose of the pi-bounded ring with its dimethylsilyl-dipropyl phenyl amido group and the two terminal chloride atoms coordinated to central metal to form a so-called constrained geometry catalyst (CGC) structure. The result of molecular mechanics (MM) calculations on compound I shows that bond lengths and bond angles from the MM calculation are comparable to the data obtained from the X-ray diffraction study. The relation of the structure of CGCs and their catalytic activity by MM calculations is also discussed.

Effects of pressure and molecular weight on the miscibility of polystyrene and cyclohexane

With the aid of Sanchez-Lacombe lattice fluid theory (SLLFT), the phase diagrams were calculated for the system cyclohexane (CH)/polystyrene (PS) with different molecular weights at different pressures. The experimental data is in reasonable agreement with SLLFT calculations. The total Gibbs interaction energy, g*(12) for different molecular weights PS at different pressures was expressed, by means of a universal relationship, as g(12)* =f(12)* + (P - P-0) nu*(12) demixing curves were then calculated at fixed (near critical) compositions of CH and PS systems for different molecular weights. The pressures of optimum miscibility obtained from the Gibbs interaction energy are close to those measured by Wolf and coworkers. Furthermore, a reasonable explanation was given for the earlier observation of Saeki et al., i.e., the phase separation temperatures of the present system increase with the increase of pressure for the low molecular weight of the polymer whereas they decrease for the higher molecular weight polymers. The effects of molecular weight, pressure, temperature and composition on the Flory Huggins interaction parameter can be described by a general equation resulting from fitting the interaction parameters by means of Sanchez-Lacombe lattice fluid theory.

Molecular dynamics simulation of sub-transition for polyethersulfone

Molecular dynamics (MD) simulations of a polyethersulfone (PES) chain are carried out in the amorphous state by using the Dreiding 2.21 force field at four temperatures. Two types of molecular motion, i.e, rotations of phenylene rings and torsions of large segments containing two oxygen atoms, two sulfur atoms, and five phenylene rings on the backbone, are simulated. The modeling results show that the successive phenylene rings should be in-phase cooperative rotations, whereas the successive large segments should be out-of-phase cooperative torsions. By calculating the diffusion coefficient for the phenylene ring rotations, it is found that this rotation contributes to the beta -transition of PES.

Molecular dynamics simulation of polystyrene-block-poly(methyl methacrylate)

Molecular dynamics is applied to the system of polystyrene-block-poly(methyl methacrylate). The simulation shows that for the block copolymer system, a layered structure, which reflects microphase separation, is obtained and this structure is stable. In order to elucidate that the formation of the layered structure is reasonable, some static properties such as the radial distribution function and the dipole moment are analyzed in some detail.

MOLECULAR-GEOMETRY AND CHAIN ENTANGLEMENT - PARAMETERS FOR THE TUBE MODEL

The tube diameter in the reptation model is the distance between a given chain segment and its nearest segment in adjacent chains. This dimension is thus related to the cross-sectional area of polymer chains and the nearest approach among chains, without effects of thermal fluctuation and steric repulsion. Prior calculated tube diameters are much larger, about 5 times, than the actual chain cross-sectional areas. This is ascribed to the local freedom required for mutual rearrangement among neighboring chain segments. This tube diameter concept seems to us to infer a relationship to the corresponding entanglement spacing. Indeed, we report here that the critical molecular weight, M(c), for the onset of entanglements is found to be M(c) = 28 A/([R2]0/M), where A is the chain cross-sectional area and [R2]0 the mean-square end-to-end distance of a freely jointed chain of molecular weight M. The new, computed relationship between the critical number of backbone atoms for entanglement and the chain cross-sectional area of polymers, N(c) = A0,44, is concordant with the cross-sectional area of polymer chains being the parameter controlling the critical entanglement number of backbone atoms of flexible polymers.

Molecular dynamics-based approach to study the anisotropic self-diffusion of molecules in porous materials with multiple cage types: Application to H-2 in losod

Van den Berg, A. W. C., Flikkema, E., Lems, S., Bromley, S. T., Jansen, J. C. (2006). Molecular dynamics-based approach to study the anisotropic self-diffusion of molecules in porous materials with multiple cage types: Application to H-2 in losod. Journal of physical chemistry b, 110 (1), 501-506. RAE2008

The HITRAN molecular spectroscopic database: Edition of 2000 including updates through 2001

This paper describes the status circa 2001, of the HITRAN compilation that comprises the public edition available through 2001. The HITRAN compilation consists of several components useful for radiative transfer calculation codes: high-resolution spectroscopic parameters of molecules in the gas phase, absorption cross-sections for molecules with very dense spectral features, aerosol refractive indices, ultraviolet line-by-line parameters and absorption cross-sections, and associated database management software. The line-by-line portion of the database contains spectroscopic parameters for 38 molecules and their isotopologues and isotopomers suitable for calculating atmospheric transmission and radiance properties. Many more molecular species are presented in the infrared cross-section data than in the previous edition, especially the chlorofluorocarbons and their replacement gases. There is now sufficient representation so that quasi-quantitative simulations can be obtained with the standard radiance codes. In addition to the description and justification of new or modified data that have been incorporated since the last edition of HITRAN (1996), future modifications are indicated for cases considered to have a significant impact on remote-sensing experiments. © 2003 Elsevier Ltd. All rights reserved.

Modelling the nano-scale phenomena in condensed matter physics via computer-based numerical simulations

A review of the atomistic modelling of the behaviour of nano-scale structures and processes via molecular dynamics (MD) simulation method of a canonical ensemble is presented. Three areas of application in condensed matter physics are considered. We focus on the adhesive and indentation properties of the solid surfaces in nano-contacts, the nucleation and growth of nano-phase metallic and semi-conducting atomic and molecular films on supporting substrates, and the nano- and multi-scale crack propagation properties of metallic lattices. A set of simulations selected from these fields are discussed, together with a brief introduction to the methodology of the MD simulation. The pertinent inter-atomic potentials that model the energetics of the metallic and semi-conducting systems are also given.

Free energy and molecular dynamics calculations for the cubic-tetragonal phase transition in zirconia

The high-temperature cubic-tetragonal phase transition of pure stoichiometric zirconia is studied by molecular dynamics (MD) simulations and within the framework of the Landau theory of phase transformations. The interatomic forces are calculated using an empirical, self-consistent, orthogonal tight-binding model, which includes atomic polarizabilities up to the quadrupolar level. A first set of standard MD calculations shows that, on increasing temperature, one particular vibrational frequency softens. The temperature evolution of the free-energy surfaces around the phase transition is then studied with a second set of calculations. These combine the thermodynamic integration technique with constrained MD simulations. The results seem to support the thesis of a second-order phase transition but with unusual, very anharmonic behavior above the transition temperature.

An investigation into the use of atomic datasets in simulations of the Ni-like gadolinium x-ray laser

The effect of differing the datasets used in the modelling of the Ni-like Gd x-ray laser (XRL) is examined through the 1.50 hydro-atomic code, EHYBRID. Two atomic datasets, including energy levels and radiative and collisional excitation rates, are used as input data for the code. It is found that the behaviour of the XRL is somewhat different than might be expected from superficial examination of the atomic data. The similarities in the gain profiles at low densities are found to have encouraging implications. in our attempts to model XRLs.

Ab initio molecular dynamics simulation of a room temperature ionic liquid

Ab initio molecular dynamics simulations have been performed for the first time on the room-temperature organic ionic liquid dimethyl imidazolium chloride [DMIM][Cl] using density functional theory. The aim is to compare the local liquid structure with both that obtained from two different classical force fields and from neutron scattering experiments. The local structure around the cation shows significant differences compared to both the classical calculations and the neutron results. In particular, and unlike in the gas-phase ion pair, chloride ions tend to be located near a ring C-H proton in a position suggesting hydrogen bonding. The results are used to suggest ways in which the classical potentials may be improved.

Surface alloying and mixing at the Mn/Fe(001) interface: Real-time photoelectron spectroscopy and modified embedded atom simulations

Structural and magnetic properties of thin Mn films on the Fe(001) surface have been investigated by a combination of photoelectron spectroscopy and computer simulation in the temperature range 300 Kless than or equal toTless than or equal to750 K. Room-temperature as deposited Mn overlayers are found to be ferromagnetic up to 2.5-monolayer (ML) coverage, with a magnetic moment parallel to that of the iron substrate. The Mn atomic moment decreases with increasing coverage, and thicker samples (4-ML and 4.5-ML coverage) are antiferromagnetic. Photoemission measurements performed while the system temperature is rising at constant rate (dT/dtsimilar to0.5 K/s) detect the first signs of Mn-Fe interdiffusion at T=450 K, and reveal a broad temperature range (610 Kless than or equal toTless than or equal to680 K) in which the interface appears to be stable. Interdiffusion resumes at Tgreater than or equal to680 K. Molecular dynamics and Monte Carlo simulations allow us to attribute the stability plateau at 610 Kless than or equal toTless than or equal to680 K to the formation of a single-layer MnFe surface alloy with a 2x2 unit cell and a checkerboard distribution of Mn and Fe atoms. X-ray-absorption spectroscopy and analysis of the dichroic signal show that the alloy has a ferromagnetic spin structure, collinear with that of the substrate. The magnetic moments of Mn and Fe atoms in the alloy are estimated to be 0.8mu(B) and 1.1mu(B), respectively.

Kinetics of phase transformations in a model with metastable fluid-fluid separation: A molecular dynamics study

By molecular dynamics (MD) simulations we study the crystallization process in a model system whose particles interact by a spherical pair potential with a narrow and deep attractive well adjacent to a hard repulsive core. The phase diagram of the model displays a solid-fluid equilibrium, with a metastable fluid-fluid separation. Our computations are restricted to fairly small systems (from 2592 to 10368 particles) and cover long simulation times, with constant energy trajectories extending up to 76x10(6) MD steps. By progressively reducing the system temperature below the solid-fluid line, we first observe the metastable fluid-fluid separation, occurring readily and almost reversibly upon crossing the corresponding line in the phase diagram. The nucleation of the crystal phase takes place when the system is in the two-fluid metastable region. Analysis of the temperature dependence of the nucleation time allows us to estimate directly the nucleation free energy barrier. The results are compared with the predictions of classical nucleation theory. The critical nucleus is identified, and its structure is found to be predominantly fcc. Following nucleation, the solid phase grows steadily across the system, incorporating a large number of localized and extended defects. We discuss the relaxation processes taking place both during and after the crystallization stage. The relevance of our simulation for the kinetics of protein crystallization under normal experimental conditions is discussed. (C) 2002 American Institute of Physics.

A dielectric continuum molecular dynamics method

We introduce a novel method to simulate hydrated macromolecules with a dielectric continuum representation of the surrounding solvent. In our approach, the interaction between the solvent and the molecular degrees of freedom is described by means of a polarization density free energy functional which is minimum at electrostatic equilibrium. After a pseudospectral expansion of the polarization and a discretization of the functional, we construct the equations of motion for the system based on a Car-Parrinello technique. In the limit of the adiabatic evolution of the polarization field variables, our method provides the solution of the dielectric continuum problem "on the fly," while the molecular coordinates are propagated. In this first study, we show how our dielectric continuum molecular dynamics method can be successfully applied to hydrated biomolecules, with low cost compared to free energy simulations with explicit solvent. To our knowledge, this is the first time that stable and conservative molecular dynamic simulations of solutes can be performed for a dielectric continuum model of the solvent. (C) 2001 American Institute of Physics.