In pursuit of an accurate spatial and temporal model of biomolecules at the atomistic level: a perspective on computer simulation.

Despite huge advances in the computational techniques available for simulating biomolecules at the quantum-mechanical, atomistic and coarse-grained levels, there is still a widespread perception amongst the experimental community that these calculations are highly specialist and are not generally applicable by researchers outside the theoretical community. In this article, the successes and limitations of biomolecular simulation and the further developments that are likely in the near future are discussed. A brief overview is also provided of the experimental biophysical methods that are commonly used to probe biomolecular structure and dynamics, and the accuracy of the information that can be obtained from each is compared with that from modelling. It is concluded that progress towards an accurate spatial and temporal model of biomacromolecules requires a combination of all of these biophysical techniques, both experimental and computational.

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.

Probing the oligomeric state and interaction surfaces of Fukutin-I in dilauroylphosphatidylcholine bilayers.

Fukutin-I is localised to the endoplasmic reticulum or Golgi apparatus within the cell, where it is believed to function as a glycosyltransferase. Its localisation within the cell is thought to to be mediated by the interaction of its N-terminal transmembrane domain with the lipid bilayers surrounding these compartments, each of which possesses a distinctive lipid composition. However, it remains unclear at the molecular level how the interaction between the transmembrane domains of this protein and the surrounding lipid bilayer drives its retention within these compartments. In this work, we employed chemical cross-linking and fluorescence resonance energy transfer measurements in conjunction with multiscale molecular dynamics simulations to determine the oligomeric state of the protein within dilauroylphosphatidylcholine bilayers to identify interactions between the transmembrane domains and to ascertain any role these interactions may play in protein localisation. Our studies reveal that the N-terminal transmembrane domain of Fukutin-I exists as dimer within dilauroylphosphatidylcholine bilayers and that this interaction is driven by interactions between a characteristic TXXSS motif. Furthermore residues close to the N-terminus that have previously been shown to play a key role in the clustering of lipids are shown to also play a major role in anchoring the protein in the membrane.

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.