KLHL3 is a BTB-BACK-Kelch family protein that serves as a substrate adapter in Cullin3 (Cul3) E3 ubiquitin ligase complexes. KLHL3 is highly expressed in distal nephron tubules where it is involved in the regulation of electrolyte homeostasis and blood pressure. Mutations in KLHL3 have been identified in patients with inherited hypertension disorders, and several of the disease-associated mutations are located in the presumed Cul3 binding region. Here, we report the crystal structure of a complex between the KLHL3 BTB-BACK domain dimer and two copies of an N terminal fragment of Cul3. We use isothermal titration calorimetry to directly demonstrate that several of the disease mutations in the KLHL3 BTB-BACK domains disrupt the association with Cul3. Both the BTB and BACK domains contribute to the Cul3 interaction surface, and an extended model of the dimeric CRL3 complex places the two E2 binding sites in a suprafacial arrangement with respect to the presumed substrate-binding sites.

Uracil (U) can be found in DNA as a mismatch paired either to adenine (A) or to guanine (G). Removal of U from DNA is performed by a class of enzymes known as uracil-DNA-glycosylases (UDG). Recent studies suggest that recognition of U-A and U-G mismatches by UDG takes place via an extra-helical mechanism. In this work, we use molecular dynamics simulations to analyze the structure, dynamics and open base pair kinetics of U-A base pairs relative to their natural T-A counterpart in 12 dodecamers. Our results show that the presence of U does not alter the local conformation of B-DNA. Breathing dynamics and base pair closing kinetics are only weakly dependent on the presence of U versus T, with open T-A and U-A pairs lifetimes in the nanosecond timescale. Additionally, we observed spontaneous base flipping in U-A pairs. We analyze the structure and dynamics for this event and compare the results to available crystallographic data of open base pair conformations. Our results are in agreement with both structural and kinetic data derived from NMR imino proton exchange measurements, providing the first detailed description at the molecular level of elusive events such as spontaneous base pair opening and flipping in mismatched U-A sequences in DNA. Based on these results, we propose that base pair flipping can occur spontaneously at room temperature via a 3-step mechanism with an open base pair intermediate. Implications for the molecular basis of U recognition by UDG are discussed.

The BTB domain (also known as the POZ domain) is a versatile protein-protein interaction motif that participates in a wide range of cellular functions, including transcriptional regulation, cytoskeleton dynamics, ion channel assembly and gating, and targeting proteins for ubiquitination. Several BTB domain structures have been experimentally determined, revealing a highly conserved core structure.

The use of detergents for the structural study of membrane proteins is discussed with an emphasis on practical issues relating to membrane solubilization, protein aggregation, detergent purity and detergent quantitation. Detergents are useful reagents as mimics of lipid bilayers because of their self-assembling properties, but as a result, they have complex properties in solution. It can be difficult to maintain a solubilized membrane protein in a native conformational state, and the non-specific aggregation of detergent-solubilized proteins is a common problem. Empirical "stability screens" can be helpful in choosing which detergents, and which detergent concentrations, may be optimal for a given system.

Cations play key roles in regulating G-protein-coupled receptors (GPCRs), although their mechanisms are poorly understood. Here, 19F NMR is used to delineate the effects of cations on functional states of the adenosine A2A GPCR. While Na+ reinforces an inactive ensemble and a partial-agonist stabilized state, Ca2+ and Mg2+ shift the equilibrium toward active states. Positive allosteric effects of divalent cations are more pronounced with agonist and a G-protein-derived peptide. In cell membranes, divalent cations enhance both the affinity and fraction of the high affinity agonist-bound state. Molecular dynamics simulations suggest high concentrations of divalent cations bridge specific extracellular acidic residues, bringing TM5 and TM6 together at the extracellular surface and allosterically driving open the G-protein-binding cleft as shown by rigidity-transmission allostery theory. An understanding of cation allostery should enable the design of allosteric agents and enhance our understanding of GPCR regulation in the cellular milieu.

