Understanding nature's design for a nanosyringe

Molecular dynamics simulations of membrane proteins

Membrane protein structures are underrepresented in the Protein Data Bank (PDB) because of difficulties associated with expression and crystallization. As such, it is one area in which computational studies, particularly molecular dynamics (MD), can provide useful additional information. Recently, there has been substantial progress in the simulation of lipid bilayers and membrane proteins embedded within them. Initial efforts at simulating membrane proteins embedded within a lipid bilayer were relatively slow and interactive processes, but recent advances now mean that the setup and running of membrane protein simulations is somewhat more straightforward, although not without its problems. In this chapter, we outline practical methods for setting up and running MD simulations of a membrane protein embedded within a lipid bilayer and discuss methodologies that are likely to contribute future improvements.

Computer simulation study of fullerene translocation through lipid membranes

Recent toxicology studies suggest that nanosized aggregates of fullerene molecules can enter cells and alter their functions, and also cross the blood-brain barrier. However, the mechanisms by which fullerenes penetrate and disrupt cell membranes are still poorly understood. Here we use computer simulations to explore the translocation of fullerene clusters through a model lipid membrane and the effect of high fullerene concentrations on membrane properties. The fullerene molecules rapidly aggregate in water but disaggregate after entering the membrane interior. The permeation of a solid-like fullerene aggregate into the lipid bilayer is thermodynamically favoured and occurs on the microsecond timescale. High concentrations of fullerene induce changes in the structural and elastic properties of the lipid bilayer, but these are not large enough to mechanically damage the membrane. Our results suggest that mechanical damage is an unlikely mechanism for membrane disruption and fullerene toxicity.

Simulations of electrophoretic RNA transport through transmembrane carbon nanotubes

The study of interactions between carbon nanotubes and cellular components, such as membranes and biomolecules, is fundamental for the rational design of nanodevices interfacing with biological systems. In this work, we use molecular dynamics simulations to study the electrophoretic transport of RNA through carbon nanotubes embedded in membranes. Decorated and naked carbon nanotubes are inserted into a dodecane membrane and a dimyristoylphosphatidylcholine lipid bilayer, and the system is subjected to electrostatic potential differences. The transport properties of this artificial pore are determined by the structural modifications of the membrane in the vicinity of the nanotube openings and they are quantified by the nonuniform electrostatic potential maps at the entrance and inside the nanotube. The pore is used to transport electrophoretically a short RNA segment and we find that the speed of translocation exhibits an exponential dependence on the applied potential differences. The RNA is transported while undergoing a repeated stacking and unstacking process, affected by steric interactions with the membrane headgroups and by hydrophobic interaction with the walls of the nanotube. The RNA is structurally reorganized inside the nanotube, with its backbone solvated by water molecules near the axis of the tube and its bases aligned with the nanotube walls. Upon exiting the pore, the RNA interacts with the membrane headgroups and remains attached to the dodecane membrane while it is expelled into the solvent in the case of the lipid bilayer. The results of the simulations detail processes of molecular transport into cellular compartments through manufactured nanopores and they are discussed in the context of applications in biotechnology and nanomedicine.

Computer modeling in biotechnology: a partner in development

Computational modeling can be a useful partner in biotechnology, in particular, in nanodevice engineering. Such modeling guides development through nanoscale views of biomolecules and devices not available through experimental imaging methods. We illustrate the role of computational modeling, mainly of molecular dynamics, through four case studies: development of silicon bionanodevices for single molecule electrical recording, development of carbon nano-tube-biomolecular systems as in vivo sensors, development of lipoprotein nanodiscs for assays of single membrane proteins, and engineering of oxygen tolerance into the enzyme hydrogenase for photosynthetic hydrogen gas production. The four case studies show how molecular dynamics approaches were adapted to the specific technical uses through (i) multi-scale extensions, (ii) fast quantum chemical force field evaluation, (iii) coarse graining, and (iv) novel sampling methods. The adapted molecular dynamics simulations provided key information on device behavior and revealed development opportunities, arguing that the "computational microscope" is an indispensable nanoengineering tool.

Use of degradable and nondegradable nanomaterials for controlled release

