PeSTo-Carbs: Geometric Deep Learning for Prediction of Protein-Carbohydrate Binding Interfaces

The Protein Structure Transformer (PeSTo), a geometric transformer, has exhibited exceptional performance in predicting protein-protein binding interfaces and distinguishing interfaces with nucleic acids, lipids, small molecules, and ions. In this study, we introduce PeSTo-Carbs, an extension of PeSTo specifically engineered to predict protein-carbohydrate binding interfaces. We evaluate the performance of this approach using independent test sets and compare them with those of previous methods. Furthermore, we highlight the model's capability to specialize in predicting interfaces involving cyclodextrins, a biologically and pharmaceutically significant class of carbohydrates. Our method consistently achieves remarkable accuracy despite the scarcity of available structural data for cyclodextrins.

A humoral stress response protects Drosophila tissues from antimicrobial peptides

7An efficient immune system must provide protection against a broad range of pathogens without causing excessive collateral tissue damage. While immune effectors have been well characterized, we know less about the resilience mechanisms protecting the host from its own immune response. Antimicrobial peptides (AMPs) are small, cationic peptides that contribute to innate defenses by targeting negatively charged membranes of microbes. While protective against pathogens, AMPs can be cytotoxic to host cells. Here, we reveal that a family of stress-induced proteins, the Turandots, protect the Drosophila respiratory system from AMPs, increasing resilience to stress. Flies lacking Turandot genes are susceptible to environmental stresses due to AMP-induced tracheal apoptosis. Turandot proteins bind to host cell membranes and mask negatively charged phospholipids, protecting them from cationic pore-forming AMPs. Collectively, these data demonstrate that Turandot stress proteins mitigate AMP cytotoxicity to host tissues and therefore improve their efficacy.

Deep Learning-Assisted Single-Molecule Detection of Protein Post-translational Modifications with a Biological Nanopore

Protein post-translational modifications (PTMs) play a crucial role in countless biological processes, profoundly modulating protein properties on both spatial and temporal scales. Protein PTMs have also emerged as reliable biomarkers for several diseases. However, only a handful of techniques are available to accurately measure their levels, capture their complexity at a single molecule level, and characterize their multifaceted roles in health and disease. Nanopore sensing provides high sensitivity for the detection of low-abundance proteins, holding the potential to impact single-molecule proteomics and PTM detection, in particular. Here, we demonstrate the ability of a biological nanopore, the pore-forming toxin aerolysin, to detect and distinguish α-synuclein-derived peptides bearing single or multiple PTMs, namely, phosphorylation, nitration, and oxidation occurring at different positions and in various combinations. The characteristic current signatures of the α-synuclein peptide and its PTM variants could be confidently identified by using a deep learning model for signal processing. We further demonstrate that this framework can quantify α-synuclein peptides at picomolar concentrations and detect the C-terminal peptides generated by digestion of full-length α-synuclein. Collectively, our work highlights the advantage of using nanopores as a tool for simultaneous detection of multiple PTMs and facilitates their use in biomarker discovery and diagnostics.

Bicc1 ribonucleoprotein complexes specifying organ laterality are licensed by ANKS6-induced structural remodeling of associated ANKS3

Organ laterality of vertebrates is specified by accelerated asymmetric decay of Dand5 mRNA mediated by Bicaudal-C1 (Bicc1) on the left side, but whether binding of this or any other mRNA to Bicc1 can be regulated is unknown. Here, we found that a CRISPR-engineered truncation in ankyrin and sterile alpha motif (SAM)-containing 3 (ANKS3) leads to symmetric mRNA decay mediated by the Bicc1-interacting Dand5 3' UTR. AlphaFold structure predictions of protein complexes and their biochemical validation by in vitro reconstitution reveal a novel interaction of the C-terminal coiled coil domain of ANKS3 with Bicc1 that inhibits binding of target mRNAs, depending on the conformation of ANKS3 and its regulation by ANKS6. The dual regulation of RNA binding by mutually opposing structured protein domains in this multivalent protein network emerges as a novel mechanism linking associated laterality defects and possibly other ciliopathies to perturbed dynamics in Bicc1 ribonucleoparticle (RNP) formation.

Structural mechanism of mitochondrial membrane remodelling by human OPA1

Distinct morphologies of the mitochondrial network support divergent metabolic and regulatory processes that determine cell function and fate1-3. The mechanochemical GTPase optic atrophy 1 (OPA1) influences the architecture of cristae and catalyses the fusion of the mitochondrial inner membrane4,5. Despite its fundamental importance, the molecular mechanisms by which OPA1 modulates mitochondrial morphology are unclear. Here, using a combination of cellular and structural analyses, we illuminate the molecular mechanisms that are key to OPA1-dependent membrane remodelling and fusion. Human OPA1 embeds itself into cardiolipin-containing membranes through a lipid-binding paddle domain. A conserved loop within the paddle domain inserts deeply into the bilayer, further stabilizing the interactions with cardiolipin-enriched membranes. OPA1 dimerization through the paddle domain promotes the helical assembly of a flexible OPA1 lattice on the membrane, which drives mitochondrial fusion in cells. Moreover, the membrane-bending OPA1 oligomer undergoes conformational changes that pull the membrane-inserting loop out of the outer leaflet and contribute to the mechanics of membrane remodelling. Our findings provide a structural framework for understanding how human OPA1 shapes mitochondrial morphology and show us how human disease mutations compromise OPA1 functions.

PeSTo: parameter-free geometric deep learning for accurate prediction of protein binding interfaces

Proteins are essential molecular building blocks of life, responsible for most biological functions as a result of their specific molecular interactions. However, predicting their  binding  interfaces remains a challenge. In this study, we present a geometric transformer that acts directly on atomic coordinates labeled only with element names. The resulting model-the Protein Structure Transformer, PeSTo-surpasses the current state of the art in predicting protein-protein interfaces and can also predict and differentiate between interfaces involving nucleic acids, lipids, ions, and small molecules with high confidence. Its low computational cost enables processing high volumes of structural data, such as molecular dynamics ensembles allowing for the discovery of interfaces that remain otherwise inconspicuous in static experimentally solved structures. Moreover, the growing foldome provided by de novo structural predictions can be easily analyzed, providing new opportunities to uncover unexplored biology.

