Molecular dynamics simulation of nonsteroidal antiinflammatory drugs,naproxen and relafen,in a lipid bilayer membrane

Naproxen and relafen, as nonsteroidal antiinflammatory drugs, were simulated in neutral and charged forms and their effects on a lipid bilayer membrane were investigated by molecular dynamics simulation using Groningen machine for chemical simulations software (GROMACS). Simulation of 10 systems was performed, which included different dosages of the drug molecules, naproxen and Relafen, in charged and neutral forms, and a mixture of naproxen and Relafen in neutral forms. The effects of the mixture and the individual drugs' dosages on membrane properties, such as electrostatic potential, order parameter, diffusion coefficients, and hydrogen bond formation, were analyzed. Hydration of the drugs in the membrane system was investigated using radial distribution function analysis. Using fully hydrated dimyristoylphosphatidylcholine (DMPC) as a reference system, 128 lipid molecules and water molecules were simulated exclusively, and the same simulation technique was performed on 10 other systems, including drug mixtures and a DMPC membrane. Angular distributions of lipid chains of the membrane were calculated, and the effects of the drug insertion and chain orientation in the membrane were evaluated. © 2013 Wiley Periodicals, Inc.     

A molecular-dynamics study of lipid bilayers: effects of the hydrocarbon chain length on permeability

In this paper, we investigate the effects of the hydrocarbon chain length of lipid molecules on the permeation process of small molecules through lipid bilayers. We perform molecular-dynamics simulations using three kinds of lipid molecules with different chain length: dilauroylphosphatidylcholine, dimyristoylphosphatidylcholine, and dipalmiltoylphosphatidylcholine. Free-energy profiles of O2, CO, NO, and water molecules are calculated by means of the cavity insertion Widom method and the probability ratio method. We show that the lipid membrane with longer chains has a larger and wider energy barrier. The local diffusion coefficients of water across the bilayers are also calculated by the force autocorrelation function method and the velocity autocorrelation function method. The local diffusion coefficients in the bilayers are not altered significantly by the chain length. We estimate the permeability coefficients of water across the three membranes according to the solubility-diffusion model; we find that the water permeability decreases modestly with increasing chain length of the lipid molecules.     

Dynamics of polymers in a particle-based mesoscopic solvent

We study the dynamics of flexible polymer chains in solution by combining multiparticle-collision dynamics (MPCD), a mesoscale simulation method, and molecular-dynamics simulations. Polymers with and without excluded-volume interactions are considered. With an appropriate choice of the collision time step for the MPCD solvent, hydrodynamic interactions build up properly. For the center-of-mass diffusion coefficient, scaling with respect to polymer length is found to hold already for rather short chains. The center-of-mass velocity autocorrelation function displays a long-time tail which decays algebraically as (Dt)(-3/2) as a function of time t, where D is the diffusion coefficient. The analysis of the intramolecular dynamics in terms of Rouse modes yields excellent agreement between simulation data and results of the Zimm model for the mode-number dependence of the mode-amplitude correlation functions.     

Water dynamics in silica nanopores: the self-intermediate scattering functions

The dynamics of water molecules confined in approximately cylindrical silica nanopores is investigated using molecular simulation. The model systems are pores of diameter varying between 20 and 40 ? containing water at room temperature and at full hydration, prepared using grand canonical Monte Carlo simulation. Water dynamics in these systems is studied via molecular dynamics simulation. The results of the basic characterization of these systems have been reported in A. A. Milischuk and B. M. Ladanyi [J. Chem. Phys. 135, 174709 (2011)]. The main focus of the present study is the self-intermediate scattering function (ISF), F(S)(Q, t), of water hydrogens, the observable in quasi-elastic neutron scattering experiments. We investigate how F(S)(Q, t) depends on the pore diameter, the direction and magnitude of the momentum transfer Q, and the proximity of water molecules to the silica surface. We also study the contributions to F(S)(Q, t) from rotational and translational motions of water molecules and the extent of rotation-translation coupling present in F(S)(Q, t). We find that F(S)(Q, t) depends strongly on the pore diameter and that this dependence is due mainly to the contributions to the ISF from water translational motion and can be attributed to the decreased mobility of water molecules near the silica surface. The relaxation rate depends on the direction of Q and is faster for Q in the axial than in the radial direction. As the magnitude of Q increases, this difference diminishes but does not disappear. We find that its source is mainly the anisotropy in translational diffusion at low Q and in molecular reorientation at higher Q values.     