Drug-delivery devices are fundamentally important in improving the pharmacological profiles of therapeutic molecules. Nanocontrolled-release systems are attracting a lot of attention currently owing to their large surface area and their ability to target delivery to specific sites in the human body. In addition, they can penetrate the cell membrane for gene, nucleic acid and bioactive peptide/protein delivery. Representative applications of nanodrug-delivery systems include controlled-release wound dressings, controlled-release scaffolds for tissue regeneration and implantable biodegradable nanomaterial-based medical devices integrated with drug-delivery functions. We review the present status and future perspectives of various types of nanocontrolled-release systems. Although many of the well-established degradable and nondegradable controlled-release vehicles are being investigated for their processing into nanocarriers, several new emerging nanomaterials are being studied for their controlled-release properties. The release of multiple bioactive agents, each with its own kinetic profile, is becoming possible. In addition, integration of the nanocontrolled-release systems with other desirable functions to create new, cross-discipline applications can also be realized.

Multiscale Coarse-Graining of Monosaccharides

A systematic multiscale coarse-graining (MS-CG) algorithm is applied to build coarse-grained models for monosaccharides in aqueous solution. The methodology is demonstrated for the example of alpha-D-glucopyranose. The nonbonded interactions are directly derived from the force-matching approach, whereas the bonded interactions are obtained through Boltzmann statistical analyses of the underlying atomistic trajectory. The MS-CG model is shown to reproduce many structural and thermodynamic properties in the constant NPT ensemble. Although the model is derived at a single temperature, pressure, and concentration, it is shown to be reasonably transferable to other thermodynamic states. In this model, long-range interactions are effectively mapped into short-range forces with a moderate cutoff and are evaluated by table look-up. As a result, molecular dynamics employing the MS-CG model is approximately 3 orders of magnitude more efficient than its atomistic counterpart. Consequently, the model is particularly suitable for simulating carbohydrate systems at large length and long time scales. Results for an alpha-(1-->4)-d-glucan with 14 glucose units are also presented, demonstrating that the MS-CG algorithm is also applicable to the coarse-graining of other saccharide systems.

Mixed atomistic and coarse-grained molecular dynamics: simulation of a membrane-bound ion channel

The recently developed multiscale coarse-graining (MS-CG) method (Izvekov, S.; Voth, G. A. J. Phys. Chem. B 2005, 109, 2469; J. Chem. Phys. 2005, 123, 134105) is used to build a mixed all-atom and coarse-grained (AA-CG) model of the gramicidin A (gA) ion channel embedded in a dimyristoylphosphatidylcholine (DMPC) lipid bilayer and water environment. In this model, the gA peptide was described in full atomistic detail, while the lipid and water molecules were described using coarse-grained representations. The atom-CG and CG-CG interactions in the mixed AA-CG model were determined using the MS-CG method. Molecular dynamics (MD) simulations were performed using the resulting AA-CG model. The results from simulations of the AA-CG model compare very favorably to those from all-atom MD simulations of the entire system. Since the MS-CG method employs a general and systematic approach to obtain effective interactions from the underlying all-atom models, the present approach to rigorously develop mixed AA-CG models has the potential to be extended to many other systems.

Single-walled carbon nanotubes exhibit strong antimicrobial activity

We provide the first direct evidence that highly purified single-walled carbon nanotubes (SWNTs) exhibit strong antimicrobial activity. By using a pristine SWNT with a narrow diameter distribution, we demonstrate that cell membrane damage resulting from direct contact with SWNT aggregates is the likely mechanism leading to bacterial cell death. This finding may be useful in the application of SWNTs as building blocks for antimicrobial materials.

Oxygen adsorption by carbon nanotubes and its application in radiotherapy

The ability to deliver large molecules, for example nucleic acids, to cells using carbon nanotubes has been reviewed. Potential applications of functionalised nanotubes to deliver oxygen to cancer cells to enhance the effects of radiotherapy are considered.

Coarse-grained molecular dynamics simulations of membrane proteins and peptides

