Monte Carlo and molecular dynamics simulations of methane in potassium montmorillonite clay hydrates at elevated pressures and temperatures

The structure and dynamics of methane in hydrated potassium montmorillonite clay have been studied under conditions encountered in sedimentary basin and compared to those of hydrated sodium montmorillonite clay using computer simulation techniques. The simulated systems contain two molecular layers of water and followed gradients of 150 barkm-1 and 30 Kkm-1 up to a maximum burial depth of 6 km. Methane particle is coordinated to about 19 oxygen atoms, with 6 of these coming from the clay surface oxygen. Potassium ions tend to move away from the center towards the clay surface, in contrast to the behavior observed with the hydrated sodium form. The clay surface affinity for methane was found to be higher in the hydrated K-form. Methane diffusion in the two-layer hydrated K-montmorillonite increases from 0.39×10-9 m2s-1 at 280 K to 3.27×10-9 m2s-1 at 460 K compared to 0.36×10-9 m2s-1 at 280 K to 4.26×10-9 m2s-1 at 460 K in Na-montmorillonite hydrate. The distributions of the potassium ions were found to vary in the hydrates when compared to those of sodium form. Water molecules were also found to be very mobile in the potassium clay hydrates compared to sodium clay hydrates. © 2004 Elsevier Inc. All All rights reserved.

Crystal structure and physical properties of the new linear chain compound [Cu(1,2-bis(tetrazol-1-yl)ethane)3](ClO4)2

The synthesis and crystal structure of a novel one-dimensional Cu(II) compound [Cu(1,2-bis(tetrazol-1-yl)ethane)3](ClO4)2 are described. The single-crystal X-ray structure determination was carried out at 298 K. The molecular structure consists of a linear chain in which the Cu(II) ions are linked by three N4,N4' coordinating bis(tetrazole) ligands in syn conformation. The Cu(II) ions are in a Jahn-Teller distorted octahedral environment (Cu(1)-N(11)=2.034(2) Å, Cu(1)-N(21)=2.041(2) Å and Cu(1)-N(31)=2.391(2) Å). The Cu⋯Cu separations are 7.420(3) Å.

Histidine hydrogen bonding in MHC at pH 5 and pH 7 modeled by molecular docking and molecular dynamics simulations

Hydrogen bonds play important roles in maintaining the structure of proteins and in the formation of most biomolecular protein-ligand complexes. All amino acids can act as hydrogen bond donors and acceptors. Among amino acids, Histidine is unique, as it can exist in neutral or positively charged forms within the physiological pH range of 5.0 to 7.0. Histidine can thus interact with other aromatic residues as well as forming hydrogen bonds with polar and charged residues. The ability of His to exchange a proton lies at the heart of many important functional biomolecular interactions, including immunological ones. By using molecular docking and molecular dynamics simulation, we examine the influence of His protonation/deprotonation on peptide binding affinity to MHC class II proteins from locus HLA-DP. Peptide-MHC interaction underlies the adaptive cellular immune response, upon which the next generation of commercially-important vaccines will depend. Consistent with experiment, we find that peptides containing protonated His residues bind better to HLA-DP proteins than those with unprotonated His. Enhanced binding at pH 5.0 is due, in part, to additional hydrogen bonds formed between peptide His+ and DP proteins. In acidic endosomes, protein His79β is predominantly protonated. As a result, the peptide binding cleft narrows in the vicinity of His79β, which stabilizes the peptide - HLA-DP protein complex. © 2014 Bentham Science Publishers.

Investigating the role of cholesterol in the formation of non-ionic surfactant based bilayer vesicles:thermal analysis and molecular dynamics

The aim of this research was to investigate the molecular interactions occurring in the formulation of non-ionic surfactant based vesicles composed monopalmitoyl glycerol (MPG), cholesterol (Chol) and dicetyl phosphate (DCP). In the formulation of these vesicles, the thermodynamic attributes and surfactant interactions based on molecular dynamics, Langmuir monolayer studies, differential scanning calorimetry (DSC), hot stage microscopy and thermogravimetric analysis (TGA) were investigated. Initially the melting points of the components individually, and combined at a 5:4:1 MPG:Chol:DCP weight ratio, were investigated; the results show that lower (90 C) than previously reported (120-140 C) temperatures could be adopted to produce molten surfactants for the production of niosomes. This was advantageous for surfactant stability; whilst TGA studies show that the individual components were stable to above 200 C, the 5:4:1 MPG:Chol:DCP mixture show ∼2% surfactant degradation at 140 C, compared to 0.01% was measured at 90 C. Niosomes formed at this lower temperature offered comparable characteristics to vesicles prepared using higher temperatures commonly reported in literature. In the formation of niosome vesicles, cholesterol also played a key role. Langmuir monolayer studies demonstrated that intercalation of cholesterol in the monolayer did not occur in the MPG:Chol:DCP (5:4:1 weight ratio) mixture. This suggests cholesterol may support bilayer assembly, with molecular simulation studies also demonstrating that vesicles cannot be built without the addition of cholesterol, with higher concentrations of cholesterol (5:4:1 vs 5:2:1, MPG:Chol:DCP) decreasing the time required for niosome assembly. © 2013 Elsevier B.V.

