Introduction to the RSC Advances themed collection on New insights into biomolecular systems from large-scale simulations

Megan O'Mara, Sarah Rauscher and Stacey Wetmore introduce the RSC Advances themed collection on New insights into biomolecular systems from large-scale simulations.

Main-chain Macromolecular Hydrazone Photoswitches

Hydrazones-consisting of a dynamic imine bond and an acidic NH proton-have recently emerged as versatile photoswitches underpinned by their ability to form thermally bistable isomers, (Z) and (E), respectively. Herein, we introduce two photoresponsive homopolymers containing structurally different hydrazones as main-chain repeating units, synthesized via head-to-tail Acyclic Diene METathesis (ADMET) polymerization. Their key difference lies in the hydrazone design, specifically the location of the aliphatic arm connecting the rotor of the hydrazone photoswitch to the aliphatic polymer backbone. Critically, we demonstrate that their main photoresponsive property, i.e., their hydrodynamic volume, changes in opposite directions upon photoisomerization (λ=410 nm) in dilute solution. Further, the polymers-independent of the design of the individual hydrazone monomer-feature a photoswitchable glass transition temperature (Tg ) by close to 10 °C. The herein established design strategy allows to photochemically manipulate macromolecular properties by simple structural changes.

The efficacy of the analgesic GlyT2 inhibitor, ORG25543, is determined by two connected allosteric sites

Glycine Transporter 2 (GlyT2) inhibitors have shown considerable potential as analgesics for the treatment of neuropathic pain but also display considerable side effects. One potential source of side effects is irreversible inhibition. In this study, we have characterized the mechanism of ORG25543 inhibition of GlyT2 by first considering three potential ligand binding sites on GlyT2-the substrate site, the vestibule allosteric site and the lipid allosteric site. The three sites were tested using a combination of molecular dynamics simulations and analysis of the inhibition of glycine transport of a series point mutated GlyT2 using electrophysiological methods. We demonstrate that the lipid allosteric site on GlyT2 is the most likely binding site for ORG25543. We also demonstrate that cholesterol derived from the cell membrane can form specific interactions with inhibitor-bound transporters to form an allosteric network of regulatory sites. These observations will guide the future design of GlyT2 inhibitors with the objective of minimising on-target side effects and improving the therapeutic window for the treatment of patients suffering from neuropathic pain.

Extracellular vesicle lipids in cancer immunoevasion

Recent studies have revealed that cancer cell-derived extracellular vesicles (EVs) modulate immunological responses. Lipids have diverse biological functions, and are known to promote tumor malignancy. However, the immunoevasive roles of EV lipids in cancer progression remain poorly understood. Nevertheless, the study of cancer cell-derived EV lipids holds great promise for diagnostic and therapeutic interventions.

Structural basis of promiscuous substrate transport by Organic Cation Transporter 1

Organic Cation Transporter 1 (OCT1) plays a crucial role in hepatic metabolism by mediating the uptake of a range of metabolites and drugs. Genetic variations can alter the efficacy and safety of compounds transported by OCT1, such as those used for cardiovascular, oncological, and psychological indications. Despite its importance in drug pharmacokinetics, the substrate selectivity and underlying structural mechanisms of OCT1 remain poorly understood. Here, we present cryo-EM structures of full-length human OCT1 in the inward-open conformation, both ligand-free and drug-bound, indicating the basis for its broad substrate recognition. Comparison of our structures with those of outward-open OCTs provides molecular insight into the alternating access mechanism of OCTs. We observe that hydrophobic gates stabilize the inward-facing conformation, whereas charge neutralization in the binding pocket facilitates the release of cationic substrates. These findings provide a framework for understanding the structural basis of the promiscuity of drug binding and substrate translocation in OCT1.

Membrane cholesterol regulates inhibition and substrate transport by the glycine transporter, GlyT2

Membrane cholesterol binds to and modulates the function of various SLC6 neurotransmitter transporters, including stabilizing the outward-facing conformation of the dopamine and serotonin transporters. Here, we investigate how cholesterol binds to GlyT2 (SLC6A5), modulates glycine transport rate, and influences bioactive lipid inhibition of GlyT2. Bioactive lipid inhibitors are analgesics that bind to an allosteric site accessible from the extracellular solution when GlyT2 adopts an outward-facing conformation. Using molecular dynamics simulations, mutagenesis, and cholesterol depletion experiments, we show that bioactive lipid inhibition of glycine transport is modulated by the recruitment of membrane cholesterol to a binding site formed by transmembrane helices 1, 5, and 7. Recruitment involves cholesterol flipping from its membrane orientation, and insertion of the 3' hydroxyl group into the cholesterol binding cavity, close to the allosteric site. The synergy between cholesterol and allosteric inhibitors provides a novel mechanism of inhibition and a potential avenue for the development of potent GlyT2 inhibitors as alternative therapeutics for the treatment of neuropathic pain and therapeutics that target other SLC6 transporters.

The Impact of Antimicrobial Peptides on the Acinetobacter baumannii Inner Membrane Is Modulated by Lipid Polyunsaturation

