Grid warehousing of molecular dynamics protein unfolding data

With the increasing awareness of protein folding disorders, the explosion of genomic information, and the need for efficient ways to predict protein structure, protein folding and unfolding has become a central issue in molecular sciences research. Molecular dynamics computer simulations are increasingly employed to understand the folding and unfolding of proteins. Running protein unfolding simulations is computationally expensive and finding ways to enhance performance is a grid issue on its own. However, 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. This paper describes efforts to provide a grid-enabled data warehouse for protein unfolding data. We outline the challenge and present first results in the design and implementation of the data warehouse.

Synthesis and crystal structure of a novel octa-aqua bridged star-shaped Ni4K complex

A penta-nuclear. star-shaped hetero-metallic compound containing a unique Ni4KO8 core has been synthesized. The X-ray single crystal structure analysis reveals that in the complex, [K(Ni(LH)(2))(4)(OH2)(8)](Br)(ClO4)(8)center dot 11H(2)O (LH=(CH3)(2)HN+(CH2)(3)N=CHC6H4O-) the eight coordinate central K+ ion is encapsulated by four terminal [Ni(LH)(2)](2+) units through the double water bridges between K+ and each Ni(II) ions.

Synthesis, characterization and molecular structure of 1,1′ bis (diphenyl phosphino ferrocene) dioxide complex of the cis-uranyl dichloride

A 1,1' bis(diphenyl phosphino ferrocene) dioxide complex of the uranyl dichloride was synthesized and characterized by elemental analysis, H-1, P-31{H-1} NMR and X-ray diffraction methods. The structure of the compound shows that the uranium(VI) ion is surrounded by four oxygen and two chlorine atoms in an octahedral geometry. Two oxygen atoms from the bis (diphenyl phosphino ferrocene) dioxide and two chlorine atoms form a square planar arrangement. Two uranyl oxygen atoms occupy the axial positions. The bis(diphenyl phosphino ferrocene) dioxide ligand acts as a bidentate chelating ligand with a bite angle of 82.90(16)degrees around the uranyl group. The two chlorine atoms are mutually cis with a CI-U-Cl angle of 97.75(7)degrees.

Synthesis, crystal structures and magnetic properties of a phenoxo-bridged dinuclear Cu(II) complex and a dicyanamide bridged novel molecular rectangle based on it

The phenoxo-bridged dinuclear Cu-II complex [Cu2L2-(NCNCN)(2)] (1) and the dicyanamide-bridged molecular rectangle [Cu4L4{mu(1,5)-(NCNCN)(2)}]center dot(ClO4)(2)(H2O)(2) (2) were synthesized using the tridentate reduced Schiff-base ligand HL {2-[(2-dimethylamino-ethylamino) methyl] phenol}. The complexes were characterized by X-ray structural analyses and variable-temperature magnetic susceptibility measurements. Complex 2 was formed through the joining of the phenoxo-bridged dinuclear Cu2O2 cores of 1 via the mu(1,5)-bridging mode of dicyanamide. The structural properties of the Cu2O2 cores in two complexes are significantly different. The geometry of the copper ions is distorted trigonal bipyramid in 1 but is nearly square-pyramidal in 2. These differences have a marked effect on the magnetic properties of two compounds. Although both are antiferromagnetically coupled, the coupling constants (J = -185.2 and -500.9 cm(-1) for 1 and 2, respectively) differ considerably.

X‑ray crystal structure of rac-[Ru(phen)2dppz]2+ with d(ATGCAT)2 shows enantiomer orientations and water ordering

We report an atomic resolution X-ray crystal structure containing both enantiomers of rac-[Ru(phen)2dppz]2+ with the d-(ATGCAT)2 DNA duplex (phen = phenanthroline; dppz = dipyridophenazine). The first example of any enantiomeric pair crystallized with a DNA duplex shows different orientations of the Λ and Δ binding sites, separated by a clearly defined structured water monolayer. Job plots show that the same species is present in solution. Each enantiomer is bound at a TG/CA step and shows intercalation from the minor groove. One water molecule is directly located on one phenazine N atom in the Δ-enantiomer only.

Mg/Ca partitioning between aqueous solution and aragonite mineral: a molecular dynamics study

We have calculated the concentrations of Mg in the bulk and surfaces of aragonite CaCO3 in equilibrium with aqueous solution, based on molecular dynamics simulations and grand-canonical statistical mechanics. Mg is incorporated in the surfaces, in particular in the (001) terraces, rather than in the bulk of aragonite particles. However, the total Mg content in the bulk and surface of aragonite particles was found to be too small to account for the measured Mg/Ca ratios in corals. We therefore argue that most Mg in corals is either highly metastable in the aragonite lattice, or is located outside the aragonite phase of the coral skeleton, and we discuss the implications of this finding for Mg/Ca paleothermometry.