Crystal structure and magnetic behaviour of a five-coordinate iron(III) complex of pyridoxal-4-methylthiosemicarbazone

The crystal structure and magnetic properties of a penta-coordinate iron(III) complex of pyridoxal-4-methylthiosemicarbazone, [Fe(Hmthpy)Cl](CHCHSO), are reported. The synthesised ligand and the metal complex were characterised by spectroscopic methods (H NMR, IR, and mass spectroscopy), elemental analysis, and single crystal X-ray diffraction. The complex crystallises as dark brown microcrystals. The crystal data determined at 100(1) K revealed a triclinic system, space group P over(1, ¯) (Z = 2). The ONSCl geometry around the iron(III) atom is intermediate between trigonal bipyramidal and square pyramidal (t = 0.40). The temperature dependence of the magnetic susceptibility (5-300 K) is consistent with a high spin Fe(III) ion (S = 5/2) exhibiting zero-field splitting. Interpretation of these data yielded: D = 0.34(1) cm and g = 2.078(3). © 2007 Elsevier B.V. All rights reserved.

HLA-DP2 binding prediction by molecular dynamics simulations

Major histocompatibility complex (MHC) II proteins bind peptide fragments derived from pathogen antigens and present them at the cell surface for recognition by T cells. MHC proteins are divided into Class I and Class II. Human MHC Class II alleles are grouped into three loci: HLA-DP, HLA-DQ, and HLA-DR. They are involved in many autoimmune diseases. In contrast to HLA-DR and HLA-DQ proteins, the X-ray structure of the HLA-DP2 protein has been solved quite recently. In this study, we have used structure-based molecular dynamics simulation to derive a tool for rapid and accurate virtual screening for the prediction of HLA-DP2-peptide binding. A combinatorial library of 247 peptides was built using the "single amino acid substitution" approach and docked into the HLA-DP2 binding site. The complexes were simulated for 1 ns and the short range interaction energies (Lennard-Jones and Coulumb) were used as binding scores after normalization. The normalized values were collected into quantitative matrices (QMs) and their predictive abilities were validated on a large external test set. The validation shows that the best performing QM consisted of Lennard-Jones energies normalized over all positions for anchor residues only plus cross terms between anchor-residues.

Characterization of the structure of RAMP1 by mutagenesis and molecular modeling

Receptor activity modifying proteins (RAMPs) are a family of single-pass transmembrane proteins that dimerize with G-protein-coupled receptors. They may alter the ligand recognition properties of the receptors (particularly for the calcitonin receptor-like receptor, CLR). Very little structural information is available about RAMPs. Here, an ab initio model has been generated for the extracellular domain of RAMP1. The disulfide bond arrangement (Cys 27-Cys82, Cys40-Cys72, and Cys 57-Cys104) was determined by site-directed mutagenesis. The secondary structure (a-helices from residues 29-51, 60-80, and 87-100) was established from a consensus of predictive routines. Using these constraints, an assemblage of 25,000 structures was constructed and these were ranked using an all-atom statistical potential. The best 1000 conformations were energy minimized. The lowest scoring model was refined by molecular dynamics simulation. To validate our strategy, the same methods were applied to three proteins of known structure; PDB:1HP8, PDB:1V54 chain H (residues 21-85), and PDB:1T0P. When compared to the crystal structures, the models had root mean-square deviations of 3.8 Å, 4.1 Å, and 4.0 Å, respectively. The model of RAMP1 suggested that Phe93, Tyr 100, and Phe101 form a binding interface for CLR, whereas Trp74 and Phe92 may interact with ligands that bind to the CLR/RAMP1 heterodimer. © 2006 by the Biophysical Society.

Structure, dynamics, and energetics of siRNA-cationic vector complexation:a molecular dynamics study

