Mechanisms and dynamics of protein clustering on a solid surface

A methodology for discovering the mechanisms and dynamics of protein clustering on solid surfaces is presented. In situ atomic force microscopy images are quantitatively compared to Monte Carlo simulations using cluster statistics to differentiate various models. We study lysozyme adsorption on mica as a model system and find that all surface-supported clusters are mobile, not just the monomers, with diffusion constant inversely related to cluster size. The surface monomer diffusion constant is measured to be D1∼9×10-16  cm2 s-1, such a low value being difficult to measure using other techniques.

New insights on growth mechanisms of protein clusters at surfaces: an AFM and simulation study

Despite its relevance to a wide range of technological and fundamental areas, a quantitative understanding of protein surface clustering dynamics is often lacking. In inorganic crystal growth, surface clustering of adatoms is well described by diffusion-aggregation models. In such models, the statistical properties of the aggregate arrays often reveal the molecular scale aggregation processes. We investigate the potential of these theories to reveal hitherto hidden facets of protein clustering by carrying out concomitant observations of lysozyme adsorption onto mica surfaces, using atomic force microscopy. and Monte Carlo simulations of cluster nucleation and growth. We find that lysozyme clusters diffuse across the substrate at a rate that varies inversely with size. This result suggests which molecular scale mechanisms are responsible for the mobility of the proteins on the substrate. In addition the surface diffusion coefficient of the monomer can also be extracted from the comparison between experiments and simulations. While concentrating on a model system of lysozyme-on-mica, this 'proof of concept' study successfully demonstrates the potential of our approach to understand and influence more biomedically applicable protein-substrate couples.

Predicting mesh density for adaptive modelling of the global atmosphere

The shallow water equations are solved using a mesh of polygons on the sphere, which adapts infrequently to the predicted future solution. Infrequent mesh adaptation reduces the cost of adaptation and load-balancing and will thus allow for more accurate mapping on adaptation. We simulate the growth of a barotropically unstable jet adapting the mesh every 12 h. Using an adaptation criterion based largely on the gradient of the vorticity leads to a mesh with around 20 per cent of the cells of a uniform mesh that gives equivalent results. This is a similar proportion to previous studies of the same test case with mesh adaptation every 1–20 min. The prediction of the mesh density involves solving the shallow water equations on a coarse mesh in advance of the locally refined mesh in order to estimate where features requiring higher resolution will grow, decay or move to. The adaptation criterion consists of two parts: that resolved on the coarse mesh, and that which is not resolved and so is passively advected on the coarse mesh. This combination leads to a balance between resolving features controlled by the large-scale dynamics and maintaining fine-scale features.

Tropospheric planetary wave dynamics and mixture modeling: Two preferred regimes and a regime shift

Investigation of preferred structures of planetary wave dynamics is addressed using multivariate Gaussian mixture models. The number of components in the mixture is obtained using order statistics of the mixing proportions, hence avoiding previous difficulties related to sample sizes and independence issues. The method is first applied to a few low-order stochastic dynamical systems and data from a general circulation model. The method is next applied to winter daily 500-hPa heights from 1949 to 2003 over the Northern Hemisphere. A spatial clustering algorithm is first applied to the leading two principal components (PCs) and shows significant clustering. The clustering is particularly robust for the first half of the record and less for the second half. The mixture model is then used to identify the clusters. Two highly significant extratropical planetary-scale preferred structures are obtained within the first two to four EOF state space. The first pattern shows a Pacific-North American (PNA) pattern and a negative North Atlantic Oscillation (NAO), and the second pattern is nearly opposite to the first one. It is also observed that some subspaces show multivariate Gaussianity, compatible with linearity, whereas others show multivariate non-Gaussianity. The same analysis is also applied to two subperiods, before and after 1978, and shows a similar regime behavior, with a slight stronger support for the first subperiod. In addition a significant regime shift is also observed between the two periods as well as a change in the shape of the distribution. The patterns associated with the regime shifts reflect essentially a PNA pattern and an NAO pattern consistent with the observed global warming effect on climate and the observed shift in sea surface temperature around the mid-1970s.

Nested parallelism for multi-core HPC systems using Java

