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.