Molecular dynamics simulation of organic-inorganic nanocomposites: Layering behavior and interlayer structure of organoclays

Isothermal-isobaric (NPT) molecular dynamics simulation has been performed to investigate the layering behavior and structure of nanoconfined quaternary alkylammoniums in organoclays. This work is focused on systems consisting of two clay layers and a number of alkylammoniums, and involves the use of modified Dreiding force field. The simulated basal spacings of organoclays agree satisfactorily with the experimental results in the literature. The atomic density profiles in the direction normal to the clay surface indicate that the alkyl chains within the interlayer space of montmorillonite exhibit an obvious layering behavior. The headgroups of long alkyl chains are distributed within two layers close to the clay surface, whereas the distributions of methyl and methylene groups are strongly dependent on the alkyl chain length and clay layer charge. Monolayer, bilayer, and pseudo-trilayer structures are found in organoclays modified with single long alkyl chains, which are identical to the structural models based on the measured basal spacings. A pseudo-quadrilayer structure, for the first time to our knowledge, is also identified in organoclays with double long alkyl chains. In the mixture structure of paraffin-type and multilayer, alkyl chains do not lie flat within a single layer but interlace, and also jump to the next layer in pseudo-trilayer as well as next nearest layer in pseudo-quadrilayer.

Molecular dynamics simulation of the structural and dynamic properties of dioctadecyldimethyl ammoniums in organoclays

The structural and dynamic properties of dioctadecyldimethylammoniums (DODDMA) intercalated into 2:1 layered clays are investigated using isothermal-isobaric (NPT) molecular dynamics (MD) simulation. The simulated results are in reasonably good agreement with the available experimental measurements, such as X-ray diffraction (XRD), atom force microscopy (AFM), Fourier transform infrared (FTIR), and nuclear magnetic resonance (NMR) spectroscopies. The nitrogen atoms are found to be located mainly within two layers close to the clay surface whereas methylene groups form a pseudoquadrilayer structure. The results of tilt angle and order parameter show that interior two-bond segments of alkyl chains prefer an arrangement parallel to the clay surface, whereas the segments toward end groups adopt a random orientation. In addition, the alkyl chains within the layer structure lie almost parallel to the clay surface whereas those out of the layer structure are essentially perpendicular to the surface. The trans conformations are predominant in all cases although extensive gauche conformations are observed, which is in agreement with previous simulations on n-butane. Moreover, an odd-even effect in conformation distributions is observed mainly along the chains close to the head and tail groups. The diffusion constants of both nitrogen atoms and methylene groups in these nanoconfined alkyl chains increase with the temperature and methelene position toward the tail groups.

Molecular dynamics simulations of the hydrophobin SC3 at a hydrophobic/hydrophilic interface

