Evaluation of the binding ability of a novel dioxatetraazamacrocyclic receptor that contains two phenanthroline units: selective uptake of carboxylate anions

The novel dioxatetraaza macrocycle [26]phen(2)N(4)O(2), which incorporates two phenanthroline units, has been synthesized, and its acid-base behavior has been evaluated by potentiometric and H-1 NMR methods. Six protonation constants were determined, and the protonation sequence was established by NMR. The location of the fifth proton on the phen nitrogen was confirmed by X-ray determinations of the crystal structures of the receptor as bromide and chloride salts. The two compounds have the general molecular formula {(H-5[26]phen(2)N(4)O(2))X-n(H2O)(5-n)}X(n-1)(.)mH(2)O, where X = Cl, n = 3, and m = 6 or X = Br, n = 4, and m = 5.5. In the solid state, the (H-5[26]phen(2)N(4)O(2))(5+) cation adopts a "horseshoe" topology with sufficient room to encapsulate three or four halogen anions through the several N-(HX)-X-... hydrogen-bonding interactions. Two supermolecules {(H-5[26]phen(2)N(4)O(2))X-n(H2O)(5-n)}((5-n)+) form an interpenetrating dimeric species, which was also found by ESI mass spectrum. Binding studies of the protonated macrocycle with aliphatic (ox(2-), mal(2-), suc(2-), cit(3-), cta(3-)) and aromatic (bzc(-), naphc(-), anthc(-), pyrc(-), ph(2-), iph(2-), tph(2-), btc(3-)) anions were determined in water by potentiometric methods. These studies were complemented by H-1 NMR titrations in D2O of the receptor with selected anions. The H-i[26]phen(2)N(4)O(2)(i+) receptor can selectively uptake highly charged or extended aromatic carboxylate anions, such as btc(3-) and pyrc(-), in the pH ranges of 4.0-8.5 and < 4.0, respectively, from aqueous solution that contain the remaining anions as pollutants or contaminants. To obtain further insight into these structural and experimental findings, molecular dynamics (MD) simulations were carried out in water solution.

Thallium pi-cation complexation by calix[4]tubes: Tl-205 NMR and X-ray evidence

Thallium cation complexation by calix[4]tubes has been investigated by a combination of (TI)-T-205, H-1 NMR and ES MS demonstrating the solution formation of a dithallium complex in which the cations are held in the calix[4]arene cavities. In addition, the structure of the complex has been determined in the solid state revealing the cations to be held exclusively by pi-cation interactions. Furthermore, this crystal structure has been used as the basis for molecular dynamics simulations to confirm that binding of the smaller K+ cation in the calix[4]tube cryptand like array occurs via the axial route featuring a g-cation intermediate.

Development of novel brush-type chiral stationary phases based on terpenoid selectors: HPLC evaluation and theoretical investigation of enantioselective binding interactions

The terpenoid chiral selectors dehydroabietic acid, 12,14-dinitrodehydroabietic acid and friedelin have been covalently linked to silica gel yielding three chiral stationary phases CSP 1, CSP 2 and CSP 3, respectively. The enantiodiscriminating capability of each one of these phases was evaluated by HPLC with four families of chiral aromatic compounds composed of alcohols, amines, phenylalanine and tryptophan amino acid derivatives and beta-lactams. The CSP 3 phase, containing a selector with a large friedelane backbone is particularly suitable for resolving free alcohols and their derivatives bearing fluorine substituents, while CSP 2 with a dehydroabietic architecture is the only phase that efficiently discriminates 1, 1'-binaphthol atropisomers. CSP 3 also gives efficient resolution of the free amines. All three phases resolve well the racemates of N-trifluoracetyl and N-3,5-dinitrobenzoyl phenylalanine amino acid ester derivatives. Good enantioseparation of beta-lactams and N-benzoyl tryptophan amino acid derivatives was achieved on CSP 1. In order to understand the structural factors that govern the chiral molecular recognition ability of these phases, molecular dynamics simulations were carried out in the gas phase with binary diastereomeric complexes formed by the selectors of CSP 1 and CSP 2 and several amino acid derivatives. Decomposition of molecular mechanics energies shows that van der Waals interactions dominate the formation of the diastereomeric transient complexes while the electrostatic binding interactions are primarily responsible for the enantioselective binding of the (R)- and (S)-analytes. Analysis of the hydrogen bonds shows that electrostatic interactions are mainly associated with the formation of N-(HO)-O-...=C enantio selective hydrogen bonds between the amide binding sites from the selectors and the carbonyl groups of the analytes. The role of mobile phase polarity, a mixture of n-hexane and propan-2-ol in different ratios, was also evaluated through molecular dynamics simulations in explicit solvent. (c) 2006 Elsevier Ltd. All rights reserved.

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.

