Complex temporal patterns in molecular dynamics:a direct measure of the phase-space exploration by the trajectory at macroscopic time scales

Computer simulated trajectories of bulk water molecules form complex spatiotemporal structures at the picosecond time scale. This intrinsic complexity, which underlies the formation of molecular structures at longer time scales, has been quantified using a measure of statistical complexity. The method estimates the information contained in the molecular trajectory by detecting and quantifying temporal patterns present in the simulated data (velocity time series). Two types of temporal patterns are found. The first, defined by the short-time correlations corresponding to the velocity autocorrelation decay times (â‰0.1â€ps), remains asymptotically stable for time intervals longer than several tens of nanoseconds. The second is caused by previously unknown longer-time correlations (found at longer than the nanoseconds time scales) leading to a value of statistical complexity that slowly increases with time. A direct measure based on the notion of statistical complexity that describes how the trajectory explores the phase space and independent from the particular molecular signal used as the observed time series is introduced. © 2008 The American Physical Society.

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.

Controlling protein molecular dynamics:how to accelerate folding while preserving the native state

The dynamics of peptides and proteins generated by classical molecular dynamics (MD) is described by using a Markov model. The model is built by clustering the trajectory into conformational states and estimating transition probabilities between the states. Assuming that it is possible to influence the dynamics of the system by varying simulation parameters, we show how to use the Markov model to determine the parameter values that preserve the folded state of the protein and at the same time, reduce the folding time in the simulation. We investigate this by applying the method to two systems. The first system is an imaginary peptide described by given transition probabilities with a total folding time of 1 micros. We find that only small changes in the transition probabilities are needed to accelerate (or decelerate) the folding. This implies that folding times for slowly folding peptides and proteins calculated using MD cannot be meaningfully compared to experimental results. The second system is a four residue peptide valine-proline-alanine-leucine in water. We control the dynamics of the transitions by varying the temperature and the atom masses. The simulation results show that it is possible to find the combinations of parameter values that accelerate the dynamics and at the same time preserve the native state of the peptide. A method for accelerating larger systems without performing simulations for the whole folding process is outlined.

Equation of state for water in the small compressibility region

The equation of state for dense fluids has been derived within the framework of the Sutherland and Katz potential models. The equation quantitatively agrees with experimental data on the isothermal compression of water under extrapolation into the high pressure region. It establishes an explicit relationship between the thermodynamic experimental data and the effective parameters of the molecular potential.

Sensitivity of peptide conformational dynamics on clustering of a classical molecular dynamics trajectory

We investigate the sensitivity of a Markov model with states and transition probabilities obtained from clustering a molecular dynamics trajectory. We have examined a 500 ns molecular dynamics trajectory of the peptide valine-proline-alanine-leucine in explicit water. The sensitivity is quantified by varying the boundaries of the clusters and investigating the resulting variation in transition probabilities and the average transition time between states. In this way, we represent the effect of clustering using different clustering algorithms. It is found that in terms of the investigated quantities, the peptide dynamics described by the Markov model is sensitive to the clustering; in particular, the average transition times are found to vary up to 46%. Moreover, inclusion of nonphysical sparsely populated clusters can lead to serious errors of up to 814%. In the investigation, the time step used in the transition matrix is determined by the minimum time scale on which the system behaves approximately Markovian. This time step is found to be about 100 ps. It is concluded that the description of peptide dynamics with transition matrices should be performed with care, and that using standard clustering algorithms to obtain states and transition probabilities may not always produce reliable results.

The study and characterisation of water-based latices used for can coating interiors

Water-based latices, used in the production of internal liners for beer/beverage cans, were investigated using a number of analytical techniques. The epoxy-graft-acrylic polymers, used to prepare the latices, and films, produced from those latices, were also examined. It was confirmed that acrylic polymer preferentially grafts onto higher molecular weight portions of the epoxy polymer. The amount of epoxy remaining ungrafted was determined to be 80%. This figure is higher than was previously thought. Molecular weight distribution studies were carried out on the epoxy and epoxy-g-acrylic resins. A quantitative method for determining copolymer composition using GPC was evaluated. The GPC method was also used to determine polymer composition as a function of molecular weight. IR spectroscopy was used to determine the total level of acrylic modification of the polymers and NMR was used to determine the level of grafting. Particle size determinations were carried out using transmission electron microscopy and dynamic light scattering. Levels of stabilising amine greatly affected the viscosity of the latex, particle size and amount of soluble polymer but the core particle size, as determined using TEM, was unaffected. NMR spectra of the latices produced spectra only from solvents and amine modifiers. Using solid-state CP/MAS/freezing techniques spectra from the epoxy component could be observed. FT-IR spectra of the latices were obtained after special subtraction of water. The only difference between the spectra of the latices and those of the dry film were due to the presence of the solvents in the former. A distinctive morphology in the films produced from the latices was observed. This suggested that the micelle structure of the latex survives the film forming process. If insufficient acrylic is present, large epoxy domains are produced which gives rise to poor film characteristics. Casting the polymers from organic solutions failed to produce similar morphology.

