Structure and dynamics of a confined ionic liquid. topics of relevance to dye-sensitized solar cells

The behavior of a model ionic liquid (IL) confined between two flat parallel walls was studied at various interwall distances using computer simulations. The results focus both on structural and dynamical properties. Mass and charge density along the confinement axis reveal a structure of layers parallel to the walls that leads to an oscillatory profile in the electrostatic potential. Orientational correlation functions indicate that cations at the interface orient tilted with respect to the surface and that any other orientational order is lost thereafter. The diffusion coefficients of the ions exhibit a maximum as a function of the confinement distance, a behavior that results from a combination of the structure of the liquid as a whole and a faster molecular motion in the vicinity of the walls. We discuss the relevance of the present results and elaborate on topics that need further attention regarding the effects of ILs in the functioning of IL-based dye-sensitized solar cells.

Thermodynamics, structure, and dynamics in room temperature ionic liquids: The case of 1-butyl-3-methyl imidazolium hexafluorophosphate ([bmim][PF6])

A detailed investigation of the phase diagram of 1-butyl-3-methyl imidazolium hexafluorophosphate ([bmim][PF6]) is presented on the basis of a wide set of experimental data accessing thermodynamic, structural, and dynamical properties of this important room temperature ionic liquid (RTIL). The combination of quasi adiabatic, continuous calorimetry, wide angle neutron and X-ray diffraction, and quasi elastic neutron scattering allows the exploration of many novel features of this material. Thermodynamic and microscopic structural information is derived on both glassy and crystalline states and compared with results that recently appeared in the literature allowing direct information to be obtained on the existence of two crystalline phases that were not previously characterized and confirming the view that RTILs show a substantial degree of order (even in their amorphous states), which resembles the crystalline order. We highlight a strong connection between structure and dynamics, showing the existence of three temperature ranges in the glassy state across which both the spatial correlation and the dynamics change. The complex crystalline polymorphism in [bmim][PF6] also is investigated; we compare our findings with the corresponding findings for similar RTILs. These results provide a strong experimental basis for the exploration of the features of the phase diagram of RTILs and for the further study of longer alkyl chain salts.

Structural basis of substrate selectivity in the glycerol-3-phosphate: phosphate antiporter GlpT.

Major facilitators represent the largest superfamily of secondary active transporter proteins and catalyze the transport of an enormous variety of small solute molecules across biological membranes. However, individual superfamily members, although they may be architecturally similar, exhibit strict specificity toward the substrates they transport. The structural basis of this specificity is poorly understood. A member of the major facilitator superfamily is the glycerol-3-phosphate (G3P) transporter (GlpT) from the Escherichia coli inner membrane. GlpT is an antiporter that transports G3P into the cell in exchange for inorganic phosphate (Pi). By combining large-scale molecular-dynamics simulations, mutagenesis, substrate-binding affinity, and transport activity assays on GlpT, we were able to identify key amino acid residues that confer substrate specificity upon this protein. Our studies suggest that only a few amino acid residues that line the transporter lumen act as specificity determinants. Whereas R45, K80, H165, and, to a lesser extent Y38, Y42, and Y76 contribute to recognition of both free Pi and the phosphate moiety of G3P, the residues N162, Y266, and Y393 function in recognition of only the glycerol moiety of G3P. It is the latter interactions that give the transporter a higher affinity to G3P over Pi.

Salt-bridge dynamics control substrate-induced conformational change in the membrane transporter GlpT.

Active transport of substrates across cytoplasmic membranes is of great physiological, medical and pharmaceutical importance. The glycerol-3-phosphate (G3P) transporter (GlpT) of the E. coli inner membrane is a secondary active antiporter from the ubiquitous major facilitator superfamily that couples the import of G3P to the efflux of inorganic phosphate (Pi) down its concentration gradient. Integrating information from a novel combination of structural, molecular dynamics simulations and biochemical studies, we identify the residues involved directly in binding of substrate to the inward-facing conformation of GlpT, thus defining the structural basis for the substrate-specificity of this transporter. The substrate binding mechanism involves protonation of a histidine residue at the binding site. Furthermore, our data suggest that the formation and breaking of inter- and intradomain salt bridges control the conformational change of the transporter that accompanies substrate translocation across the membrane. The mechanism we propose may be a paradigm for organophosphate:phosphate antiporters.