Since its introduction in 1993, the Message Passing Interface (MPI) has become a de facto standard for writing High Performance Computing (HPC) applications on clusters and Massively Parallel Processors (MPPs). The recent emergence of multi-core processor systems presents a new challenge for established parallel programming paradigms, including those based on MPI. This paper presents a new Java messaging system called MPJ Express. Using this system, we exploit multiple levels of parallelism - messaging and threading - to improve application performance on multi-core processors. We refer to our approach as nested parallelism. This MPI-like Java library can support nested parallelism by using Java or Java OpenMP (JOMP) threads within an MPJ Express process. Practicality of this approach is assessed by porting to Java a massively parallel structure formation code from Cosmology called Gadget-2. We introduce nested parallelism in the Java version of the simulation code and report good speed-ups. To the best of our knowledge it is the first time this kind of hybrid parallelism is demonstrated in a high performance Java application. (C) 2009 Elsevier Inc. All rights reserved.

Multistep electrochemical reduction path of clusters [Os3(CO)10(α-diimine)]: comparison of electrochemical and photochemical Os-Os(α-diimine) bond cleavage

Electrochemical reduction of the triangular clusters [Os-3(CO)(10)(alpha-dimine)] (alpha-dimine = 2,2'-bipyridine (bpy), 2,2'-bipyrimidine (bpym)) and [Os-3(CO)(10)(mu-bpym) ReBr(CO)(3)] produces primarily the corresponding radical anions. Their stability is strongly determined by the pi acceptor ability of the reducible alpha-dimine ligand, which decreases in the order mu-bpym > bpym >> bpy. Along this series, increasing delocalisation of the odd electron density in the radical anion over the Os(alpha-dimine) chelate ring causes weakening of the axial (CO)(4)Os-Os(CO)(2)(alpha-dimine) bond and its facile cleavage for alpha-diimine = bpy. In contrast, the cluster radical anion is inherently stable for the bridging bpym ligand, the strongest pi-acceptor in the studied series. In the absence of the partial delocalisation of the unpaired electron over the Re( bpym) chelate bond, the Os-3-core of the radical anion remains intact only at low temperatures. Subsequent one-electron reduction of [Os-3(CO)(10)(bpym)](center dot-) at T = 223 K gives the open-triosmium core (= Os-3*) dianion, [Os-3*(CO)(10)(bpym)](2-). Its oxidation leads to the recovery of parent [Os-3(CO)(10)( bpym)]. At room temperature, [Os-3*( CO)(10)(bpym)](2-) is formed along a two-electron (ECE) reduction path. The chemical step (C) results in the formation of an open- core radical anion that is directly reducible at the cathodic potential of the parent cluster in the second electrochemical (E) step. In weakly coordinating tetrahydrofuran, [Os-3*(CO)(10)( bpym)](2-) rapidly attacks yet non- reduced parent cluster molecules, producing the relatively stable open- core dimer [Os-3*(CO)(10)(bpym)](2)(2-) featuring two open- triangle cluster moieties connected with an ( bpym) Os - Os( bpym) bond. In butyronitrile, [Os-3*( CO)(10)(bpym)](2-) is stabilised by the solvent and the dimer [Os-3*(CO)(10)(bpym)](2)(2-) is then mainly formed by reoxidation of the dianion on reverse potential scan. The more reactive cluster [Os-3(CO)(10)(bpy)] follows the same reduction path, as supported by spectroelectrochemical results and additional valuable evidence obtained from cyclic voltammetric scans. The ultimate process in the reduction mechanism is fragmentation of the cluster core triggered by the reduction of the dimer [Os-3*(CO)(10)(alpha- diimine)](2)(2-). The products formed are [Os-2(CO)(8)](2-) and {Os(CO)(2)(alpha- diimine)}(2). The latter dinuclear fragments constitute a linear polymeric chain [Os( CO)(2)(alpha-dimine)] n that is further reducible at the alpha-dimine ligands. For alpha-dimine = bpy, the charged polymer is capable of reducing carbon dioxide. The electrochemical opening of the triosmium core in the [Os-3( CO)(10)(alpha-dimine)] clusters exhibits several common features with their photochemistry. The same Os-alpha-dimine bond dissociates in both cases but the intimate mechanisms are different.

Heterosite effects in novel heteronuclear clusters [OS2Ru(CO)(11)(PPh3)] and Os2Ru(CO)(10)(2-acetylpyridine-N-isopropylimine)]