Molecular dynamics (MD) simulations provide a valuable approach to the dynamics, structure, and stability of membrane-protein systems. Coarse-grained (CG) models, in which small groups of atoms are treated as single particles, enable extended (>100 ns) timescales to be addressed. In this study, we explore how CG-MD methods that have been developed for detergents and lipids may be extended to membrane proteins. In particular, CG-MD simulations of a number of membrane peptides and proteins are used to characterize their interactions with lipid bilayers. CG-MD is used to simulate the insertion of synthetic model membrane peptides (WALPs and LS3) into a lipid (PC) bilayer. WALP peptides insert in a transmembrane orientation, whilst the LS3 peptide adopts an interfacial location, both in agreement with experimental biophysical data. This approach is extended to a transmembrane fragment of the Vpu protein from HIV-1, and to the coat protein from fd phage. Again, simulated protein/membrane interactions are in good agreement with solid state NMR data for these proteins. CG-MD has also been applied to an M3-M4 fragment from the CFTR protein. Simulations of CFTR M3-M4 in a detergent micelle reveal formation of an alpha-helical hairpin, consistent with a variety of biophysical data. In an I231D mutant, the M3-M4 hairpin is additionally stabilized via an inter-helix Q207/D231 interaction. Finally, CG-MD simulations are extended to a more complex membrane protein, the bacterial sugar transporter LacY. Comparison of a 200 ns CG-MD simulation of LacY in a DPPC bilayer with a 50 ns atomistic simulation of the same protein in a DMPC bilayer shows that the two methods yield comparable predictions of lipid-protein interactions. Taken together, these results demonstrate the utility of CG-MD simulations for studies of membrane/protein interactions.

Bioengineered flagella protein nanotubes with cysteine loops: self-assembly and manipulation in an optical trap

An E. coli flagellin protein, termed FliTrx, was investigated for use as a novel form of self-assembling protein nanotube. This protein was genetically engineered to display constrained peptide loops with a series of different thiol, cationic, anionic, and imidazole functional groups. "Cys-loop" thiol variants consisting of 6 and 12 cysteine residues were isolated in the form of disulfide-linked nanotube bundles, a novel nanomaterial. Bundles were characterized by fluorescence microscopy, transmission electron microscopy, and optical trapping.

SFG studies on interactions between antimicrobial peptides and supported lipid bilayers

The mode of action of antimicrobial peptides (AMPs) in disrupting cell membrane bilayers is of fundamental importance in understanding the efficiency of different AMPs, which is crucial to design antibiotics with improved properties. Recent developments in the field of sum frequency generation (SFG) vibrational spectroscopy have made it a powerful and unique biophysical technique in investigating the interactions between AMPs and a single substrate supported planar lipid bilayer. We will review some of the recent progress in applying SFG to study membrane lipid bilayers and discuss how SFG can provide novel information such as real-time bilayer structure change and AMP orientation during AMP-lipid bilayer interactions in a very biologically relevant manner. Several examples of applying SFG to monitor such interactions between AMPs and a dipalmitoyl phosphatidylglycerol (DPPG) bilayer are presented. Different modes of actions are observed for melittin, tachyplesin I, d-magainin 2, MSI-843, and a synthetic antibacterial oligomer, demonstrating that SFG is very effective in the study of AMPs and AMP-lipid bilayer interactions.

Synthetic ion channels and pores (2004-2005)

This critical review covers synthetic ion channels and pores created between January 2004 and December 2005 comprehensively. The discussion of a rich collection of structural motifs may particularly appeal to organic, biological, supramolecular and polymer chemists. Functions addressed include ion selectivity and molecular recognition, as well as responsiveness to light, heat, voltage and membrane composition. The practical applications involved concern certain topics in medicinal chemistry (antibiotics, drug delivery), catalysis and sensing. An introduction to principles and methods is provided for the non-specialist; some new sources of inspiration from fields beyond chemistry are highlighted.

Observing a molecular knife at work

Sum frequency generation (SFG) vibrational spectroscopy has been employed to study the molecular interactions between a single substrate supported lipid bilayer and an amphiphilic antibiotic compound 1, with a design based on the common structural motif of natural antimicrobial peptides. The interfacial sensitivity of SFG allows real-time in situ monitoring of ordering changes in both leaflets of the bilayer and orientation of 1 simultaneously. A critical concentration of about 0.8 microg/mL of 1 is found, above which the inner leaflet of the bilayer is significantly perturbed. This concentration corresponds well to the minimum inhibition concentration of 1 that is obtained from bacterial experiments. Orientation of 1 in the bilayer is shown to be perpendicular to the bilayer surface, in agreement with simulation results. SFG can be developed into a very informative technique for studying the cell membrane and the interactions of membrane-active molecules.