Structural motifs of cholesterol nanoparticles

The growth sequence of gas-phase cholesterol clusters (Ch(N)) with up to N=36 molecules has been investigated by atomistic simulation based on an empirical force field model. The results of long annealings from high temperature show that the geometric motifs characterizing the structure of pure cholesterol crystals already appear in nanometric aggregates. In all clusters molecules tend to align along a common direction. For cluster sizes above the smallest ones, dispersion interactions among the hydrocarbon body and tails of cholesterol cooperate with hydrogen bonding to give rise to a bilayer structure. Analysis of snapshots from the annealing shows that the condensation of hydrogen bonds into a connected network of rings and chains is an important step in the self-organization of cholesterol clusters. The effect of solvation on the equilibrium properties of medium-size aggregates is investigated by short molecular dynamics simulations for the N=30 and N=40 clusters in water at near ambient conditions and in supercritical carbon dioxide at T=400 K.

Molecular Modeling on Inhibitor Complexes and Active-Site Dynamics of Cytochrome P450 C17, a Target for Prostate Cancer Therapy

A molecular model for the P450 enzyme cytochrome P450 C17 (CYP17) is presented based on sequence alignments of multiple template structures and homology modeling. This enzyme plays a central role in the biosynthesis of testosterone and is emerging as a major target in prostate cancer, with the recently developed inhibitor abiraterone currently in advanced clinical trials. The model is described in detail, together with its validation, by providing structural explanations to available site-directed mutagenesis data. The CYP17 molecule in this model is in the form of a triangular prism, with an edge of similar to 55 angstrom and a thickness of similar to 37 angstrom. It is predominantly helical, comprising 13 alpha helices interspersed by six 3(10) helices and 11 beta-sheets. Multinanosecond molecular dynamics simulations in explicit solvent have been carried out, and principal components analysis has been used to reveal the details of dynamics around the active site. Coarse-grained methods have also been used to verify low-frequency motions, which have been correlated with active-site gating. The work also describes the results of docking synthetic inhibitors, including the drug abiraterone and the natural substrate pregnenolone, in the CYP17 active site together with molecular dynamics simulations on the complexes. (C) 2010 Elsevier Ltd. All rights reserved.

Transonic Aeroelastic Stability Predictions Under the Influence of Structural Variability

Flutter prediction as currently practiced is almost always deterministic in nature, based on a single structural model that is assumed to represent a fleet of aircraft. However, it is also recognized that there can be significant structural variability, even for different flights of the same aircraft. The safety factor used for flutter clearance is in part meant to account for this variability. Simulation tools can, however, represent the consequences of structural variability in the flutter predictions, providing extra information that could be useful in planning physical tests and assessing risk. The main problem arising for this type of calculation when using high-fidelity tools based on computational fluid dynamics is the computational cost. The current paper uses an eigenvalue-based stability method together with Euler-level aerodynamics and different methods for propagating structural variability to stability predictions. The propagation methods are Monte Carlo, perturbation, and interval analysis. The feasibility of this type of analysis is demonstrated. Results are presented for the Goland wing and a generic fighter configuration.

Spontaneous nucleation of structural defects in inhomogeneous ion chains