The Gram-negative pathogen Acinetobacter baumannii is a primary contributor to nosocomial multi-drug-resistant (MDR) infections. To combat the rise of MDR infections, novel features of A. baumannii need to be considered for the development of new treatment options. One such feature is the preferential scavenging of exogenous lipids, including host-derived polyunsaturated fatty acids (PUFAs), for membrane phospholipid synthesis. These alterations in membrane composition impact both the lipid chemistry and the membrane biophysical properties. In this work we examine how antimicrobial peptides (AMPs) interact with the inner membranes of A. baumannii in the presence and absence of polyunsaturated phospholipids. Using coarse-grained molecular dynamics simulations of complex A. baumannii inner membrane models derived from lipidomes of bacteria grown in the presence and absence of PUFAs, we examine the impact of the adsorption of four prototypical AMPs (CAMEL, LL-37, pexiganan, and magainin-2) on the membrane biophysical properties. Our simulations reveal that the impact of AMP adsorption on the membrane biophysical properties was dependent on both the membrane composition and the specific AMP involved. Both lipid headgroup charge and tail unsaturation played important roles in driving the interactions that occurred both within the membrane and between the membrane and AMPs. The changes to the membrane biophysical properties also showed a complex relationship with the AMP's physical properties, such as AMP charge, chain length, and charge-to-mass ratio. Cumulatively, this work highlights the importance of studying AMPs using a complex membrane environment and provides insights into the mechanistic action of AMPs in polyunsaturated lipid-rich bacterial membranes.

Visible Light Reactive Single-Chain Nanoparticles

We introduce a single-chain nanoparticle (SCNP) system capable of catalyzing the photooxidation of nonpolar alkenes up to three times more efficiently than an equivalent small-molecule photosensitizer at an identical concentration. Specifically, we construct a polymer chain constituted of poly(ethylene glycol) methyl ether methacrylate and glycidyl methacrylate which we compact via multifunctional thiol-epoxide ligation and functionalize with Rose Bengal (RB) in a one pot reaction, affording SCNPs with a hydrophilic shell and hydrophobic photocatalytic regions. Photooxidation of the internal alkene in oleic acid proceeds under green light. RB confined within the SCNP is three times more effective for nonpolar alkenes than free RB in solution, which we hypothesize is due to the spatial proximity of the photosensitizing units to the substrate in the hydrophobic region. Our approach demonstrates that SCNP based catalysts can afford enhanced photocatalysis via confinement effects in a homogeneous reaction environment.

Protonophoric and mitochondrial uncoupling activity of aryl-carbamate substituted fatty acids

Aryl-urea substituted fatty acids are protonophores and mitochondrial uncouplers that utilise a urea-based synthetic anion transport moiety to carry out the protonophoric cycle. Herein we show that replacement of the urea group with carbamate, a functional group not previously reported to possess anion transport activity, produces analogues that retain the activity of their urea counterparts. Thus, the aryl-carbamate substituted fatty acids uncouple oxidative phosphorylation and inhibit ATP production by collapsing the mitochondrial proton gradient. Proton transport proceeds via self-assembly of the deprotonated aryl-carbamates into membrane permeable dimeric species, formed by intermolecular binding of the carboxylate group to the carbamate moiety. These results highlight the anion transport capacity of the carbamate functional group.

Dynamics of the Acinetobacter baumannii inner membrane under exogenous polyunsaturated fatty acid stress

Exogenous polyunsaturated fatty acids (PUFAs) are readily incorporated into the synthesis pathways of A. baumannii membrane phospholipids, where they contribute to reduced bacterial fitness and increased antimicrobial susceptibility. Here we examine the impact of PUFA membrane modification on membrane organisation and biophysical properties using coarse grained MARTINI simulations of chemically representative membrane models developed from mass-spectrometry datasets of an untreated, arachidonic acid (AA) treated and docosahexaenoic acid (DHA) treated A. baumannii membranes. Enzymatic integration of AA or DHA into phospholipids of the A. baumannii membrane resulted in modulation of membrane biophysical properties. Membrane thickness decreased slightly following PUFA treatment, concomitant with changes in the lateral area per lipid of each lipid headgroup class. PUFA treatment resulted in a decrease in membrane ordering and an increase in lipid lateral diffusion. Changes in lateral membrane organisation were observed in the PUFA treated membranes, with a concurrent increase in ordered cardiolipin domains and disordered PUFA-containing domains. Notably, separation between ordered and disordered domains was enhanced and was more pronounced for DHA relative to AA, providing a possible mechanism for greater antimicrobial action of DHA relative to AA observed experimentally. Furthermore, the membrane active antimicrobial, pentamidine, preferentially adsorbs to cardiolipin domains of the A. baumannii model membranes. This interaction, and membrane penetration of pentamidine, was enhanced following PUFA treatment. Cumulatively, this work explores the wide-ranging effects of PUFA incorporation on the A. baumannii membrane and provides a molecular basis for bacterial inner membrane disruption by PUFAs.

Hydroxy Groups Enhance [2]Rotaxane Anion Binding Selectivity

We report the synthesis of two [2]rotaxanes containing an interlocked three dimensional binding cavity formed from a pyridinium bis(amide) axle component containing two phenol donors, and an isophthalamide based macrocycle. In the competitive solvent mixture 1 : 1 CDCl₃ : CD₃OD, one of the receptors exhibits a much higher selectivity preference for chloride than an analogous rotaxane without the hydroxy groups. X‐ray crystal structures reveal the chloride anion guest encapsulated within the interlocked binding cavity, though not all of the hydrogen bond donors are utilised. Computational semi‐empirical simulations indicate that secondary intermolecular interactions occur between the axle hydroxy hydrogen bond donors and the [2]rotaxane macrocycle components, contributing to a more preorganised binding pocket, which may be responsible for the observed enhanced selectivity.