The separation of secondary oil-water dispersions in particulate beds

The literature relating to haze formation, methods of separation, coalescence mechanisms, and models by which droplets <100 μm are collected, coalesced and transferred, have been reviewed with particular reference to particulate bed coalescers. The separation of secondary oil-water dispersions was studied experimentally using packed beds of monosized glass ballotini particles. The variables investigated were superficial velocity, bed depth, particle size, and the phase ratio and drop size distribution of inlet secondary dispersion. A modified pump loop was used to generate secondary dispersions of toluene or Clairsol 350 in water with phase ratios between 0.5-6.0 v/v%.Inlet drop size distributions were determined using a Malvern Particle Size Analyser;effluent, coalesced droplets were sized by photography. Single phase flow pressure drop data were correlated by means of a Carman-Kozeny type equation. Correlations were obtained relating single and two phase pressure drops, as (ΔP2/μc)/ΔP1/μd) = kp Ua Lb dcc dpd Cine A flow equation was derived to correlate the two phase pressure drop data as, ΔP2/(ρcU2) = 8.64*107 [dc/D]-0.27 [L/D]0.71 [dp/D]-0.17 [NRe]1.5 [e1]-0.14 [Cin]0.26  In a comparison between functions to characterise the inlet drop size distributions a modification of the Weibull function provided the best fit of experimental data. The general mean drop diameter was correlated by: q_p q_p p_q /β      Γ ((q-3/β) +1) d qp = d fr  .α        Γ ((P-3/β +1 The measured and predicted mean inlet drop diameters agreed within ±15%. Secondary dispersion separation depends largely upon drop capture within a bed. A theoretical analysis of drop capture mechanisms in this work indicated that indirect interception and London-van der Waal's mechanisms predominate. Mathematical models of dispersed phase concentration m the bed were developed by considering drop motion to be analogous to molecular diffusion.The number of possible channels in a bed was predicted from a model in which the pores comprised randomly-interconnected passage-ways between adjacent packing elements and axial flow occured in cylinders on an equilateral triangular pitch. An expression was derived for length of service channels in a queuing system leading to the prediction of filter coefficients. The insight provided into the mechanisms of drop collection and travel, and the correlations of operating parameters, should assist design of industrial particulate bed coalescers.

Solubility enhancement of poorly water soluble drugs using liposome technology

The aim of this work is to investigate the various parameters that could control the encapsulation of lipophilic drugs and investigate the influence of the physical properties of poorly water-soluble drugs on bilayer loading. Initial work investigated on the solubilisation of ibuprofen, a model insoluble drug. Drug loading was assessed using HPLC and UV spectrophotometric analysis. Preliminary studies focused on the influence of bilayer composition on drug loading to obtain an optimum cholesterol concentration. This was followed up by studies investigating the effect of longer alkyl chain lipids, unsaturated alkyl chain lipids and charged lipids. The studies also focused on the effects of pH of the hydration medium and addition of the single chain surfactant a-tocopherol. The work was followed up by investigation of a range of insoluble drugs including flurbiprofen, indomethacin, sulindac, mefenamic acid, lignocaine and progesterone to investigate the influence of drugs properties and functional group on liposomal loading. The results show that no defined trend could be obtained linking the drug loading to the different drug properties including molecular weight, log P and other drug specific characteristics. However, the presence of the oppositely charged lipids improved the encapsulation of all the drugs investigated with a similar effect obtained with the substitution of the longer chain lipids. The addition of the single chain surfactant a-tocopherol resulted in enhancement of drug loading and possibly is governed by the log P of the drug candidate. Environmental scanning-electron microscopy (ESEM) was used to dynamically follow the changes in liposome morphology in real time during dehydration thereby providing a alternative assay of liposome formulation and stability. The ESEM analysis clearly demonstrated ibuprofen incorporation enhanced the stability of PC:Chol liposomes.

The influence of the chemical nature of dispersed phase on stability in oil-in-water emulsions

Some of the problems arising from the inherent instability of emulsions are discussed. Aspects of emulsion stability are described and particular attention is given to the influence of the chemical nature of the dispersed phase on adsorbed film structure and stability, Emulsion stability has been measured by a photomicrographic technique. Electrophoresis, interfacial tension and droplet rest-time data were also obtained. Emulsions were prepared using a range of oils, including aliphatic and aromatic hydrocarbons, dispersed In a solution of sodium dodecyl sulphate. In some cases a small amount of alkane or alkanol was incorporated into the oil phase. In general the findings agree with the classical view that the stability of oil-in-water emulsions is favoured by a closely packed interfacial film and appreciable electric charge on the droplets. The inclusion of non-ionic alcohol leads to enhanced stability, presumably owing to the formation of a "mixed" interfacial film which is more closely packed and probably more coherent than that of the anionic surfactant alone. In some instances differences in stability cannot he accounted for simply by differences in interfacial adsorption or droplet charge. Alternative explanations are discussed and it is postulated that the coarsening of emulsions may occur not only hy coalescence but also through the migration of oil from small droplets to larger ones by molecular diffusion. The viability of using the coalescence rates of droplets at a plane interface as a guide to emulsion stability has been researched. The construction of a suitable apparatus and the development of a standard testing procedure are described. Coalescence-time distributions may be correlated by equations similar to those presented by other workers, or by an analysis based upon the log-normal function. Stability parameters for a range of oils are discussed in terms of differences in film drainage and the natl1re of the interfacial film. Despite some broad correlations there is generally poor agreement between droplet and emulsion stabilities. It is concluded that hydrodynamic factors largely determine droplet stability in the systems studied. Consequently droplet rest-time measurements do not provide a sensible indication of emulsion stability,

Molecular regulation of iron uptake in pseudomonas aeruginosa

The microbial demand for iron is often met by the elaboration of siderophores into the surrounding medium and expression of cognate outer membrane receptors for the ferric siderophore complexes. Conditions of iron limitation, such as those encountered in vivo, cause Pseudomonas aeruginosa to express two high-affinity iron-uptake systems based on pyoverdin and pyochelin. These systems will operate both in the organism's natural habitat, soil and water, where the solubility of iron at neutral pH is extremely low, and in the human host where the availability of free iron is too low to sustain bacterial growth due to the iron-binding glycoproteins transferrin and lactoferrin. Cross-feeding and radiolabelled iron uptake experiments demonstrated that pyoverdin biosynthesis and uptake were highly heterogeneous amongst P.aeruginosa strains, that growth either in the presence of pyoverdin or pyochelin resulted in induction of specific IROMPs, and that induction of iron uptake is siderophore-specific. The P.aeruginosa Tn5 mutant PH1 is deficient in ferripyoverdin uptake and resistant to pyocin Sa, suggesting that the site of interaction of pyocin Sa is a ferripyoverdin receptor. Additional Tn5 mutants appeared to exploit different strategies to achieve pyocin Sa-resistance, involving modifications in expression of pyoverdin-mediated iron uptake, indicating that complex regulatory systems exist to enable these organisms to compete effectively for iron. Modulation of expression of IROMPs prompted a study of the mechanism of uptake of a semi-synthetic C(7) α-formamido substituted cephalosporin BRL 41897A. Sensitivity to this agent correlated with expression of the 75 kDa ferri-pyochelin receptor and demonstrated the potential of high-affinity iron uptake systems for targeting of novel antibiotics. Studies with ferri-pyoverdin uptake-deficient mutant PH1 indicated that expression of outer membrane protein G (OprG), which is usually expressed under iron-rich conditions and repressed under iron-deficient conditions, was perturbed. Attempts were made to clone the oprG gene using a degenerate probe based on the N-terminal amino acid sequence. A strongly hybridising HindIll restriction fragment was cloned and sequenced, but failed to reveal an open reading frame correspondmg to OprG. However, there appears to be good evidence that a part of the gene codmg for the hydrophilic membrane-associated ATP-binding component of a hitherto uncharacterised periplasmic- binding-protein-dependent transport system has been isolated. The full organisation and sequence of the operon, and substrate for this putative transport system, are yet: to be elucidated,

Ionic and molecular diffusion in cementitious materials

The work described in this thesis is an attempt to provide improved understanding of the effects of several factors affecting diffusion in hydrated cement pastes and to aid the prediction of ionic diffusion processes in cement-based materials. Effect of pore structure on diffusion was examined by means of comparative diffusion studies of quaternary ammonium ions with different ionic radii. Diffusivities of these ions in hydrated pastes of ordinary portland cement with or without addition of fly ash were determined by a quasi-steady state technique. The restriction of the pore geometry on diffusion was evaluated from the change of diffusivity in response to the change of ionic radius. The pastes were prepared at three water-cement ratios, 0.35, 0.50 and 0.65. Attempts were made to study the effect of surface charge or the electrochemical double layer at the pore/solution interface on ionic diffusion. An approach was to evaluate the zeta potentials of hydrated cement pastes through streaming potential measurements. Another approach was the comparative studies of the diffusion kinetics of chloride and dissolved oxygen in hydrated pastes of ordinary portland cement with addition of 0 and 20% fly ash. An electrochemical technique for the determination of oxygen diffusivity was also developed. Non-steady state diffusion of sodium potassium, chloride and hydroxyl ions in hydrated ordinary portland cement paste of water-cement ratio 0.5 was studied with the aid of computer-modelling. The kinetics of both diffusion and ionic binding were considered for the characterization of the concentration profiles by Fick's first and second laws. The effect of the electrostatic interactions between ions on the overall diffusion rates was also considered. A general model concerning the prediction of ionic diffusion processes in cement-based materials has been proposed.

The role of lipid geometry in designing liposomes for the solubilisation of poorly water soluble drugs

Liposomes are well recognised for their ability to improve the delivery of a range of drugs. More commonly they are applied for the delivery of water-soluble drugs, but given their structural attributes they can also be employed as solubilising agents for low solubility drugs as well as drug targeting agents. To further explore the potential of liposomes as solubilising agents, we have investigated the role of bilayer packaging in promoting drug solubilisation in liposome bilayers. The effect of alkyl chain length and symmetry was investigated to consider if using 'mis-matched' phospholipids could be used to create 'voids' within the bilayers, and enhance bilayer loading capacity. Lipid packing was investigated using Langmuir studies, which demonstrated that increasing the alkyl chain length enhanced lipid packing, with condensed monolayer forming, whilst asymmetric lipids formed less condensed monolayers. However this more open packing did not translate into improved drug loading, with the longer chain, condensed bilayers formed from long-chain, saturated lipids offering higher drug loading capacity. These studies demonstrate that liposomes formulated from longer chain, saturated lipids offer enhanced solubilisation capacity. However the molecular size, rather than lipophilicity, of the drug to be incorporated was also a key factor dominating bilayer incorporation efficiency.

Ion interactions with the carbon nanotube surface in aqueous solutions:understanding the molecular mechanisms

We study the molecular mechanisms of alkali halide ion interactions with the single-wall carbon nanotube surface in water by means of fully atomistic molecular dynamics simulations. We focus on the basic physical-chemical principles of ion–nanotube interactions in aqueous solutions and discuss them in light of recent experimental findings on selective ion effects on carbon nanotubes.

Water-peptide dynamics during conformational transitions

Transitions between metastable conformations of a dipeptide are investigated using classical molecular dynamics simulation with explicit water molecules. The distribution of the surrounding water at different moments before the transitions and the dynamical correlations of water with the peptide's configurational motions indicate that the water molecules represent an integral part of the molecular system during the conformational changes, in contrast to the metastable periods when water and peptide dynamics are essentially decoupled.

Identification and molecular mechanisms of the rapid tonicity-induced relocalization of the aquaporin 4 channel

The aquaporin family of integral membrane proteins is comprised of channels that mediate cellular water flow. Aquaporin 4 (AQP4) is highly expressed in the glial cells of the central nervous system and facilitates the osmotically-driven pathological brain swelling associated with stroke and traumatic brain injury. Here we show that AQP4 cell surface expression can be rapidly and reversibly regulated in response to changes of tonicity in primary cortical rat astrocytes and in transfected HEK293 cells. The translocation mechanism involves protein kinase A (PKA) activation, influx of extracellular calcium and activation of calmodulin. We identify five putative PKA phosphorylation sites and use site-directed mutagenesis to show that only phosphorylation at one of these sites, serine- 276, is necessary for the translocation response. We discuss our findings in the context of the identification of new therapeutic approaches to treating brain oedema.