The design and synthesis of safe and efficient nonviral vectors for gene delivery has attracted significant attention in recent years. Previous experiments have revealed that the charge density of a polycation (the carrier) plays a crucial role in complexation and the release of the gene from the complex in the cytosol. In this work, we adopt an atomistic molecular dynamics simulation approach to study the complexation of short strand duplex RNA with six cationic carrier systems of varying charge and surface topology. The simulations reveal detailed molecular-level pictures of the structures and dynamics of the RNA-polycation complexes. Estimates for the binding free energy indicate that electrostatic contributions are dominant followed by van der Waals interactions. The binding free energy between the 8(+)polymers and the RNA is found to be larger than that of the 4(+)polymers, in general agreement with previously published data. Because reliable binding free energies provide an effective index of the ability of the polycationic carrier to bind the nucleic acid and also carry implications for the process of gene release within the cytosol, these novel simulations have the potential to provide us with a much better understanding of key mechanistic aspects of gene-polycation complexation and thereby advance the rational design of nonviral gene delivery systems.

Flexibility of short-strand RNA in aqueous solution as revealed by molecular dynamics simulation:are A′-RNA and A-RNA distinct conformational structures?

We use molecular dynamics simulations to compare the conformational structure and dynamics of a 21-base pair RNA sequence initially constructed according to the canonical A-RNA and A'-RNA forms in the presence of counterions and explicit water. Our study aims to add a dynamical perspective to the solid-state structural information that has been derived from X-ray data for these two characteristic forms of RNA. Analysis of the three main structural descriptors commonly used to differentiate between the two forms of RNA namely major groove width, inclination and the number of base pairs in a helical twist over a 30 ns simulation period reveals a flexible structure in aqueous solution with fluctuations in the values of these structural parameters encompassing the range between the two crystal forms and more. This provides evidence to suggest that the identification of distinct A-RNA and A'-RNA structures, while relevant in the crystalline form, may not be generally relevant in the context of RNA in the aqueous phase. The apparent structural flexibility observed in our simulations is likely to bear ramifications for the interactions of RNA with biological molecules (e.g. proteins) and non-biological molecules (e.g. non-viral gene delivery vectors). © CSIRO 2009.

Structure and dynamics of multiple cationic vectors-siRNA complexation by all-atomic molecular dynamics simulations

Understanding the molecular mechanism of gene condensation is a key component to rationalizing gene delivery phenomena, including functional properties such as the stability of the gene-vector complex and the intracellular release of the gene. In this work, we adopt an atomistic molecular dynamics simulation approach to study the complexation of short strand duplex RNA with four cationic carrier systems of varying charge and surface topology at different charge ratios. At lower charge ratios, polymers bind quite effectively to siRNA, while at high charge ratios, the complexes are saturated and there are free polymers that are unable to associate with RNA. We also observed reduced fluctuations in RNA structures when complexed with multiple polymers in solution as compared to both free siRNA in water and the single polymer complexes. These novel simulations provide a much better understanding of key mechanistic aspects of gene-polycation complexation and thereby advance progress toward rational design of nonviral gene delivery systems.

A molecular dynamics model of the atomic structure of dysprosium alumino-phosphate glass