The mechanism of water diffusion in narrow carbon nanotubes

Carbon nanotubes show exceptional physical properties that render them promising candidates as building blocks for nanostructured materials. Many ambitious applications, ranging from gene therapy to membrane separations, require the delivery of fluids, in particular aqueous solutions, through the interior of carbon nanotubes. To foster these and other applications, it is necessary to understand the thermodynamic and transport properties of water confined within long narrow carbon nanotubes. Previous theoretical work considered either short carbon nanotubes or short periods of time. By conducting molecular dynamics simulations in the microcanonical ensemble for water confined in infinitely long carbon nanotubes of diameter 1.08 nm, we show here that confined water molecules diffuse through a fast ballistic motion mechanism for up to 500 ps at room temperature. By comparing the results obtained for the diffusion of water to those obtained for the diffusion of a reference Lennard-Jones fluid, we prove here that long-lasting hydrogen bonds are responsible for the ballistic diffusion of water clusters in narrow carbon nanotubes, as opposed to spatial mismatches between pore-fluid and fluid-fluid attractive interactions which, as shown previously by others, are responsible for the concerted motion of simple fluids in molecular sieves. Additionally we prove here for the first time that, despite the narrow diameter of the carbon nanotubes considered which may suggest the existence of single-file diffusion, when the trajectories of confined water are studied at time scales in excess of 500 ps, a Fickian-type diffusion mechanism prevails. Our results are important for designing nano fluidic apparatuses to develop, for example, novel drug-delivery devices.

Probing the Kinetics of Membrane-Mediated Helix Folding

The kinetics of peptide-membrane association have been studied previously using stopped-flow tryptophan fluorescence; however, such experiments do not directly report the coil-to-helix transition process, which is a hallmark of peptide-membrane interaction. Herein, we report a new method for directly assessing the kinetics of the helix formation accompanied by the peptide-membrane association. This method is based on the technique of fluorescence resonance energy transfer (FRET) and an amino acid FRET pair, p-cyano-L-phenylalanine and tryptophan. To demonstrate the utility of this method, we have studied the membrane-mediated helix folding dynamics of a mutant of magainin 2, an antibiotic peptide found in the skin of the African clawed frog, Xenopus laevis. Our results indicate that the coil-to-helix transition occurs during the binding of the peptide to the lipid vesicle (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1-palmitoyl-2-oleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)], 3:1, wt/wt) but prior to the full insertion of the peptide into the hydrophobic region of the lipid bilayers.

Modelling of proteins in membranes

This review describes some recent theories and simulations of mesoscopic and microscopic models of lipid membranes with embedded or attached proteins. We summarize results supporting our understanding of phenomena for which the activities of proteins in membranes are expected to be significantly affected by the lipid environment. Theoretical predictions are pointed out, and compared to experimental findings, if available. Among others, the following phenomena are discussed: interactions of interfacially adsorbed peptides, pore-forming amphipathic peptides, adsorption of charged proteins onto oppositely charged lipid membranes, lipid-induced tilting of proteins embedded in lipid bilayers, protein-induced bilayer deformations, protein insertion and assembly, and lipid-controlled functioning of membrane proteins.

Insertion and assembly of membrane proteins via simulation

Interactions of lipids are central to the folding and stability of membrane proteins. Coarse-grained molecular dynamics simulations have been used to reveal the mechanisms of self-assembly of protein/membrane and protein/detergent complexes for representatives of two classes of membrane protein, namely, glycophorin (a simple alpha-helical bundle) and OmpA (a beta-barrel). The accuracy of the coarse-grained simulations is established via comparison with the equivalent atomistic simulations of self-assembly of protein/detergent micelles. The simulation of OmpA/bilayer self-assembly reveals how a folded outer membrane protein can be inserted in a bilayer. The glycophorin/bilayer simulation supports the two-state model of membrane folding, in which transmembrane helix insertion precedes dimer self-assembly within a bilayer. The simulations also suggest that a dynamic equilibrium exists between the glycophorin helix monomer and dimer within a bilayer. The simulated glycophorin helix dimer is remarkably close in structure to that revealed by NMR. Thus, coarse-grained methods may help to define mechanisms of membrane protein (re)folding and will prove suitable for simulation of larger scale dynamic rearrangements of biological membranes.