Pore structure controls stability and molecular flux in engineered protein cages

Protein cages are a common architectural motif used by living organisms to compartmentalize and control biochemical reactions. While engineered protein cages have featured in the construction of nanoreactors and synthetic organelles, relatively little is known about the underlying molecular parameters that govern stability and flux through their pores. In this work, we systematically designed 24 variants of the Thermotoga maritima encapsulin cage, featuring pores of different sizes and charges. Twelve pore variants were successfully assembled and purified, including eight designs with exceptional thermal stability. While negatively charged mutations were better tolerated, we were able to form stable assemblies covering a full range of pore sizes and charges, as observed in seven new cryo-EM structures at 2.5- to 3.6-Å resolution. Molecular dynamics simulations and stopped-flow experiments revealed the importance of considering both pore size and charge, together with flexibility and rate-determining steps, when designing protein cages for controlling molecular flux.

Conserved FimH mutations in the global Escherichia coli ST131 multi-drug resistant lineage weaken interdomain interactions and alter adhesin function

The binding of the type 1 fimbrial adhesin FimH to mannosylated receptors is allosterically regulated to enhance the fitness of uropathogenic Escherichia coli (UPEC) during urinary tract infection (UTI). Mutations in the two FimH domains (pilin and lectin) located outside the mannose binding pocket have been shown to influence mannose binding affinity, yet the details of the allostery mechanism are not fully elucidated. Here we characterised different FimH conformational states (termed low-affinity tense and high-affinity relaxed conformations) of natural FimH variants using molecular dynamics (MD) simulation techniques and report key structural dynamics differences between them. The clinically dominant FimH30 variant from the pandemic multidrug resistant E. coli ST131 lineage contains an R166H mutation that weakens FimH interdomain interactions and allows enhanced mannose interactions with pre-existing high-affinity relaxed conformations. When expressed in an isogenic ST131 strain background, FimH30 mediated high human cell adhesion and invasion, and enhanced biofilm formation over other variants. Collectively, our computational and experimental findings support a model of FimH protein allostery that is mediated by shifts in the pre-existing conformational equilibrium of FimH, additional to the sequential step-wise process of structural perturbations transmitted from one site to another within the protein. Importantly, it is the first study to shed light into how natural mutations in a clinically dominant FimH variant influence the protein’s conformational landscape optimising its function for ST131 fitness at intestinal and extraintestinal niches.

Simple synthetic route to an enzyme-inspired transesterification catalyst

Self-assembling transesterification catalyst inspired by the catalytic triad.

Effect of the Force Field on Molecular Dynamics Simulations of the Multidrug Efflux Protein P-Glycoprotein

Molecular dynamics (MD) simulations have been used extensively to study P-glycoprotein (P-gp), a flexible multidrug transporter that is a key player in the development of multidrug resistance to chemotherapeutics. A substantial body of literature has grown from simulation studies that have employed various simulation conditions and parameters, including AMBER, CHARMM, OPLS, GROMOS, and coarse-grained force fields, drawing conclusions from simulations spanning hundreds of nanoseconds. Each force field is typically parametrized and validated on different data and observables, usually of small molecules and peptides; there have been few comparisons of force field performance on large protein-membrane systems. Here we compare the conformational ensembles of P-gp embedded in a POPC/cholesterol bilayer generated over 500 ns of replicate simulation with five force fields from popular biomolecular families: AMBER 99SB-ILDN, CHARMM 36, OPLS-AA/L, GROMOS 54A7, and MARTINI. We find considerable differences among the ensembles with little conformational overlap, although they correspond to similar extents to structural data obtained from electron paramagnetic resonance and cross-linking studies. Moreover, each trajectory was still sampling new conformations at a high rate after 500 ns of simulation, suggesting the need for more sampling. This work highlights the need to consider known limitations of the force field used (e.g., biases toward certain secondary structures) and the simulation itself (e.g., whether sufficient sampling has been achieved) when interpreting accumulated results of simulation studies of P-gp and other transport proteins.

The structural basis of bacterial manganese import

Metal ions are essential for all forms of life. In prokaryotes, ATP-binding cassette (ABC) permeases serve as the primary import pathway for many micronutrients including the first-row transition metal manganese. However, the structural features of ionic metal transporting ABC permeases have remained undefined. Here, we present the crystal structure of the manganese transporter PsaBC from Streptococcus pneumoniae in an open-inward conformation. The type II transporter has a tightly closed transmembrane channel due to "extracellular gating" residues that prevent water permeation or ion reflux. Below these residues, the channel contains a hitherto unreported metal coordination site, which is essential for manganese translocation. Mutagenesis of the extracellular gate perturbs manganese uptake, while coordination site mutagenesis abolishes import. These structural features are highly conserved in metal-specific ABC transporters and are represented throughout the kingdoms of life. Collectively, our results define the structure of PsaBC and reveal the features required for divalent cation transport.

The role of plasmalogens, Forssman lipids, and sphingolipid hydroxylation in modulating the biophysical properties of the epithelial plasma membrane