A new synthetic route towards the mixed-metal cluster [OS2Ru(CO)(12)] is described together with the syntheses of its PPh3 and iPr-AcPy (iPr-AcPy = 2-acetylpyridine-N-isopropylimine) derivatives. The molecular structures of the novel clusters [Os2Ru(CO)(11)(PPh3)] and [Os2Ru(CO)(10)(iPr-AcPy)] were determined on the basis of crystalline solid solutions of the Os2Ru and corresponding Os-3 species. The structures reveal that coordination of the Lewis bases occurs exclusively at the ruthenium site of [Os2Ru(CO)(12)], which is in agreement with density functional theory (DFT) calculations on several structural isomers of these compounds. According to the time-dependent DFT results, the lowest optically accessible excited state of [Os2Ru(CO)(10)(iPr-AcPy)] has a prevailing sigma(Ru-Os-2)pi*(iPr-AcPy) character, with a partial sigma sigma*(Ru-Os-2) contribution. In weakly coordinating 2-chlorobutane, the excited state has a lifetime tau = 10.4 +/- 1.2 ps and produces biradicals considerably faster than observed for [Os-3(CO)10(iPr-AcPy) (tau = 25.3 +/- 0.7ps)]. In coordinating acetonitrile, the excited state of [Os2Ru(CO)(10)(iPr-AcPy)] decays mono-exponentially with a lifetime tau = 2.1 +/- 0.2 ps. In contrast to [Os-3(CO)(10)(iPr-AcPy)] that forms biradicals as the main primary photoproduct even in strongly coordinating solvents, zwitterion formation from the solvated lowest excited state is observed for the heterometallic cluster. This is concluded from time-resolved absorption studies in the microsecond time domain. Due to the lower tendency of the coordinatively unsaturated Ru+(CO)(2)(iPr-AcPy-/0) moiety to bind a Lewis base, the heteronuclear biradical and zwitterionic photoproducts live significantly shorter than their triosmium counterparts. The influence of the weaker Os-2-Ru(iPr-AcPy) bond on the redox reactivity is clearly reflected in very reactive radical anions formed upon electrochemical reduction of [Os2Ru(CO)(10)(iPr-AcPy)]. The dimer [-OS(CO)(4)-Os(CO)(4)-Ru(CO)(2)(iPr-AcPy)](2)(2-) is the only IR-detectable intermediate reduction product. The dinuclear complex [Os-2(CO)(8)](2-) and insoluble [Ru(CO)(2)(iPr-AcPy)](n), are the ultimate reduction products, proving fragmentation of the OS2Ru core.

Low-lying excited states and primary photoproducts of [Os-3(CO)(10)-(s-cis-L)] (L = cyclohexa-1,3-diene, buta-1,3-diene)] clusters studied by picosecond time-resolved UV/Vis and IR spectroscopy and by density functional theory