The role of molecular modeling in bionanotechnology

Molecular modeling is advocated here as a key methodology for research and development in bionanotechnology. Molecular modeling provides nanoscale images at atomic and even electronic resolution, predicts the nanoscale interaction of unfamiliar combinations of biological and inorganic materials, and evaluates strategies for redesigning biopolymers for nanotechnological uses. The methodology is illustrated in this paper through reviewing three case studies. The first one involves the use of single-walled carbon nanotubes as biomedical sensors where a computationally efficient, yet accurate, description of the influence of biomolecules on nanotube electronic properties through nanotube-biomolecule interactions was developed; this development furnishes the ability to test nanotube electronic properties in realistic biological environments. The second case study involves the use of nanopores manufactured into electronic nanodevices based on silicon compounds for single molecule electrical recording, in particular, for DNA sequencing. Here, modeling combining classical molecular dynamics, material science and device physics, described the interaction of biopolymers, e.g., DNA, with silicon nitrate and silicon oxide pores, furnished accurate dynamic images of pore translocation processes, and predicted signals. The third case study involves the development of nanoscale lipid bilayers for the study of embedded membrane proteins and cholesterol. Molecular modeling tested scaffold proteins, redesigned apolipoproteins found in mammalian plasma that hold the discoidal membranes in the proper shape, and predicted the assembly as well as final structure of the nanodiscs. In entirely new technological areas such as bionanotechnology, qualitative concepts, pictures and suggestions are sorely needed; these three case studies document that molecular modeling can serve a critical role in this respect, even though it may still fall short on quantitative precision.

Membrane Activity of Biomimetic Facially Amphiphilic Antibiotics

Membranes are a central feature of all biological systems, and their ability to control many cellular processes is critically important. As a result, a better understanding of how molecules bind to and select between biological membranes is an active area of research. Antimicrobial host defense peptides are known to be membrane-active and, in many cases, exhibit discrimination between prokaryotic and eukaryotic cells. The design of synthetic molecules that capture the biological activity of these natural peptides has been shown. In this report, the interaction between our biomimetic structures and different biological membranes is reported using both model vesicle and in vitro bacterial cell experiments. Compound 1 induces 12% leakage at 20 microg/mL against phosphatidylglycerol (PG)-phosphatidylethanolamine (PE) vesicles vs only 3% leakage at 200 microg/mL against phosphatidyl-L-serine (PS)-phosphatidylcholine (PC) vesicles. Similarly, a 40% reduction in fluorescence is measured in lipid movement experiments for PG-PE compared to 10% for PS-PC at 600 s. A 30 degrees C increase in the phase transition of stearoyl-oleoyl-phosphatidylserine is observed in the presence of 1. These results show that lipid composition is more important for selectivity than overall net charge. Additionally, the overall concentration of a given lipid is another important factor. An effort is made to connect model vesicle studies with in vitro data and naturally occurring lipid compositions.

Molecular Dynamics Simulations of Surfactant Self-Organization at a Solid−Liquid Interface

Self-organization of aqueous surfactants at a planar graphite-like surface is studied by means of coarse-grain molecular dynamics simulations. The nonionic surfactant, n-alkyl poly(ethylene oxide), and water are both represented by coarse-grain models while an implicit representation is used for the graphite surface. The observed morphology of the aggregated surfactants depends on the alkyl chain length. Surfactants with a short chain form a monolayer on the graphite surface with a thickness roughly equal to that of the alkane tail. On the other hand, longer-tail surfactants form continuous hemicylinders on the surface with diameter approximately 5.0 +/- 0.5 nm, in good agreement with experimental AFM data.

Carbon nanotubes for the delivery of therapeutic molecules