A coarse-grain model of the epithelial plasma membrane was developed from high-resolution lipidomic data and simulated using the MARTINI force field to characterize its biophysical properties. Plasmalogen lipids, Forssman glycosphingolipids, and hydroxylated Forssman glycosphingolipids and sphingomyelin were systematically added to determine their structural effects. Plasmalogen lipids have a minimal effect on the overall biophysical properties of the epithelial plasma membrane. In line with the hypothesized role of Forssman lipids in the epithelial apical membrane, the introduction of Forssman lipids initiates the formation of glycosphingolipid-rich nanoscale lipid domains, which also include phosphatidylethanolamine (PE), sphingomyelin (SM), and cholesterol (CHOL). This decreases the lateral diffusion in the extracellular leaflet, as well as the area per lipid of domain forming lipids, most notably PE. Finally, hydroxylation of the Forssman glycosphingolipids and sphingomyelin further modulates the lateral organization of the membrane. Through comparison to the previously studied average and neuronal plasma membranes, the impact of membrane lipid composition on membrane properties was characterized. Overall, this study furthers our understanding of the biophysical properties of complex membranes and the impact of lipid diversity in modulating membrane properties.

The allosteric inhibition of glycine transporter 2 by bioactive lipid analgesics is controlled by penetration into a deep lipid cavity

The role of lipids in modulating membrane protein function is an emerging and rapidly growing area of research. The rational design of lipids that target membrane proteins for the treatment of pathological conditions is a novel extension in this field and provides a step forward in our understanding of membrane transporters. Bioactive lipids show considerable promise as analgesics for the treatment of chronic pain and bind to a high-affinity allosteric-binding site on the human glycine transporter 2 (GlyT2 or SLC6A5). Here, we use a combination of medicinal chemistry, electrophysiology, and computational modeling to develop a rational structure-activity relationship for lipid inhibitors and demonstrate the key role of the lipid tail interactions for GlyT2 inhibition. Specifically, we examine how lipid inhibitor head group stereochemistry, tail length, and double-bond position promote enhanced inhibition. Overall, the l-stereoisomer is generally a better inhibitor than the d-stereoisomer, longer tail length correlates with greater potency, and the position of the double bond influences the activity of the inhibitor. We propose that the binding of the lipid inhibitor deep into the allosteric-binding pocket is critical for inhibition. Furthermore, this provides insight into the mechanism of inhibition of GlyT2 and highlights how lipids can modulate the activity of membrane proteins by binding to cavities between helices. The principles identified in this work have broader implications for the development of a larger class of compounds that could target SLC6 transporters for disease treatment.

Investigating the lipid fingerprint of SLC6 neurotransmitter transporters: a comparison of dDAT, hDAT, hSERT, and GlyT2

The local lipid annulus, or "fingerprint", of four SLC6 transporters (dDAT, hDAT, hSERT, and GlyT2) embedded in a complex neuronal membrane were compared and characterised using molecular dynamics. Our analysis included the development of new tools to improve membrane leaflet detection and the analysis of leaflet-dependent properties. Overall, the lipid fingerprints of the four transporters are comprised of similar lipids when grouped by headgroup or tail saturation. The enrichment and depletion of specific lipids, including sites of cholesterol contacts, varies between transporters. The subtle differences in lipid fingerprints results in varying membrane biophysical properties near the transporter. Our results highlight that the lipid-fingerprint of SLC6 transporters in complex membranes is highly dependent on membrane composition. Our results further characterize how the presence and identity of membrane proteins affects the complex interplay of lipid-protein interactions, influencing the local lipid environment and membrane biophysical properties.

Aryl urea substituted fatty acids: a new class of protonophoric mitochondrial uncoupler that utilises a synthetic anion transporter

Respiring mitochondria establish a proton gradient across the mitochondrial inner membrane (MIM) that is used to generate ATP. Protein-independent mitochondrial uncouplers collapse the proton gradient and disrupt ATP production by shuttling protons back across the MIM in a protonophoric cycle. Continued cycling relies on the formation of MIM-permeable anionic species that can return to the intermembrane space after deprotonation in the mitochondrial matrix. Previously described protonophores contain acidic groups that are part of delocalised π-systems that provide large surfaces for charge delocalisation and facilitate anion permeation across the MIM. Here we present a new class of protonophoric uncoupler based on aryl-urea substituted fatty acids in which an acidic group and a π-system are separated by a long alkyl chain. The aryl-urea group in these molecules acts as a synthetic anion receptor that forms intermolecular hydrogen bonds with the fatty acid carboxylate after deprotonation. Dispersal of the negative charge across the aryl-urea system produces lipophilic dimeric complexes that can permeate the MIM and facilitate repeated cycling. Substitution of the aryl-urea group with lipophilic electron withdrawing groups is critical to complex lipophilicity and uncoupling activity. The aryl-urea substituted fatty acids represent the first biological example of mitochondrial uncoupling mediated by the interaction of a fatty acid and an anion receptor moiety, via self-assembly.

High resolution crystal structure of a KRAS promoter G-quadruplex reveals a dimer with extensive poly-A π-stacking interactions for small-molecule recognition