Structure and functionality of a multimeric human COQ7:COQ9 complex

Coenzyme Q (CoQ) is a redox-active lipid essential for core metabolic pathways and antioxidant defense. CoQ is synthesized upon the mitochondrial inner membrane by an ill-defined "complex Q" metabolon. Here, we present structure-function analyses of a lipid-, substrate-, and NADH-bound complex comprising two complex Q subunits: the hydroxylase COQ7 and the lipid-binding protein COQ9. We reveal that COQ7 adopts a ferritin-like fold with a hydrophobic channel whose substrate-binding capacity is enhanced by COQ9. Using molecular dynamics, we further show that two COQ7:COQ9 heterodimers form a curved tetramer that deforms the membrane, potentially opening a pathway for the CoQ intermediates to translocate from the bilayer to the proteins' lipid-binding sites. Two such tetramers assemble into a soluble octamer with a pseudo-bilayer of lipids captured within. Together, these observations indicate that COQ7 and COQ9 cooperate to access hydrophobic precursors within the membrane and coordinate subsequent synthesis steps toward producing CoQ.

Biological nanopores for single-molecule sensing

Evolution has found countless ways to transport material across cells and cellular compartments separated by membranes. Protein assemblies are the cornerstone for the formation of channels and pores that enable this regulated passage of molecules in and out of cells, contributing to maintaining most of the fundamental processes that sustain living organisms. As in several other occasions, we have borrowed from the natural properties of these biological systems to push technology forward and have been able to hijack these nano-scale proteinaceous pores to learn about the physical and chemical features of molecules passing through them. Today, a large repertoire of biological pores is exploited as molecular sensors for characterizing biomolecules that are relevant for the advancement of life sciences and application to medicine. Although the technology has quickly matured to enable nucleic acid sensing with transformative implications for genomics, biological pores stand as some of the most promising candidates to drive the next developments in single-molecule proteomics.

Visualization, Interactive Handling and Simulation of Molecules in Commodity Augmented Reality in Web Browsers Using moleculARweb’s Virtual Modeling Kits

