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.

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.

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.

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.

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.

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.

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.

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.

Loss, Gain and Altered Function of GlyR α2 Subunit Mutations in Neurodevelopmental Disorders.

Glycine receptors (GlyRs) containing the α2 subunit govern cell fate, neuronal migration and synaptogenesis in the developing cortex and spinal cord. Rare missense variants and microdeletions in the X-linked GlyR α2 subunit gene (GLRA2) have been associated with human autism spectrum disorder (ASD), where they typically cause a loss-of-function via protein truncation, reduced cell-surface trafficking and/or reduced glycine sensitivity (e.g., GLRA2Δex8-9 and extracellular domain variants p.N109S and p.R126Q). However, the GlyR α2 missense variant p.R323L in the intracellular M3-M4 domain results in a gain-of-function characterized by slower synaptic decay times, longer duration active periods and increases in channel conductance. This study reports the functional characterization of four missense variants in GLRA2 associated with ASD or developmental disorders (p.V-22L, p.N38K, p.K213E, p.T269M) using a combination of bioinformatics, molecular dynamics simulations, cellular models of GlyR trafficking and electrophysiology in artificial synapses. The GlyR α2V-22L variant resulted in altered predicted signal peptide cleavage and a reduction in cell-surface expression, suggestive of a partial loss-of-function. Similarly, GlyR α2N38K homomers showed reduced cell-surface expression, a reduced affinity for glycine and a reduced magnitude of IPSCs in artificial synapses. By contrast, GlyR α2K213E homomers showed a slight reduction in cell-surface expression, but IPSCs were larger, with faster rise/decay times, suggesting a gain-of-function. Lastly, GlyR α2T269M homomers exhibited a high glycine sensitivity accompanied by a substantial leak current, suggestive of an altered function that could dramatically enhance glycinergic signaling. These results may explain the heterogeneity of clinical phenotypes associated with GLRA2 mutations and reveal that missense variants can result in a loss, gain or alteration of GlyR α2 function. In turn, these GlyR α2 missense variants are likely to either negatively or positively deregulate cortical progenitor homeostasis and neuronal migration in the developing brain, leading to changes in cognition, learning, and memory.

Lipid-mediated antimicrobial resistance: a phantom menace or a new hope?

The proposition of a post-antimicrobial era is all the more realistic with the continued rise of antimicrobial resistance. The development of new antimicrobials is failing to counter the ever-increasing rates of bacterial antimicrobial resistance. This necessitates novel antimicrobials and drug targets. The bacterial cell membrane is an essential and highly conserved cellular component in bacteria and acts as the primary barrier for entry of antimicrobials into the cell. Although previously under-exploited as an antimicrobial target, the bacterial cell membrane is attractive for the development of novel antimicrobials due to its importance in pathogen viability. Bacterial cell membranes are diverse assemblies of macromolecules built around a central lipid bilayer core. This lipid bilayer governs the overall membrane biophysical properties and function of its membrane-embedded proteins. This mini-review will outline the mechanisms by which the bacterial membrane causes and controls resistance, with a focus on alterations in the membrane lipid composition, chemical modification of constituent lipids, and the efflux of antimicrobials by membrane-embedded efflux systems. Thorough insight into the interplay between membrane-active antimicrobials and lipid-mediated resistance is needed to enable the rational development of new antimicrobials. In particular, the union of computational approaches and experimental techniques for the development of innovative and efficacious membrane-active antimicrobials is explored.

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 CDCl3 : CD3 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.

Heteroleptic Tripalladium(II) Cages.

There is a concerted attempt to develop self-assembled metallo-cages of greater structural complexity, and heteroleptic PdII cages are emerging as prime candidates in these efforts. Most of these are dinuclear: few examples of higher nuclearity have been reported. We demonstrate here a robust method for the formation of tripalladium(II) cages from the 2 : 3 : 3 combination of a tritopic ligand, PdII , and a selection of ditopic ligands of the correct size and geometry.