Aberrant KRAS signaling is a driver of many cancers and yet remains an elusive target for drug therapy. The nuclease hypersensitive element of the KRAS promoter has been reported to form secondary DNA structures called G-quadruplexes (G4s) which may play important roles in regulating KRAS expression, and has spurred interest in structural elucidation studies of the KRAS G-quadruplexes. Here, we report the first high-resolution crystal structure (1.6 Å) of a KRAS G-quadruplex as a 5'-head-to-head dimer with extensive poly-A π-stacking interactions observed across the dimer. Molecular dynamics simulations confirmed that the poly-A π-stacking interactions are also maintained in the G4 monomers. Docking and molecular dynamics simulations with two G4 ligands that display high stabilization of the KRAS G4 indicated the poly-A loop was a binding site for these ligands in addition to the 5'-G-tetrad. Given sequence and structural variability in the loop regions provide the opportunity for small-molecule targeting of specific G4s, we envisage this high-resolution crystal structure for the KRAS G-quadruplex will aid in the rational design of ligands to selectively target KRAS.

Understanding the Link between Lipid Diversity and the Biophysical Properties of the Neuronal Plasma Membrane

Cell membranes contain incredible diversity in the chemical structures of their individual lipid species and the ratios in which these lipids are combined to make membranes. Nevertheless, our current understanding of how each of these components affects the properties of the cell membrane remains elusive, in part due to the difficulties in studying the dynamics of membranes at high spatiotemporal resolution. In this work, we use coarse-grained molecular dynamics simulations to investigate how individual lipid species contribute to the biophysical properties of the neuronal plasma membrane. We progress through eight membranes of increasing chemical complexity, ranging from a simple POPC/CHOL membrane to a previously published neuronal plasma membrane [Ingólfsson, H. I., et al. (2017) Biophys. J. 113 (10), 2271-2280] containing 49 distinct lipid species. Our results show how subtle chemical changes can affect the properties of the membrane and highlight the lipid species that give the neuronal plasma membrane its unique biophysical properties. This work has potential far-reaching implications for furthering our understanding of cell membranes.

Cross-linking, DEER-spectroscopy and molecular dynamics confirm the inward facing state of P-glycoprotein in a lipid membrane

The drug efflux pump P-glycoprotein (P-gp) displays a complex transport mechanism involving multiple drug binding sites and two centres for nucleotide hydrolysis. Elucidating the molecular mechanism of transport remains elusive and the availability of P-gp structures in distinct natural and ligand trapped conformations will accelerate our understanding. The present investigation sought to provide biochemical data to validate specific features of these structures; with particular focus on the transmembrane domain that provides the transport conduit. Hence our focus was on transmembrane helices six and twelve (TM6/TM12), which are believed to participate in drug binding, as they line the central transport conduit and provide a direct link to the catalytic centres. A series of P-gp mutants were generated with a single cysteine in both TM6 and TM12 to facilitate measurement of inter-helical distances using cross-linking and DEER strategies. Experimental results were compared to published structures per se and those refined by MD simulations. This analysis revealed that the refined inward-facing murine structure (4M1M) of P-gp provides a good representation of the proximity, topography and relative motions of TM6 and TM12 in reconstituted human P-gp.

Comparing Nonbonded Metal Ion Models in the Divalent Cation Binding Protein PsaA

Divalent metal cations are essential for many biological processes; however, accurately modeling divalent metal ions has proved a significant challenge for molecular dynamics force fields. Here we show that the choice of ion model influences the observed dynamics in PsaA, a metal binding protein from Streptococcus pneumoniae. We conduct extensive unbiased simulations and free energy calculations of PsaA bound to its cognate ligand Mn²⁺ and inhibitory ligand Zn²⁺ using three nonbonded ion models: a 12-6 model, a 12-6-4 model, and a multisite model. The observed coordination geometries and metal binding dynamics are sensitive to the choice of ion model, with the most dramatic differences observed in free energy calculations of ion release. We show that the conformational ensemble of Mn-bound PsaA is more similar to the crystallographic metal bound open state. This work extends the current model of PsaA metal binding and provides a framework for the rationalization of experimentally determined metal binding behavior. Our findings support the use of the 12-6-4 ion model for further simulations of divalent cation binding proteins.

The Fats of Life: Using Computational Chemistry to Characterise the Eukaryotic Cell Membrane

Our current knowledge of the structural dynamics and complexity of lipid bilayers is still developing. Computational techniques, especially molecular dynamics simulations, have increased our understanding significantly as they allow us to model functions that cannot currently be experimentally resolved. Here we review available computational tools and techniques, the role of the major lipid species, insights gained into lipid bilayer structure and function from molecular dynamics simulations, and recent progress towards the computational modelling of the physiological complexity of eukaryotic lipid bilayers.

Identification of an allosteric binding site on the human glycine transporter, GlyT2, for bioactive lipid analgesics

The treatment of chronic pain is poorly managed by current analgesics, and there is a need for new classes of drugs. We recently developed a series of bioactive lipids that inhibit the human glycine transporter GlyT2 (SLC6A5) and provide analgesia in animal models of pain. Here, we have used functional analysis of mutant transporters combined with molecular dynamics simulations of lipid-transporter interactions to understand how these bioactive lipids interact with GlyT2. This study identifies a novel extracellular allosteric modulator site formed by a crevice between transmembrane domains 5, 7, and 8, and extracellular loop 4 of GlyT2. Knowledge of this site could be exploited further in the development of drugs to treat pain, and to identify other allosteric modulators of the SLC6 family of transporters.