Functionalised carbon nanotubes (f-CNTs) are emerging as new tools in the field of nanobiotechnology and nanomedicine. This is because they can be easily manipulated and modified by encapsulation with biopolymers or by covalent linking of solubilising groups to the external walls and tips. The possibility of incorporating f-CNTs into biological systems has opened the way to the exploration of their potential applications in biology and medicinal chemistry. Within the different fields of applications (i.e., biosensors, composite materials, molecular electronics), one use of CNTs is as new carrier systems for the delivery of therapeutic molecules. Research discussed in this review is focused on recent advances in the development of CNT technology for the delivery of drugs, antigens and genes.

Molecular dynamics simulations of salicylate effects on the micro- and mesoscopic properties of a dipalmitoylphosphatidylcholine bilayer

Salicylate, an amphiphilic molecule and a popular member of the nonsteroidal anti-inflammatory drug family, is known to affect hearing through reduction of the electromechanical coupling in the outer hair cells of the ear. This reduction of electromotility by salicylate has been widely studied, but the molecular mechanism of the phenomenon is still unknown. In this study, we investigated one aspect of salicylate's action, namely the perturbation of electrical and mechanical membrane properties by salicylate in the absence of cytoskeletal or membrane-bound motor proteins such as prestin. In particular, we simulated the interaction of salicylate with a dipalmitoylphosphatidylcholine (DPPC) bilayer via atomically detailed molecular dynamics simulations to observe the effect of salicylate on the microscopic and mesoscopic properties of the bilayer. The results demonstrate that salicylate interacts with the bilayer by associating at the water-DPPC interface in a nearly perpendicular orientation and penetrating more deeply into the bilayer than either sodium or chloride. This association has several affects on the membrane properties. First, binding of salicylate to the membrane displaces chloride from the bilayer-water interface. Second, salicylate influences the electrostatic potential and dielectric properties of the bilayer, with significant changes at the water-lipid bilayer interface. Third, salicylate association results in structural changes, including decreased headgroup area per lipid and increased lipid tail order. However, salicylate does not significantly alter the mechanical properties of the DPPC bilayer; bulk compressibility, area compressibility, and bending modulus were only perturbed by small, statistically insignificant amounts by the presence of salicylate. The observations from these simulations are in qualitative agreement with experimental data and support the conclusion that salicylate influences the electrical but not the mechanical properties of DPPC membranes.

Applications of carbon nanotubes in drug delivery

The development of new and efficient drug delivery systems is of fundamental importance to improve the pharmacological profiles of many classes of therapeutic molecules. Many different types of drug delivery systems are currently available. Within the family of nanomaterials, carbon nanotubes (CNT) have emerged as a new alternative and efficient tool for transporting and translocating therapeutic molecules. CNT can be functionalised with bioactive peptides, proteins, nucleic acids and drugs, and used to deliver their cargos to cells and organs. Because functionalised CNT display low toxicity and are not immunogenic, such systems hold great potential in the field of nanobiotechnology and nanomedicine.

Ion transport through membrane-spanning nanopores studied by molecular dynamics simulations and continuum electrostatics calculations

Narrow hydrophobic regions are a common feature of biological channels, with possible roles in ion-channel gating. We study the principles that govern ion transport through narrow hydrophobic membrane pores by molecular dynamics simulation of model membranes formed of hexagonally packed carbon nanotubes. We focus on the factors that determine the energetics of ion translocation through such nonpolar nanopores and compare the resulting free-energy barriers for pores with different diameters corresponding to the gating regions in closed and open forms of potassium channels. Our model system also allows us to compare the results from molecular dynamics simulations directly to continuum electrostatics calculations. Both simulations and continuum calculations show that subnanometer wide pores pose a huge free-energy barrier for ions, but a small increase in the pore diameter to approximately 1 nm nearly eliminates that barrier. We also find that in those wider channels the ion mobility is comparable to that in the bulk phase. By calculating local electrostatic potentials, we show that the long range Coulomb interactions of ions are strongly screened in the wide water-filled channels. Whereas continuum calculations capture the overall energetics reasonably well, the local water structure, which is not accounted for in this model, leads to interesting effects such as the preference of hydrated ions to move along the pore wall rather than through the center of the pore.

Mechanics of receptor-mediated endocytosis