Structure and dynamics of a fluid phase bilayer on a solid support as observed by a molecular dynamics computer simulation

Simulations of a 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine lipid bilayer interacting with a solid surface of hydroxylated nanoporous amorphous silica have been carried out over a range of lipid-solid substrate distances. The porous solid surface allowed the water layer to dynamically adjust its thickness, maintaining equal pressures above and below the membrane bilayer. Qualitative estimates of the force between the surfaces leads to an estimated lipid-silicon distance in very good agreement with the results of neutron scattering experiments. Detailed analysis of the simulation at the separation suggested by experiment shows that for this type of solid support the water layer between surfaces is very narrow, consisting only of bound waters hydrating the lipid head groups and hydrophilic silica surface. The reduced hydration, however, has only minor effects on the head group hydration, the orientation of water molecules at the interface, and the membrane dipole potential. Whereas these structural properties were not sensitive to the presence of the solid substrate, the calculated diffusion coefficient for translation of the lipid molecules was altered significantly by the silica surface.     

A method for DNA detection in a microchannel: fluid dynamics phenomena and optimization of microchannel structure

A simple DNA diagnosis method using microfluidics has been developed which requires simple and straightforward procedures such as injection of sample and probe DNA solutions. This method takes advantage of the highly accurate control of fluids in microchannels, and is superior to DNA microarray diagnosis methods due to its simplicity, highly quantitative determination, and high-sensitivity. The method is capable of detecting DNA hybridization for molecules as small as a 20 mer. This suggests the difference in microfluidic behavior between single strand DNA (ssDNA) and double stranded DNA (dsDNA). In this work, influence of both the inertial force exerted on DNA molecules and the diffusion of DNA molecules was investigated. Based on the determination of these parameters for both ssDNA and dsDNA by experiments, a numerical model describing the phenomena in the microchannel was designed. Computational simulation results using this model were in good agreement with previously reported experimental results. The simulation results showed that appropriate selection of the analysis point and the design of microchannel structure are important to bring out the diffusion and inertial force effects suitably and increase the sensitivity of the detection of DNA hybridization, that is, the analytical performance of the microfluidic DNA chip.     

Surface response methodology for the study of supported membrane formation

We report on the investigations of the formation of the tethered lipid bilayer by vesicle deposition on amine-functionalized surfaces. The tethered bilayer was created by the deposition of egg-PC vesicles containing 1,2-distearoyl-sn-glycero-3-phosphoethanolamine-poly-(ethyleneglycol)-N-hydroxysuccinimide as anchoring molecules on an amine-coated surface. This approach is an easy route for the formation of a biomimetic-supported membrane. A Doelhert experimental design was applied to determine the conditions leading to the formation of a continuous and defect-free tethered bilayer on different surfaces (gold and glass). Doehlert designs allow modeling of the experimental responses by second-order polynomial equations as a function of experimental factors. Four factors expected to influence bilayer formation were studied: the lipid concentration in the vesicle suspension, the mass percentage of anchoring molecules in the vesicles, the contact time between the vesicles and the surface, and the resting time of the membrane after buffer rinse. The optimization of the membrane preparation parameters was achieved by monitoring lipid assembly formation using surface plasmon resonance spectroscopy on gold and by fluorescence recovery after photobleaching on glass. Three characteristic responses were systematically measured: the bilayer thickness, the lipid diffusion coefficient, and the lipid mobile fraction. The simultaneous inspection of the three characteristics revealed that a restricted experimental domain leads to properties that are in accordance with a bilayer presence. The factors of this domain are a lipid concentration from 0.1 to 1 mg/mL, 4-8% of anchoring molecules in the vesicles, 1-4 h of contact time between vesicles and surface, and 21-24 h of resting time after buffer rinse. Under these conditions, a membrane having a lipid mass per surface between 545 +/- 5 and 590 +/- 10 ng/cm2, a diffusion coefficient of between 2.5 +/- 0.3 x 10(-8) and 3.60 +/- 0.5 x 10(-8) cm2/s, and a mobile fraction between 94 +/- 2 and 99 +/- 1% was formed. These findings were confirmed by atomic force microscopy observations, which showed the presence of a continuous and homogeneous bilayer in the determined experimental domain. This formation procedure presents many advantages; it provides an easily obtainable biomimetic membrane model for proteins studies and offers a versatile tethered bilayer because it can be adapted easily to various types of supports.     