Mutation p.R356Q in the Collybistin Phosphoinositide Binding Site Is Associated With Mild Intellectual Disability

The recruitment of inhibitory GABA_A receptors to neuronal synapses requires a complex interplay between receptors, neuroligins, the scaffolding protein gephyrin and the GDP-GTP exchange factor collybistin (CB). Collybistin is regulated by protein-protein interactions at the N-terminal SH3 domain, which can bind neuroligins 2/4 and the GABA_AR α2 subunit. Collybistin also harbors a RhoGEF domain which mediates interactions with gephyrin and catalyzes GDP-GTP exchange on Cdc42. Lastly, collybistin has a pleckstrin homology (PH) domain, which binds phosphoinositides, such as phosphatidylinositol 3-phosphate (PI3P/PtdIns3P) and phosphatidylinositol 4-monophosphate (PI4P/PtdIns4P). PI3P located in early/sorting endosomes has recently been shown to regulate the postsynaptic clustering of gephyrin and GABA_A receptors and consequently the strength of inhibitory synapses in cultured hippocampal neurons. This process is disrupted by mutations in the collybistin gene (ARHGEF9), which cause X-linked intellectual disability (XLID) by a variety of mechanisms converging on disrupted gephyrin and GABA_A receptor clustering at central synapses. Here we report a novel missense mutation (chrX:62875607C>T, p.R356Q) in ARHGEF9 that affects one of the two paired arginine residues in the PH domain that were predicted to be vital for binding phosphoinositides. Functional assays revealed that recombinant collybistin CB3_SH3- ^R356Q was deficient in PI3P binding and was not able to translocate EGFP-gephyrin to submembrane microaggregates in an in vitro clustering assay. Expression of the PI3P-binding mutants CB3_SH3- ^R356Q and CB3_SH3- ^R356N/R357N in cultured hippocampal neurones revealed that the mutant proteins did not accumulate at inhibitory synapses, but instead resulted in a clear decrease in the overall number of synaptic gephyrin clusters compared to controls. Molecular dynamics simulations suggest that the p.R356Q substitution influences PI3P binding by altering the range of structural conformations adopted by collybistin. Taken together, these results suggest that the p.R356Q mutation in ARHGEF9 is the underlying cause of XLID in the probands, disrupting gephyrin clustering at inhibitory GABAergic synapses via loss of collybistin PH domain phosphoinositide binding.

A potential new, stable state of the E-cadherin strand-swapped dimer in solution

E-cadherin is a transmembrane glycoprotein that facilitates inter-cellular adhesion in the epithelium. The ectodomain of the native structure is comprised of five repeated immunoglobulin-like domains. All E-cadherin crystal structures show the protein in one of three alternative conformations: a monomer, a strand-swapped trans homodimer and the so-called X-dimer, which is proposed to be a kinetic intermediate to forming the strand-swapped trans homodimer. However, previous studies have indicated that even once the trans strand-swapped dimer is formed, the complex is highly dynamic and the E-cadherin monomers may reorient relative to each other. Here, molecular dynamics simulations have been used to investigate the stability and conformational flexibility of the human E-cadherin trans strand-swapped dimer. In four independent, 100 ns simulations, the dimer moved away from the starting structure and converged to a previously unreported structure, which we call the Y-dimer. The Y-dimer was present for over 90% of the combined simulation time, suggesting that it represents a stable conformation of the E-cadherin dimer in solution. The Y-dimer conformation is stabilised by interactions present in both the trans strand-swapped dimer and X-dimer crystal structures, as well as additional interactions not found in any E-cadherin dimer crystal structures. The Y-dimer represents a previously unreported, stable conformation of the human E-cadherin trans strand-swapped dimer and suggests that the available crystal structures do not fully capture the conformations that the human E-cadherin trans homodimer adopts in solution.

Hydrogen bond-Driven Self-Assembly between Amidinium Cations and Carboxylate Anions: A Combined Molecular Dynamics, NMR Spectroscopy, and Single Crystal X-ray Diffraction Study

A combination of molecular dynamics (MD), NMR spectroscopy, and single crystal X-ray diffraction (SCXRD) techniques was used to probe the self-assembly of para- and meta-bis(amidinium) compounds with para-, meta-, and ortho-dicarboxylates. Good concordance was observed between the MD and experimental results. In DMSO solution, the systems form several rapidly exchanging assemblies, in part because a range of hydrogen bonding interactions is possible between the amidinium and carboxylate moieties. Upon crystallization, the majority of the systems form 1D supramolecular polymers, which are held together by short N-H⋅⋅⋅O hydrogen bonds.

Synthetically controlling dendrimer flexibility improves delivery of large plasmid DNA

Tools for editing the genome and epigenome have revolutionised the field of molecular biology and represent a new frontier in targeted therapeutic intervention. Although efficiencies and specificities of genome editing technologies have improved with the development of TALEs and CRISPR platforms, intracellular delivery of these larger constructs still remains a challenge using existing delivery agents. Viral vectors, including lentiviruses and adeno-associated viruses, as well as some non-viral strategies, such as cationic polymers and liposomes, are limited by packaging capacity, poor delivery, toxicity, and immunogenicity. We report a highly controlled synthetic strategy to engineer a flexible dendritic polymer using click chemistry to overcome the aforementioned delivery challenges associated with genome engineering technologies. Using a systematic approach, we demonstrate that high transfection efficiencies and packaging capacity can be achieved using this non-viral delivery methodology to deliver zinc fingers, TALEs and CRISPR/dCas9 platforms.

