Ligand diffusion on protein surface observed in molecular dynamics simulation

The process of binding of small ligands to dihydrofolate reductase protein has been investigated using all-atom molecular dynamics simulations. The existence of a mechanism that facilitates the search of the binding site by the ligand is demonstrated. The mechanism consists of ligand diffusing on the protein’s surface. It has been discussed in the literature before, but has not been explicitly confirmed for realistic molecular systems. The strength of this nonspecific binding is roughly estimated and found to be essential for the binding kinetics.

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.

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.

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.

21-residue peptide's dynamics at and between elementary structural transitions

Elementary conformational changes of the backbone of a 21-residue peptide A5(A3RA)3A are studied using molecular dynamics simulations in explicit water. The processes of the conformational transitions and the regimes of stationary fluctuations between them are investigated using minimal perturbations of the system. The perturbations consist of a few degrees rotation of the velocity of one of the systems' atoms and keep the system on the same energy surface. It is found that (i) the system dynamics is insignificantly changed by the perturbations in the regimes between the transitions; (ii) it is very sensitive to the perturbations just before the transitions that prevents the peptide from making the transitions; and (iii) the perturbation of any atom of the system, including distant water molecules is equally effective in preventing the transition. The latter implies strongly collective dynamics of the peptide and water during the transitions.

Characterisation and surface-profiling techniques for composite particles produced by dry powder coating in pharmaceutical drug delivery

The production of composite particles using dry powder coating is a one-step, environmentally friendly, process for the fabrication of particles with targeted properties and favourable functionalities. Diverse functionalities, such flowability enhancement, content uniformity, and dissolution, can be developed from dry particle coating. In this review, we discuss the particle functionalities that can be tailored and the selection of characterisation techniques relevant to understanding their molecular basis. We address key features in the powder blend sampling process and explore the relevant characterisation techniques, focussing on the functionality delivered by dry coating and on surface profiling that explores the dynamics and surface characteristics of the composite blends. Dry particle coating is a solvent- and heat-free process that can be used to develop functionalised particles. However, assessment of the resultant functionality requires careful selection of sensitive analytical techniques that can distinguish particle surface changes within nano and/or micrometre ranges.

Self-assembly of trehalose molecules on a lysozyme surface: the broken glass hypothesis

To help understand how sugar interactions with proteins stabilise biomolecular structures, we compare the three main hypotheses for the phenomenon with the results of long molecular dynamics simulations on lysozyme in aqueous trehalose solution (0.75 M). We show that the water replacement and water entrapment hypotheses need not be mutually exclusive, because the trehalose molecules assemble in distinctive clusters on the surface of the protein. The flexibility of the protein backbone is reduced under the sugar patches supporting earlier findings that link reduced flexibility of the protein with its higher stability. The results explain the apparent contradiction between different experimental and theoretical results for trehalose effects on proteins.

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.

Deflection of jets discharged into a reservoir with a free surface

Deflections of jets discharged into a reservoir with a free surface are investigated numerically. The jets are known to deflect towards either side of the free surface or the bottom, whose direction is not determined uniquely in some experimental conditions, i.e. there are multiple stable states realizable in the same condition. The origin of the multiple stable states is explored by utilizing homotopy transformations in which the top boundary of the reservoir is transformed from a rigid to a free boundary and also the location of the outlet throat is continuously moved from mid-height to the top. We depicted bifurcation diagrams of the flow compiling the data of numerical simulations, from which we identified the origin as an imperfect pitchfork bifurcation, and obtained an insight into the mechanism for the direction to be determined. The parameter region where such multiple stable states are possible is also delimited. © 2011 The Japan Society of Fluid Mechanics and IOP Publishing Ltd.

A decade of silicone hydrogel development:surface properties, mechanical properties, and ocular compatibility

Since the initial launch of silicone hydrogel lenses, there has been a considerable broadening in the range of available commercial material properties. The very mobile silicon–oxygen bonds convey distinctive surface and mechanical properties on silicone hydrogels, in which advantages of enhanced oxygen permeability, reduced protein deposition, and modest frictional interaction are balanced by increased lipid and elastic response. There are now some 15 silicone hydrogel material variants available to practitioners; arguably, the changes that have taken place have been strongly influenced by feedback based on clinical experience. Water content is one of the most influential properties, and the decade has seen a progressive rise from lotrafilcon-A (24%) to efrofilcon-A (74%). Moduli have decreased over the same period from 1.4 to 0.3 MPa, but not solely as a result of changes in water content. Surface properties do not correlate directly with water content, and ingenious approaches have been used to achieve desirable improvements (e.g., greater lubricity and lower contact angle hysteresis). This is demonstrated by comparing the hysteresis value of the earliest (lotrafilcon-A, >40°) and most recent (delefilcon-A, <10°) coated silicone hydrogels. Although wettability is important, it is not of itself a good predictor of ocular response because this involves a much wider range of physicochemical and biochemical factors. The interference of the lens with ocular dynamics is complex leading separately to tissue–material interactions involving anterior and posterior lens surfaces. The biochemical consequences of these interactions may hold the key to a greater understanding of ocular incompatibility and end of day discomfort.

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.

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.

The effect of contact lens wear on dynamic ocular surface temperature

Aim: To determine the dynamic emitted temperature changes of the anterior eye during and immediately after wearing different materials and modalities of soft contact lenses. Method: A dynamic, non-contact infrared camera (Thermo-Tracer TH7102MX, NEC San-ei) was used to record the ocular surface temperature (OST) in 48 subjects (mean age 21.7 ± 1.9 years) wearing: lotrafilcon-A contact lenses on a daily wear (LDW; n = 8) or continuous wear (LCW; n = 8) basis; balafilcon-A contact lenses on a daily wear (BDW; n = 8) or continuous wear (BCW; n = 8) basis; etafilcon-A contact lenses on a daily disposable regimen (EDW; n = 8); and no lenses (controls; n = 8). OST was measured continuously five times, for 8 s after a blink, following a minimum of 2 h wear and immediately following lens removal. Absolute temperature, changes in temperature post-blink and the dynamics of temperature changes were calculated. Results: OST immediately following contact lens wear was significantly greater compared to non-lens wearers (37.1 ± 1.7 °C versus 35.0 ± 1.1 °C; p < 0.005), predominantly in the LCW group (38.6 ± 1.0 °C; p < 0.0001). Lens surface temperature was highly correlated (r = 0.97) to, but lower than OST (by -0.62 ± 0.3 °C). There was no difference with modality of wear (DW 37.5 ± 1.6 °C versus CW 37.8 ± 1.9 °C; p = 0.63), but significant differences were found between etafilcon A and silicone hydrogel lens materials (35.3 ± 1.1 °C versus 37.5 ± 1.5 °C; p < 0.0005). Ocular surface cooling following a blink was not significantly affected by contact lens wear with (p = 0.07) or without (p = 0.47) lenses in situ. Conclusions: Ocular surface temperature is greater with hydrogel and greater still with silicone hydrogel contact lenses in situ, regardless of modality of wear. The effect is likely to be due to the thermal transmission properties of a contact lens. © 2004 British Contact Lens Association. Published by Elsevier Ltd. All rights reserved.

Molecular dynamics simulation of brittle fracture in silicon

The fracture process involves converting potential energy from a strained body into surface energy, thermal energy, and the energy needed to create lattice defects. In dynamic fracture, energy is also initially converted into kinetic energy. This paper uses molecular dynamics (MD) to simulate brittle frcture in silicon and determine how energy is converted from potential energy (strain energy) into other forms.

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.