Syma Khalid

Coarse-grained MD simulations of membrane protein-bilayer self-assembly

Complete determination of a membrane protein structure requires knowledge of the protein position within the lipid bilayer. As the number of determined structures of membrane proteins increases so does the need for computational methods which predict their position in the lipid bilayer. Here we present a coarse-grained molecular dynamics approach to lipid bilayer self-assembly around membrane proteins. We demonstrate that this method can be used to predict accurately the protein position in the bilayer for membrane proteins with a range of different sizes and architectures.

Intramolecular DNA coiling mediated by metallo-supramolecular cylinders: Differential binding of P and M helical enantiomers

We have designed a synthetic tetracationic metallo-supramolecular cylinder that targets the major groove of DNA with a binding constant in excess of 107 M−1 and induces DNA bending and intramolecular coiling. The two enantiomers of the helical molecule bind differently to DNA and have different structural effects. We report the characterization of the interactions by a range of biophysical techniques. The M helical cylinder binds to the major groove and induces dramatic intramolecular coiling. The DNA bending is less dramatic for the P enantiomer.

Molecular Dynamics Simulations of Phosphatidylcholine Membranes: A Comparative Force Field Study

Molecular dynamics simulations provide a route to studying the dynamics of lipid bilayers at atomistic or near atomistic resolution. Over the past 10 years or so, molecular dynamics simulations have become an established part of the biophysicists tool kit for the study of model biological membranes. As simulation time scales move from tens to hundreds of nanoseconds and beyond, it is timely to re-evaluate the accuracy of simulation models. We describe a comparative analysis of five freely available force fields that are commonly used to model lipid bilayers. We focus our analysis on 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. We show that some bilayer properties have a pronounced force field dependence, while others are less sensitive. In general, we find strengths and weaknesses, with respect to experimental data, in all of the force fields we have studied. We do, however, find some combinations of simulation and force field parameters that should be avoided when simulating DPPC and POPC membranes. We anticipate that the results presented for some of the membrane properties will guide future improvements for several force fields studied in this work.

Electroporation of the E. coli and S. Aureus membranes: molecular dynamics simulations of complex bacterial membranes.

Bacterial membranes are complex organelles composed of a variety of lipid types. The differences in their composition are a key factor in determining their relative permeabilities. The success of antibacterial agents depends upon their interaction with bacterial membranes, yet little is known about the molecular-level interactions within membranes of different bacterial species. To address this, we have performed molecular dynamics simulations of two bacterial membranes: the outer membrane of E. coli and the cell membrane of S. aureus . We have retained the chemical complexity of the membranes by considering the details of their lipidic components. We identify the extended network of lipid-lipid interactions that stabilize the membranes. Our simulations of electroporation show that the S. aureus cell membrane is less resistant to poration than the E. coli outer membrane. The mechanisms of poration for the two membranes have subtle differences; for the E. coli outer membrane, relative differences in mobilities of the lipids of both leaflets are key in the process of poration.

Outer membrane protein G: Engineering a quiet pore for biosensing

Bacterial outer membrane porins have a robust β-barrel structure and therefore show potential for use as stochastic sensors based on single-molecule detection. The monomeric porin OmpG is especially attractive compared with multisubunit proteins because appropriate modifications of the pore can be easily achieved by mutagenesis. However, the gating of OmpG causes transient current blockades in single-channel recordings that would interfere with analyte detection. To eliminate this spontaneous gating activity, we used molecular dynamics simulations to identify regions of OmpG implicated in the gating. Based on our findings, two approaches were used to enhance the stability of the open conformation by site-directed mutagenesis. First, the mobility of loop 6 was reduced by introducing a disulfide bond between the extracellular ends of strands β12 and β13. Second, the interstrand hydrogen bonding between strands β11 and β12 was optimized by deletion of residue D215. The OmpG porin with both stabilizing mutations exhibited a 95% reduction in gating activity. We used this mutant for the detection of adenosine diphosphate at the single-molecule level, after equipping the porin with a cyclodextrin molecular adapter, thereby demonstrating its potential for use in stochastic sensing applications.

Self-Assembly of a Simple Membrane Protein: Coarse-Grained Molecular Dynamics Simulations of the Influenza M2 Channel

The transmembrane (TM) domain of the M2 channel protein from influenza A is a homotetrameric bundle of alpha-helices and provides a model system for computational approaches to self-assembly of membrane proteins. Coarse-grained molecular dynamics (CG-MD) simulations have been used to explore partitioning into a membrane of M2 TM helices during bilayer self-assembly from lipids. CG-MD is also used to explore tetramerization of preinserted M2 TM helices. The M2 helix monomer adopts a membrane spanning orientation in a lipid (DPPC) bilayer. Multiple extended CG-MD simulations (5 x 5 micros) were used to study the tetramerization of inserted M2 helices. The resultant tetramers were evaluated in terms of the most populated conformations and the dynamics of their interconversion. This analysis reveals that the M2 tetramer has 2x rotationally symmetrical packing of the helices. The helices form a left-handed bundle, with a helix tilt angle of approximately 16 degrees. The M2 helix bundle generated by CG-MD was converted to an atomistic model. Simulations of this model reveal that the bundles stability depends on the assumed protonation state of the H37 side chains. These simulations alongside comparison with recent x-ray (3BKD) and NMR (2RLF) structures of the M2 bundle suggest that the model yielded by CG-MD may correspond to a closed state of the channel.

DNA and lipid bilayers: self-assembly and insertion