Hydrophobic Shielding Drives Catalysis of Hydride Transfer in a Family of F420H2-Dependent Enzymes

A family of flavin/deazaflavin-dependent oxidoreductases (FDORs) from mycobacteria has been recently characterized and found to play a variety of catalytic roles, including the activation of prodrugs such as the candidate anti-tuberculosis drug pretomanid (PA-824). However, our understanding of the catalytic mechanism used by these enzymes is relatively limited. To address this, we have used a combination of quantum mechanics and molecular dynamics calculations to study the catalytic mechanism of the activation of pretomanid by the deazaflavin-dependent nitroreductase (Ddn) from Mycobacterium tuberculosis. The preferred pathway involves an initial hydride transfer step from the deprotonated cofactor (i.e., F₄₂₀H^-), with subsequent protonation, before a series of spontaneous intramolecular reactions to form the final reactive nitrogen species. The most likely proton source is a hydroxonium ion within the solvent accessible active site. Intriguingly, catalysis of the rate-determining hydride transfer step is aided by three tyrosine residues that form a hydrophobic barrier around the active site that, upon reaction, is then disrupted to allow increased water accessibility to facilitate the subsequent proton transfer step. The catalytic mechanism we propose is consistent with previous experimental observations of the Ddn enzyme and will inform the design of improved prodrugs in the future.

Rv2074 is a novel F420 H2 -dependent biliverdin reductase in Mycobacterium tuberculosis

Bilirubin is a potent antioxidant that is produced from the reduction of the heme degradation product biliverdin. In mammalian cells and Cyanobacteria, NADH/NADPH-dependent biliverdin reductases (BVRs) of the Rossmann-fold have been shown to catalyze this reaction. Here, we describe the characterization of Rv2074 from Mycobacterium tuberculosis, which belongs to a structurally and mechanistically distinct family of F420 H2 -dependent BVRs (F-BVRs) that are exclusively found in Actinobacteria. We have solved the crystal structure of Rv2074 bound to its cofactor, F420 , and used this alongside molecular dynamics simulations, site-directed mutagenesis and NMR spectroscopy to elucidate its catalytic mechanism. The production of bilirubin by Rv2074 could exploit the anti-oxidative properties of bilirubin and contribute to the range of immuno-evasive mechanisms that have evolved in M. tuberculosis to allow persistent infection.

Structural and dynamic perspectives on the promiscuous transport activity of P-glycoprotein

The multidrug transporter P-glycoprotein (P-gp) is expressed in the blood-brain barrier endothelium where it effluxes a range of drug substrates, preventing their accumulation within the brain. P-gp has been studied extensively for 40 years because of its crucial role in the absorption, distribution, metabolism and elimination of a range of pharmaceutical compounds. Despite this, many aspects of the structure-function mechanism of P-gp are unresolved. Here we review the emerging role of molecular dynamics simulation techniques in our understanding of the membrane-embedded conformation of P-gp. We discuss its conformational plasticity in the presence and absence of ATP, and recent efforts to characterize the drug binding sites and uptake pathways.

Molecular basis for the interaction of the mammalian amino acid transporters B0AT1 and B0AT3 with their ancillary protein collectrin

Many solute carrier 6 (SLC6) family transporters require ancillary subunits to modify their expression and activity. The main apical membrane neutral amino acid transporters in mouse intestine and kidney, B(0)AT1 and B(0)AT3, require the ancillary protein collectrin or ACE2 for plasma membrane expression. Expression and activity of SLC6 neurotransmitter transporters are modulated by interaction with syntaxin 1A. Utilizing monocarboxylate-B(0)AT1/3 fusion constructs, we discovered that collectrin is also necessary for B(0)AT1 and B(0)AT3 catalytic function. Syntaxin 1A and syntaxin 3 inhibit the membrane expression of B(0)AT1 by competing with collectrin for access. A mutagenesis screening approach identified residues on trans-membrane domains 1α, 5, and 7 on one face of B(0)AT3 as a key region involved in interaction with collectrin. Mutant analysis established residues that were involved in collectrin-dependent functions as follows: plasma membrane expression of B(0)AT3, catalytic activation, or both. These results identify a potential binding site for collectrin and other SLC6 ancillary proteins.

Characterizing the conformational dynamics of metal-free PsaA using molecular dynamics simulations and electron paramagnetic resonance spectroscopy

Prokaryotic metal-ion receptor proteins, or solute-binding proteins, facilitate the acquisition of metal ions from the extracellular environment. Pneumococcal surface antigen A (PsaA) is the primary Mn(2+)-recruiting protein of the human pathogen Streptococcus pneumoniae and is essential for its in vivo colonization and virulence. The recently reported high-resolution structures of metal-free and metal-bound PsaA have provided the first insights into the mechanism of PsaA-facilitated metal binding. However, the conformational dynamics of metal-free PsaA in solution remain unknown. Here, we use continuous wave electron paramagnetic resonance (EPR) spectroscopy and molecular dynamics (MD) simulations to study the relative flexibility of the structural domains in metal-free PsaA and its distribution of conformations in solution. The results show that the crystal structure of metal-free PsaA is a good representation of the dominant conformation in solution, but the protein also samples structurally distinct conformations that are not captured by the crystal structure. Further, these results suggest that the metal binding site is both larger and more solvent exposed than indicated by the metal-free crystal structure. Collectively, this study provides atomic-resolution insight into the conformational dynamics of PsaA prior to metal binding and lays the groundwork for future EPR and MD based studies of PsaA in solution.