Molecular dynamics (MD) has been used to identify the relative distribution of dysprosium in the phosphate glass DyAl0.30P3.05O9.62. The MD model has been compared directly with experimental data obtained from neutron diffraction to enable a detailed comparison beyond the total structure factor level. The MD simulation gives Dy ... Dy correlations at 3.80(5) and 6.40(5) angstrom with relative coordination numbers of 0.8(1) and 7.3(5), thus providing evidence of minority rare-earth clustering within these glasses. The nearest neighbour Dy-O peak occurs at 2.30 angstrom with each Dy atom having on average 5.8 nearest neighbour oxygen atoms. The MD simulation is consistent with the phosphate network model based on interlinked PO4 tetrahedra where the addition of network modifiers Dy3+ depolymerizes the phosphate network through the breakage of P-(O)-P bonds whilst leaving the tetrahedral units intact. The role of aluminium within the network has been taken into explicit account, and A1 is found to be predominantly (78 tetrahedrally coordinated. In fact all four A1 bonds are found to be to P (via an oxygen atom) with negligible amounts of Al-O-Dy bonds present. This provides an important insight into the role of Al additives in improving the mechanical properties of these glasses.

The structure and dynamics of 2-dimensional fluids in swelling clays

The interlayer pores of swelling 2:1 clays provide an ideal 2-dimensional environment in which to study confined fluids. In this paper we discuss our understanding of the structure and dynamics of interlayer fluid species in expanded clays, based primarily on the outcome of recent molecular modelling and neutron scattering studies. Counterion solvation is compared with that measured in bulk solutions, and at a local level the cation-oxygen coordination is found to be remarkably similar in these two environments. However, for the monovalent ions the contribution to the first coordination shell from the clay surfaces increases with counterion radius. This gives rise to inner-sphere (surface) complexes in the case of potassium and caesium. In this context, the location of the negative clay surface charge (i.e. arising from octahedral or tetrahedral substitution) is also found to be of major importance. Divalent cations, such as calcium, eagerly solvate to form outer-sphere complexes. These complexes are able to pin adjacent clay layers together, and thereby prevent colloidal swelling. Confined water molecules form hydrogen bonds to each other and to the clays' surfaces. In this way their local environment relaxes to close to the bulk water structure within two molecular layers of the clay surface. Finally, we discuss the way in which the simple organic molecules methane, methanol and ethylene glycol behave in the interlayer region of hydrated clays. Quasi-elastic neutron scattering of isotopically labelled interlayer CH 3OD and (CH2OD)2 in deuterated clay allows us to measure the diffusion of the CH3- and CH2-groups in both clay and liquid environments. We find that in both the one-layer methanol solvates and the two-layer glycol solvates the diffusion of the most mobile organic molecules is close to that in the bulk solution.

Structure and physical properties of [µ-Tris(1,4-bis(tetrazol-1-yl)butane-N4,N4‘')iron(II)] Bis(hexafluorophosphate), a new Fe(II) spin-crossover compound with a three-dimensional threefold interlocked crystal lattice

[μ-Tris(1,4-bis(tetrazol-1-yl)butane-N4,N4‘)iron(II)] bis(hexafluorophosphate), [Fe(btzb)3](PF6)2, crystallizes in a three-dimensional 3-fold interlocked structure featuring a sharp two-step spin-crossover behavior. The spin conversion takes place between 164 and 182 K showing a discontinuity at about T1/2 = 174 K and a hysteresis of about 4 K between T1/2 and the low-spin state. The spin transition has been independently followed by magnetic susceptibility measurements, 57Fe-Mössbauer spectroscopy, and variable temperature far and midrange FTIR spectroscopy. The title compound crystallizes in the trigonal space group P30¯(No. 147) with a unit cell content of one formula unit plus a small amount of disordered solvent. The lattice parameters were determined by X-ray diffraction at several temperatures between 100 and 300 K. Complete crystal structures were resolved for 9 of these temperatures between 100 (only low spin, LS) and 300 K (only high spin, HS), Z = 1 [Fe(btzb)3](PF  6)2:  300 K (HS), a = 11.258(6) Å, c = 8.948(6) Å, V = 982.2(10) Å3; 100 K (LS), a = 10.989(3) Å, c = 8.702(2) Å, V = 910.1(4) Å3. The molecular structure consists of octahedral coordinated iron(II) centers bridged by six N4,N4‘ coordinating bis(tetrazole) ligands to form three 3-dimensional networks. Each of these three networks is symmetry related and interpenetrates each other within a unit cell to form the interlocked structure. The Fe−N bond lengths change between 1.993(1) Å at 100 K in the LS state and 2.193(2) Å at 300 K in the HS state. The nearest Fe separation is along the c-axis and identical with the lattice parameter c.

Molecular dynamics simulation of methane in sodium montmorillonite clay hydrates at elevated pressures and temperatures

Computer simulation has been used to study the structure and dynamics of methane in hydrated sodium montmorillonite clays under conditions encountered in sedimentary basins. Systems containing approximately one, two, three and four molecular layers of water have followed gradients of 150 bar km-1 and 30Kkm-1, to a maximum burial depth of 6 km (900 bar and 460 K). Methane is coordinated to approximately 19 oxygen atoms, of which typically 6 are provided by the clay surface. Only in the three- and four-layer hydrates is methane able to leave the clay surface. Diffusion depends strongly on the porosity (water content) and burial depth: self-diffusion coefficients are in the range 0.12 × 10-9m2s-1 for water and 0.04 × 10−9m2s−1 < D < 8.64 × 10−9m2s−1 for methane. Bearing in mind that porosity decreases with burial depth, it is estimated that maximum diffusion occurs at around 3 km. This is in good agreement with the known location of methane reservoirs in sedimentary basins.