PURA and GLUT1: Sweet partners for brain health

PURA, also known as Pur-alpha, is an evolutionarily conserved DNA/RNA-binding protein crucial for various cellular processes, including DNA replication, transcriptional regulation, and translational control. Comprising three PUR domains, it engages with nucleic acids and has a role in protein-protein interactions. The manifestation of PURA syndrome, arising from mutations in the PURA gene, presents neurologically with developmental delay, hypotonia, and seizures. In our prior work from 2018, we highlighted the unique case of a PURA patient displaying hypoglycorrhachia, suggesting a potential association with GLUT1 dysfunction in this syndrome. In this current study, we expand the patient cohort with PURA mutations exhibiting hypoglycorrhachia and aim to unravel the molecular basis of this phenomenon. We established an in vitro model in HeLa cells to modulate PURA expression and investigated GLUT1 function and expression. Our findings indicate that PURA levels directly impact glucose uptake through the functioning of GLUT1, without influencing significantly GLUT1 expression. Moreover, our study reveals evidence for a possible physical interaction between PURA and GLUT1, demonstrated by colocalization and co-immunoprecipitation of both proteins. Computational analyses, employing molecular dynamics, further corroborates these findings, demonstrating that PURA:GLUT1 interactions are plausible, and that the stability of the complex is altered when PURA is truncated and/or mutated. In conclusion, our results suggest that PURA plays a pivotal role in driving the function of GLUT1 for glucose uptake, potentially forming a regulatory complex. Additional investigations are warranted to elucidate the precise mechanisms governing this complex and its significance in ensuring proper GLUT1 function.

Myomerger Induces Membrane Hemifusion and Regulates Fusion Pore Expansion

Membrane fusion is a crucial mechanism in a wide variety of important events in cell biology from viral infection to exocytosis. However, despite many efforts and much progress, cell-cell fusion has remained elusive to our understanding. Along the life of the fusion pore, large conformational changes take place from the initial lipid bilayer bending, passing through the hemifusion intermediates, and ending with the formation of the first nascent fusion pore. In this sense, computer simulations are an ideal technique for describing such complex lipid remodeling at the molecular level. In this work, we studied the role played by the muscle-specific membrane protein Myomerger during the formation of the fusion pore. We have conducted μs length atomistic and coarse-grained molecular dynamics, together with free-energy calculations using ad hoc collective variables. Our results show that Myomerger favors the hemifusion diaphragm-stalk transition, reduces the nucleation-expansion energy difference, and promotes the formation of nonenlarging fusion pores.

Selection of a novel cell-internalizing RNA aptamer specific for CD22 antigen in B cell acute lymphoblastic leukemia

Despite improvements in B cell acute lymphoblastic leukemia (B-ALL) treatment, a significant number of patients experience relapse of the disease, resulting in poor prognosis and high mortality. One of the drawbacks of current B-ALL treatments is the high toxicity associated with the non-specificity of chemotherapeutic drugs. Targeted therapy is an appealing strategy to treat B-ALL to mitigate these toxic off-target effects. One such target is the B cell surface protein CD22. The restricted expression of CD22 on the B-cell lineage and its ligand-induced internalizing properties make it an attractive target in cases of B cell malignancies. To target B-ALL and the CD22 protein, we performed cell internalization SELEX (Systematic Evolution of Ligands by EXponential enrichment) followed by molecular docking to identify internalizing aptamers specific for B-ALL cells that bind the CD22 cell-surface receptor. We identified two RNA aptamers, B-ALL1 and B-ALL2, that target human malignant B cells, with B-ALL1 the first documented RNA aptamer interacting with the CD22 antigen. These B-ALL-specific aptamers represent an important first step toward developing novel targeted therapies for B cell malignancy treatments.

The SARS-CoV-2 mutation landscape is shaped before replication starts

Mutation landscapes and signatures have been thoroughly studied in SARS-CoV-2. Here, we analyse those patterns and link their changes to the viral replication tissue in the respiratory tract. Surprisingly, a substantial difference in those patterns is observed in samples from vaccinated patients. Hence, we propose a model to explain where those mutations could originate during the replication cycle.