Spectroelectrochemical study of complexes [Mo(CO)2(h3-allyl)(adiimine)(NCS)] (a-diimine ¼ bis(2,6-dimethylphenyl)- acenaphthenequinonediimine and 2,20-bipyridine) exhibiting different molecular structure and redox reactivity

The redox properties and reactivity of [Mo(CO)2(η3-allyl)(α-diimine)(NCS)] (α-diimine = bis(2,6-dimethylphenyl)-acenaphthenequinonediimine (2,6-xylyl-BIAN) and 2,2′-bipyridine (bpy)) were studied using cyclic voltammetry and IR/UV–Vis spectroelectrochemistry. [Mo(CO)2(η3-allyl)(2,6-xylyl-BIAN)(NCS)] was shown by X-ray crystallography to have an asymmetric (B-type) conformation. The extended aromatic system of the strong π-acceptor 2,6-xylyl-BIAN ligand stabilises the primary 1e−-reduced radical anion, [Mo(CO)2(η3-allyl)(2,6-xylyl-BIAN•−)(NCS)]−, that can be reduced further to give the solvento anion [Mo(CO)2(η3-allyl)(2,6-xylyl-BIAN)(THF)]−. The initial reduction of [Mo(CO)2(η3-allyl)(bpy)(NCS)] in THF at ambient temperature results in the formation of [Mo(CO)2(η3-allyl)(bpy)]2 by reaction of the remaining parent complex with [Mo(CO)2(η3-allyl)(bpy)]− produced by dissociation of NCS− from [Mo(CO)2(η3-allyl)(bpy•−)(NCS)]−. Further reduction of the dimer [Mo(CO)2(η3-allyl)(bpy)]2 restores [Mo(CO)2(η3-allyl)(bpy)]−. In PrCN at 183 K, [Mo(CO)2(η3-allyl)(2,6-xylyl-BIAN•−)(NCS)]− converts slowly to 2e−-reduced [Mo(CO)2(η3-allyl)(2,6-xylyl-BIAN)(PrCN)]− and free NCS−. At room temperature, the reduction path in PrCN involves mainly the dimer [Mo(CO)2(η3-allyl)(bpy)]2; however, the detailed course of the reduction within the spectroelectrochemical cell is complicated and involves a mixture of several unassigned products. Finally, it has been shown that the five-coordinate anion [Mo(CO)2(η3-allyl)(bpy)]− promotes in THF reduction of CO2 to CO and formate via the formation of the intermediate [Mo(CO)2(η3-allyl)(bpy)(O2CH)] and its subsequent reduction.

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.

Structure and viscoelasticity of an electrorheological fluid in oscillatory shear: computer simulation investigation

The three-dimensional molecular dynamics simulation method has been used to study the dynamic responses of an electrorheological (ER) fluid in oscillatory shear. The structure and related viscoelastic behaviour of the fluid are found to be sensitive to the amplitude of the strain. With the increase of the strain amplitude, the structure formed by the particles changes from isolated columns to sheet-like structures which may be perpendicular or parallel to the oscillating direction. Along with the structure evolution, the field-induced moduli decrease significantly with an increase in strain amplitude. The viscoelastic behaviour of the structures obtained in the cases of different strain amplitudes was examined in the linear response regime and an evident structure dependence of the moduli was found. The reason for this lies in the anisotropy of the arrangement of the particles in these structures. Short-range interactions between the particles cannot be neglected in determining the viscoelastic behaviour of ER fluids at small strain amplitude, especially for parallel sheets. The simulation results were compared with available experimental data and good agreement was reached for most of them.

On the structure and dynamics of Indian monsoon depressions

ERA-Interim reanalysis data from the past 35 years have been used with a newly-developed feature tracking algorithm to identify Indian monsoon depressions originating in or near the Bay of Bengal. These were then rotated, centralised and combined to give a fully three-dimensional 106-depression composite structure – a considerably larger sample than any previous detailed study on monsoon depressions and their structure. Many known features of depression structure are confirmed, particularly the existence of a maximum to the southwest of the centre in rainfall and other fields, and a westward axial tilt in others. Additionally, the depressions are found to have significant asymmetry due to the presence of the Himalayas; a bimodal mid-tropospheric potential vorticity core; a separation into thermally cold- (~–1.5K) and neutral- (~0K) cores near the surface with distinct properties; and that the centre has very large CAPE and very small CIN. Variability as a function of background state has also been explored, with land/coast/sea, diurnal, ENSO, active/break and Indian Ocean Dipole contrasts considered. Depressions are found to be markedly stronger during the active phase of the monsoon, as well as during La Niña. Depressions on land are shown to be more intense and more tightly constrained to the central axis. A detailed schematic diagram of a vertical cross-section through a composite depression is also presented, showing its inherent asymmetric structure.

