Daniel A. Holdbrook

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.

Conformational dynamics and membrane interactions of the E. coli outer membrane protein FecA: A molecular dynamics simulation study

The TonB-dependent transporters mediate high-affinity binding and active transport of a variety of substrates across the outer membrane of Escherichia coli. The substrates transported by these proteins are large, scarce nutrients that are unable to gain entry into the cell by passive diffusion across the complex, asymmetric bilayer that constitutes the outer membrane. Experimental studies have identified loop regions that are essential for the correct functioning of these proteins. A number of these loops have been implicated in ligand binding. We report the first simulations of an E. coli outer membrane protein in an asymmetric model membrane that incorporates lipopolysaccharide (LPS) molecules. Comparative simulations of the apo and holo forms of the TonB-dependent transporter FecA in different membrane models enable us to identify the nature of the LPS-protein interactions and determine how these interactions impact upon the conformational dynamics of this protein. In particular, our simulations provide molecular-level insights into the influence of the environment and ligand on the dynamics of the functionally important loops of FecA. In addition, we provide insights into the nature of the protein-ligand interactions and ligand induced conformational change in FecA.

South-east Asian Zika virus strain linked to cluster of cases in Singapore, August 2016.

Zika virus (ZIKV) is an ongoing global public health emergency with 70 countries and territories reporting evidence of ZIKV transmission since 2015. On 27 August 2016, Singapore reported its first case of local ZIKV transmission and identified an ongoing cluster. Here, we report the genome sequences of ZIKV strains from two cases and find through phylogenetic analysis that these strains form an earlier branch distinct from the recent large outbreak in the Americas.

Stability and membrane orientation of the fukutin transmembrane domain: a combined multiscale molecular dynamics and circular dichroism study.

The N-terminal domain of fukutin-I has been implicated in the localization of the protein in the endoplasmic reticulum and Golgi Apparatus. It has been proposed to mediate this through its interaction with the thinner lipid bilayers found in these compartments. Here we have employed multiscale molecular dynamics simulations and circular dichroism spectroscopy to explore the structure, stability, and orientation of the short 36-residue N-terminus of fukutin-I (FK1TMD) in lipids with differing tail lengths. Our results show that FK1TMD adopts a stable helical conformation in phosphatidylcholine lipids when oriented with its principal axis perpendicular to the bilayer plane. The stability of the helix is largely insensitive to the lipid tail length, preventing hydrophobic mismatch by virtue of its mobility and ability to tilt within the lipid bilayers. This suggests that changes in FK1TMD tilt in response to bilayer properties may be implicated in the regulation of its trafficking. Coarse-grained simulations of the complex Golgi membrane suggest the N-terminal domain may induce the formation of microdomains in the surrounding membrane through its preferential interaction with 1,2-dipalmitoyl-sn-glycero-3-phosphatidylinositol 4,5-bisphosphate lipids.

The NorM MATE Transporter from N. gonorrhoeae: Insights into Drug and Ion Binding from Atomistic Molecular Dynamics Simulations

The multidrug and toxic compound extrusion transporters extrude a wide variety of substrates out of both mammalian and bacterial cells via the electrochemical gradient of protons and cations across the membrane. The substrates transported by these proteins include toxic metabolites and antimicrobial drugs. These proteins contribute to multidrug resistance in both mammalian and bacterial cells and are therefore extremely important from a biomedical perspective. Although specific residues of the protein are known to be responsible for the extrusion of solutes, mechanistic details and indeed structures of all the conformational states remain elusive. Here, we report the first, to our knowledge, simulation study of the recently resolved x-ray structure of the multidrug and toxic compound extrusion transporter, NorM from Neisseria gonorrhoeae (NorM_NG). Multiple, atomistic simulations of the unbound and bound forms of NorM in a phospholipid lipid bilayer allow us to identify the nature of the drug-protein/ion-protein interactions, and secondly determine how these interactions contribute to the conformational rearrangements of the protein. In particular, we identify the molecular rearrangements that occur to enable the Na(+) ion to enter the cation-binding cavity even in the presence of a bound drug molecule. These include side chain flipping of a key residue, GLU-261 from pointing toward the central cavity to pointing toward the cation binding side when bound to a Na(+) ion. Our simulations also provide support for cation binding in the drug-bound and apo states of NorM_NG.

Stability and membrane interactions of an autotransport protein: MD simulations of the Hia translocator domain in a complex membrane environment

Hia is a trimeric autotransporter found in the outer membrane of Haemphilus influenzae. The X-ray structure of Hia translocator domain revealed each monomer to consist of an α-helix connected via a loop to a 4-stranded β-sheet, thus the topology of the trimeric translocator domain is a 12-stranded β-barrel containing 3 α-helices that protrude from the mouth of the β-barrel into the extracellular medium. Molecular dynamics simulations of the Hia monomer and trimer have been employed to explore the interactions between the helices, β-barrel and connecting loops that may contribute to the stability of the trimer. In simulations of the Hia monomer we show that the central α-helix may stabilise the fold of the 4-stranded β-sheet. In simulations of the Hia trimer, a H-bond network involving residues in the β-barrel, α-helices and loops has been identified as providing stability for the trimeric arrangement of the monomers. Glutamine residues located in the loops connecting the α-helices to the β-barrel are orientated in a triangular arrangement such that each forms 2 hydrogen bonds to each of the corresponding glutamines in the other loops. In the absence of the loops, the β-barrel becomes distorted. Simulations show that while the trimeric translocator domain β-barrel is inherently flexible, it is unlikely to accommodate the passenger domain in a folded conformation. Simulations of Hia in an asymmetric model of the outer membrane have revealed membrane-protein interactions that anchor the protein within its native membrane environment.