α-Synuclein is required for sperm exocytosis at a post-fusion stage

The sperm acrosome is a large dense-core granule whose contents are secreted by regulated exocytosis at fertilization through the opening of numerous fusion pores between the acrosomal and plasma membranes. In other cells, the nascent pore generated when the membrane surrounding a secretory vesicle fuses with the plasma membrane may have different fates. In sperm, pore dilation leads to the vesiculation and release of these membranes, together with the granule contents. α-Synuclein is a small cytosolic protein claimed to exhibit different roles in exocytic pathways in neurons and neuroendocrine cells. Here, we scrutinized its function in human sperm. Western blot revealed the presence of α-synuclein and indirect immunofluorescence its localization to the acrosomal domain of human sperm. Despite its small size, the protein was retained following permeabilization of the plasma membrane with streptolysin O. α-Synuclein was required for acrosomal release, as demonstrated by the inability of an inducer to elicit exocytosis when permeabilized human sperm were loaded with inhibitory antibodies to human α-synuclein. The antibodies halted calcium-induced secretion when introduced after the acrosome docked to the cell membrane. Two functional assays, fluorescence and transmission electron microscopies revealed that the stabilization of open fusion pores was responsible for the secretion blockage. Interestingly, synaptobrevin was insensitive to neurotoxin cleavage at this point, an indication of its engagement in cis SNARE complexes. The very existence of such complexes during AE reflects a new paradigm. Recombinant α-synuclein rescued the inhibitory effects of the anti-α-synuclein antibodies and of a chimeric Rab3A-22A protein that also inhibits AE after fusion pore opening. We applied restrained molecular dynamics simulations to compare the energy cost of expanding a nascent fusion pore between two model membranes and found it higher in the absence than in the presence of α-synuclein. Hence, our results suggest that α-synuclein is essential for expanding fusion pores.

Intrinsic Disorder in α-Synuclein Regulates the Exocytotic Fusion Pore Transition

Today, it is widely accepted that intrinsic disorder is strongly related to the cell cycle, during mitosis, differentiation, and apoptosis. Of particular interest are hybrid proteins possessing both structured and unstructured domains that are critical in human health and disease, such as α-synuclein. In this work, we describe how α-synuclein interacts with the nascent fusion pore as it evolves toward expansion. We unveil the key role played by its intrinsically disordered region as a thermodynamic regulator of the nucleation-expansion energy barrier. By analyzing a truncated variant of α-synuclein that lacks the disordered region, we find that the landscape of protein interactions with PIP₂ and POPS lipids is highly altered, ultimately increasing the energy cost for the fusion pore to transit from nucleation to expansion. We conclude that the intrinsically disordered region in full-length α-synuclein recognizes and allocates pivotal protein:lipid interactions during membrane remodeling in the first stages of the fusion pore.

A lung-specific mutational signature enables inference of viral and bacterial respiratory niche

Exposure to different mutagens leaves distinct mutational patterns that can allow inference of pathogen replication niches. We therefore investigated whether SARS-CoV-2 mutational spectra might show lineage-specific differences, dependent on the dominant site(s) of replication and onwards transmission, and could therefore rapidly infer virulence of emergent variants of concern (VOCs). Through mutational spectrum analysis, we found a significant reduction in G>T mutations in the Omicron variant, which replicates in the upper respiratory tract (URT), compared to other lineages, which replicate in both the URT and lower respiratory tract (LRT). Mutational analysis of other viruses and bacteria indicates a robust, generalizable association of high G>T mutations with replication within the LRT. Monitoring G>T mutation rates over time, we found early separation of Omicron from Beta, Gamma and Delta, while mutational patterns in Alpha varied consistent with changes in transmission source as social restrictions were lifted. Mutational spectra may be a powerful tool to infer niches of established and emergent pathogens.

Enhanced Expansion and Reduced Kiss-and-Run Events in Fusion Pores Steered by Synaptotagmin-1 C2B Domains