Structure of single-wall peptide nanotubes: in situ flow aligning X-ray diffraction

The structure of single wall peptide nanotubes is presented for the model surfactant-like peptide A6K. Capillary flow alignment of a sample in the nematic phase at high concentration in water leads to oriented X-ray diffraction patterns. Analysis of these, accompanied by molecular dynamics simulations, suggests the favourable self-assembly of antiparallel peptide dimers into beta-sheet ribbons that wrap helically to form the nanotube wall.

Ferric and ferrous iron in nitroso-myoglobin: computer simulations of stable and metastable states and their infrared spectra

The binding of NO to iron is involved in the biological function of many heme proteins. Contrary to ligands like CO and O-2, which only bind to ferrous (Fe-II) iron, NO binds to both ferrous and ferric (Fe-II) iron. In a particular protein, the natural oxidation state can therefore be expected to be tailored to the required function. Herein, we present an ob initio potential-energy surface for ferric iron interacting with NO. This potential-energy surface exhibits three minima corresponding to eta'-NO coordination (the global minimum), eta(1)-ON coordination and eta(2) coordination. This contrasts with the potential-energy surface for Fe-II-NO, which ex- hibits only two minima (the eta(2) coordination mode for Fe-II is a transition state, not a minimum). In addition, the binding energies of NO are substantially larger for Fe-III than for Fe-II. We have performed molecular dynamics simulations for NO bound to ferric myoglobin (Mb(III)) and compare these with results obtained for Mb(II). Over the duration of our simulations (1.5 ns), all three binding modes are found to be stable at 200 K and transiently stable at 300 K, with eventual transformation to the eta(1)-NO global-minimum conformation. We discuss the implication of these results related to studies of rebinding processes in myoglobin.

Potential energy surface and moleculardynamics of MbNO: existence of an unsuspected FeON minimum

Ligands such as CO, O2, or NO are involved in the biological function of myoglobin. Here we investigate the energetics and dynamics of NO interacting with the Fe(II) heme group in native myoglobin using ab initio and molecular dynamics simulations. At the global minimum of the ab initio potential energy surface (PES), the binding energy of 23.4 kcal/mol and the Fe-NO structure compare well with the experimental results. Interestingly, the PES is found to exhibit two minima: There exists a metastable, linear Fe-O-N minimum in addition to the known, bent Fe-N-O global minimum conformation. Moreover, the T-shaped configuration is found to be a saddle point, in contrast to the corresponding minimum for NO interacting with Fe(III). To use the ab initio results for finite temperature molecular dynamics simulations, an analytical function was fitted to represent the Fe-NO interaction. The simulations show that the secondary minimum is dynamically stable up to 250 K and has a lifetime of several hundred picoseconds at 300 K. The difference in the topology of the heme-NO PES from that assumed previously (one deep, single Fe-NO minimum) suggests that it is important to use the full PES for a quantitative understanding of this system. Why the metastable state has not been observed in the many spectroscopic studies of myoglobin interacting with NO is discussed, and possible approaches to finding it are outlined.

Whither tube theory: From believing to measuring

This short contribution examines the difficulties that have not yet been fully overcome in the many developments made from the simplest (and original) tube model for entangled polymers. It is concluded that many more length scales have to be considered sequentially when deriving a continuum rheological model from molecular considerations than have been considered in the past. In particular, most unresolved issues of the tube theory are related to the length scales of tube diameter, and molecular dynamics simulations is the perfect route to resolve them. The power of molecular simulations is illustrated by two examples: stress contributions from bonded and non-bonded interaction, and the inter-chain coupling, which is usually neglected in the tube theory.

Transition metal clusters on graphite

The topic of this work is 3d transition metals deposited on graphite. Spin-polarised density-functional calculations are used to obtain the magnetic moments of deposited adatoms and dimers. Interatomic potentials are also deduced. These are used in molecular dynamics simulations to study cluster formation and to investigate cluster morphology.

Adsorption from aqueous solutions on opened carbon nanotubes: organic compounds speed up delivery of water from inside

We report the results of first systematic studies of organic adsorption from aqueous solutions onto relatively long single walled carbon nanotubes (four tubes, in initial and oxidised forms). Using molecular dynamics simulations (GROMACS package) we discuss the behaviour of tube-water as well as tube-adsorbate systems, for three different adsorbates (benzene, phenol and paracetamol).

Self-assembly of Fmoc-tetrapeptides based on the RGDS cell adhesion motif