Site of Cholesterol Oxidation Impacts Its Localization and Domain Formation in the Neuronal Plasma Membrane.

Cholesterol is integral to the structure of mammalian cell membranes. Oxidation of cholesterol alters how it behaves in the membrane and influences the membrane biophysical properties. Elevated levels of oxidized cholesterol are associated with neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, multiple sclerosis, and Huntington's disease. Previous work has investigated the impact of oxidized cholesterol in the context of simple model membrane systems. However, there is a growing body of literature that shows that complex membranes possessing physiological phospholipid distributions have different properties from those of binary or trinary model membranes. In the current work, the impact of oxidized cholesterol on the biophysical properties of a complex neuronal plasma membrane is investigated using coarse-grained Martini molecular dynamics simulations. Comparison of the native neuronal membrane to neuronal membranes containing 10% tail-oxidized or 10% head-oxidized cholesterol shows that the site of oxidization changes the behavior of the oxidized cholesterol in the membrane. Furthermore, species-specific domain formation is observed between each oxidized cholesterol and minor lipid classes. Although both tail-oxidized and head-oxidized cholesterols modulate the biophysical properties of the membrane, smaller changes are observed in the complex neuronal membrane than seen in the previous work on simple binary or trinary model membranes. This work highlights the presence of compensatory effects of lipid diversity in the complex neuronal membrane. Overall, this study improves our molecular-level understanding of the effects of oxidized cholesterol on the properties of neuronal tissue and emphasizes the importance of studying membranes with realistic lipid compositions.

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.

A Unique Sequence Is Essential for Efficient Multidrug Efflux Function of the MtrD Protein of Neisseria gonorrhoeae.

Antimicrobial resistance in Neisseria gonorrhoeae has reached an alarming level, severely impacting the effective treatment of gonorrhea. Belonging to the resistance-nodulation-cell division (RND) superfamily of efflux transporters, the MtrD membrane protein of N. gonorrhoeae provides resistance to a broad range of antimicrobial compounds. A unique feature of MtrD is an 11-residue sequence (from N917 to P927 [N917-P927]) that connects transmembrane helices (TMS) 9 and 10; this sequence is not present in homologous RND proteins. This study explores the structural and functional roles of the N917-P927 region by means of mutant analysis and molecular dynamics simulations. We show that N917-P927 plays a key role in modulating substrate access to the binding cleft and influences the overall orientation of the protein within the inner membrane necessary for optimal functioning. Removal of N917-P927 significantly reduced MtrD-mediated resistance to a range of antimicrobials and mutations of three single amino acids impacted MtrD-mediated multidrug resistance. Furthermore, molecular dynamics simulations showed deletion of N917-P927 in MtrD may dysregulate access of the substrate to the binding cleft and closure of the substrate-binding pocket during the transport cycle. These findings indicate that N917-P927 is a key region for interacting with the inner membrane, conceivably influencing substrate capture from the membrane-periplasm interface and thus is essential for full multidrug resistance capacity of MtrD. IMPORTANCE The historical sexually transmitted infection gonorrhea continues to be a major public health concern with an estimated global annual incidence of 86.9 million cases. N. gonorrhoeae has been identified by the World Health Organization as one of the 12 antimicrobial-resistant bacterial species that poses the greatest risk to human health. As the major efflux pump in gonococci, the MtrD transporter contributes to the cell envelope barrier in this organism and pumps antimicrobials from the periplasm and inner membrane, resulting in resistance. This study demonstrates that a unique region of the MtrD protein that connects TMS 9 and TMS 10 forms a structure that may interact with the inner membrane positioning TMS 9 and stabilizing the protein facilitating substrate capture from the inner membrane-periplasm interface. Analysis of mutants of this region identified that it was essential for MtrD-mediated multidrug resistance. Characterization of the structure and function of this unique local region of MtrD has implications for drug efflux mechanisms used by related proteins and is important knowledge for development of antibiotics that bypass efflux.

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.