The fusion pore controls the release of exocytotic vesicle contents through a precise orchestration of lipids from the fusing membranes and proteins. There is a major lipid reorganization during the different stages in life of the fusion pore (membrane fusion, nucleation, and expansion) that can be scrutinized thermodynamically. In this work, using umbrella sampling simulations we describe the expansion of the fusion pore. We have calculated free energy profiles to drive a nascent, just nucleated, fusion pore to its expanded configuration. We have quantified the effects on the free energy of one and two Synaptotagmin-1 C2B domains in the cytosolic space. We show that C2B domains cumulatively reduce the cost for expansion, favoring the system to evolve toward full fusion. Finally, by conducting thousands of unbiased molecular dynamics simulations, we show that C2B domains significantly decrease the probability of kiss-and-run events.

Synaptotagmin-1 C2B domains cooperatively stabilize the fusion stalk via a master-servant mechanism

Synaptotagmin-1 is a low-affinity Ca2+ sensor that triggers synchronous vesicle fusion. It contains two similar C2 domains (C2A and C2B) that cooperate in membrane binding, being the C2B domain the...

Artificial stabilization of the fusion pore by intra-organelle styrene-maleic acid copolymers

Styrene-maleic acid copolymers have become an advantageous detergent-free alternative for membrane protein isolation. Since their discovery, experimental membrane protein extraction and purification by keeping intact their lipid environment has become significantly easier. With the aim of identifying new applications of these interesting copolymers, their molecular binding and functioning mechanisms have recently become intense objects of study. In this work, we describe the use of styrene-maleic acid copolymers as an artificial tool to stabilize the fusion pore. We show that when these copolymers circumscribe the water channel that defines the fusion pore, they keep it from shrinking and closing. We describe how only intra-organelle copolymers have stabilizing capabilities while extra-organelle ones have negligible or even contrary effects on the fusion pore life-time.

Order–disorder skewness in alpha-synuclein: a key mechanism to recognize membrane curvature

Currently, membrane curvature is understood as an active mechanism to control cells spatial organization and activity.

A connection between reversible tyrosine phosphorylation and SNARE complex disassembly activity of N-ethylmaleimide-sensitive factor unveiled by the phosphomimetic mutant N-ethylmaleimide-sensitive factor-Y83E

N-ethylmaleimide-sensitive factor (NSF) disassembles fusion-incompetent cis soluble-NSF attachment protein receptor (SNARE) complexes making monomeric SNAREs available for subsequent trans pairing and fusion. In most cells the activity of NSF is constitutive, but in Jurkat cells and sperm it is repressed by tyrosine phosphorylation; the phosphomimetic mutant NSF–Y83E inhibits secretion in the former. The questions addressed here are if and how the NSF mutant influences the configuration of the SNARE complex. Our model is human sperm, where the initiation of exocytosis (acrosome reaction (AR)) de-represses the activity of NSF through protein tyrosine phosphatase 1B (PTP1B)-mediated dephosphorylation. We developed a fluorescence microscopy-based method to show that capacitation increased, and challenging with an AR inducer decreased, the number of cells with tyrosine-phosphorylated PTP1B substrates in the acrosomal domain. Results from bioinformatic and biochemical approaches using purified recombinant proteins revealed that NSF–Y83E bound PTP1B and thereupon inhibited its catalytic activity. Mutant NSF introduced into streptolysin O-permeabilized sperm impaired cis SNARE complex disassembly, blocking the AR; subsequent addition of PTP1B rescued exocytosis. We propose that NSF–Y83E prevents endogenous PTP1B from dephosphorylating sperm NSF, thus maintaining NSF’s activity in a repressed mode and the SNARE complex unable to dissociate. The contribution of this paper to the sperm biology field is the detection of PTP1B substrates, one of them likely being NSF, whose tyrosine phosphorylation status varies during capacitation and the AR. The contribution of this paper to the membrane traffic field is to have generated direct evidence that explains the dominant-negative role of the phosphomimetic mutant NSF–Y83E.

Grab recruitment by Rab27A-Rabphilin3a triggers Rab3A activation in human sperm exocytosis