Hydrophobins are small (similar to 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/ hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily P-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in a-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.

Molecular dynamics study of MscL interactions with a curved lipid bilayer

Mechanosensitivity is a ubiquitous sensory mechanism found in living organisms. The simplest known mechanotransducing mechanism is found in bacteria in the form of the mechanosensitive membrane channel of large conductance, MscL. This channel has been studied extensively using a variety of methods at a functional and structural level. The channel is gated by membrane tension in the lipid bilayer alone. It serves as a safety valve protecting bacterial cells against hypoosmotic shock. MscL of Escherichia coli embedded in bilayers composed of asymmetric amounts of single-tailed and double-tailed lipids has been shown to gate spontaneously, even in the absence of membrane tension. To gain insight into the effect of the lipid membrane composition and geometry on MscL structure, a fully solvated, all-atom model of MscL in a stress-free curved bilayer composed of double- and single-tailed lipids was studied using a 9.5-ns molecular dynamics simulation. The bilayer was modeled as a domed structure accommodating the asymmetric composition of the monolayers. During the course of the simulation a spontaneous restructuring of the periplasmic loops occurred, leading to interactions between one of the loops and phospholipid headgroups. Previous experimental studies of the role of the loops agree with the observation that opening starts with a restructuring of the periplasmic loop, suggesting an effect of the curved bilayer. Because of limited resources, only one simulation of the large system was performed. However, the results obtained suggest that through the geometry and composition of the bilayer the protein structure can be affected even on short timescales.

Simulation of glassy state relaxations in polymers: A static analysis of methyl group and methoxy group rotation in poly(vinyl methyl ether)

We show that the simple quasi-static technique, also called the adiabatic mapping technique, can be used to determine the energetics of rotation of methyl and methoxy groups in amorphous poly(vinyl methyl ether) even though the latter process is too slow to be amenable to direct molecular dynamics simulation. For the methyl group rotation, we find that the mean and standard deviation of the simulated rotational barrier heights agree well with experimental data from quasi-elastic neutron scattering. In the case of the methoxy groups we find that just 4% of the groups contribute more than 90% of the observed dielectric relaxation strength. The groups which make the most contribution are those which, by virtue of their particular conformation and local environment, have two alternative positions of similar energy.

Complexes of 1-thia-4,7,10-triazacyclododecane ([12]aneN(3)S) - Synthesis, structure and stability constants

The 12-membered macrocyclic ligand 1-thia-4,7, 10-triazacyclododecane ([12]aneN(3)S) has been synthesised, although upon crystallization from acetonitrile a product in which carbon dioxide had added to one secondary amine in the macrocyclic ring (H[12]aneN(3)SCO(2). H2O) was isolated and subsequently characterised by X-ray crystallography. The protonation constants for [12]aneN(3)S and stability constants with Zn(II), Pb(II), Cd(II) and Cu(II) have been determined either potentiometrically or spectrophotometrically in aqueous solution, and compared with those measured or reported for the ligands 1-oxa-4,7,10-triazacyclododecane ([12]aneN(3)O) and 1,4,7,10-tetraazacyclododecane ([12]aneN(4)). The magnitudes of the stability constants are consistent with trends observed previously for macrocyclic ligands as secondary amine donors are replaced with oxygen and thioether donors although the stability constant for the [Hg([12]aneN(4))](2+) complex has been estimated from an NMR experiment to be at least three orders of magnitude larger than reported previously. Zinc(II), mercury(II), lead(II), copper(II) and nickel(II) complexes of [12]aneN(3)S have been isolated and characterised by X-ray crystallography. In the case of copper(II), two complexes [Cu([12]aneN(3)S)(H2O)](ClO4)(2) and [Cu-2([12]aneN(3)S)(2)(OH)(2)](ClO4)(2) were isolated, depending on the conditions employed. Molecular mechanics calculations have been employed to investigate the relative metal ion size preferences of the [3333], asym-[2424] and sym-[2424] conformation isomers. The calculations predict that the asym-[2424] conformer is most stable for M-N bond lengths in the range 2.00-2.25 Angstrom whilst for the larger metal ions the [3333] conformer is dominant. The disorder seen in the structure of the [Zn([12]aneN(3)S)(NO3)](+) complex is also explained by the calculations. (C) 1999 Elsevier Science Ltd. All rights reserved.

The Ni-II, Hg-II and Cu-II complexes of 12-membered-ring mixed-donor macrocycles

The structures of diaqua(1,7-dioxa-4-thia-10-azacyclododecane)nickel dinitrate, [Ni(C8H17NO2S)(H2O)(2)](NO3)(2), (I), bis(nitrato-O,O')(1,4,7-trioxa-10-azacyclododecane)mercury, [Hg(NO3)(2)(C8H17NO3)], (II), and aqua(nitrato-O)(1-oxa-4,7,10-triazacyclododecane)copper nitrate, [Cu(NO3)(C8H19N3O)(H2O)]NO3, (III), reveal each macrocycle binding in a tetradentate manner. The conformations of the ligands in (I) and (III) are the same and distinct from that identified for (II). These differences are in agreement with molecular-mechanics predictions of ligand conformation as a function of metal-ion size.

Interfacial interactions and structure of organic-inorganic nanohybrids

Understanding the interfacial interactions and structure is important to better design and application of organic-inorganic nanohybrids. This paper presents our recent molecular dynamic studies on organoclays and polymer nanocomposites, including the layering behavior of organoclays, structural and dynamic properties of dioctadecyldimethyl ammoniums in organoclays, and interfacial interactions and structure of polyurethane nanocomposites. The results demonstrate that the layering behaviors of organoclays are closely related to the chain length of quaternary alkyl ammoniums and cation exchangeable capacity of clays. In addition to typical layered structures such as monolayer, bilayer and pseudo-trilayer, a pseudo-quadrilayer structure was also observed in organoclays modified with dioctadecyldimethyl ammoniums (DODDMA). In such a structure, alkyl chains do not lie flat within a single layer but interlace, and also jump to the next layer or even the next nearest layer. Moreover, the diffusion constants of nitrogen and methylene atoms increase with the temperature and methelene towards the tail groups. For polyurethane nanocomposite, the van der Waals interaction between apolar alkyl chains and soft segments of polyurethane predominates the interactions between organoclay and polyurethane. Different from most bulk polyurethane systems, there is no distinct phase-separated structure for the polyurethane.

N-methylation of diamino-substituted macrocyclic complexes: Intermolecular reactions

Two N-based isomeric copper(II) complexes of the macrocycle trans-6,13-dimethyl-6,13-bis(dimethylamino)1,4,8,11-tetraazacyclotetradecane (L(3)) have been synthesized and characterised spectroscopically and structurally: alpha-[CuL(3)(OH2)(2)]Cl-2, monoclinic, space group C2/m, a = 12.908(4), b = 12.433(2), c = 7.330(2) Angstrom, beta = 105.87(2)degrees, Z = 2; beta-[CuL(3)(OClO3)(2)]. 2H(2)O, monoclinic, space group P2(1)/c, a = 9.708(3), b = 9.686(3), c = 14.202(4) Angstrom, beta = 106.17(1)degrees, Z = 2. The two isomers exhibit very similar co-ordination spheres but significantly different visible electronic maxima. This difference is attributed to an intramolecular N ... H contact between the pendant dimethylamino group and an adjacent secondary amine H atom.

Structural and electron self-exchange rate variations in isomeric (hexaamine)cobalt(III/II) complexes

The syntheses and characterisation of the new macrocyclic hexaamine trans-(5(S),7(S),12(R),14(R)-tetramethyl)-1,4,8,11-tetraazacyclotetradecane-6,13-diamine (L-6) and its Co-III complex are reported. The X-ray crystal structural analyses of [CoL6]Cl-2(ClO4) [monoclinic, space group C2/c, a = 16.468(3) Angstrom, b = 9.7156(7) Angstrom, c = 15.070(3) Angstrom, beta = 119.431(8)degrees, Z = 4] and the closely related cis-diamino-substituted macrocyclic complex [CoL2](ClO4)(3) . 2H(2)O (L-2 = cis-6,13-dimethyl-1,4,8,11-tetraazacyclotetradecane-6,13-diamine) [orthorhombic, space group Pna2(1), a = 16.8220(8) Angstrom, b = 10.416(2) Angstrom, c = 14.219(3) Angstrom, Z = 4] reveal significant variations in the observed Co-N bond lengths and coordination geometries, which may be attributed to the trans or cis disposition of the pendent primary amines. The Co-III/II self-exchange electron transfer rate constants for these and other closely related hexaamines have been determined, and variations of some 2 orders of magnitude are found between pairs of trans and cis isomeric Co-III complexes.

Diverse solid-state and solution structures within a series of hexaamine dicopper(II) complexes

Mono- and dicopper(II) complexes of a series of potentially bridging hexaamine ligands have been prepared and characterized in the solid state by X-ray crystallography. The crystal structures of the following Cu-II complexes are reported: [Cu(HL3)](ClO4)(3), C11H31Cl3CuN6O12, monoclinic, P2(1)/n, a = 8.294(2) Angstrom, b = 18.364(3) Angstrom, c = 15.674(3) Angstrom, beta = 94.73(2)degrees, Z = 4; {[Cu-2(L-4)(CO3)](2)}(ClO4)(4). 4H(2)O, C40H100Cl4Cu4N12O26, triclinic, P (1) over bar, a = 9.4888(8) Angstrom, b=13.353(1) Angstrom,. c = 15.329(1) Angstrom, alpha = 111.250(7)degrees, beta = 90.068(8)degrees, gamma = 105.081(8)degrees, Z=1; [Cu-2(L-5)(OH2)(2)](ClO4)(4), C(13)H(36)Cl(4)Cu(2)Z(6)O(18), monoclinic, P2(1)/c, a = 7.225(2) Angstrom. b = 8.5555(5) Angstrom, c = 23.134(8) Angstrom, beta = 92.37(1)degrees, Z = 2; [Cu-2(L-6)(OH2)(2)](ClO4)(4). 3H(2)O, C14H44Cl4Cu2N6O21, monoclinic, P2(1)/a, a = 15.204(5) Angstrom, b = 7.6810(7) Angstrom, c = 29.370(1) Angstrom, beta = 100.42(2)degrees, Z = 4. Solution spectroscopic properties of the bimetallic complexes indicate that significant conformational changes occur upon dissolution, and this has been probed with EPR spectroscopy and molecular mechanics calculations.

Crown ether appended cyclam receptors for cationic guests

A series of crown ether appended macrocyclic amines has been prepared comprising benzo-12-crown-4, benzo-15-crown-5, or benzo-18-crown-6 attached to a diamino-substituted cyclam. The Co-III complexes of these three receptors have been prepared and characterized spectroscopically and structurally. Crystal structures of each receptor in complex with an alkali metal ion and structures of the benzo-12-crown-4 and benzo-15-crown-5-receptors without guest ions are reported. 2D NMR and molecular mechanics modeling have been used to examine conformational variations upon guest ion complexation. Addition of cations to these receptors results in an appreciable anodic shift in the Co-III:II 11 redox potential, even in aqueous solution, but little cation selectivity is observed. Evidence for complex formation has been corroborated by Na-23 and Li-7 NMR spectroscopy and electrospray mass spectrometry.

Synthesis and structural properties of patellamide a derivatives and their copper(II) compounds

The synthesis, characterization and copper(II) coordination chemistry of three new cyclic peptide ligands, PatJ(1) (cyclo-(Ile -Thr- (Gly)Thz-lle-Thr(Gly)Thz)), PatJ(2) (cyclo-(Ile-Thr(Gly)Thz-(D)-Ile-Thr-(Gly)Thz)), and PatL (cyclo-(Ile-Ser-(Gly)Thz-Ile-Ser(Gly)Thz)) are reported. All of these cyclic peptides and PatN (cyclo-(Ile-Ser(Gly)Thz-Ile-Thr-(Gly)Thz)) are derivatives of patellamide A and have a [24]azacrown-8 macrocyclic structure. All four synthetic cyclic peptides have two thiazole rings but, in contrast to patellamide A, no oxazoline rings. The molecular structure of PatJ1, determined by X-ray crystallography, has a saddle conformation with two close-to-co-parallel thiazole rings, very similar to the geometry of patellamide D. The two coordination sites of PatJ1 with thiazole-N and amide-N donors are each well preorganized for transition metal ion binding. The coordination of copper(II) was monitored by UV/Vis spectroscopy, and this reveals various (meta)stable mono- and dinuclear copper(II) complexes whose stoichiometry was confirmed by mass spectra. Two types of dinuclear copper(II) complexes, [Cu-2(H4L)(OH2)(n)](2+) (n = 6, 8) and [Cu-2(H4L)(OH2)(n)] (n=4, 6; L=PatN, PatL, PatJ1, PatJ2) have been identified and analyzed structurally by EPR spectroscopy and a combination of spectra simulations and molecular mechanics calculations (MM-EPR). The four structures are similar to each other and have a saddle conformation, that is, derived from the crystal structure of PatJ(1) by a twist of the two thiozole rings. The small but significant structural differences are characterized by the EPR simulations.

The influence of cis/trans isomerism on the physical properties of a cyano-bridged dinuclear mixed valence complex

The cyano-bridged complexes cis-[L14CoIIINCFeII(CN)5]– and cis-[L14CoIIINCFeIII(CN)5] (L14= 6-methyl-1,4,8,11-tetraazacyclotetradecan-6-amine) are prepared and characterised spectroscopically, electrochemically and structurally: Na{cis-[L14CoIIINCFeII(CN)5]}·9H2O, monoclinic space group P21/c, a= 14.758(3), b= 10.496(1), c= 19.359(3) , = 92.00(2)°, Z= 4; cis-[L14CoIIINCFeIII(CN)5]·4H2O, orthorhombic space group P212121, a= 9.492(1), b= 14.709(2), c= 18.760(3) , Z= 4. In both complexes, the pendant amine is cis to the bridging cyanide ligand. An analysis of the metal-to-metal charge transfer (MMCT) transition in these systems with Hush theory has been carried out. This has revealed that the change in the configuration of the macrocycle both decreases the redox isomer energy difference (E1/2) and increases the reorganisational energy () of the cis-[L14CoIIINCFeII(CN)5]– complex with respect to the trans-[L14CoIIINCFeII(CN)5]– complex, the result being that both isomers display an MMCT transition of similar energy. The variation in redox isomer energy differences of the configurational isomers has been related to strain energy differences by molecular mechanics analysis of the [CoL14Cl]2+/+ precursor complexes.