Brownian Motion of Lipid Molecules

Lateral and transmembrane diffusions in lipid molecules are discussed in terms of translational and rotational motions of a long cylinder. A procedure is suggested to deduce the dimensions (length and radius) of lipid molecules using experimental diffusion or viscosity data. It is shown that the hydrodynamic equation for transmembrane diffusion predicts a slower motion than for lateral diffusion, in qualitative agreement with experimental results. Further, the consequences of coupling transmembrane motion with rotational motion are discussed.     

Diffusion in the presence of periodically spaced permeable membranes

The diffusion of molecules in biological tissues and some other microheterogeneous systems is affected by the presence of permeable barriers. This leads to the slowdown of diffusion at long times as compared to barrier-free diffusion. At short times the effect of barriers is weak. In consequence, the diffusion coefficient D(t) decreases as a function of time. We derive an exact solution for the Laplace transform of D(t) for diffusion in a space separated into layers by equally spaced, parallel identical planes of arbitrary permeability. Additionally, we give an approximation to D(t) which is reasonably accurate over the whole range of the partition permeability from zero (the case of isolated layers) to infinity (the case of no barriers).     

Molecular level investigation of organization in ternary lipid bilayer: a computational approach

The differential organization of lipid components in a multicomponent membrane leads to formation of domains having diverse composition and size. Cholesterol and glycosphingolipids are known to be important components of such lateral assembly. We report here the ordering of cholesterol around ganglioside GM1 and the nature of the cluster from an all-atom simulation of a ternary lipid system. The results are compared with a binary bilayer and a pure phospholipid bilayer. The difference in molecular rearrangements in ternary and binary lipid mixture shows the role of GM1 in the rearrangement of cholesterol. Calculation of the radial distribution function, rotational reorientation, and residence time analysis of cholesterol shows that cholesterol is preferentially accumulating near gangliosides, while the lateral translational motion, rotational diffusion, and order parameter of phospholipids characterize the amount of rigidity imparted on the phospholipid bilayer.     

Molecular dynamics simulations of ternary membrane mixture: phosphatidylcholine, phosphatidic acid, and cholesterol

A ternary mixture of 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl phosphatidic acid (POPA), and cholesterol (CHOL) works effectively for a functional conformation of nicotinic acetylcholine receptor (nAChR) that can undergo agonist-induced conformation changes, but POPC alone can stabilize only a desensitized state of nAChR. To gain insights into the lipid mixture that has strong impact to nAChR functions, we performed more than 50 ns all atom molecular dynamic (MD) simulations at 303 K on a fully hydrated bilayer consisting of 240 POPC, 80 POPA, and 80 CHOL (3:1:1). The MD simulation revealed various interactions between different types of molecular pairs that ultimately regulated lipid organization. The heterogeneous interactions among three different constituents resulted in a broad spectrum of lipid properties, including extensive distributions of average area per lipid and varied lipid ordering as a function of lipid closeness to CHOL. Higher percentage of POPA than POPC had close association with CHOL, which coincided with relatively higher ordering of POPA molecules in their acyl chains near lipid head groups. Lower fraction of gauche dihedrals was also found in the same region of POPA. Although the CHOL molecules had the effects on the enhancement of surrounding lipid order, relatively low lipid order parameters and high fraction of gauche bonds were observed in the ternary mixture. Collectively, these results suggest that the dynamical structure of the ternary system could be determinant for a functional nAChR.     