The protein elastin imparts extensibility, elastic recoil, and resilience to tissues including arterial walls, skin, lung alveoli, and the uterus. Elastin and elastin-like peptides are hydrophobic, disordered, and undergo liquid-liquid phase separation upon self-assembly. Despite extensive study, the structure of elastin remains controversial. We use molecular dynamics simulations on a massive scale to elucidate the structural ensemble of aggregated elastin-like peptides. Consistent with the entropic nature of elastic recoil, the aggregated state is stabilized by the hydrophobic effect. However, self-assembly does not entail formation of a hydrophobic core. The polypeptide backbone forms transient, sparse hydrogen-bonded turns and remains significantly hydrated even as self-assembly triples the extent of non-polar side chain contacts. Individual chains in the assembly approach a maximally-disordered, melt-like state which may be called the liquid state of proteins. These findings resolve long-standing controversies regarding elastin structure and function and afford insight into the phase separation of disordered proteins.

The Saccharomyces cerevisiae respiratory supercomplex factor 2 (Rcf2) is a 224-residue protein located in the mitochondrial inner membrane where it is involved in the formation of supercomplexes composed of cytochrome bc1 and cytochrome c oxidase. We previously demonstrated that Rcf2 forms a dimer in dodecylphosphocholine micelles, and here we report the solution NMR structure of this Rcf2 dimer. Each Rcf2 monomer has two soluble α helices and five putative transmembrane (TM) α helices, including an unexpectedly charged TM helix at the C terminus, which mediates dimer formation. The NOE contacts indicate the presence of inter-monomer salt bridges and hydrogen bonds at the dimer interface, which stabilize the Rcf2 dimer structure. Moreover, NMR chemical shift change mapping upon lipid titrations as well as molecular dynamics analysis reveal possible structural changes upon embedding Rcf2 into a native lipid environment. Our results contribute to the understanding of respiratory supercomplex formation and regulation.

Saposin B (also known as cerebroside sulfate activator or CSAct) is a small non-enzymatic glycoprotein required for the breakdown of cerebroside sulfates (sulfatides) in lysosomes. Saposin B contains three intramolecular disulfide bridges, exists as a dimer and is remarkably heat, protease, and pH stable. We have expressed the protein in a thioredoxin reductase deficient strain of Escherichia coli and purified the protein by heat treatment, followed by ion-exchange, gel filtration, and hydrophobic interaction chromatographies. The protein is properly folded as judged by the observed disulfide bond topology, the hydrogen-deuterium exchange rate, and the level of stimulation of sulfatide hydrolysis by arylsulfatase A. Crystals of human saposin B were grown by vapor diffusion and diffract to a resolution of 2.2A. Despite obtaining only merohedrally twinned P3(1) native crystals, an untwined seleomethionine-substituted crystal belonging to space group P3(1)21 was also grown. The three-dimensional structure of saposin B protein will provide insights into how this 79 amino acid protein is able to solubilize relatively large membrane-bound lipid ligands.

Acid sphingomyelinase (ASM) is a lysosomal phosphodiesterase that catalyzes the hydrolysis of sphingomyelin to produce ceramide and phosphocholine. While other lysosomal sphingolipid hydrolases require a saposin activator protein for full activity, the ASM polypeptide incorporates a built-in N-terminal saposin domain and does not require an external activator protein. Here, we report the crystal structure of human ASM and describe the organization of the three main regions of the enzyme: the N-terminal saposin domain, the proline-rich connector, and the catalytic domain. The saposin domain is tightly associated along an edge of the large, bowl-shaped catalytic domain and adopts an open form that exposes a hydrophobic concave surface approximately 30Å from the catalytic center. The calculated electrostatic potential of the enzyme is electropositive at the acidic pH of the lysosome, consistent with the strict requirement for the presence of acidic lipids in target membranes. Docking studies indicate that sphingomyelin binds with the ceramide-phosphate group positioned at the binuclear zinc center and molecular dynamic simulations indicate that the intrinsic flexibility of the saposin domain is important for monomer-dimer exchange and for membrane interactions. Overall, ASM uses a combination of electrostatic and hydrophobic interactions to cause local disruptions of target bilayers in order to bring the lipid headgroup to the catalytic center in a membrane-bound reaction.