Structural defects in ion crystals can be formed during a linear quench of the transverse trapping frequency across the mechanical instability from a linear chain to a zigzag structure. The density of defects after the sweep can be conveniently described by the Kibble-Zurek mechanism (KZM). In particular, the number of kinks in the zigzag ordering can be derived from a time-dependent Ginzburg-Landau equation for the order parameter, here the zigzag transverse size, under the assumption that the ions are continuously laser cooled. In a linear Paul trap, the transition becomes inhomogeneous, since the charge density is larger in the center and more rarefied at the edges. During the linear quench, the mechanical instability is first crossed in the center of the chain, and a front, at which the mechanical instability is crossed during the quench, is identified that propagates along the chain from the center to the edges. If the velocity of this front is smaller than the sound velocity, the dynamics become adiabatic even in the thermodynamic limit and no defect is produced. Otherwise, the nucleation of kinks is reduced with respect to the case in which the charges are homogeneously distributed, leading to a new scaling of the density of kinks with the quenching rate. The analytical predictions are verified numerically by integrating the Langevin equations of motion of the ions, in the presence of a time-dependent transverse confinement. We argue that the non-equilibrium dynamics of an ion chain in a Paul trap constitutes an ideal scenario to test the inhomogeneous extension of the KZM, which lacks experimental evidence to date.

Structural Defects in Ion Chains by Quenching the External Potential: The Inhomogeneous Kibble-Zurek Mechanism

The nonequilibrium dynamics of an ion chain in a highly anisotropic trap is studied when the transverse trap frequency is quenched across the value at which the chain undergoes a continuous phase transition from a linear to a zigzag structure. Within Landau theory, an equation for the order parameter, corresponding to the transverse size of the zigzag structure, is determined when the vibrational motion is damped via laser cooling. The number of structural defects produced during a linear quench of the transverse trapping frequency is predicted and verified numerically. It is shown to obey the scaling predicted by the Kibble-Zurek mechanism, when extended to take into account the spatial inhomogeneities of the ion chain in a linear Paul trap.

Binding modes of diketo-acid inhibitors of HIV-1 integrase: A comparative molecular dynamics simulation study

HIV-1 integrase (IN) has become an attractive target since drug resistance against HIV-1 reverse transcriptase (RT) and protease (PR) has appeared. Diketo acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN: however the action mechanism is not well understood. Here, to study the inhibition mechanism of DKAs we performed 10 ns comparative molecular dynamics simulations on HIV-1 IN bound with three most representative DMA inhibitors: Shionogi inhibitor, S-1360 and two Merck inhibitors L-731,988 and L-708,906. Our simulations show that the acidic part of S-1360 formed salt bridge and cation-pi interactions with Lys159. In addition, the catalytic Glu152 in S-1360 was pushed away from the active site to form an ion-pair interaction with Arg199. The Merck inhibitors can maintain either one or both of these ion-pair interaction features. The difference in potencies of the DMA inhibitors is thus attributed to the different binding modes at the catalytic site. Such structural information at atomic level, not only demonstrates the action modes of DMA inhibitors but also provides a novel starting point for structural-based design of HIV-1 IN inhibitors.

Molecular dynamics studies of the STAT3 homodimer:DNA complex: relationships between STAT3 mutations and protein-DNA recognition.

Signal Transducers and Activators of Transcription (STAT) proteins are a group of latent cytoplasmic transcription factors involved in cytokine signaling. STAT3 is a member of the STAT family and is expressed at elevated levels in a large number of diverse human cancers and is now a validated target for anticancer drug discovery.. Understanding the dynamics of the STAT3 dimer interface, accounting for both protein-DNA and protein-protein interactions, with respect to the dynamics of the latent unphosphorylated STAT3 monomer, is important for designing potential small-molecule inhibitors of the activated dimer. Molecular dynamics (MD) simulations have been used to study the activated STAT3 homodimer:DNA complex and the latent unphosphorylated STAT3 monomer in an explicit water environment. Analysis of the data obtained from MD simulations over a 50 ns time frame has suggested how the transcription factor interacts with DNA, the nature of the conformational changes, and ways in which function may be affected. Examination of the dimer interface, focusing on the protein-DNA interactions, including involvement of water molecules, has revealed the key residues contributing to the recognition events involved in STAT3 protein-DNA interactions. This has shown that the majority of mutations in the DNA-binding domain are found at the protein-DNA interface. These mutations have been mapped in detail and related to specific protein-DNA contacts. Their structural stability is described, together with an analysis of the model as a starting-point for the discovery of novel small-molecule STAT3 inhibitors.