Most viruses and bioparticles endocytosed by cells have characteristic sizes in the range of tens to hundreds of nanometers. The process of viruses entering and leaving animal cells is mediated by the binding interaction between ligand molecules on the viral capid and their receptor molecules on the cell membrane. How does the size of a bioparticle affect receptor-mediated endocytosis? Here, we study how a cell membrane containing diffusive mobile receptors wraps around a ligand-coated cylindrical or spherical particle. It is shown that particles in the size range of tens to hundreds of nanometers can enter or exit cells via wrapping even in the absence of clathrin or caveolin coats, and an optimal particles size exists for the smallest wrapping time. This model can also be extended to include the effect of clathrin coat. The results seem to show broad agreement with experimental observations.

Interfacial folding and membrane insertion of designed peptides studied by molecular dynamics simulations

The mechanism of interfacial folding and membrane insertion of designed peptides is explored by using an implicit membrane generalized Born model and replica-exchange molecular dynamics. Folding/insertion simulations initiated from fully extended peptide conformations in the aqueous phase, at least 28 A away from the membrane interface, demonstrate a general mechanism for structure formation and insertion (when it occurs). The predominately hydrophobic peptides from the synthetic WALP and TMX series first become localized at the membrane-solvent interface where they form significant helical secondary structure via a helix-turn-helix motif that inserts the central hydrophobic residues into the membrane interior, and then fluctuations occur that provide a persistent helical structure throughout the peptide and it inserts with its N-terminal end moving across the membrane. More specifically, we observed that: (i) the WALP peptides (WALP16, WALP19, and WALP23) spontaneously insert in the membrane as just noted; (ii) TMX-1 also inserts spontaneously after a similar mechanism and forms a transmembrane helix with a population of approximately 50% at 300 K; and (iii) TMX-3 does not insert, but exists in a fluctuating membrane interface-bound form. These findings are in excellent agreement with available experimental data and demonstrate the potential for new implicit solvent/membrane models together with advanced simulation protocols to guide experimental programs in exploring the nature and mechanism of membrane-associated folding and insertion of biologically important peptides.

Controllable Water Channel Gating of Nanometer Dimensions

The dynamics of water molecules in a single-walled carbon nanotube (SWNT) under continuous deformations was studied with molecular dynamics simulations. The flux and occupancy remain almost fixed within a deformation of 2.0 A but decrease sharply for a further deformation of 0.6 A. The nanopore is an excellent on-off gate that is both effectively resistant to deformation noises and sensitive to available signals. Biological water channels are expected to share this advantage due to similar wavelike water distributions. The minimal external force required for triggering an open-close transition falls within the working range of many available experimental facilities, which provides the possibility of developing SWNT-based nanoscale devices.

Lipid bilayer perturbations around a transmembrane nanotube: a coarse grain molecular dynamics study

The perturbations induced in a lipid bilayer by the presence of a transmembrane nanotube are investigated using coarse grained molecular dynamics. Meniscus formation by the lipids and tilting of the nanotube occur in response to hydrophobic mismatch, although these two effects do not compensate completely for the total mismatch. The lipid head-to-tail vector field is examined and shows strong ordering in the membrane plane regardless of the nanotube length. Molecular layering at the lipid-nanotube interface is reported. This study extends previous theoretical approaches to a more realistic setting.

Structure and dynamics of model pore insertion into a membrane

A cylindrical transmembrane molecule is constructed by linking hydrophobic sites selected from a coarse grain model. The resulting hollow tube assembly serves as a representation of a transmembrane channel, pore, or a carbon nanotube. The interactions of a coarse grain di-myristoyl-phosphatidyl-choline hydrated bilayer with both a purely hydrophobic tube and a tube with hydrophilic caps are studied. The hydrophobic tube rotates in the membrane and becomes blocked by lipid tails after a few tens of nanoseconds. The hydrophilic sites of the capped tube stabilize it by anchoring the tube in the lipid headgroup/water interfacial region of each membrane leaflet. The capped tube remains free of lipid tails. The capped tube spontaneously conducts coarse grain water sites; the free-energy profile of this process is calculated using three different methods and is compared to the barrier for water permeation through the lipid bilayer. Spontaneous tube insertion into an undisturbed lipid bilayer is also studied, which we reported briefly in a previous publication. The hydrophobic tube submerges into the membrane core in a carpetlike manner. The capped tube laterally fuses with the closest leaflet, and then, after plunging into the membrane interior, rotates to assume a transbilayer orientation. Two lipids become trapped at the end of the tube as it penetrates the membrane. The hydrophilic headgroups of these lipids associate with the lower tube cap and assist the tube in crossing the interior of the membrane. When the rotation is complete these lipids detach from the tube caps and fuse with the lower leaflet lipids.