Analysis of Decomposition for Structure I Methane Hydrate by Molecular Dynamics Simulation

Under multi-nodes of temperatures and pressures, microscopic decomposition mechanisms of structure I methane hydrate in contact with bulk water molecules have been studied through LAMMPS software by molecular dynamics simulation. Simulation system consists of 482 methane molecules in hydrate and 3027 randomly distributed bulk water molecules. Through analyses of simulation results, decomposition number of hydrate cages, density of methane molecules, radial distribution function for oxygen atoms, mean square displacement and coefficient of diffusion of methane molecules have been studied. A significant result shows that structure I methane hydrate decomposes from hydrate-bulk water interface to hydrate interior. As temperature rises and pressure drops, the stabilization of hydrate will weaken, decomposition extent will go deep, and mean square displacement and coefficient of diffusion of methane molecules will increase. The studies can provide important meanings for the microscopic decomposition mechanisms analyses of methane hydrate.     

Multiscale simulation reveals multiple pathways for H2 and O2 transport in a [NiFe]-hydrogenase

Hydrogenases are enzymes that catalyze the reversible conversion of hydrogen molecules to protons and electrons. The mechanism by which the gas molecules reach the active site is important for understanding the function of the enzyme and may play a role in the selectivity for hydrogen over inhibitor molecules. Here, we develop a general multiscale molecular simulation approach for the calculation of diffusion rates and determination of pathways by which substrate or inhibitor gases can reach the protein active site. Combining kinetic data from both equilibrium simulations and enhanced sampling, we construct a master equation describing the movement of gas molecules within the enzyme. We find that the time-dependent gas population of the active site can be fit to the same phenomenological rate law used to interpret experiments, with corresponding diffusion rates in very good agreement with experimental data. However, in contrast to the conventional picture, in which the gases follow a well-defined hydrophobic tunnel, we find that there is a diverse network of accessible pathways by which the gas molecules can reach the active site. The previously identified tunnel accounts for only about 60% of the total flux. Our results suggest that the dramatic decrease in the diffusion rate for mutations involving the residue Val74 could be in part due to the narrowing of the passage Val74-Arg476, immediately adjacent to the binding site, explaining why mutations of Leu122 had only a negligible effect in experiment. Our method is not specific to the [NiFe]-hydrogenase and should be generally applicable to the transport of small molecules in proteins.     

Heterogeneous and anomalous diffusion inside lipid tubules

Self-assembled lipid tubules with crystalline bilayer walls are promising candidates for controlled drug delivery vehicles on the basis of their ability to release preloaded biological molecules in a sustained manner. While a previous study has shown that the release rate of protein molecules from lipid tubules depends on the associated molecular mass, suggesting that the pertinent diffusion follows the well-known Stokes-Einstein relationship, only a few attempts have been made toward investigating the details of molecular diffusion in the tubule interior. Herein, we have characterized the diffusion rates of several molecules encapsulated in lipid tubules formed by 1,2-bis(10,12-tricosadiynoyl)-sn-glycero-3-phosphocholine (DC8,9PC) using the techniques of fluorescence recovery after photobleaching (FRAP) and fluorescence correlation spectroscopy (FCS). Our results show that the mobility of these molecules depends not only on their positions in the DC8,9PC tubules but also on their respective concentrations. While the former indicates that the interior of the DC8,9PC tubules is heterogeneous in terms of diffusion, the latter further highlights the possibility of engineering specific conditions for achieving sustained release of a "drug molecule" over a targeted period of time. In addition, our FCS results indicate that the molecular diffusions inside the crystalline bilayer walls of the DC8,9PC tubules strongly deviate from the normal, stochastic processes, with features characterizing not only anomalous subdiffusions but also motions that are superdiffusive in nature.     

Motion of single DNA molecules at a liquid-solid interface as revealed by variable-angle evanescent-field microscopy