Self-assembly in aqueous solution has been investigated for two Fmoc [Fmoc ¼ N-(fluorenyl)-9-methoxycarbonyl] tetrapeptides comprising the RGDS cell adhesion motif from fibronectin or the scrambled sequence GRDS. The hydrophobic Fmoc unit confers amphiphilicity on the molecules, and introduces aromatic stacking interactions. Circular dichroism and FTIR spectroscopy show that the self-assembly of both peptides at low concentration is dominated by interactions among Fmoc units, although Fmoc-GRDS shows b-sheet features, at lower concentration than Fmoc-RGDS. Fibre X-ray diffraction indicates b-sheet formation by both peptides at sufficiently high concentration. Strong alignment effects are revealed by linear dichroism experiments for Fmoc-GRDS. Cryo-TEM and smallangle X-ray scattering (SAXS) reveal that both samples form fibrils with a diameter of approximately 10 nm. Both Fmoc-tetrapeptides form self-supporting hydrogels at sufficiently high concentration. Dynamic shear rheometry enabled measurements of the moduli for the Fmoc-GRDS hydrogel, however syneresis was observed for the Fmoc-RGDS hydrogel which was significantly less stable to shear. Molecular dynamics computer simulations were carried out considering parallel and antiparallel b-sheet configurations of systems containing 7 and 21 molecules of Fmoc-RGDS or Fmoc-GRDS, the results being analyzed in terms of both intermolecular structural parameters and energy contributions.

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 Å.

Towards data warehousing and mining of protein unfolding simulation data

OBJECTIVES: The prediction of protein structure and the precise understanding of protein folding and unfolding processes remains one of the greatest challenges in structural biology and bioinformatics. Computer simulations based on molecular dynamics (MD) are at the forefront of the effort to gain a deeper understanding of these complex processes. Currently, these MD simulations are usually on the order of tens of nanoseconds, generate a large amount of conformational data and are computationally expensive. More and more groups run such simulations and generate a myriad of data, which raises new challenges in managing and analyzing these data. Because the vast range of proteins researchers want to study and simulate, the computational effort needed to generate data, the large data volumes involved, and the different types of analyses scientists need to perform, it is desirable to provide a public repository allowing researchers to pool and share protein unfolding data. METHODS: To adequately organize, manage, and analyze the data generated by unfolding simulation studies, we designed a data warehouse system that is embedded in a grid environment to facilitate the seamless sharing of available computer resources and thus enable many groups to share complex molecular dynamics simulations on a more regular basis. RESULTS: To gain insight into the conformational fluctuations and stability of the monomeric forms of the amyloidogenic protein transthyretin (TTR), molecular dynamics unfolding simulations of the monomer of human TTR have been conducted. Trajectory data and meta-data of the wild-type (WT) protein and the highly amyloidogenic variant L55P-TTR represent the test case for the data warehouse. CONCLUSIONS: Web and grid services, especially pre-defined data mining services that can run on or 'near' the data repository of the data warehouse, are likely to play a pivotal role in the analysis of molecular dynamics unfolding data.

Relaxation of surface tension in the free-surface boundary layers of simple Lennard-Jones liquids

In this paper we use molecular dynamics to answer a classical question: how does the surface tension on a liquid/gas interface appear? After defining surface tension from the first principles and performing several consistency checks, we perform a dynamic experiment with a single simple liquid nanodroplet. At time zero, we remove all molecules of the interfacial layer of molecules, creating a fresh bare interface with the bulk arrangement of molecules. After that the system evolves towards equilibrium, and the expected surface tension is re-established. We found that the system relaxation consists of three distinct stages. First, the mechanical balance is quickly re-established. During this process the notion of surface tension is meaningless. In the second stage, the surface tension equilibrates, and the density profile broadens to a value which we call “intrinsic” interfacial width. During the third stage, the density profile continues to broaden due to capillary wave excitations, which does not however affect the surface tension.We have observed this scenario for monatomic Lennard-Jones (LJ) liquid as well as for binary LJ mixtures at different temperatures, monitoring a wide range of physical observables.

Relaxation of surface tension in the liquid-solid interfaces of Lennard-Jones liquids

We have established the surface tension relaxation time in the liquid-solid interfaces of Lennard-Jones (LJ) liquids by means of direct measurements in molecular dynamics (MD) simulations. The main result is that the relaxation time is found to be almost independent of the molecular structures and viscosity of the liquids (at seventy-fold change) used in our study and lies in such a range that in slow hydrodynamic motion the interfaces are expected to be at equilibrium. The implications of our results for the modelling of dynamic wetting processes and interpretation of dynamic contact angle data are discussed.