Sperm must undergo the regulated exocytosis of its dense core granule (the acrosome reaction, AR) to fertilize the egg. We have previously described that Rabs3 and 27 are organized in a RabGEF cascade within the signaling pathway elicited by exocytosis stimuli in human sperm. Here, we report the identity and the role of two molecules that link these secretory Rabs in the RabGEF cascade: Rabphilin3a and GRAB. Like Rab3 and Rab27, GRAB and Rabphilin3a are present, localize to the acrosomal region and are required for calcium-triggered exocytosis in human sperm. Sequestration of either protein with specific antibodies introduced into streptolysin O-permeabilized sperm impairs the activation of Rab3 in the acrosomal region elicited by calcium, but not that of Rab27. Biochemical and functional assays indicate that Rabphilin3a behaves as a Rab27 effector during the AR and that GRAB exhibits GEF activity toward Rab3A. Recombinant, active Rab27A pulls down Rabphilin3a and GRAB from human sperm extracts. Conversely, immobilized Rabphilin3a recruits Rab27 and GRAB; the latter promotes Rab3A activation. The enzymatic activity of GRAB toward Rab3A was also suggested by in silico and in vitro assays with purified proteins. In summary, we describe here a signaling module where Rab27A-GTP interacts with Rabphilin3a, which in turn recruits a guanine nucleotide-exchange activity toward Rab3A. This is the first description of the interaction of Rabphilin3a with a GEF. Because the machinery that drives exocytosis is highly conserved, it is tempting to hypothesize that the RabGEF cascade unveiled here might be part of the molecular mechanisms that drive exocytosis in other secretory systems.

Transmembrane domain dimerization induces cholesterol rafts in curved lipid bilayers

Are the dimerization of transmembrane (TM) domains and the reorganization of the lipid bilayer two independent events?

Human Adipose-Derived Mesenchymal Stem/Stromal Cells Handling Protocols. Lipid Droplets and Proteins Double-Staining

Human Adipose-derived mesenchymal stem/stromal cells (hASCs) are of great interest because of their potential for therapeutic approaches. The method described here covers every single step necessary for hASCs isolation from subcutaneous abdominal adipose tissue, multicolor phenotyping by flow cytometry, and quantitative determination of adipogenic differentiation status by means of lipid droplets (LDs) accumulation, and Western blot analysis. Moreover, to simultaneously analyze both LDs accumulation and cellular proteins localization by fluorescence microscopy, we combined Oil Red O (ORO) staining with immunofluorescence detection. For LDs quantification we wrote a program for automatic ORO-stained digital image processing implemented in Octave, a freely available software package. Our method is based on the use of the traditional low cost neutral lipids dye ORO, which can be imaged both by bright-field and fluorescence microscopy. The utilization of ORO instead of other more expensive lipid-specific dyes, together with the fact that the whole method has been designed employing cost-effective culture reagents (standard culture medium and serum), makes it affordable for tight-budget research laboratories. These may be replaced, if necessary or desired, by defined xeno-free reagents for clinical research and applications.

Bending Lipid Bilayers: A Closed-Form Collective Variable for Effective Free-Energy Landscapes in Quantitative Biology

Curvature-related processes are of major importance during protein-membrane interactions. The illusive simplicity of membrane reshaping masks a complex molecular process crucial for a wide range of biological functions like fusion, endo- and exocytosis, cell division, cytokinesis, and autophagy. To date, no functional expression of a reaction coordinate capable of biasing molecular dynamics simulations to produce membrane curvature has been reported. This represents a major drawback given that the adequate identification of proper collective variables to enhance sampling is fundamental for restrained dynamics techniques. In this work, we present a closed-form equation of a collective variable that induces bending in lipid bilayers in a controlled manner, allowing for straightforward calculation of free energy landscapes of important curvature-related events, using standard methods such as umbrella sampling and metadynamics. As a direct application of the collective variable, we calculate the bending free energies of a ternary lipid bilayer in the presence and the absence of a Bin/Amphiphysin/Rvs domain with an N-terminal amphipathic helix (N-BAR), a well-known peripheral membrane protein that induces curvature.

On the role of residue phosphorylation in 14-3-3 partners: AANAT as a case study