Membrane-less organelles in cells are large, dynamic protein/protein or protein/RNA assemblies that have been reported in some cases to have liquid droplet properties. However, the molecular interactions underlying the recruitment of components are not well understood. Herein, we study how the ability to form higher-order assemblies influences the recruitment of the speckle-type POZ protein (SPOP) to nuclear speckles. SPOP, a cullin-3-RING ubiquitin ligase (CRL3) substrate adaptor, self-associates into higher-order oligomers; that is, the number of monomers in an oligomer is broadly distributed and can be large. While wild-type SPOP localizes to liquid nuclear speckles, self-association-deficient SPOP mutants have a diffuse distribution in the nucleus. SPOP oligomerizes through its BTB and BACK domains. We show that BTB-mediated SPOP dimers form linear oligomers via BACK domain dimerization, and we determine the concentration-dependent populations of the resulting oligomeric species. Higher-order oligomerization of SPOP stimulates CRL3(SPOP) ubiquitination efficiency for its physiological substrate Gli3, suggesting that nuclear speckles are hotspots of ubiquitination. Dynamic, higher-order protein self-association may be a general mechanism to concentrate functional components in membrane-less cellular bodies.

As part of the mitochondrial respiratory chain, cytochrome c oxidase utilizes the energy produced by the reduction of O2 to water to fuel vectorial proton transport. The mechanism coupling proton pumping to redox chemistry is unknown. Recent advances have provided evidence that each of the four observable transitions in the complex catalytic cycle consists of a similar sequence of events. However, the physico-chemical basis underlying this recurring sequence has not been identified. We identify this recurring pattern based on a comprehensive model of the catalytic cycle derived from the analysis of oxygen chemistry and available experimental evidence. The catalytic cycle involves the periodic repetition of a sequence of three states differing in the spatial distribution of charge in the active site: [0|1], [1|0], and [1|1], where the total charge of heme a and the binuclear center appears on the left and on the right, respectively. This sequence recurs four times per turnover despite differences in the redox chemistry. This model leads to a simple, robust, and reproducible sequence of electron and proton transfer steps and rationalizes the pumping mechanism in terms of electrostatic coupling of proton translocation to redox chemistry. Continuum electrostatic calculations support the proposed mechanism and suggest an electrostatic origin for the decoupled and inactive phenotypes of ionic mutants in the principal proton-uptake pathway.

The intracellular transport of cholesterol is subject to tight regulation. The structure of the lysosomal integral membrane protein type 2 (LIMP-2, also known as SCARB2) reveals a large cavity that traverses the molecule and resembles the cavity in SR-B1 that mediates lipid transfer. The detection of cholesterol within the LIMP-2 structure and the formation of cholesterol-like inclusions in LIMP-2 knockout mice suggested the possibility that LIMP2 transports cholesterol in lysosomes. We present results of molecular modeling, crosslinking studies, microscale thermophoresis and cell-based assays that support a role of LIMP-2 in cholesterol transport. We show that the cavity in the luminal domain of LIMP-2 can bind and deliver exogenous cholesterol to the lysosomal membrane and later to lipid droplets. Depletion of LIMP-2 alters SREBP-2-mediated cholesterol regulation, as well as LDL-receptor levels. Our data indicate that LIMP-2 operates in parallel with Niemann Pick (NPC)-proteins, mediating a slower mode of lysosomal cholesterol export.

Ral Guanine Nucleotide Dissociation Stimulator Like 1 (RGL1) is a RAS effector protein that activates Ral GTPase by stimulating nucleotide exchange. Most structures of RAS-effector complexes are for the HRAS isoform; relatively few KRAS-effector structures have been solved, even though KRAS mutations are more frequent in human cancers. We determined crystal structures of KRAS/RGL1-RAS-association (RA) domain complexes and characterized the interaction in solution using nuclear magnetic resonance spectroscopy, size-exclusion chromatography combined with multi-angle light scattering and biolayer interferometry. We report structures of wild-type KRAS and the oncogenic G12V mutant in complex with the RA domain of RGL1 at < 2 Å resolution. KRASWT/RGL1-RA crystallized as a 1:1 heterodimer, whilst KRASG12V/RGL1-RA crystallized as a heterotetrameric structure in which RGL1-RA dimerized via domain-swapping the C-terminal beta-strand. Solution data indicated that KRASWT and KRASG12V in complex with RGL1-RA both exist predominantly as 1:1 dimers, while tetramerization occurs through very slow association. Through detailed structural analyses, the distance and angle between RAS α1 helix and RBD/RA α1 helix were found to differ significantly among RAS and RBD/RA complexes. The KRAS/RGL1-RA structures possess some of the largest α1RAS/α1Effector distances (21.7-22.2 Å), whereas the corresponding distances in previously reported RAS/RAF complexes are significantly shorter (15.2-17.7 Å). Contact map analysis identified unique structural signatures involving contacts between the β1-β2 loop of RA and the α1 helix of RAS, clearly distinguishing the KRAS/RGL1-RA (and other RAS/RA complexes) from RAS/RBD complexes. These results demonstrate that RAS effectors employ an assortment of finely-tuned docking surfaces to achieve optimal interactions with RAS.