Computational studies of membrane channels

The determination of the structure of several members of the K+ channel and aquaporin family represents a unique opportunity to explain the mechanism of these biomolecular systems. With their ability to go beyond static structures, molecular dynamics simulations offer a unique route for relating functional properties to membrane channel structure. The recent progress in this area is reviewed.

Concentration effects of volatile anesthetics on the properties of model membranes: a coarse-grain approach

To gain insights into the molecular level mechanism of drug action at the membrane site, we have carried out extensive molecular dynamics simulations of a model membrane in the presence of a volatile anesthetic using a coarse-grain model. Six different anesthetic (halothane)/lipid (dimyristoylphosphatidylcholine) ratios have been investigated, going beyond the low doses typical of medical applications. The volatile anesthetics were introduced into a preassembled fully hydrated 512-molecule lipid bilayer and each of the molecular dynamics simulations were carried out at ambient conditions, using the NPT ensemble. The area per lipid increases monotonically with the halothane concentration and the lamellar spacing decreases, whereas the lipid bilayer thickness shows no appreciable differences and only a slight increase upon addition of halothane. The density profiles of the anesthetic molecules display a bimodal distribution along the membrane normal with maxima located close to the lipid-water interface region. We have studied how halothane molecules fluctuate between the two maxima of the bimodal distribution and we observed a different mechanism at low and high anesthetic concentrations. Through the investigation of the reorientational motions of the lipid tails, we found that the anesthetic molecules increase the segmental order of the lipids close to the membrane surface.

Asymmetric ABC-Triblock Copolymer Membranes Induce a Directed Insertion of Membrane Proteins

Asymmetric molecules and materials provide an important basis for the organization and function of biological systems. It is well known that, for example, the inner and outer leaflets of biological membranes are strictly asymmetric with respect to lipid composition and distribution. This plays a crucial role for many membrane-related processes like carrier-mediated transport or insertion and orientation of integral membrane proteins. Most artificial membrane systems are, however, symmetric with respect to their midplane and membrane proteins are incorporated with random orientation. Here we describe a new approach to induce a directed insertion of membrane proteins into asymmetric membranes formed by amphiphilic ABC triblock copolymers with two chemically different water-soluble blocks A and C. In a comparative study we have reconstituted His-tag labeled Aquaporin 0 in lipid, ABA block copolymer, and ABC block copolymer vesicles. Immunolabeling, colorimetric, and fluorescence studies clearly show that a preferential orientation of the protein is only observed in the asymmetric ABC triblock copolymer membranes.

Microscopic Kinetics of DNA Translocation through synthetic nanopores

We have previously demonstrated that a nanometer-diameter pore in a nanometer-thick metal-oxide-semiconductor-compatible membrane can be used as a molecular sensor for detecting DNA. The prospects for using this type of device for sequencing DNA are avidly being pursued. The key attribute of the sensor is the electric field-induced (voltage-driven) translocation of the DNA molecule in an electrolytic solution across the membrane through the nanopore. To complement ongoing experimental studies developing such pores and measuring signals in response to the presence of DNA, we conducted molecular dynamics simulations of DNA translocation through the nanopore. A typical simulated system included a patch of a silicon nitride membrane dividing water solution of potassium chloride into two compartments connected by the nanopore. External electrical fields induced capturing of the DNA molecules by the pore from the solution and subsequent translocation. Molecular dynamics simulations suggest that 20-basepair segments of double-stranded DNA can transit a nanopore of 2.2 x 2.6 nm(2) cross section in a few microseconds at typical electrical fields. Hydrophobic interactions between DNA bases and the pore surface can slow down translocation of single-stranded DNA and might favor unzipping of double-stranded DNA inside the pore. DNA occluding the pore mouth blocks the electrolytic current through the pore; these current blockades were found to have the same magnitude as the blockade observed when DNA transits the pore. The feasibility of using molecular dynamics simulations to relate the level of the blocked ionic current to the sequence of DNA was investigated.