E. Jayne Wallace

The interaction of C60 and its derivatives with a lipid bilayer via molecular dynamics simulations

Coarse-grained molecular dynamics simulations have been used to explore the interactions of C(60) and its derivatives with lipid bilayers. Pristine C(60) partitions into the bilayer core, whilst C(60)(OH)(20) experiences a central energetic barrier to permeation across the bilayer. For intermediate levels of derivatization, e.g. C(60)(OH)(10), this central barrier is smaller and there is an energetic well at the bilayer/water interface, thus promoting entry into cells via bilayer permeation whilst maintaining solubility in water.

A multiscale simulation study of carbon nanotube interactions with designed amphiphilic peptide helices

The dispersion and manipulation of carbon nanotubes (CNTs) are of great importance if we are to utilise the unique properties of CNTs in a range of biological, electrical and mechanical applications. Recently, a designed amphiphilic peptide helix termed nano-1 has been shown to solubilise CNTs in aqueous solution. Furthermore, the peptide is capable of assembling these coated tubes into fibres. We use a multiscale molecular dynamics approach to study the adsorption profile of nano-1 on a CNT surface. We find that nano-1 interacts with a CNT in a preferred orientation, such that its hydrophobic surface is in contact with the tube. The adsorption profile is unchanged upon increasing the number of peptides on the CNT. Interestingly, when few peptides are adsorbed onto the CNT surface we find that the secondary structure of the peptide is unstable. However, the helical secondary structure is stabilised upon increasing the number of peptides on the CNT surface. This study sheds light on the adsorption of peptides on CNTs, and may be exploitable to enhance the selective solubilisation and manipulation of CNTs.

Designing a hydrophobic barrier within biomimetic nanopores.

Nanopores in membranes have a range of potential applications. Biomimetic design of nanopores aims to mimic key functions of biological pores within a stable template structure. Molecular dynamics simulations have been used to test whether a simple β-barrel protein nanopore can be modified to incorporate a hydrophobic barrier to permeation. Simulations have been used to evaluate functional properties of such nanopores, using water flux as a proxy for ionic conductance. The behavior of these model pores has been characterized as a function of pore size and of the hydrophobicity of the amino acid side chains lining the narrow central constriction of the pore. Potential of mean force calculations have been used to calculate free energy landscapes for water and for ion permeation in selected models. These studies demonstrate that a hydrophobic barrier can indeed be designed into a β-barrel protein nanopore, and that the height of the barrier can be adjusted by modifying the number of consecutive rings of hydrophobic side chains. A hydrophobic barrier prevents both water and ion permeation even though the pore is sterically unoccluded. These results both provide insights into the nature of hydrophobic gating in biological pores and channels, and furthermore demonstrate that simple design features may be computationally transplanted into β-barrel membrane proteins to generate functionally complex nanopores.

Functional Annotation of Ion Channel Structures by Molecular Simulation.

Summary Ion channels play key roles in cell membranes, and recent advances are yielding an increasing number of structures. However, their functional relevance is often unclear and better tools are required for their functional annotation. In sub-nanometer pores such as ion channels, hydrophobic gating has been shown to promote dewetting to produce a functionally closed (i.e., non-conductive) state. Using the serotonin receptor (5-HT3R) structure as an example, we demonstrate the use of molecular dynamics to aid the functional annotation of channel structures via simulation of the behavior of water within the pore. Three increasingly complex simulation analyses are described: water equilibrium densities; single-ion free-energy profiles; and computational electrophysiology. All three approaches correctly predict the 5-HT3R crystal structure to represent a functionally closed (i.e., non-conductive) state. We also illustrate the application of water equilibrium density simulations to annotate different conformational states of a glycine receptor.

DNA sequencing with MspA: Molecular Dynamics simulations reveal free-energy differences between sequencing and non-sequencing mutants.

MspA has been identified as a promising candidate protein as a component of a nanopore-based DNA-sequencing device. However the wildtype protein must be engineered to incorporate all of the features desirable for an accurate and efficient device. In the present study we have utilized atomistic molecular dynamics to perform umbrella-sampling calculations to calculate the potential of mean force (PMF) profiles for translocation of the four DNA nucleotides through MspA. We show there is an energetic barrier to translocation of individual nucleotides through a mutant that closely resembles the wildtype protein, but not through a mutant engineered for the purpose of sequencing. Crucially we are able to quantify the change in free energy for mutating key residues. Thus providing a quantitative characterisation of the energetic impact of individual amino acid sidechains on nucleotide translocation through the pore of MspA.

Free-energy calculations reveal the subtle differences in the interactions of DNA bases with α-hemolysin.

Next generation DNA sequencing methods that utilize protein nanopores have the potential to revolutionize this area of biotechnology. While the technique is underpinned by simple physics, the wild-type protein pores do not have all of the desired properties for efficient and accurate DNA sequencing. Much of the research efforts have focused on protein nanopores, such as α-hemolysin from Staphylococcus aureus. However, the speed of DNA translocation has historically been an issue, hampered in part by incomplete knowledge of the energetics of translocation. Here we have utilized atomistic molecular dynamics simulations of nucleotide fragments in order to calculate the potential of mean force (PMF) through α-hemolysin. Our results reveal specific regions within the pore that play a key role in the interaction with DNA. In particular, charged residues such as D127 and K131 provide stabilizing interactions with the anionic DNA and therefore are likely to reduce the speed of translocation. These regions provide rational targets for pore optimization. Furthermore, we show that the energetic contributions to the protein-DNA interactions are a complex combination of electrostatics and short-range interactions, often mediated by water molecules.

Electric-Field-Driven Translocation of ssDNA through Hydrophobic Nanopores

The accurate sequencing of DNA using nanopores requires control over the speed of DNA translocation through the pores and also of the DNA conformation. Our studies show that ssDNA translocates through hourglass-shaped pores with hydrophobic constriction regions when an electric field is applied. The constriction provides a barrier to translocation and thereby slows down DNA movement through the pore compared with pores without the constriction. We show that ssDNA moves through these hydrophobic pores in an extended conformation and therefore does not form undesirable secondary structures that may affect the accuracy of partial current blockages for DNA sequencing.