A variable-angle total-internal-reflection fluorescence microscope (VATIRFM) capable of providing a large range of incident angles was constructed for imaging single DNA molecule dynamics at a solid/liquid interface. An algorithm using a public-domain image-processing program, ImageJ, was developed for single-molecule counting. The experimental counts at various incident angles with different evanescent-field layer (EFL) thicknesses are affected by molecular diffusion. The dynamics of molecules near the surface and the observed counts in the VATIRFM are elucidated using a limited one-dimensional random-walk diffusion model. The simulation fits well with the experimental counting results. Further analysis using the simulation reveals the details of single-molecule motion. One implication is that the measured intensities cannot be used directly to determine the distances of molecules from the surface, though the majority of fluorescence does come from the EFL. Another implication is that rather than providing molecular concentrations within EFL the experimental counting results depict the distance-dependent dynamics of molecules near the surface. Thus, the VATIRFM could be a powerful technique to study the surface repulsion/attraction of molecules within a few hundred nanometers of the surface. Further studies show that molecules at low ionic strengths experience electrostatic repulsion at distances much further away from the surface than the calculated thickness of the electrical double layer.     

Multi-dimensional pulsed field gradient magic angle spinning NMR experiments on membranes

The benefits of gradient techniques in the study of lipid membranes are demonstrated on a sample of 1-palmitoyl-2-oleoyl-sn-glycero-3 phosphocholine (POPC) liposomes embedded with ibuprofen. Most techniques from gradient NMR spectroscopy on solution samples are directly applicable to membrane samples subjected to magic angle spinning (MAS). Gradient-enhanced homo- and heteronuclear chemical shift correlation techniques were used to make resonance assignments. Gradient NOESY experiments provide insight into the location and dynamics of lipids, ibuprofen and water. Application of gradients not only reduces experiment time but also the t(1) noise in the multi-dimensional spectra. Diffusion measurements with pulsed field gradients characterize lateral movements of lipid and drug molecules in membranes. The theoretical framework for data analysis of MAS diffusion experiments on randomly oriented multilamellar liposomes is presented.     

Molecular dynamics simulation of the nematic liquid crystal phase in the presence of an intense magnetic field

The influence of an intense external field on the dynamics of the nematic liquid crystal phase is investigated using a molecular dynamics simulation for the Gay-Berne nematogen under isobaric-isothermal conditions. The molecular dynamics as a function of the second-rank orientational order parameter P<2> for a system consisting of a nematic liquid crystal in the presence of an intense magnetic field is compared with that of a similar system without the field. The translational motion of molecules is determined as a function of the translational diffusion coefficient tensor and the anisotropy and compared with the values predicted theoretically. The rotational dynamics of molecules is analyzed using the first- and the second-rank orientational time correlation functions. The translational diffusion coefficient parallel with respect to the director is constrained by the intense field, although the perpendicular one is decreased as the P<2> is increased, just as it is in the system without the field. However, no essential effect of the strong magnetic field is observed in the rotational molecular dynamics. Further, the rotational diffusion coefficient parallel with respect to the director obtained from the first-rank orientational time correlation function in the simulation is qualitatively in agreement with that in the real nematic liquid crystalline molecules. The P<2> dependence of the rotational diffusion coefficient for the system with the intense magnetic field shows a tendency similar to that for the system without the field.     

Detecting phase transitions in phosphatidylcholine vesicles by Raman microscopy and self-modeling curve resolution

The study of phospholipid phase transitions is important for understanding drug- and protein-membrane interactions as well as other phenomena such as trans-membrane diffusion and vesicle fusion. A temperature-controlled stage on a confocal Raman microscope has allowed phase transitions in optically trapped phospholipid vesicles to be monitored. Raman spectra were acquired and analyzed using self-modeling curve resolution, a multivariate statistical analysis technique. This method revealed the subtle spectral changes indicative of sub- and pretransitions and main transitions in vesicles composed of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC). The Raman scattering results were compared to differential scanning calorimetry (DSC) experiments and found to be in good agreement. This method of observing lipid phase transition profiles requires little sample preparation and a minimal amount of lipid (相似文献