The recent elucidation of atomistic structures of Cl- channel CFTR provides opportunities for understanding the molecular basis of cystic fibrosis. Despite having been activated through phosphorylation and provided with ATP ligands, several near-atomistic cryo-EM structures of CFTR are in a closed state, as inferred from the lack of a continuous passage through a hydrophobic bottleneck region located in the extracellular portion of the pore. Here, we present repeated, microsecond-long molecular dynamics simulations of human CFTR solvated in a lipid bilayer and aqueous NaCl. At equilibrium, Cl- ions enter the channel through a lateral intracellular portal and bind to two distinct cationic sites inside the channel pore but do not traverse the narrow, de-wetted bottleneck. Simulations conducted in the presence of a strong hyperpolarizing electric field led to spontaneous Cl- translocation events through the bottleneck region of the channel, suggesting that the protein relaxed to a functionally open state. Conformational changes of small magnitude involving transmembrane helices 1 and 6 preceded ion permeation through diverging exit routes at the extracellular end of the pore. The pore bottleneck undergoes wetting prior to Cl- translocation, suggesting that it acts as a hydrophobic gate. Although permeating Cl- ions remain mostly hydrated, partial dehydration occurs at the binding sites and in the bottleneck. The observed Cl- pathway is largely consistent with the loci of mutations that alter channel conductance, anion binding, and ion selectivity, supporting the model of the open state of CFTR obtained in the present study.

The ability of enzymes to distinguish between fatty acyl groups can involve molecular measuring devices termed hydrocarbon rulers, but the molecular basis for acyl-chain recognition in any membrane-bound enzyme remains to be defined. PagP is an outer membrane acyltransferase that helps pathogenic bacteria to evade the host immune response by transferring a palmitate chain from a phospholipid to lipid A (endotoxin). PagP can distinguish lipid acyl chains that differ by a single methylene unit, indicating that the enzyme possesses a remarkably precise hydrocarbon ruler. We present the 1.9 A crystal structure of PagP, an eight-stranded beta-barrel with an unexpected interior hydrophobic pocket that is occupied by a single detergent molecule. The buried detergent is oriented normal to the presumed plane of the membrane, whereas the PagP beta-barrel axis is tilted by approximately 25 degrees. Acyl group specificity is modulated by mutation of Gly88 lining the bottom of the hydrophobic pocket, thus confirming the hydrocarbon ruler mechanism for palmitate recognition. A striking structural similarity between PagP and the lipocalins suggests an evolutionary link between these proteins.

Glioblastoma is one of the most challenging forms of cancer to treat. Here we describe a computational platform that integrates the analysis of copy number variations and somatic mutations and unravels the landscape of in-frame gene fusions in glioblastoma. We found mutations with loss of heterozygosity in LZTR1, encoding an adaptor of CUL3-containing E3 ligase complexes. Mutations and deletions disrupt LZTR1 function, which restrains the self renewal and growth of glioma spheres that retain stem cell features. Loss-of-function mutations in CTNND2 target a neural-specific gene and are associated with the transformation of glioma cells along the very aggressive mesenchymal phenotype. We also report recurrent translocations that fuse the coding sequence of EGFR to several partners, with EGFR-SEPT14 being the most frequent functional gene fusion in human glioblastoma. EGFR-SEPT14 fusions activate STAT3 signaling and confer mitogen independence and sensitivity to EGFR inhibition. These results provide insights into the pathogenesis of glioblastoma and highlight new targets for therapeutic intervention.