Mode-specific chemisorption of CH4 on Pt{110}-(1 x 2) explored by first-principles molecular dynamics

The chemisorption of CH4 on Pt{110}-(1 x 2) has been studied by vibrational analysis of the reaction pathway defined by the potential energy surface and, in time reversal, by first-principles molecular dynamics simulations of CH4 associative desorption, with the electronic structure treated explicitly using density functional theory. We find that the symmetric stretch vibration ν1 is strongly coupled to the reaction coordinate; our results therefore provide a firm theoretical basis for recently reported state-resolved reactivity measurements, which show that excitation of the ν1 normal mode is the most efficient way to enhance the reaction probability

Bond-selective energy redistribution in the chemisorption of CH3D and CD3H on Pt{1 1 0}-(1 2): a first-principles molecular dynamics study

We have investigated the chemisorption of CH3D and CD3H on Pt{11 0}-(1 2) by performing first-principles molecular dynamics simulations of the recombinative desorption of CH3D (from adsorbed methyl and deuterium) and of CD3H (from adsorbed trideuteromethyl and hydrogen). Vibrational analysis of the symmetry adapted internal coordinates of the desorbing molecules shows that excitation of the single C– D (C–H) bond in the parent molecule is strongly correlated with energy excess in the reaction coordinate. The results of the molecular dynamics simulations are consistent with observed mode- and bond-specific reactivity measurements for chemisorption of methane and its isotopomers on platinum and nickel surfaces.

Crystallographic structure and substrate-binding interactions of the molybdate-binding protein of the phytopathogen Xanthomonas axonopodis pv. citri

In Xanthomonas axonopodis pv. citri (Xac or X citri), the modA gene codes for a periplasmic protein (ModA) that is capable of binding molybdate and tungstate as part of the ABC-type transporter required for the uptake of micronutrients. In this study, we report the crystallographic structure of the Xac ModA protein with bound molybdate. The Xac ModA structure is similar to orthologs with known three-dimensional structures and consists of two nearly symmetrical domains separated by a hinge region where the oxyanion-binding site lies. Phylogenetic analysis of different ModA orthologs based on sequence alignments revealed three groups of molybdate-binding proteins: bacterial phytopathogens, enterobacteria and soil bacteria. Even though the ModA orthologs are segregated into different groups, the ligand-binding hydrogen bonds are mostly conserved, except for Archaeglobus fulgidus ModA. A detailed discussion of hydrophobic interactions in the active site is presented and two new residues, Ala(38) and Ser(151), are shown to be part of the ligand-binding pocket. (c) 2007 Elsevier B.V All rights reserved.

The crystal structure of the leptospiral hypothetical protein LIC12922 reveals homology with the periplasmic chaperone SurA

Leptospirosis is a world spread zoonosis caused by members of the genus Leptospira. Although leptospires were identified as the causal agent of leptospirosis almost 100 years ago, little is known about their biology, which hinders the development of new treatment and prevention strategies. One of the several aspects of the leptospiral biology not yet elucidated is the process by which outer membrane proteins (OMPs) traverse the periplasm and are inserted into the outer membrane. The crystal structure determination of the conserved hypothetical protein LIC12922 from Leptospira interrogans revealed a two domain protein homologous to the Escherichia coli periplasmic chaperone SurA. The LIC12922 NC-domain is structurally related to the chaperone modules of E. coli SurA and trigger factor, whereas the parvulin domain is devoid of peptidyl prolyl cis-trans isomerase activity. Phylogenetic analyses suggest a relationship between LIC12922 and the chaperones PrsA, PpiD and SurA. Based on our structural and evolutionary analyses, we postulate that LIC12922 is a periplasmic chaperone involved in OMPs biogenesis in Leptospira spp. Since LIC12922 homologs were identified in all spirochetal genomes sequenced to date, this assumption may have implications for the OMPs biogenesis studies not only in leptospires but in the entire Phylum Spirochaetes. (C) 2010 Elsevier Inc. All rights reserved.

Electronic properties of liquid hydrogen fluoride: A sequential quantum mechanical/Born-Oppenheimer molecular dynamics approach

The electronic properties of liquid hydrogen fluoride (HF) were investigated by carrying out sequential quantum mechanics/Born-Oppenheimer molecular dynamics. The structure of the liquid is in good agreement with recent experimental information. Emphasis was placed on the analysis of polarisation effects, dynamic polarisability and electronic excitations in liquid HF. Our results indicate an increase in liquid phase of the dipole moment (similar to 0.5 D) and isotropic polarisability (5%) relative to their gas-phase values. Our best estimate for the first vertical excitation energy in liquid HF indicates a blue-shift of 0.4 +/- 0.2 eV relative to that of the gas-phase monomer (10.4 eV). (C) 2010 Elsevier B.V. All rights reserved.