Dynamics of Crowded Vesicles: Local and Global Responses to Membrane Composition.

The bacterial cell envelope is composed of a mixture of different lipids and proteins, making it an inherently complex organelle. The interactions between integral membrane proteins and lipids are crucial for their respective spatial localization within bacterial cells. We have employed microsecond timescale coarse-grained molecular dynamics simulations of vesicles of varying sizes and with a range of protein and lipid compositions, and used novel approaches to measure both local and global system dynamics, the latter based on spherical harmonics analysis. Our results suggest that both hydrophobic mismatch, enhanced by embedded membrane proteins, and curvature based sorting, due to different modes of undulation, may drive assembly in vesicular systems. Interestingly, the modes of undulation of the vesicles were found to be altered by the specific protein and lipid composition of the vesicle. Strikingly, lipid dynamics were shown to be coupled to proteins up to 6 nm from their surface, a substantially larger distance than has previously been observed, resulting in multi-layered annular rings enriched with particular types of phospholipid. Such large protein-lipid complexes may provide a mechanism for long-range communication. Given the complexity of bacterial membranes, our results suggest that subtle changes in lipid composition may have major implications for lipid and protein sorting under a curvature-based membrane-sorting model.

The membranes of Gram-negative bacteria: progress in molecular modelling and simulation

Molecular modelling and simulations have been employed to study the membranes of Gram-negative bacteria for over 20 years. Proteins native to these membranes, as well as antimicrobial peptides and drug molecules have been studied using molecular dynamics simulations in simple models of membranes, usually only comprising one lipid species. Thus, traditionally, the simulations have reflected the majority of in vitro membrane experimental setups, enabling observations from the latter to be rationalized at the molecular level. In the last few years, the sophistication and complexity of membrane models have improved considerably, such that the heterogeneity of the lipid and protein composition of the membranes can now be considered both at the atomistic and coarse-grain levels of granularity. Importantly this means relevant biology is now being retained in the models, thereby linking the in silico and in vivo scenarios. We discuss recent progress in simulations of proteins in simple lipid bilayers, more complex membrane models and finally describe some efforts to overcome timescale limitations of atomistic molecular dynamics simulations of bacterial membranes.

Revisiting the interaction between the chaperone Skp and lipopolysaccharide

The bacterial outer membrane comprises two main classes of components, lipids and membrane proteins. These nonsoluble compounds are conveyed across the aqueous periplasm along specific molecular transport routes: the lipid lipopolysaccharide (LPS) is shuttled by the Lpt system, whereas outer membrane proteins (Omps) are transported by chaperones, including the periplasmic Skp. In this study, we revisit the specificity of the chaperone-lipid interaction of Skp and LPS. High-resolution NMR spectroscopy measurements indicate that LPS interacts with Skp nonspecifically, accompanied by destabilization of the Skp trimer and similar to denaturation by the nonnatural detergent lauryldimethylamine-N-oxide (LDAO). Bioinformatic analysis of amino acid conservation, structural analysis of LPS-binding proteins, and MD simulations further confirm the absence of a specific LPS binding site on Skp, making a biological relevance of the interaction unlikely. Instead, our analysis reveals a highly conserved salt-bridge network, which likely has a role for Skp function.

Multiscale molecular dynamics simulation approaches to the structure and dynamics of viruses

Viral pathogens are a significant source of human morbidity and mortality, and have a major impact on societies and economies around the world. One of the challenges inherent in targeting these pathogens with drugs is the tight integration of the viral life cycle with the hosts cellular machinery. However, the reliance of the virus on the host cell replication machinery is also an opportunity for therapeutic targeting, as successful entry- and exit-inhibitors have demonstrated. An understanding of the extracellular and intracellular structure and dynamics of the virion - as well as of the entry and exit pathways in host and vector cells - is therefore crucial to the advancement of novel antivirals. In recent years, advances in computing architecture and algorithms have begun to allow us to use simulations to study the structure and dynamics of viral ultrastructures at various stages of their life cycle in atomistic or near-atomistic detail. In this review, we outline specific challenges and solutions that have emerged to allow for structurally detailed modelling of viruses in silico. We focus on the history and state of the art of atomistic and coarse-grained approaches to simulate the dynamics of the large, macromolecular structures associated with viral infection, and on their usefulness in explaining and expanding upon experimental data. We discuss the types of interactions that need to be modeled to describe major components of the virus particle and advances in modelling techniques that allow for the treatment of these systems, highlighting recent key simulation studies.