Combined picosecond transient absorption and time-resolved infrared studies were performed, aimed at characterising low-lying excited states of the cluster [Os-3(CO)(10)(s-cis-L)] (L= cyclohexa-1,3-diene, 1) and monitoring the formation of its photoproducts. Theoretical (DFT and TD-DFT) calculations on the closely related cluster with L=buta-1,3-diene (2') have revealed that the low-lying electronic transitions of these [Os-3(CO)(10)(s-cis-1,3-diene)] clusters have a predominant sigma(core)pi*(CO) character. From the lowest sigmapi* excited state, cluster 1 undergoes fast Os-Os(1,3-diene) bond cleavage (tau=3.3 ps) resulting in the formation of a coordinatively unsaturated primary photoproduct (1a) with a single CO bridge. A new insight into the structure of the transient has been obtained by DFT calculations. The cleaved Os-Os(1,3-diene) bond is bridged by the donor 1,3-diene ligand, compensating for the electron deficiency at the neighbouring Os centre. Because of the unequal distribution of the electron density in transient la, a second CO bridge is formed in 20 ps in the photoproduct [Os-3(CO)(8)(mu-CO)(2)- (cyclohexa-1,3-diene)] (1b). The latter compound, absorbing strongly around 630 nm, mainly regenerates the parent cluster with a lifetime of about 100 ns in hexane. Its structure, as suggested by the DFT calculations, again contains the 1,3-diene ligand coordinated in a bridging fashion. Photoproduct 1b can therefore be assigned as a high-energy coordination isomer of the parent cluster with all Os-Os bonds bridged.

An MPI-based implementation of intelligent agents on clusters

Recent research in multi-agent systems incorporate fault tolerance concepts, but does not explore the extension and implementation of such ideas for large scale parallel computing systems. The work reported in this paper investigates a swarm array computing approach, namely 'Intelligent Agents'. A task to be executed on a parallel computing system is decomposed to sub-tasks and mapped onto agents that traverse an abstracted hardware layer. The agents intercommunicate across processors to share information during the event of a predicted core/processor failure and for successfully completing the task. The feasibility of the approach is validated by implementation of a parallel reduction algorithm on a computer cluster using the Message Passing Interface.

Molecular dynamics simulations of threadlike cetyltrimethylammonium chloride micelles: effects of sodium chloride and sodium salicylate salts

We use atomistic molecular dynamics simulations to probe the effects of added sodium chloride (NaCl) and sodium salicylate (NaSal) salts on the spherical-to-threadlike micelle shape transition in aqueous solutions of cetyltrimethylammonium chloride (CTAC) surfactants. Long threadlike micelles are found to be unstable and break into spherical micelles at low concentrations or NaCl, but remain stable for 20 ns above a threshold value of [NaCl] approximate to 3.0 M, which is about 2.5 times larger than the experimental salt concentration at which the transition between spherical and rodlike micelles occurs. The chloride counterions associate weakly oil the surface of the CTAC micelles with the degree of counterion dissociation decreasing slightly with increasing [NaCl] on spherical micelles, but dropping significantly on the threadlike micelles tit high [NaCl]. This effect indicates that the electrolyte ions drive the micellar shape transition by screening the electrostatic repulsions between the micellar headgroups, The aromatic salicylate counterions, on the other hand, penetrate inside the micelle with their hydrophilic groups staying in the surfactant headgroup region and the hydrophobic groups partially embedded into the hydrophobic core of the micelle. The strong association of the salicylate ions with the surfactant headgroups leads to dense packing of the surfactant molecules, which effectively reduces the surface area per surfactant, and increases intramicellar ordering of the surfactant headgroups, favoring the formation of long threadlike micelles. Simulation predictions of the geometric and electrostatic properties of the spherical and threadlike micelles are in good agreement with experiments.

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.

Supramolecular structures of enzyme clusters

The structural characterization of subtilisin mesoscale clusters, which were previously shown to induce supramolecular order in biocatalytic self-assembly of Fmocdipeptides, was carried out by synchrotron small-angle X-ray, dynamic, and static light scattering measurements. Subtilisin molecules self-assemble to form supramolecular structures in phosphate buffer solutions. Structural arrangement of subtilisin clusters at 55 degrees Centigrade was found to vary systematically with increasing enzyme concentration. Static light scattering measurements showed the cluster structure to be consistent with a fractal-like arrangement, with fractal dimension varying from 1.8 to 2.6 with increasing concentration for low to moderate enzyme concentrations. This was followed by a structural transition around the enzyme concentration of 0.5 mg mL-1 to more compact structures with significantly slower relaxation dynamics, as evidenced by dynamic light scattering measurements. These concentration-dependent supramolecular enzyme clusters provide tunable templates for biocatalytic self-assembly.

High resolution Northern Hemisphere wintertime mid-latitude dynamics during the Last Glacial Maximum

Hourly winter weather of the Last Glacial Maximum (LGM) is simulated using the Community Climate Model version 3 (CCM3) on a globally resolved T170 (75 km) grid. Results are compared to a longer LGM climatological run with the same boundary conditions and monthly saves. Hourly-scale animations are used to enhance interpretations. The purpose of the study is to explore whether additional insights into ice age conditions can be gleaned by going beyond the standard employment of monthly average model statistics to infer ice age weather and climate. Results for both LGM runs indicate a decrease in North Atlantic and increase in North Pacific cyclogenesis. Storm trajectories react to the mechanical forcing of the Laurentide Ice Sheet, with Pacific storms tracking over middle Alaska and northern Canada, terminating in the Labrador Sea. This result is coincident with other model results in also showing a significant reduction in Greenland wintertime precipitation – a response supported by ice core evidence. Higher-temporal resolution puts in sharper focus the close tracking of Pacific storms along the west coast of North America. This response is consistent with increased poleward heat transport in the LGM climatological run and could help explain “early” glacial warming inferred in this region from proxy climate records. Additional analyses shows a large increase in central Asian surface gustiness that support observational inferences that upper-level winds associated with Asian- Pacific storms transported Asian dust to Greenland during the LGM.