moleculARweb (https://molecularweb.epfl.ch) began as a website for education and outreach in chemistry and structural biology through augmented reality (AR) content that runs in the web browsers of regular devices like smartphones, tablets, and computers. Here we present two evolutions of moleculARweb’s Virtual Modeling Kits (VMK), tools where users can build and view molecules, and explore their mechanics, in 3D AR by handling the molecules in full 3D with custom-printed cube markers (VMK 2.0) or by moving around a simulated scene with mouse or touch gestures (VMK 3.0). Upon simulation the molecules experience visually realistic torsions, clashes, and hydrogen-bonding interactions that the user can manually switch on and off to explore their effects. Moreover, by manually tuning a fictitious temperature the users can accelerate conformational transitions or 'freeze' specific conformations for careful inspection in 3D. Even some phase transitions and separations can be simulated. We here showcase these and other features of the new VMKs connecting them to possible specific applications to teaching and self-learning of concepts from general, organic, biological and physical chemistry; and in assisting with small tasks in molecular modelling for research. Last, in a short discussion section we overview what future developments are needed for the 'dream tool' for the future of chemistry education and work.

Prolyl endopeptidase-like is a (thio)esterase involved in mitochondrial respiratory chain function

Deficiency of the serine hydrolase prolyl endopeptidase-like (PREPL) causes a recessive metabolic disorder characterized by neonatal hypotonia, feeding difficulties, and growth hormone deficiency. The pathophysiology of PREPL deficiency and the physiological substrates of PREPL remain largely unknown. In this study, we connect PREPL with mitochondrial gene expression and oxidative phosphorylation by analyzing its protein interactors. We demonstrate that the long PREPLL isoform localizes to mitochondria, whereas PREPLS remains cytosolic. Prepl KO mice showed reduced mitochondrial complex activities and disrupted mitochondrial gene expression. Furthermore, mitochondrial ultrastructure was abnormal in a PREPL-deficient patient and Prepl KO mice. In addition, we reveal that PREPL has (thio)esterase activity and inhibition of PREPL by Palmostatin M suggests a depalmitoylating function. We subsequently determined the crystal structure of PREPL, thereby providing insight into the mechanism of action. Taken together, PREPL is a (thio)esterase rather than a peptidase and PREPLL is involved in mitochondrial homeostasis.

BET1 variants establish impaired vesicular transport as a cause for muscular dystrophy with epilepsy

BET1 is required, together with its SNARE complex partners GOSR2, SEC22b, and Syntaxin-5 for fusion of endoplasmic reticulum-derived vesicles with the ER-Golgi intermediate compartment (ERGIC) and the cis-Golgi. Here, we report three individuals, from two families, with severe congenital muscular dystrophy (CMD) and biallelic variants in BET1 (P1 p.(Asp68His)/p.(Ala45Valfs*2); P2 and P3 homozygous p.(Ile51Ser)). Due to aberrant splicing and frameshifting, the variants in P1 result in low BET1 protein levels and impaired ER-to-Golgi transport. Since in silico modeling suggested that p.(Ile51Ser) interferes with binding to interaction partners other than SNARE complex subunits, we set off and identified novel BET1 interaction partners with low affinity for p.(Ile51Ser) BET1 protein compared to wild-type, among them ERGIC-53. The BET1/ERGIC-53 interaction was validated by endogenous co-immunoprecipitation with both proteins colocalizing to the ERGIC compartment. Mislocalization of ERGIC-53 was observed in P1 and P2's derived fibroblasts; while in the p.(Ile51Ser) P2 fibroblasts specifically, mutant BET1 was also mislocalized along with ERGIC-53. Thus, we establish BET1 as a novel CMD/epilepsy gene and confirm the emerging role of ER/Golgi SNAREs in CMD.

S-acylation controls SARS-CoV-2 membrane lipid organization and enhances infectivity

SARS-CoV-2 virions are surrounded by a lipid bilayer that contains membrane proteins such as spike, responsible for target-cell binding and virus fusion. We found that during SARS-CoV-2 infection, spike becomes lipid modified, through the sequential action of the S-acyltransferases ZDHHC20 and 9. Particularly striking is the rapid acylation of spike on 10 cytosolic cysteines within the ER and Golgi. Using a combination of computational, lipidomics, and biochemical approaches, we show that this massive lipidation controls spike biogenesis and degradation, and drives the formation of localized ordered cholesterol and sphingolipid-rich lipid nanodomains in the early Golgi, where viral budding occurs. Finally, S-acylation of spike allows the formation of viruses with enhanced fusion capacity. Our study points toward S-acylating enzymes and lipid biosynthesis enzymes as novel therapeutic anti-viral targets.

Investigating Crosstalk Among PTMs Provides Novel Insight Into the Structural Basis Underlying the Differential Effects of Nt17 PTMs on Mutant Httex1 Aggregation

Post-translational modifications (PTMs) within the first 17 amino acids (Nt17) of the Huntingtin protein (Htt) have been shown to inhibit the aggregation and attenuate the toxicity of mutant Htt proteins in vitro and in various models of Huntington’s disease. Here, we expand on these studies by investigating the effect of methionine eight oxidation (oxM8) and its crosstalk with lysine 6 acetylation (AcK6) or threonine 3 phosphorylation (pT3) on the aggregation of mutant Httex1 (mHttex1). We show that M8 oxidation delays but does not inhibit the aggregation and has no effect on the final morphologies of mHttex1aggregates. The presence of both oxM8 and AcK6 resulted in dramatic inhibition of Httex1 fibrillization. Circular dichroism spectroscopy and molecular dynamics simulation studies show that PTMs that lower the mHttex1 aggregation rate (oxM8, AcK6/oxM8, pT3, pT3/oxM8, and pS13) result in increased population of a short N-terminal helix (first eight residues) in Nt17 or decreased abundance of other helical forms, including long helix and short C-terminal helix. PTMs that did not alter the aggregation rate (AcK6) of mHttex1 exhibit a similar distribution of helical conformation as the unmodified peptides. These results show that the relative abundance of N- vs. C-terminal helical conformations and long helices, rather than the overall helicity of Nt17, better explains the effect of different Nt17 PTMs on mHttex1; thus, explaining the lack of correlation between the effect of PTMs on the overall helicity of Nt17 and mHttex1 aggregation in vitro. Taken together, our results provide novel structural insight into the differential effects of single PTMs and crosstalk between different PTMs in regulating mHttex1 aggregation.

CLoNe: automated clustering based on local density neighborhoods for application to biomolecular structural ensembles

Motivation

Proteins are intrinsically dynamic entities. Flexibility sampling methods, such as molecular dynamics or those arising from integrative modeling strategies, are now commonplace and enable the study of molecular conformational landscapes in many contexts. Resulting structural ensembles increase in size as technological and algorithmic advancements take place, making their analysis increasingly demanding. In this regard, cluster analysis remains a go-to approach for their classification. However, many state-of-the-art algorithms are restricted to specific cluster properties. Combined with tedious parameter fine-tuning, cluster analysis of protein structural ensembles suffers from the lack of a generally applicable and easy to use clustering scheme.

Results

We present CLoNe, an original Python-based clustering scheme that builds on the Density Peaks algorithm of Rodriguez and Laio. CLoNe relies on a probabilistic analysis of local density distributions derived from nearest neighbors to find relevant clusters regardless of cluster shape, size, distribution and amount. We show its capabilities on many toy datasets with properties otherwise dividing state-of-the-art approaches and improves on the original algorithm in key aspects. Applied to structural ensembles, CLoNe was able to extract meaningful conformations from membrane binding events and ligand-binding pocket opening as well as identify dominant dimerization motifs or inter-domain organization. CLoNe additionally saves clusters as individual trajectories for further analysis and provides scripts for automated use with molecular visualization software.

Availability and implementation

www.epfl.ch/labs/lbm/resources, github.com/LBM-EPFL/CLoNe.

Supplementary information

Supplementary data are available at Bioinformatics online.

Palmitoylated acyl protein thioesterase APT2 deforms membranes to extract substrate acyl chains

Many biochemical reactions require controlled recruitment of proteins to membranes. This is largely regulated by posttranslational modifications. A frequent one is S-acylation, which consists of the addition of acyl chains and can be reversed by poorly understood acyl protein thioesterases (APTs). Using a panel of computational and experimental approaches, we dissect the mode of action of the major cellular thioesterase APT2 (LYPLA2). We show that soluble APT2 is vulnerable to proteasomal degradation, from which membrane binding protects it. Interaction with membranes requires three consecutive steps: electrostatic attraction, insertion of a hydrophobic loop and S-acylation by the palmitoyltransferases ZDHHC3 or ZDHHC7. Once bound, APT2 is predicted to deform the lipid bilayer to extract the acyl chain bound to its substrate and capture it in a hydrophobic pocket to allow hydrolysis. This molecular understanding of APT2 paves the way to understand the dynamics of APT2-mediated deacylation of substrates throughout the endomembrane system.

Extensive tissue-specific expression variation and novel regulators underlying circadian behavior

Natural genetic variation affects circadian rhythms across the evolutionary tree, but the underlying molecular mechanisms are poorly understood. We investigated population-level, molecular circadian clock variation by generating >700 tissue-specific transcriptomes of Drosophila melanogaster (w1118 ) and 141 Drosophila Genetic Reference Panel (DGRP) lines. This comprehensive circadian gene expression atlas contains >1700 cycling genes including previously unknown central circadian clock components and tissue-specific regulators. Furthermore, >30% of DGRP lines exhibited aberrant circadian gene expression, revealing abundant genetic variation-mediated, intertissue circadian expression desynchrony. Genetic analysis of one line with the strongest deviating circadian expression uncovered a novel cry mutation that, as shown by protein structural modeling and brain immunohistochemistry, disrupts the light-driven flavin adenine dinucleotide cofactor photoreduction, providing in vivo support for the importance of this conserved photoentrainment mechanism. Together, our study revealed pervasive tissue-specific circadian expression variation with genetic variants acting upon tissue-specific regulatory networks to generate local gene expression oscillations.

The Human RNA Helicase DDX21 Presents a Dimerization Interface Necessary for Helicase Activity

Members of the DEAD-box helicase family are involved in all fundamental processes of RNA metabolism, and as such, their malfunction is associated with various diseases. Currently, whether and how oligomerization impacts their biochemical and biological functions is not well understood. In this work, we show that DDX21, a human DEAD-box helicase with RNA G-quadruplex resolving activity, is dimeric and that its oligomerization state influences its helicase activity. Solution small-angle X-ray scattering (SAXS) analysis uncovers a flexible multi-domain protein with a central dimerization domain. While the Arg/Gly rich C termini, rather than dimerization, are key to maintaining high affinity for RNA substrates, in vitro helicase assays indicate that an intact dimer is essential for both DDX21 ATP-dependent double-stranded RNA unwinding and ATP-independent G-quadruplex remodeling activities. Our results suggest that oligomerization plays a key role in regulating RNA DEAD-box helicase activity.

•

The human RNA DEAD-box helicase DDX21 is dimeric

•

DDX21 dimerization is mediated by a hydrophobic central core domain

•

SAXS-based modeling reveals that DDX21 is intrinsically flexible

•

Dimerization and C-terminal domains mediate G-quadruplex and dsRNA unwinding

Aerolysin nanopores decode digital information stored in tailored macromolecular analytes

Digital data storage is a growing need for our society and finding alternative solutions than those based on silicon or magnetic tapes is a challenge in the era of "big data." The recent development of polymers that can store information at the molecular level has opened up new opportunities for ultrahigh density data storage, long-term archival, anticounterfeiting systems, and molecular cryptography. However, synthetic informational polymers are so far only deciphered by tandem mass spectrometry. In comparison, nanopore technology can be faster, cheaper, nondestructive and provide detection at the single-molecule level; moreover, it can be massively parallelized and miniaturized in portable devices. Here, we demonstrate the ability of engineered aerolysin nanopores to accurately read, with single-bit resolution, the digital information encoded in tailored informational polymers alone and in mixed samples, without compromising information density. These findings open promising possibilities to develop writing-reading technologies to process digital data using a biological-inspired platform.

Structural basis for recognition of bacterial cell wall teichoic acid by pseudo-symmetric SH3b-like repeats of a viral peptidoglycan hydrolase

Endolysins are bacteriophage-encoded peptidoglycan hydrolases targeting the cell wall of host bacteria via their cell wall-binding domains (CBDs). The molecular basis for selective recognition of surface carbohydrate ligands by CBDs remains elusive. Here, we describe, in atomic detail, the interaction between the Listeria phage endolysin domain CBD500 and its cell wall teichoic acid (WTA) ligands. We show that 3′O-acetylated GlcNAc residues integrated into the WTA polymer chain are the key epitope recognized by a CBD binding cavity located at the interface of tandem copies of beta-barrel, pseudo-symmetric SH3b-like repeats. This cavity consists of multiple aromatic residues making extensive interactions with two GlcNAc acetyl groups via hydrogen bonds and van der Waals contacts, while permitting the docking of the diastereomorphic ligands. Our multidisciplinary approach tackled an extremely challenging protein–glycopolymer complex and delineated a previously unknown recognition mechanism by which a phage endolysin specifically recognizes and targets WTA, suggesting an adaptable model for regulation of endolysin specificity.

Combining genetic, biochemical and computational approaches, we elucidated the molecular mechanisms underlying the recognition of Listeria wall teichoic acid by bacteriophage-encoded SH3b repeats.

Molecular Bases of the Membrane Association Mechanism Potentiating Antibiotic Resistance by New Delhi Metallo-β-lactamase 1

Resistance to last-resort carbapenem antibiotics is an increasing threat to human health, as it critically limits therapeutic options. Metallo-β-lactamases (MBLs) are the largest family of carbapenemases, enzymes that inactivate these drugs. Among MBLs, New Delhi metallo-β-lactamase 1 (NDM-1) has experienced the fastest and largest worldwide dissemination. This success has been attributed to the fact that NDM-1 is a lipidated protein anchored to the outer membrane of bacteria, while all other MBLs are soluble periplasmic enzymes. By means of a combined experimental and computational approach, we show that NDM-1 interacts with the surface of bacterial membranes in a stable, defined conformation, in which the active site is not occluded by the bilayer. Although the lipidation is required for a long-lasting interaction, the globular domain of NDM-1 is tuned to interact specifically with the outer bacterial membrane. In contrast, this affinity is not observed for VIM-2, a natively soluble MBL. Finally, we identify key residues involved in the membrane interaction with NDM-1, which constitute potential targets for developing therapeutic strategies able to combat resistance granted by this enzyme.

The Polyglutamine Expansion at the N-Terminal of Huntingtin Protein Modulates the Dynamic Configuration and Phosphorylation of the C-Terminal HEAT Domain

The polyQ expansion in huntingtin protein (HTT) is the prime cause of Huntington's disease (HD). The recent cryoelectron microscopy (cryo-EM) structure of HTT-HAP40 complex provided the structural information on its HEAT-repeat domains. Here, we present analyses of the impact of polyQ length on the structure and function of HTT via an integrative structural and biochemical approach. The cryo-EM analysis of normal (Q23) and disease (Q78) type HTTs shows that the structures of apo HTTs significantly differ from the structure of HTT in a HAP40 complex and that the polyQ expansion induces global structural changes in the relative movements among the HTT domains. In addition, we show that the polyQ expansion alters the phosphorylation pattern across HTT and that Ser2116 phosphorylation in turn affects the global structure and function of HTT. These results provide a molecular basis for the effect of the polyQ segment on HTT structure and activity, which may be important for HTT pathology.

State-of-the-art web services for de novo protein structure prediction

Residue coevolution estimations coupled to machine learning methods are revolutionizing the ability of protein structure prediction approaches to model proteins that lack clear homologous templates in the Protein Data Bank (PDB). This has been patent in the last round of the Critical Assessment of Structure Prediction (CASP), which presented several very good models for the hardest targets. Unfortunately, literature reporting on these advances often lacks digests tailored to lay end users; moreover, some of the top-ranking predictors do not provide webservers that can be used by nonexperts. How can then end users benefit from these advances and correctly interpret the predicted models? Here we review the web resources that biologists can use today to take advantage of these state-of-the-art methods in their research, including not only the best de novo modeling servers but also datasets of models precomputed by experts for structurally uncharacterized protein families. We highlight their features, advantages and pitfalls for predicting structures of proteins without clear templates. We present a broad number of applications that span from driving forward biochemical investigations that lack experimental structures to actually assisting experimental structure determination in X-ray diffraction, cryo-EM and other forms of integrative modeling. We also discuss issues that must be considered by users yet still require further developments, such as global and residue-wise model quality estimates and sources of residue coevolution other than monomeric tertiary structure.

Ligand Binding to the Collagen VI Receptor Triggers a Talin-to-RhoA Switch that Regulates Receptor Endocytosis

Capillary morphogenesis gene 2 (CMG2/ANTXR2) is a cell surface receptor for both collagen VI and anthrax toxin. Biallelic loss-of-function mutations in CMG2 lead to a severe condition, hyaline fibromatosis syndrome (HFS). We have here dissected a network of dynamic interactions between CMG2 and various actin interactors and regulators, describing a different behavior from other extracellular matrix receptors. CMG2 binds talin, and thereby the actin cytoskeleton, only in its ligand-free state. Extracellular ligand binding leads to src-dependent talin release and recruitment of the actin cytoskeleton regulator RhoA and its effectors. These sequential interactions of CMG2 are necessary for the control of oriented cell division during fish development. Finally, we demonstrate that effective switching between talin and RhoA binding is required for the intracellular degradation of collagen VI in human fibroblasts, which explains why HFS mutations in the cytoskeleton-binding domain lead to dysregulation of extracellular matrix homeostasis.

Will Cryo-Electron Microscopy Shift the Current Paradigm in Protein Structure Prediction?

Protein dynamics is undoubtedly a pervasive ingredient in all biological functions. However, structural biology has been strongly driven by a static-centered view of protein architecture. We argue that the recent advances of cryo-electron microscopy (EM) have the potential to more broadly explore the conformational landscapes of protein complexes and therefore will enhance our ability to predict the diverse conformations of tertiary and quaternary protein structures that are functionally relevant in physiological conditions.

Active Site-Induced Evolutionary Constraints Follow Fold Polarity Principles in Soluble Globular Enzymes

A recent analysis of evolutionary rates in >500 globular soluble enzymes revealed pervasive conservation gradients toward catalytic residues. By looking at amino acid preference profiles rather than evolutionary rates in the same data set, we quantified the effects of active sites on site-specific constraints for physicochemical traits. We found that conservation gradients respond to constraints for polarity, hydrophobicity, flexibility, rigidity and structure in ways consistent with fold polarity principles; while sites far from active sites seem to experience no physicochemical constraint, rather being highly variable and favoring amino acids of low metabolic cost. Globally, our results highlight that amino acid variation contains finer information about protein structure than usually regarded in evolutionary models, and that this information is retrievable automatically with simple fits. We propose that analyses of the kind presented here incorporated into models of protein evolution should allow for better description of the physical chemistry that underlies molecular evolution.

An experiment-informed signal transduction model for the role of the Staphylococcus aureus MecR1 protein in β-lactam resistance

The treatment of hospital- and community-associated infections by methicillin-resistant Staphylococcus aureus (MRSA) is a perpetual challenge. This Gram-positive bacterium is resistant specifically to β-lactam antibiotics, and generally to many other antibacterial agents. Its resistance mechanisms to β-lactam antibiotics are activated only when the bacterium encounters a β-lactam. This activation is regulated by the transmembrane sensor/signal transducer proteins BlaR1 and MecR1. Neither the transmembrane/metalloprotease domain, nor the complete MecR1 and BlaR1 proteins, are isolatable for mechanistic study. Here we propose a model for full-length MecR1 based on homology modeling, residue coevolution data, a new extensive experimental mapping of transmembrane topology, partial structures, molecular simulations, and available NMR data. Our model defines the metalloprotease domain as a hydrophilic transmembrane chamber effectively sealed by the apo-sensor domain. It proposes that the amphipathic helices inserted into the gluzincin domain constitute the route for transmission of the β-lactam-binding event in the extracellular sensor domain, to the intracellular and membrane-embedded zinc-containing active site. From here, we discuss possible routes for subsequent activation of proteolytic action. This study provides the first coherent model of the structure of MecR1, opening routes for future functional investigations on how β-lactam binding culminates in the proteolytic degradation of MecI.

Modelling structures in cryo-EM maps

Recent advances in structure determination of sub-cellular structures using cryo-electron microscopy and tomography have enabled us to understand their architecture in a more detailed manner and gain insight into their function. The choice of approach to use for atomic model building, fitting, refinement and validation in the 3D map resulting from these experiments depends primarily on the resolution of the map and the prior information on the corresponding model. Here, we survey some of such methods and approaches and highlight their uses in specific recent examples.

KAP1 is an antiparallel dimer with a functional asymmetry

This study reveals the architecture of human KAP1 by integrating molecular modeling with small-angle X-ray scattering and single-molecule experiments. KAP1 dimers feature a structural asymmetry at the C-terminal domains that has functional implications for recruitment of HP1.

KAP1 (KRAB domain–associated protein 1) plays a fundamental role in regulating gene expression in mammalian cells by recruiting different transcription factors and altering the chromatin state. In doing so, KAP1 acts both as a platform for macromolecular interactions and as an E3 small ubiquitin modifier ligase. This work sheds light on the overall organization of the full-length protein combining solution scattering data, integrative modeling, and single-molecule experiments. We show that KAP1 is an elongated antiparallel dimer with an asymmetry at the C-terminal domains. This conformation is consistent with the finding that the Really Interesting New Gene (RING) domain contributes to KAP1 auto-SUMOylation. Importantly, this intrinsic asymmetry has key functional implications for the KAP1 network of interactions, as the heterochromatin protein 1 (HP1) occupies only one of the two putative HP1 binding sites on the KAP1 dimer, resulting in an unexpected stoichiometry, even in the context of chromatin fibers.

A further leap of improvement in tertiary structure prediction in CASP13 prompts new routes for future assessments

We present our assessment of tertiary structure predictions for hard targets in Critical Assessment of Structure Prediction round 13 (CASP13). The analysis includes (a) assignment and discussion of best models through scores-aided visual inspection of models for each evaluation unit (EU); (b) ranking of predictors resulting from this evaluation and from global scores; and (c) evaluation of progress, state of the art, and current limitations of protein structure prediction. We witness a sizable improvement in tertiary structure prediction building on the progress observed from CASP11 to CASP12, with (a) top models reaching backbone RMSD <3 å for several EUs of size <150 residues, contributed by many groups; (b) at least one model that roughly captures global topology for all EUs, probably unprecedented in this track of CASP; and (c) even quite good models for full, unsplit targets. Better structure predictions are brought about mainly by improved residue-residue contact predictions, and since this CASP also by distance predictions, achieved through state-of-the-art machine learning methods which also progressed to work with slightly shallower alignments compared to CASP12. As we reach a new realm of tertiary structure prediction quality, new directions are proposed and explored for future CASPs: (a) dropping splitting into EUs, (b) rethinking difficulty metrics probably in terms of contact and distance predictions, (c) assessing also side chains for models of high backbone accuracy, and (d) assessing residue-wise and possibly residue-residue quality estimates.

Abstract P250: Unveiling Binding Pocket Structure Of Mas Receptor And Its Interaction With Angiotensin-(1-7)

Substantial amount of preclinical data points out Mas receptor as a promising pharmacological target to treat cardiovascular diseases. Here, we aimed to characterize Mas binding pocket and its interaction with Angiotensin-(Ang)-1-7, in order to provide insights for potential structure-based drug discovery. Homology modeling based on AT1 and AT2 crystal structures, combined with peptide docking (Schrodinger TM ) and molecular dynamics simulations (Gromacs) were employed to generate hypothesis about crucial interactions in protein-ligand complex. Proposed model was further validated in vitro using HEK cells transiently transfected with wild-type Mas or mutant Mas (Arg245Ala). Mas-receptor activation was assessed by FOXO1 dephosphorylation at Ser256 after 5 minutes incubation with Ang-1-7 (10 -7 M). Model of Mas/Ang-(1-7) complex revealed electrostatic interaction between Ang-(1-7) C-terminus and Mas Arg245. Interestingly, while Ang-1-7 incubation reduced FOXO1 phosphorylation by 35 % in HEK-Mas wild type (1.000 vs. 0.654 A.U., p<0.01, n=8), this effect was abolished in HEK-Mas mutant (1.000 vs. 0.990 A.U., n.s. , n=5). In sum, Arg245 of Mas is essential for the receptor activation, which is in accordance with A-779 antagonism (C-terminus modification of Ang-(1-7)), endorsing our proposed structure model. Further elucidation of the ligand-receptor interactions may drive rational drug discovery and lead to development of improved strategies for Mas targeting.

Nucleo-cytosolic Shuttling of ARGONAUTE1 Prompts a Revised Model of the Plant MicroRNA Pathway

Unlike in metazoans, plant microRNAs (miRNAs) undergo stepwise nuclear maturation before engaging cytosolic, sequence-complementary transcripts in association with the silencing effector protein ARGONAUTE1 (AGO1). Since their discovery, how and under which form plant miRNAs translocate to the cytosol has remained unclear, as has their sub-cellular AGO1 loading site(s). Here, we show that the N termini of all plant AGO1s contain a nuclear-localization (NLS) and nuclear-export signal (NES) that, in Arabidopsis thaliana (At), enables AtAGO1 nucleo-cytosolic shuttling in a Leptomycin-B-inhibited manner, diagnostic of CRM1(EXPO1)/NES-dependent nuclear export. Nuclear-only AtAGO1 contains the same 2'O-methylated miRNA cohorts as its nucleo-cytosolic counterpart, but it preferentially interacts with the miRNA loading chaperone HSP90. Furthermore, mature miRNA translocation and miRNA-mediated silencing both require AtAGO1 nucleo-cytosolic shuttling. These findings lead us to propose a substantially revised view of the plant miRNA pathway in which miRNAs are matured, methylated, loaded into AGO1 in the nucleus, and exported to the cytosol as AGO1:miRNA complexes in a CRM1(EXPO1)/NES-dependent manner.

Functional assays for the assessment of the pathogenicity of variants of GOSR2, an ER-to-Golgi SNARE involved in progressive myoclonus epilepsies

Progressive myoclonus epilepsies (PMEs) are inherited disorders characterized by myoclonus, generalized tonic-clonic seizures, and ataxia. One of the genes that is associated with PME is the ER-to-Golgi Q_b-SNARE GOSR2, which forms a SNARE complex with syntaxin-5, Bet1 and Sec22b. Most PME patients are homo­zygous for a p.Gly144Trp mutation and develop similar clinical presentations. Recently, a patient who was compound heterozygous for p.Gly144Trp and a previously unseen p.Lys164del mutation was identified. Because this patient presented with a milder disease phenotype, we hypothesized that the p.Lys164del mutation may be less severe compared to p.Gly144Trp. To characterize the effect of the p.Gly144Trp and p.Lys164del mutations, both of which are present in the SNARE motif of GOSR2, we examined the corresponding mutations in the yeast ortholog Bos1. Yeasts expressing the orthologous mutants in Bos1 showed impaired growth, suggesting a partial loss of function, which was more severe for the Bos1 p.Gly176Trp mutation. Using anisotropy and gel filtration, we report that Bos1 p.Gly176Trp and p.Arg196del are capable of complex formation, but with partly reduced activity. Molecular dynamics (MD) simulations showed that the hydrophobic core, which triggers SNARE complex formation, is compromised due to the glycine-to-tryptophan substitution in both GOSR2 and Bos1. In contrast, the deletion of residue p.Lys164 (or p.Arg196del in Bos1) interferes with the formation of hydrogen bonds between GOSR2 and syntaxin-5. Despite these perturbations, all SNARE complexes stayed intact during longer simulations. Thus, our data suggest that the milder course of disease in compound heterozygous PME is due to less severe impairment of the SNARE function.

Summary: Mutations in the Q_b-SNARE GOSR2 cause progressive myoclonus epilepsies. The authors report the effect of two mutations on SNARE function to investigate their correlation with progression and severity of disease.

Conserved Lipid and Small-Molecule Modulation of COQ8 Reveals Regulation of the Ancient Kinase-like UbiB Family

Human COQ8A (ADCK3) and Saccharomyces cerevisiae Coq8p (collectively COQ8) are UbiB family proteins essential for mitochondrial coenzyme Q (CoQ) biosynthesis. However, the biochemical activity of COQ8 and its direct role in CoQ production remain unclear, in part due to lack of known endogenous regulators of COQ8 function and of effective small molecules for probing its activity in vivo. Here, we demonstrate that COQ8 possesses evolutionarily conserved ATPase activity that is activated by binding to membranes containing cardiolipin and by phenolic compounds that resemble CoQ pathway intermediates. We further create an analog-sensitive version of Coq8p and reveal that acute chemical inhibition of its endogenous activity in yeast is sufficient to cause respiratory deficiency concomitant with CoQ depletion. Collectively, this work defines lipid and small-molecule modulators of an ancient family of atypical kinase-like proteins and establishes a chemical genetic system for further exploring the mechanistic role of COQ8 in CoQ biosynthesis.

Recruitment of the amyloid precursor protein by γ-secretase at the synaptic plasma membrane

Γ-secretase is a membrane-embedded protease that cleaves single transmembrane helical domains of various integral membrane proteins. The amyloid precursor protein (APP) is an important substrate due to its pathological relevance to Alzheimer's disease. The mechanism of the cleavage of APP by γ-secretase that leads to accumulation of Alzheimer's disease causing amyloid-β (Aβ) is still unknown. Coarse-grained molecular dynamics simulations in this study reveal initial lipids raft formation near the catalytic site of γ-secretase as well as changes in dynamic behavior of γ-secretase once interacting with APP. The results suggest a precursor of the APP binding mode and hint at conformational changes of γ-secretase in the nicastrin (NCT) domain upon APP binding.

Structural, physicochemical and dynamic features conserved within the aerolysin pore-forming toxin family

Aerolysin is the founding member of a major class of β-pore-forming toxins (β-PFTs) found throughout all kingdoms of life. PFTs are cytotoxic proteins produced as soluble monomers, which oligomerize at the membrane of target host cells forming pores that may lead to osmotic lysis and cell death. Besides their role in microbial infection, they have become interesting for their potential as biotechnological sensors and delivery systems. Using an approach that integrates bioinformatics with molecular modeling and simulation, we looked for conserved features across this large toxin family. The cell surface-binding domains present high variability within the family to provide membrane receptor specificity. On the contrary, the novel concentric double β-barrel structure found in aerolysin is highly conserved in terms of sequence, structure and conformational dynamics, which likely contribute to preserve a common transition mechanism from the prepore to the mature pore within the family.Our results point to the key role of several amino acids in the conformational changes needed for oligomerization and further pore formation, such as Y221, W227, P248, Q263 and L277, which we propose are involved in the release of the stem loop and the two adjacent β-strands to form the transmembrane β-barrel.

Molecular Effects of Concentrated Solutes on Protein Hydration, Dynamics, and Electrostatics

Most studies of protein structure and function are performed in dilute conditions, but proteins typically experience high solute concentrations in their physiological scenarios and biotechnological applications. High solute concentrations have well-known effects on coarse protein traits like stability, diffusion, and shape, but likely also perturb other traits through finer effects pertinent at the residue and atomic levels. Here, NMR and molecular dynamics investigations on ubiquitin disclose variable interactions with concentrated solutes that lead to localized perturbations of the protein's surface, hydration, electrostatics, and dynamics, all dependent on solute size and chemical properties. Most strikingly, small polar uncharged molecules are sticky on the protein surface, whereas charged small molecules are not, but the latter still perturb the internal protein electrostatics as they diffuse nearby. Meanwhile, interactions with macromolecular crowders are favored mainly through hydrophobic, but not through polar, surface patches. All the tested small solutes strongly slow down water exchange at the protein surface, whereas macromolecular crowders do not exert such strong perturbation. Finally, molecular dynamics simulations predict that unspecific interactions slow down microsecond- to millisecond-timescale protein dynamics despite having only mild effects on pico- to nanosecond fluctuations as corroborated by NMR. We discuss our results in the light of recent advances in understanding proteins inside living cells, focusing on the physical chemistry of quinary structure and cellular organization, and we reinforce the idea that proteins should be studied in native-like media to achieve a faithful description of their function.

Signal Sensing and Transduction by Histidine Kinases as Unveiled through Studies on a Temperature Sensor

Histidine kinases (HK) are the sensory proteins of two-component systems, responsible for a large fraction of bacterial responses to stimuli and environmental changes. Prototypical HKs are membrane-bound proteins that phosphorylate cognate response regulator proteins in the cytoplasm upon signal detection in the membrane or periplasm. HKs stand as potential drug targets but also constitute fascinating systems for studying proteins at work, specifically regarding the chemistry and mechanics of signal detection, transduction through the membrane, and regulation of catalytic outputs. In this Account, we focus on Bacillus subtilis DesK, a membrane-bound HK part of a two-component system that maintains appropriate membrane fluidity at low growth temperatures. Unlike most HKs, DesK has no extracytoplasmic signal-sensing domains; instead, sensing is carried out by 10 transmembrane helices (coming from two protomers) arranged in an unknown structure. The fifth transmembrane helix from each protomer connects, without any of the intermediate domains found in other HKs, into the dimerization and histidine phosphotransfer (DHp) domain located in the cytoplasm, which is followed by the ATP-binding domains (ABD). Throughout the years, genetic, biochemical, structural, and computational studies on wild-type, mutant, and truncated versions of DesK allowed us to dissect several aspects of DesK's functioning, pushing forward a more general understanding of its own structure/function relationships as well as those of other HKs. We have shown that the sensing mechanism is rooted in temperature-dependent membrane properties, most likely a combination of thickness, fluidity, and water permeability, and we have proposed possible mechanisms by which DesK senses these properties and transduces the signals. X-ray structures and computational models have revealed structural features of TM and cytoplasmic regions in DesK's kinase- and phosphatase-competent states. Biochemical and genetic experiments and molecular simulations further showed that reversible formation of a two-helix coiled coil in the fifth TM segment and the N-terminus of the cytoplasmic domain is essential for the sensing and signal transduction mechanisms. Together with other structural and functional works, the emerging picture suggests that diverse HKs possess distinct sensing and transduction mechanisms but share as rather general features (i) a symmetric phosphatase state and an asymmetric kinase state and (ii) similar functional outputs on the conserved DHp and ABD domains, achieved through different mechanisms that depend on the nature of the initial signal. We here advance (iii) an important role for TM prolines in transducing the initial signals to the cytoplasmic coiled coils, based on simulations of DesK's TM helices and our previous work on a related HK, PhoQ. Lastly, evidence for DesK, PhoQ, BvgS, and DctB HKs shows that (iv) overall catalytic output is tuned by a delicate balance between hydration potentials, coiled coil stability, and exposure of hydrophobic surface patches at their cytoplasmic coiled coils and at the N-terminal and C-terminal sides of their TM helices. This balance is so delicate that small perturbations, either physiological signals or induced by mutations, lead to large remodeling of the underlying conformational landscape achieving clear-cut changes in catalytic output, mirroring the required response speed of these systems for proper biological function.

Transcriptional regulatory logic of the diurnal cycle in the mouse liver

Many organisms exhibit temporal rhythms in gene expression that propel diurnal cycles in physiology. In the liver of mammals, these rhythms are controlled by transcription-translation feedback loops of the core circadian clock and by feeding-fasting cycles. To better understand the regulatory interplay between the circadian clock and feeding rhythms, we mapped DNase I hypersensitive sites (DHSs) in the mouse liver during a diurnal cycle. The intensity of DNase I cleavages cycled at a substantial fraction of all DHSs, suggesting that DHSs harbor regulatory elements that control rhythmic transcription. Using chromatin immunoprecipitation followed by DNA sequencing (ChIP-seq), we found that hypersensitivity cycled in phase with RNA polymerase II (Pol II) loading and H3K27ac histone marks. We then combined the DHSs with temporal Pol II profiles in wild-type (WT) and Bmal1-/- livers to computationally identify transcription factors through which the core clock and feeding-fasting cycles control diurnal rhythms in transcription. While a similar number of mRNAs accumulated rhythmically in Bmal1-/- compared to WT livers, the amplitudes in Bmal1-/- were generally lower. The residual rhythms in Bmal1-/- reflected transcriptional regulators mediating feeding-fasting responses as well as responses to rhythmic systemic signals. Finally, the analysis of DNase I cuts at nucleotide resolution showed dynamically changing footprints consistent with dynamic binding of CLOCK:BMAL1 complexes. Structural modeling suggested that these footprints are driven by a transient heterotetramer binding configuration at peak activity. Together, our temporal DNase I mappings allowed us to decipher the global regulation of diurnal transcription rhythms in the mouse liver.

Protein post-translational modifications: In silico prediction tools and molecular modeling

Post-translational modifications (PTMs) occur in almost all proteins and play an important role in numerous biological processes by significantly affecting proteins' structure and dynamics. Several computational approaches have been developed to study PTMs (e.g., phosphorylation, sumoylation or palmitoylation) showing the importance of these techniques in predicting modified sites that can be further investigated with experimental approaches. In this review, we summarize some of the available online platforms and their contribution in the study of PTMs. Moreover, we discuss the emerging capabilities of molecular modeling and simulation that are able to complement these bioinformatics methods, providing deeper molecular insights into the biological function of post-translational modified proteins.

Disentangling constraints using viability evolution principles in integrative modeling of macromolecular assemblies

Predicting the structure of large molecular assemblies remains a challenging task in structural biology when using integrative modeling approaches. One of the main issues stems from the treatment of heterogeneous experimental data used to predict the architecture of native complexes. We propose a new method, applied here for the first time to a set of symmetrical complexes, based on evolutionary computation that treats every available experimental input independently, bypassing the need to balance weight components assigned to aggregated fitness functions during optimization.

A Biologically-validated HCV E1E2 Heterodimer Structural Model

The design of vaccine strategies and the development of drugs targeting the early stages of Hepatitis C virus (HCV) infection are hampered by the lack of structural information about its surface glycoproteins E1 and E2, the two constituents of HCV entry machinery. Despite the recent crystal resolution of limited versions of both proteins in truncated form, a complete picture of the E1E2 complex is still missing. Here we combined deep computational analysis of E1E2 secondary, tertiary and quaternary structure with functional and immunological mutational analysis across E1E2 in order to propose an in silico model for the ectodomain of the E1E2 heterodimer. Our model describes E1-E2 ectodomain dimerization interfaces, provides a structural explanation of E1 and E2 immunogenicity and sheds light on the molecular processes and disulfide bridges isomerization underlying the conformational changes required for fusion. Comprehensive alanine mutational analysis across 553 residues of E1E2 also resulted in identifying the epitope maps of diverse mAbs and the disulfide connectivity underlying E1E2 native conformation. The predicted structure unveils E1 and E2 structures in complex, thus representing a step towards the rational design of immunogens and drugs inhibiting HCV entry.