DNA–lipid complexes are of biomedical importance as delivery vectors for gene therapy. To gain insight into the interactions of DNA with zwitterionic and cationic (dimyristoyltrimethylammonium propane (DMTAP)) lipids, we have used coarse-grained molecular dynamics simulations to study the self-assembly of DPPC and DPPC/DMTAP lipid bilayers in the presence of a DNA dodecamer. We observed the spontaneous formation of lipid bilayers from initial systems containing randomly placed lipids, water–counterions and DNA. In both the DPPC and DPPC/DMTAP simulations, the DNA molecule is located at the water–lipid headgroup interface, lying approximately parallel to the plane of the bilayer. We have also calculated the potential of mean force for transferring a DNA dodecamer through a DPPC/DMTAP bilayer. A high energetic barrier to DNA insertion into the hydrophobic core of the bilayer is observed. The DNA adopts a transmembrane orientation only in this region. Local bilayer deformation in the vicinity of the DNA molecule is observed, largely as a result of the DNA–DMTAP headgroup attraction.

On the Ability of PAMAM Dendrimers and Dendrimer/DNA Aggregates To Penetrate POPC Model Biomembranes

Poly(amido amine) (PAMAM) dendrimers have previously been shown, as cationic condensing agents of DNA, to have high potential for nonviral gene delivery. This study addresses two key issues for gene delivery: the interaction of the biomembrane with (i) the condensing agent (the cationic PAMAM dendrimer) and (ii) the corresponding dendrimer/DNA aggregate. Using in situ null ellipsometry and neutron reflection, parallel experiments were carried out involving dendrimers of generations 2 (G2), 4 (G4), and 6 (G6). The study demonstrates that free dendrimers of all three generations were able to traverse supported palmitoyloleoylphosphatidylcholine (POPC) bilayers deposited on silica surfaces. The model biomembranes were elevated from the solid surfaces upon dendrimer penetration, which offers a promising new way to generate more realistic model biomembranes where the contact with the supporting surface is reduced and where aqueous cavities are present beneath the bilayer. The largest dendrimer (G6) induced partial bilayer destruction directly upon penetration, whereas the smaller dendrimers (G2 and G4) leave the bilayer intact, so we propose that lower generation dendrimers have greater potential as transfection mediators. In addition to the experimental observations, coarse-grained simulations on the interaction between generation 3 (G3) dendrimers and POPC bilayers were performed in the absence and presence of a bilayer-supporting negatively charged surface that emulates the support. The simulations demonstrate that G3 is transported across free-standing POPC bilayers by direct penetration and not by endocytosis. The penetrability was, however, reduced in the presence of a surface, indicating that the membrane transport observed experimentally was not driven solely by the surface. The experimental reflection techniques were also applied to dendrimer/DNA aggregates of charge ratio = 0.5, and while G2/DNA and G4/DNA aggregates interact with POPC bilayers, G6/DNA displays no such interaction. These results indicate that, in contrast to free dendrimer molecules, dendrimer/DNA aggregates of low charge ratios are not able to traverse a membrane by direct penetration.

A positively charged channel within the Smc1/Smc3 hinge required for sister chromatid cohesion

Cohesins structural maintenance of chromosome 1 (Smc1) and Smc3 are rod‐shaped proteins with 50‐nm long intra‐molecular coiled‐coil arms with a heterodimerization domain at one end and an ABC‐like nucleotide‐binding domain (NBD) at the other. Heterodimerization creates V‐shaped molecules with a hinge at their centre. Inter‐connection of NBDs by Scc1 creates a tripartite ring within which, it is proposed, sister DNAs are entrapped. To investigate whether cohesins hinge functions as a possible DNA entry gate, we solved the crystal structure of the hinge from Mus musculus, which like its bacterial counterpart is characterized by a pseudo symmetric heterodimeric torus containing a small channel that is positively charged. Mutations in yeast Smc1 and Smc3 that together neutralize the channels charge have little effect on dimerization or association with chromosomes, but are nevertheless lethal. Our finding that neutralization reduces acetylation of Smc3, which normally occurs during replication and is essential for cohesion, suggests that the positively charged channel is involved in a major conformational change during S phase.

Enantiomeric resolution of supramolecular helicates with different surface topographies

The enantiomeric resolution of an extended range of di-metallo supramolecular triple-helical molecules are reported. The ligands for all complexes are symmetric with two units containing an aryl group linked via an imine bond to a pyridine. Alkyl substituents have been attached in different positions on the ligand backbone. Previous work on the parent compound, whose molecular formula is [Fe(2)(C(25)H(20)N(4))(3)]Cl4, showed that it could be resolved into enantiomerically pure solutions using cellulose and 20 mM aqueous sodium chloride. In this work a range of mobile phases have been investigated to see if the separation and speed of elution could be increased and the amount of NaCl co-eluted with the compounds decreased. Methanol, ethanol and acetonitrile were considered, together with aqueous NaCl : organic mixtures. Effective separation was most often achieved when using 90% acetonitrile : 10% 20 mM NaCl (aq) w/v, which gives scope for scaling up to incorporate the use of HPLC. The overall most efficient (i.e. fastest) separation was generally achieved where the cellulose column was packed with 20 mM NaCl (aq) and the column first eluted with 100% acetonitrile, then with 75% ethanol : 25% 20 mM NaCl (aq) until the M enantiomer had fully eluted and finally with 90% acetonitrile : 10% 20 mM NaCl (aq) until the P enantiomer had been collected. The sequence of eluents ensured minimum NaCl accompanying the enantiomers and minimum total solvent being required to elute the enantiomers, especially the second one, from the column. No helicate with a methyl group on the imine bond could be resolved and methyl groups on the pyridine rings also have an adverse effect on resolution.