Identification of Possible Binding Sites for Morphine and Nicardipine on the Multidrug Transporter P-Glycoprotein Using Umbrella Sampling Techniques

The multidrug transporter P-glycoprotein (P-gp) is central to the development of multidrug resistance in cancer. While residues essential for transport and binding have been identified, the location, composition, and specificity of potential drug binding sites are uncertain. Here molecular dynamics simulations are used to calculate the free energy profile for the binding of morphine and nicardipine to P-gp. We show that morphine and nicardipine primarily interact with key residues implicated in binding and transport from mutational studies, binding at different but overlapping sites within the transmembrane pore. Their permeation pathways were distinct but involved overlapping sets of residues. The results indicate that the binding location and permeation pathways of morphine and nicardipine are not well separated and cannot be considered as unique. This has important implications for our understanding of substrate uptake and transport by P-gp. Our results are independent of the choice of starting structure and consistent with a range of experimental studies.

Dysregulation of transition metal ion homeostasis is the molecular basis for cadmium toxicity in Streptococcus pneumoniae

Cadmium is a transition metal ion that is highly toxic in biological systems. Although relatively rare in the Earth’s crust, anthropogenic release of cadmium since industrialization has increased biogeochemical cycling and the abundance of the ion in the biosphere. Despite this, the molecular basis of its toxicity remains unclear. Here we combine metal-accumulation assays, high-resolution structural data and biochemical analyses to show that cadmium toxicity, in Streptococcus pneumoniae, occurs via perturbation of first row transition metal ion homeostasis. We show that cadmium uptake reduces the millimolar cellular accumulation of manganese and zinc, and thereby increases sensitivity to oxidative stress. Despite this, high cellular concentrations of cadmium (~17 mM) are tolerated, with negligible impact on growth or sensitivity to oxidative stress, when manganese and glutathione are abundant. Collectively, this work provides insight into the molecular basis of cadmium toxicity in prokaryotes, and the connection between cadmium accumulation and oxidative stress.

The molecular basis for the high toxicity of cadmium is unclear. Here, Begg et al. use the bacterium Streptococcus pneumoniae as a model system, and show that cadmium uptake increases sensitivity to oxidative stress by reducing intracellular concentrations of manganese and zinc through different mechanisms.

Mechanism of activation of protein kinase JAK2 by the growth hormone receptor

Signaling from JAK (Janus kinase) protein kinases to STAT (signal transducers and activators of transcription) transcription factors is key to many aspects of biology and medicine, yet the mechanism by which cytokine receptors initiate signaling is enigmatic. We present a complete mechanistic model for activation of receptor-bound JAK2, based on an archetypal cytokine receptor, the growth hormone receptor. For this, we used fluorescence resonance energy transfer to monitor positioning of the JAK2 binding motif in the receptor dimer, substitution of the receptor extracellular domains with Jun zippers to control the position of its transmembrane (TM) helices, atomistic modeling of TM helix movements, and docking of the crystal structures of the JAK2 kinase and its inhibitory pseudokinase domain with an opposing kinase-pseudokinase domain pair. Activation of the receptor dimer induced a separation of its JAK2 binding motifs, driven by a ligand-induced transition from a parallel TM helix pair to a left-handed crossover arrangement. This separation leads to removal of the pseudokinase domain from the kinase domain of the partner JAK2 and pairing of the two kinase domains, facilitating trans-activation. This model may well generalize to other class I cytokine receptors.

The central cavity of ABCB1 undergoes alternating access during ATP hydrolysis

Understanding the process that underlies multidrug recognition and efflux by P-glycoprotein (ABCB1) remains a key biological challenge. Structural data have recently become available for the murine and Caenorhabditis elegans homologues of ABCB1; however all structures were obtained in the absence of nucleotide. A feature of these structures was the presence of a central cavity that is inaccessible from the extracellular face of the protein. To determine the conformational dynamics of this region several residues in transmembrane helices TM6 (331, 343 and 354) and TM12 (980) were mutated to cysteine. Based upon structural predictions, these residues are proposed to line, or reside proximal to, the central cavity. The mutant isoforms were labelled with a paramagnetic probe enabling the application of EPR spectroscopic methods. Power saturation EPR spectra were recorded in the presence of hydrophobic (O2 ) or hydrophilic (NiEDDA) quenching agents to study the local environment of each residue. ABCB1 was trapped in both its nucleotide-bound and post-hydrolytic conformations and EPR spectra were again recorded in the presence and absence of quenching agents. The EPR line shapes provide information on the movements of these residues within TM6 and TM12 during ATP hydrolysis. Rationalization of the data with molecular dynamic simulations indicates that the cavity is converted to a configuration open to the aqueous phase following nucleotide binding, thereby suggesting alternating access to the cavity on opposite sides of the membrane during translocation.