Twenty years ago, a novel concept in protein structural biology was discovered: the intrinsically disordered regions (IDRs). These regions remain largely unstructured under native conditions and the more are studied, more properties are attributed to them. Possibly, one of the most important is their ability to conform a new type of protein-protein interaction. Besides the classical domain-to-domain interactions, IDRs follow a 'fly-casting' model including 'induced folding'. Unfortunately, it is only possible to experimentally explore initial and final states. However, the complete movie of conformational changes of protein regions and their characterization can be addressed by in silico experiments. Here, we simulate the binding of two proteins to describe how the phosphorylation of a single residue modulates the entire process. 14-3-3 protein family is considered a master regulator of phosphorylated proteins and from a modern point-of-view, protein phosphorylation is a three component system, with writers (kinases), erasers (phosphatases) and readers. This later biological role is attributed to the 14-3-3 protein family. Our molecular dynamics results show that phosphorylation of the key residue Thr31 in a partner of 14-3-3, the aralkylamine N-acetyltransferase, releases the fly-casting mechanism during binding. On the other hand, the non-phosphorylation of the same residue traps the proteins, systematically and repeatedly driving the simulations into wrong protein-protein conformations.

Molecular-level insight into the binding of arginine to a zwitterionic Langmuir monolayer

Arginine molecules bind to a DPPC monolayer, altering the interfacial electrostatic potential and the lateral mobility of the lipids, while having little effect on the compression isotherm of the monolayer.

Study on the interaction of artificial and natural food colorants with human serum albumin: A computational point of view

Due to the high amount of artificial food colorants present in infants' diets, their adverse effects have been of major concern among the literature. Artificial food colorants have been suggested to affect children's behavior, being hyperactivity the most common disorder. In this study we compare binding affinities of a group of artificial colorants (sunset yellow, quinoline yellow, carmoisine, allura red and tartrazine) and their natural industrial equivalents (carminic acid, curcumin, peonidin-3-glucoside, cyanidin-3-glucoside) to human serum albumin (HSA) by a docking approach and further refinement through atomistic molecular dynamics simulations. Due to the protein-ligand conformational interface complexity, we used collective variable driven molecular dynamics to refine docking predictions and to score them according to a hydrogen-bond criterion. With this protocol, we were able to rank ligand affinities to HSA and to compare between the studied natural and artificial food additives. Our results show that the five artificial colorants studied bind better to HSA than their equivalent natural options, in terms of their H-bonding network, supporting the hypothesis of their potential risk to human health.

Interaction of glycine, lysine, proline and histidine with dipalmitoylphosphatidylcholine lipid bilayers: a theoretical and experimental study

The interaction of unblocked glycine, lysine, proline and histidine with a DPPC bilayer was assessed using extensive atomistic molecular dynamics simulations.

H-bond refinement for electron transfer membrane-bound protein-protein complexes: cytochrome c oxidase and cytochrome c552

In this study we propose a protocol to evaluate membrane-bound cytochrome c oxidase-cytochrome c552 docking candidates. An initial rigid docking algorithm generates docking poses of the cytochrome c oxidase-cytochrome c552, candidates are then aggregated into a 512-DPPC membrane model and solvated in explicit solvent. Molecular dynamic simulations are performed to induce conformational changes to membrane-bound protein complexes. Lastly each protein-protein complex is optimized in terms of its hydrogen bond network, evaluated energetically and ranked. The protocol is directly applicable to other membrane-protein complexes, such as protein-ligand systems.

H-bond network optimization in protein-protein complexes: are all-atom force field scores enough?

Structural prediction of protein-protein complexes given the structures of the two interacting compounds in their unbound state is a key problem in biophysics. In addition to the problem of sampling of near-native orientations, one of the modeling main difficulties is to discriminate true from false positives. Here, we present a hierarchical protocol for docking refinement able to discriminate near native poses from a group of docking candidates. The main idea is to combine an efficient sampling of the full system hydrogen bond network and side chains, together with an all-atom force field and a surface generalized born implicit solvent. We tested our method on a set of twenty two complexes containing a near-native solution within the top 100 docking poses, obtaining a near native solution as the top pose in 70% of the cases. We show that all atom force fields optimized H-bond networks do improve significantly state of the art scoring functions.