Phase Behaviour, Interactions, and Structural Studies of (Amines+Ionic Liquids) Binary Mixtures

We present a study on the phase equilibrium behaviour of binary mixtures containing two 1-alkyl-3-methylimidazolium bis{(trifluoromethyl)sulfonyl}imide-based ionic liquids, [Cnmim] [NTf2] (n=2 and 4), mixed with diethylamine or triethylamine as a function of temperature and composition using different experimental techniques. Based on this work, two systems showing an LCST and one system with a possible hourglass shape are measured. Their phase behaviours are then correlated and predicted by using Flory–Huggins equations and the UNIQUAC method implemented in Aspen. The potential of the COSMO-RS methodology to predict the phase equilibria was also tested for the binary systems studied. However, this methodology is unable to predict the trends obtained experimentally, limiting its use for systems involving amines in ionic liquids. The liquid-state structure of the binary mixture ([C2mim] [NTf2]+diethylamine) is also investigated by molecular dynamics simulation and neutron diffraction. Finally, the absorption of gaseous ethane by the ([C2mim][NTf2]+diethylamine) binary mixture is determined and compared with that observed in the pure solvents.

Dynamics of the Indonesian seas circulation. Part I – The influence of bottom topography on temperature and salinity distributions

The influence of bottom topography on the distribution of temperature and salinity in the Indonesian seas region has been studied with a high-resolution model based on the Princeton Ocean Model. One of the distinctive properties of the model is an adequate reproduction of all major topographic features in the region by the model bottom relief. The three major routes of flow of Pacific water through the region have been identified. The western route follows the flow of North Pacific Water through the Sulawesi Sea, Makassar Strait, Flores Sea, and Banda Sea. This is the main branch of the Indonesian Throughflow. The eastern routes follow the flow of South Pacific water through the eastern Indonesian seas. This water enters the region either through the Halmahera Sea or by flowing to the north around Halmahera Island into the Morotai Basin and then into the Maluku Sea. A deep southward flow of South Pacific Water fills the Seram Sea below 1200 m through the Lifamatola Passage. As it enters the Seram Sea, this overflow turns eastward at depths greater than 2000 m, then upwells in the eastern part of the Seram Sea before returning westward at ~1500-2000 m. The flow continues westward across the Seram Sea, spreading to greater depths before entering the Banda Sea at the Buru-Mangole passage. It is this water that shapes the temperature and salinity of the deep Banda Sea. Topographic elevations break the Indonesian seas region down into separate basins. The difference in the distributions of potential temperature, ?, and salinity, S, in adjacent basins is primarily due to specific properties of advection of ? and S across a topographic rise. By and large, the topographic rise blocks deep flow between basins whereas water shallower than the depth of the rise is free to flow between basins. To understand this process, the structure of simulated fields of temperature and salinity has been analyzed. To identify a range of advected ? or S, special sections over the sills with isotherms or isohalines and isotachs of normal velocity have been considered. Following this approach the impact of various topographic rises on the distribution of ? and S has been identified. There are no substantial structural changes of potential temperature and salinity distributions between seasons, though values of some parameters of temperature and salinity distributions, e.g., magnitudes of maxima and minima, can change. It is shown that the main structure of the observed distributions of temperature and salinity is satisfactorily reproduced by the model throughout the entire domain.

Plume characteristics and dynamics of UV and IR laser-desorbed oligonucleotides

Laser desorption of dye-tagged oligonucleotides was studied using laser-induced fluorescence imaging. Desorption with ultra violet (UV) and infra-red (IR) lasers resulted in forward directed plumes of molecules. In the case of UV desorption, the initial shot desorbed approximately seven-fold more material than subsequent shots. In contrast, the initial shot in IR desorption resulted in the ejection of less material compared to subsequent shots and these plumes had a component directed along the path of the laser. Thermal equilibrium of the molecules in the plume was achieved after approximately 25 µs with a spread in molecular temperature which was described by a modified Maxwell-Boltzmann equation.