Invading bacteria from the Neisseriaceae, Acinetobacteriaceae, Bordetellaceae and Moraxellaceae families express the conserved outer-membrane zinc transporter zinc-uptake component D (ZnuD) to overcome nutritional restriction imposed by the host organism during infection. Here we demonstrate that ZnuD is required for efficient systemic infections by the causative agent of bacterial meningitis, Neisseria meningitidis, in a mouse model. We also combine X-ray crystallography and molecular dynamics simulations to gain insight into the mechanism of zinc recognition and transport across the bacterial outer-membrane by ZnuD. Because ZnuD is also considered a promising vaccine candidate against N. meningitidis, we use several ZnuD structural intermediates to map potential antigenic epitopes, and propose a mechanism by which ZnuD can maintain high sequence conservation yet avoid immune recognition by altering the conformation of surface-exposed loops.

The topological similarity of voltage-gated proton channels (H(V)1s) to the voltage-sensing domain (VSD) of other voltage-gated ion channels raises the central question of whether H(V)1s have a similar structure. We present the construction and validation of a homology model of the human H(V)1 (hH(V)1). Multiple structural alignment was used to construct structural models of the open (proton-conducting) state of hH(V)1 by exploiting the homology of hH(V)1 with VSDs of K(+) and Na(+) channels of known three-dimensional structure. The comparative assessment of structural stability of the homology models and their VSD templates was performed using massively repeated molecular dynamics simulations in which the proteins were allowed to relax from their initial conformation in an explicit membrane mimetic. The analysis of structural deviations from the initial conformation based on up to 125 repeats of 100-ns simulations for each system reveals structural features consistently retained in the homology models and leads to a consensus structural model for hH(V)1 in which well-defined external and internal salt-bridge networks stabilize the open state. The structural and electrostatic properties of this open-state model are compatible with proton translocation and offer an explanation for the reversal of charge selectivity in neutral mutants of Asp(112). Furthermore, these structural properties are consistent with experimental accessibility data, providing a valuable basis for further structural and functional studies of hH(V)1. Each Arg residue in the S4 helix of hH(V)1 was replaced by His to test accessibility using Zn(2+) as a probe. The two outermost Arg residues in S4 were accessible to external solution, whereas the innermost one was accessible only to the internal solution. Both modeling and experimental data indicate that in the open state, Arg(211), the third Arg residue in the S4 helix in hH(V)1, remains accessible to the internal solution and is located near the charge transfer center, Phe(150).

The synthesis of exopolysaccharides as biofilm matrix components by pathogens is a crucial factor for chronic infections and antibiotic resistance. Many periplasmic proteins involved in polymer processing and secretion in Gram-negative synthase dependent exopolysaccharide biosynthetic systems have been individually characterized. The operons responsible for the production of PNAG, alginate, cellulose and the Pel polysaccharide each contain a gene that encodes an outer membrane associated tetratricopeptide repeat (TPR) domain containing protein. While the TPR domain has been shown to bind other periplasmic proteins, the functional consequences of these interactions for the polymer remain poorly understood. Herein, we show that the C-terminal TPR region of PgaA interacts with the de-N-acetylase domain of PgaB, and increases its deacetylase activity. Additionally, we found that when the two proteins form a complex, the glycoside hydrolase activity of PgaB is also increased. To better understand structure-function relationships we determined the crystal structure of a stable TPR module, which has a conserved groove formed by three repeat motifs. Tryptophan quenching, mass spectrometry analysis and molecular dynamics simulation studies suggest that the crystallized TPR module can bind PNAG/dPNAG via its electronegative groove on the concave surface, and potentially guide the polymer through the periplasm towards the porin for export. Our results suggest a scaffolding role for the TPR domain that combines PNAG/dPNAG translocation with the modulation of its chemical structure by PgaB.

EPDR1, a member of the ependymin-related protein family, is a relatively uncharacterized protein found in the lysosomes and secretomes of most vertebrates. Despite having roles in human disease and health, the molecular functions of EPDR1 remain unknown. Here, we present crystal structures of human EPDR1 and reveal that the protein adopts a fold previously seen only in bacterial proteins related to the LolA lipoprotein transporter. EPDR1 forms a homodimer with an overall shape resembling a half-shell with two non-overlapping hydrophobic grooves on the flat side of the hemisphere. EPDR1 can interact with membranes that contain negatively charged lipids, including BMP and GM1, and we suggest that EPDR1 may function as a lysosomal activator protein or a lipid transporter. A phylogenetic analysis reveals that the fold is more widely distributed than previously suspected, with representatives identified in all branches of cellular life.