Structural biology in cellulo: Minding the gap between conceptualization and realization

Recent technological advances have deepened our perception of cellular structure. However, most structural data doesn't originate from intact cells, limiting our understanding of cellular processes. Here, we discuss current and future developments that will bring us towards a structural picture of the cell. Electron cryotomography is the standard bearer, with its ability to provide in cellulo snapshots. Single-particle electron microscopy (of purified biomolecules and of complex mixtures) and covalent crosslinking combined with mass spectrometry also have significant roles to play, as do artificial intelligence algorithms in their many guises. To integrate these multiple approaches, data curation and standardisation will be critical - as is the need to expand efforts beyond our current protein-centric view to the other (macro)molecules that sustain life.

Delineating organizational principles of the endogenous L-A virus by cryo-EM and computational analysis of native cell extracts

The high abundance of most viruses in infected host cells benefits their structural characterization. However, endogenous viruses are present in low copy numbers and are therefore challenging to investigate. Here, we retrieve cell extracts enriched with an endogenous virus, the yeast L-A virus. The determined cryo-EM structure discloses capsid-stabilizing cation-π stacking, widespread across viruses and within the Totiviridae, and an interplay of non-covalent interactions from ten distinct capsomere interfaces. The capsid-embedded mRNA decapping active site trench is supported by a constricting movement of two flexible opposite-facing loops. tRNA-loaded polysomes and other biomacromolecules, presumably mRNA, are found in virus proximity within the cell extract. Mature viruses participate in larger viral communities resembling their rare in-cell equivalents in terms of size, composition, and inter-virus distances. Our results collectively describe a 3D-architecture of a viral milieu, opening the door to cell-extract-based high-resolution structural virology.

Protein-Rich Rafts in Hybrid Polymer/Lipid Giant Unilamellar Vesicles

Considerable attention has been dedicated to lipid rafts due to their importance in numerous cell functions such as membrane trafficking, polarization, and signaling. Next to studies in living cells, artificial micrometer-sized vesicles with a minimal set of components are established as a major tool to understand the phase separation dynamics and their intimate interplay with membrane proteins. In parallel, mixtures of phospholipids and certain amphiphilic polymers simultaneously offer an interface for proteins and mimic this segregation behavior, presenting a tangible synthetic alternative for fundamental studies and bottom-up design of cellular mimics. However, the simultaneous insertion of complex and sensitive membrane proteins is experimentally challenging and thus far has been largely limited to natural lipids. Here, we present the co-reconstitution of the proton pump bo3 oxidase and the proton consumer ATP synthase in hybrid polymer/lipid giant unilamellar vesicles (GUVs) via fusion/electroformation. Variations of the current method allow for tailored reconstitution protocols and control of the vesicle size. In particular, mixing of protein-free and protein-functionalized nanosized vesicles in the electroformation film results in larger GUVs, while separate reconstitution of the respiratory enzymes enables higher ATP synthesis rates. Furthermore, protein labeling provides a synthetic mechanism for phase separation and protein sequestration, mimicking lipid- and protein-mediated domain formation in nature. The latter means opens further possibilities for re-enacting phenomena like supercomplex assembly or symmetry breaking and enriches the toolbox of bottom-up synthetic biology.

Effect of Molecular Dynamics and Internal Water Contact on the Photophysical Properties of Red pH-Sensitive Proteins

The pH dependence of the absorption and (time-resolved) fluorescence of two red-shifted fluorescent proteins, mCardinal and mNeptune, was investigated. Decay-associated spectra were measured following fluorescence excitation at 470 nm in PBS buffer with a pH that ranged from 5.5 to 8.0. The fluorescence of both proteins shows two different decay components. mCardinal exhibits an increase in the long-lived fluorescence component with acidification from 1.34 ns at pH 8.0 to 1.62 ns at pH 5.5. An additional fast decay component with 0.64 ns at pH 8.0 up to 1.1 ns at pH 5.5 was found to be blue-shifted compared to the long-lived component. The fluorescence lifetime of mNeptune is insensitive to pH. DAS of mCardinal were simulated assuming a coupled two-level system to describe the 1S state of the chromophore within two different conformations of the protein. MD simulations were conducted to correlate the experimentally observed pH-induced change in the lifetime in mCardinal with its molecular properties. While the chromophores of both protein variants are stabilized by the same number of hydrogen bonds, it was found that the chromophore in mCardinal exhibits more water contacts compared to mNeptune. In mCardinal, interaction between the chromophore and Glu-145 is reduced as compared to mNeptune, but interaction with Thr-147 which is Ser-147 in mNeptune is stronger in mCardinal. Therefore, the dynamics of the excited-state proton transfer (ESPT) might be different in mCardinal and mNeptune. The pH dependency of ESPT is suggested as a key mechanism for pH sensitivity.

Cryo-EM of a heterogeneous biochemical fraction elucidates multiple protein complexes from a multicellular thermophilic eukaryote

Biomolecular complexes and their interactions govern cellular structure and function. Understanding their architecture is a prerequisite for dissecting the cell's inner workings, but their higher-order assembly is often transient and challenging for structural analysis. Here, we performed cryo-EM on a single, highly heterogeneous biochemical fraction derived from Chaetomium thermophilum cell extracts to visualize the biomolecular content of the multicellular eukaryote. After cryo-EM single-particle image processing, results showed that a simultaneous three-dimensional structural characterization of multiple chemically diverse biomacromolecules is feasible. Namely, the thermophilic, eukaryotic complexes of (a) ATP citrate-lyase, (b) Hsp90, (c) 20S proteasome, (d) Hsp60 and (e) UDP-glucose pyrophosphorylase were characterized. In total, all five complexes have been structurally dissected in a thermophilic eukaryote in a total imaged sample area of 190.64 μm2, and two, in particular, 20S proteasome and Hsp60, exhibit side-chain resolution features. The C. thermophilum Hsp60 near-atomic model was resolved at 3.46 Å (FSC = 0.143) and shows a hinge-like conformational change of its equatorial domain, highly similar to the one previously shown for its bacterial orthologue, GroEL. This work demonstrates that cryo-EM of cell extracts will greatly accelerate the structural analysis of cellular complexes and provide unprecedented opportunities to annotate architectures of biomolecules in a holistic approach.

The structural principles underlying molybdenum insertase complex assembly

Within the cell, the trace element molybdenum (Mo) is only biologically active when complexed either within the nitrogenase-specific FeMo cofactor or within the molybdenum cofactor (Moco). Moco consists of an organic part, called molybdopterin (MPT) and an inorganic part, that is, the Mo-center. The enzyme which catalyzes the Mo-center formation is the molybdenum insertase (Mo-insertase). Mo-insertases consist of two functional domains called G- and E-domain. The G-domain catalyzes the formation of adenylated MPT (MPT-AMP), which is the substrate for the E-domain, that catalyzes the actual molybdate insertion reaction. Though the functions of E- and G-domain have been elucidated to great structural and mechanistic detail, their combined function is poorly characterized. In this work, we describe a structural model of the eukaryotic Mo-insertase Cnx1 complex that was generated based on cross-linking mass spectrometry combined with computational modeling. We revealed Cnx1 to form an asymmetric hexameric complex which allows the E- and G-domain active sites to align in a catalytic productive orientation toward each other.

Enabling cryo-EM density interpretation from yeast native cell extracts by proteomics data and AlphaFold structures

In the cellular context, proteins participate in communities to perform their function. The detection and identification of these communities as well as in-community interactions has long been the subject of investigation, mainly through proteomics analysis with mass spectrometry. With the advent of cryogenic electron microscopy and the "resolution revolution," their visualization has recently been made possible, even in complex, native samples. The advances in both fields have resulted in the generation of large amounts of data, whose analysis requires advanced computation, often employing machine learning approaches to reach the desired outcome. In this work, we first performed a robust proteomics analysis of mass spectrometry (MS) data derived from a yeast native cell extract and used this information to identify protein communities and inter-protein interactions. Cryo-EM analysis of the cell extract provided a reconstruction of a biomolecule at medium resolution (∼8 Å (FSC = 0.143)). Utilizing MS-derived proteomics data and systematic fitting of AlphaFold-predicted atomic models, this density was assigned to the 2.6 MDa complex of yeast fatty acid synthase. Our proposed workflow identifies protein complexes in native cell extracts from Saccharomyces cerevisiae by combining proteomics, cryo-EM, and AI-guided protein structure prediction.

Cryo-EM structure of a plant photosystem II supercomplex with light-harvesting protein Lhcb8 and α-tocopherol

The heart of oxygenic photosynthesis is the water-splitting photosystem II (PSII), which forms supercomplexes with a variable amount of peripheral trimeric light-harvesting complexes (LHCII). Our knowledge of the structure of green plant PSII supercomplex is based on findings obtained from several representatives of green algae and flowering plants; however, data from a non-flowering plant are currently missing. Here we report a cryo-electron microscopy structure of PSII supercomplex from spruce, a representative of non-flowering land plants, at 2.8 Å resolution. Compared with flowering plants, PSII supercomplex in spruce contains an additional Ycf12 subunit, Lhcb4 protein is replaced by Lhcb8, and trimeric LHCII is present as a homotrimer of Lhcb1. Unexpectedly, we have found α-tocopherol (α-Toc)/α-tocopherolquinone (α-TQ) at the boundary between the LHCII trimer and the inner antenna CP43. The molecule of α-Toc/α-TQ is located close to chlorophyll a614 of one of the Lhcb1 proteins and its chromanol/quinone head is exposed to the thylakoid lumen. The position of α-Toc in PSII supercomplex makes it an ideal candidate for the sensor of excessive light, as α-Toc can be oxidized to α-TQ by high-light-induced singlet oxygen at low lumenal pH. The molecule of α-TQ appears to shift slightly into the PSII supercomplex, which could trigger important structure-functional modifications in PSII supercomplex. Inspection of the previously reported cryo-electron microscopy maps of PSII supercomplexes indicates that α-Toc/α-TQ can be present at the same site also in PSII supercomplexes from flowering plants, but its identification in the previous studies has been hindered by insufficient resolution.

Structural assessment of the full-length wild-type tumor suppressor protein p53 by mass spectrometry-guided computational modeling

The tetrameric tumor suppressor p53 represents a great challenge for 3D-structural analysis due to its high degree of intrinsic disorder (ca. 40%). We aim to shed light on the structural and functional roles of p53's C-terminal region in full-length, wild-type human p53 tetramer and their importance for DNA binding. For this, we employed complementary techniques of structural mass spectrometry (MS) in an integrated approach with computational modeling. Our results show no major conformational differences in p53 between DNA-bound and DNA-free states, but reveal a substantial compaction of p53's C-terminal region. This supports the proposed mechanism of unspecific DNA binding to the C-terminal region of p53 prior to transcription initiation by specific DNA binding to the core domain of p53. The synergies between complementary structural MS techniques and computational modeling as pursued in our integrative approach is envisioned to serve as general strategy for studying intrinsically disordered proteins (IDPs) and intrinsically disordered region (IDRs).

Structural analysis of an endogenous 4-megadalton succinyl-CoA-generating metabolon

The oxoglutarate dehydrogenase complex (OGDHc) participates in the tricarboxylic acid cycle and, in a multi-step reaction, decarboxylates α-ketoglutarate, transfers succinyl to CoA, and reduces NAD+. Due to its pivotal role in metabolism, OGDHc enzymatic components have been studied in isolation; however, their interactions within the endogenous OGDHc remain elusive. Here, we discern the organization of a thermophilic, eukaryotic, native OGDHc in its active state. By combining biochemical, biophysical, and bioinformatic methods, we resolve its composition, 3D architecture, and molecular function at 3.35 Å resolution. We further report the high-resolution cryo-EM structure of the OGDHc core (E2o), which displays various structural adaptations. These include hydrogen bonding patterns confining interactions of OGDHc participating enzymes (E1o-E2o-E3), electrostatic tunneling that drives inter-subunit communication, and the presence of a flexible subunit (E3BPo), connecting E2o and E3. This multi-scale analysis of a succinyl-CoA-producing native cell extract provides a blueprint for structure-function studies of complex mixtures of medical and biotechnological value.

Solubilization, purification, and characterization of the hexameric form of phosphatidylserine synthase from Candida albicans

Phosphatidylserine (PS) synthase from Candida albicans, encoded by the CHO1 gene, has been identified as a potential drug target for new antifungals against systemic candidiasis. Rational drug design or small molecule screening are effective ways to identify specific inhibitors of Cho1, but both will be facilitated by protein purification. Due to the transmembrane nature of Cho1, methods were needed to solubilize and purify the native form of Cho1. Here, we used six non-ionic detergents and three styrene maleic acids (SMAs) to solubilize an HA-tagged Cho1 protein from the total microsomal fractions. Blue native PAGE (BN-PAGE) and immunoblot analysis revealed a single band corresponding to Cho1 in all detergent-solubilized fractions, while two bands were present in the SMA2000-solubilized fraction. Our enzymatic assay suggests that digitonin- or DDM-solubilized enzyme has the most PS synthase activity. Pull-downs of HA-tagged Cho1 in the digitonin-solubilized fraction reveal an apparent MW of Cho1 consistent with a hexamer. Furthermore, negative-staining electron microscopy analysis and AlphaFold2 structure prediction modeling suggest the hexamer is composed of a trimer of dimers. We purified Cho1 protein to near-homogeneity as a hexamer using affinity chromatography and TEV protease treatment, and optimized Cho1 enzyme activity for manganese and detergent concentrations, temperature (24°C), and pH (8.0). The purified Cho1 has a K_m for its substrate CDP-diacylglycerol of 72.20 μM with a V_max of 0.079 nmol/(μg*min) while exhibiting a sigmoidal kinetic curve for its other substrate serine, indicating cooperative binding. Purified hexameric Cho1 can potentially be used in downstream structure determination and small drug screening.

Novel exotic alleles of EARLY FLOWERING 3 determine plant development in barley

EARLY FLOWERING 3 (ELF3) is an important regulator of various physiological and developmental processes and hence may serve to improve plant adaptation which will be substantial for future plant breeding. To expand the limited knowledge on barley ELF3 in determining agronomic traits, we conducted field studies with heterogeneous inbred families (HIFs) derived from selected lines of the wild barley nested association mapping population HEB-25. During two growing seasons, phenotypes of nearly isogenic HIF sister lines, segregating for exotic and cultivated alleles at the ELF3 locus, were compared for ten developmental and yield-related traits. We determine novel exotic ELF3 alleles and show that HIF lines, carrying the exotic ELF3 allele, accelerated plant development compared to the cultivated ELF3 allele, depending on the genetic background. Remarkably, the most extreme effects on phenology could be attributed to one exotic ELF3 allele differing from the cultivated Barke ELF3 allele in only one SNP. This SNP causes an amino acid substitution (W669G), which predictively has an impact on the protein structure of ELF3, thereby possibly affecting phase separation behaviour and nano-compartment formation of ELF3 and, potentially, also affecting its local cellular interactions causing significant trait differences between HIF sister lines.

Encapsulation of cannabidiol in oil-in-water nanoemulsions and nanoemulsion-filled hydrogels: A structure and biological assessment study

HYPOTHESIS

Lipophilic cannabidiol can be solubilized in oil-in water nanoemulsions, which can then be impregnated into chitosan hydrogels forming another colloidal system that will facilitate cannabidiol's release. The delivery from both systems was compared, alongside structural and biological studies, to clarify the effect of the two carriers' structure on the release and toxicity of the systems.

EXPERIMENTS

Oil-in-water nanoemulsions (NEs) and the respective nanoemulsion-filled chitosan hydrogels (NE/HGs) were formulated as carriers of cannabidiol (CBD). Size, polydispersity and stability of the NEs were evaluated and then membrane dynamics, shape and structure of both systems were investigated with EPR spin probing, SAXS and microscopy. Biocompatibility of the colloidal delivery systems was evaluated through cytotoxicity tests over normal human skin fibroblasts. An ex vivo permeation protocol using porcine ear skin was implemented to assess the release of CBD and its penetration through the skin.

FINDINGS

Incorporation of the NEs in chitosan hydrogels does not significantly affect their structural properties as evidenced through SAXS, EPR and confocal microscopy. These findings indicate the successful development of a novel nanocarrier that preserves the NE structure with the CBD remaining encapsulated in the oil core while providing new rheological properties advantageous over NEs. Moreover, NE/HGs proved to be more efficient as a carrier for the release of CBD. Cell viability assessment revealed high biocompatibility of the proposed colloids.

Cracking the code of cellular protein-protein interactions: Alphafold and whole-cell crosslinking to the rescue

Integration of experimental and computational methods is crucial to better understanding protein-protein interactions (PPIs), ideally in their cellular context. In their recent work, Rappsilber and colleagues (O'Reilly et al, 2023) identified bacterial PPIs using an array of approaches. They combined whole-cell crosslinking, co-fractionation mass spectrometry, and open-source data mining with artificial intelligence (AI)-based structure prediction of PPIs in the well-studied organism Bacillus subtilis. This innovative approach reveals architectural knowledge for in-cell PPIs that are often lost upon cell lysis, making it applicable to genetically intractable organisms such as pathogenic bacteria.

IRAA: A statistical tool for investigating a protein-protein interaction interface from multiple structures

Understanding protein-protein interactions (PPIs) is fundamental to infer how different molecular systems work. A major component to model molecular recognition is the buried surface area (BSA), that is, the area that becomes inaccessible to solvent upon complex formation. To date, many attempts tried to connect BSA to molecular recognition principles, and in particular, to the underlying binding affinity. However, the most popular approach to calculate BSA is to use a single (or in some cases few) bound structures, consequently neglecting a wealth of structural information of the interacting proteins derived from ensembles corresponding to their unbound and bound states. Moreover, the most popular method inherently assumes the component proteins to bind as rigid entities. To address the above shortcomings, we developed a Monte Carlo method-based Interface Residue Assessment Algorithm (IRAA), to calculate a combined distribution of BSA for a given complex. Further, we apply our algorithm to human ACE2 and SARS-CoV-2 Spike protein complex, a system of prime importance. Results show a much broader distribution of BSA compared to that obtained from only the bound structure or structures and extended residue members of the interface with implications to the underlying biomolecular recognition. We derive that specific interface residues of ACE2 and of S-protein are consistently highly flexible, whereas other residues systematically show minor conformational variations. In effect, IRAA facilitates the use of all available structural data for any biomolecular complex of interest, extracting quantitative parameters with statistical significance, thereby providing a deeper biophysical understanding of the molecular system under investigation.

Cryo-Electron Microscopy Snapshots of Eukaryotic Membrane Proteins in Native Lipid-Bilayer Nanodiscs

New technologies for purifying membrane-bound protein complexes in combination with cryo-electron microscopy (EM) have recently allowed the exploration of such complexes under near-native conditions. In particular, polymer-encapsulated nanodiscs enable the study of membrane proteins at high resolution while retaining protein-protein and protein-lipid interactions within a lipid bilayer. However, this powerful technology has not been exploited to address the important question of how endogenous─as opposed to overexpressed─membrane proteins are organized within a lipid environment. In this work, we demonstrate that biochemical enrichment protocols for native membrane-protein complexes from Chaetomium thermophilum in combination with polymer-based lipid-bilayer nanodiscs provide a substantial improvement in the quality of recovered endogenous membrane-protein complexes. Mass spectrometry results revealed ∼1123 proteins, while multiple 2D class averages and two 3D reconstructions from cryo-EM data furnished prominent structural signatures. This integrated methodological approach to enriching endogenous membrane-protein complexes provides unprecedented opportunities for a deeper understanding of eukaryotic membrane proteomes.

Cryo-EM structure of the SEA complex

The SEA complex (SEAC) is a growth regulator that acts as a GTPase-activating protein (GAP) towards Gtr1, a Rag GTPase that relays nutrient status to the Target of Rapamycin Complex 1 (TORC1) in yeast1. Functionally, the SEAC has been divided into two subcomplexes: SEACIT, which has GAP activity and inhibits TORC1, and SEACAT, which regulates SEACIT2. This system is conserved in mammals: the GATOR complex, consisting of GATOR1 (SEACIT) and GATOR2 (SEACAT), transmits amino acid3 and glucose4 signals to mTORC1. Despite its importance, the structure of SEAC/GATOR, and thus molecular understanding of its function, is lacking. Here, we solve the cryo-EM structure of the native eight-subunit SEAC. The SEAC has a modular structure in which a COPII-like cage corresponding to SEACAT binds two flexible wings, which correspond to SEACIT. The wings are tethered to the core via Sea3, which forms part of both modules. The GAP mechanism of GATOR1 is conserved in SEACIT, and GAP activity is unaffected by SEACAT in vitro. In vivo, the wings are essential for recruitment of the SEAC to the vacuole, primarily via the EGO complex. Our results indicate that rather than being a direct inhibitor of SEACIT, SEACAT acts as a scaffold for the binding of TORC1 regulators.

Filling the Gap with Long n-Alkanes: Incorporation of C20 and C30 into Phospholipid Membranes

Investigating how hydrophobic molecules mix with phospholipid bilayers and how they affect membrane properties is commonplace in biophysics. Despite this, a molecular-level empirical description of a membrane model as simple as a phospholipid bilayer with long linear hydrophobic chains incorporated is still missing. Here, we present an unprecedented molecular characterization of the incorporation of two long n-alkanes, n-eicosane (C20) and n-triacontane (C30) with 20 and 30 carbons, respectively, in phosphatidylcholine (PC) bilayers using a combination of experimental techniques (²H NMR, ³¹P NMR, ¹H-¹³C dipolar recoupling solid-state NMR, X-ray scattering, and cryogenic electron microscopy) and atomistic molecular dynamics (MD) simulations. At low hydration, deuterated C20 and C30 yield ²H NMR spectra evidencing anisotropic-motion, which demonstrates their miscibility in PC membranes up to a critical alkane-to-acyl-chain volume fraction, ϕ_c. The acquired ²H NMR spectra of C20 and C30 have notably different lineshapes. At low alkane volume fractions below ϕ_c, CHARMM36 MD simulations predict such ²H NMR spectra qualitatively and thus enable an atomistic-level interpretation of the spectra. Above ϕ_c, the ²H NMR lineshapes become characteristic of motions in the intermediate-regime that, together with the MD simulation results, suggest the onset of immiscibility between the alkane molecules and the acyl chains. For all the systems investigated, the phospholipid molecular structure is unperturbed by the presence of the alkanes. However, at conditions of excess hydration and at surprisingly low alkane fractions below ϕ_c, a peak characteristic of isotropic motion is observed in both the ²H spectra of the alkanes and ³¹P spectra of the phospholipids, strongly indicating that the incorporation of the alkanes induces a reduction on the average radius of the lipid vesicles.

An Electrostatically-steered Conformational Selection Mechanism Promotes SARS-CoV-2 Spike Protein Variation

After two years since the outbreak, the COVID-19 pandemic remains a global public health emergency. SARS-CoV-2 variants with substitutions on the spike (S) protein emerge increasing the risk of immune evasion and cross-species transmission. Here, we analyzed the evolution of the S protein as recorded in 276,712 samples collected before the start of vaccination efforts. Our analysis shows that most variants destabilize the S protein trimer, increase its conformational heterogeneity and improve the odds of the recognition by the host cell receptor. Most frequent substitutions promote overall hydrophobicity by replacing charged amino acids, reducing stabilizing local interactions in the unbound S protein trimer. Moreover, our results identify "forbidden" regions that rarely show any sequence variation, and which are related to conformational changes occurring upon fusion. These results are significant for understanding the structure and function of SARS-CoV-2 related proteins which is a critical step in vaccine development and epidemiological surveillance.

Cryo-EM structure of a tetrameric photosystem I from Chroococcidiopsis TS-821, a thermophilic, unicellular, non-heterocyst-forming cyanobacterium

Photosystem I (PSI) is one of two photosystems involved in oxygenic photosynthesis. PSI of cyanobacteria exists in monomeric, trimeric, and tetrameric forms, in contrast to the strictly monomeric form of PSI in plants and algae. The tetrameric organization raises questions about its structural, physiological, and evolutionary significance. Here we report the ∼3.72 Å resolution cryo-electron microscopy structure of tetrameric PSI from the thermophilic, unicellular cyanobacterium Chroococcidiopsis sp. TS-821. The structure resolves 44 subunits and 448 cofactor molecules. We conclude that the tetramer is arranged via two different interfaces resulting from a dimer-of-dimers organization. The localization of chlorophyll molecules permits an excitation energy pathway within and between adjacent monomers. Bioinformatics analysis reveals conserved regions in the PsaL subunit that correlate with the oligomeric state. Tetrameric PSI may function as a key evolutionary step between the trimeric and monomeric forms of PSI organization in photosynthetic organisms.

Cryo-EM snapshots of a native lysate provide structural insights into a metabolon-embedded transacetylase reaction

Found across all kingdoms of life, 2-keto acid dehydrogenase complexes possess prominent metabolic roles and form major regulatory sites. Although their component structures are known, their higher-order organization is highly heterogeneous, not only across species or tissues but also even within a single cell. Here, we report a cryo-EM structure of the fully active Chaetomium thermophilum pyruvate dehydrogenase complex (PDHc) core scaffold at 3.85 Å resolution (FSC = 0.143) from native cell extracts. By combining cryo-EM with macromolecular docking and molecular dynamics simulations, we resolve all PDHc core scaffold interfaces and dissect the residing transacetylase reaction. Electrostatics attract the lipoyl domain to the transacetylase active site and stabilize the coenzyme A, while apolar interactions position the lipoate in its binding cleft. Our results have direct implications on the structural determinants of the transacetylase reaction and the role of flexible regions in the context of the overall 10 MDa PDHc metabolon architecture.

Cross-Linking Mass Spectrometry for Investigating Protein Conformations and Protein-Protein Interactions─A Method for All Seasons

Mass spectrometry (MS) has become one of the key technologies of structural biology. In this review, the contributions of chemical cross-linking combined with mass spectrometry (XL-MS) for studying three-dimensional structures of proteins and for investigating protein-protein interactions are outlined. We summarize the most important cross-linking reagents, software tools, and XL-MS workflows and highlight prominent examples for characterizing proteins, their assemblies, and interaction networks in vitro and in vivo. Computational modeling plays a crucial role in deriving 3D-structural information from XL-MS data. Integrating XL-MS with other techniques of structural biology, such as cryo-electron microscopy, has been successful in addressing biological questions that to date could not be answered. XL-MS is therefore expected to play an increasingly important role in structural biology in the future.

En route to dynamic life processes by SNARE-mediated fusion of polymer and hybrid membranes

A variety of artificial cells springs from the functionalization of liposomes with proteins. However, these models suffer from low durability without repair and replenishment mechanisms, which can be partly addressed by replacing the lipids with polymers. Yet natural membranes are also dynamically remodeled in multiple cellular processes. Here, we show that synthetic amphiphile membranes also undergo fusion, mediated by the protein machinery for synaptic secretion. We integrated fusogenic SNAREs in polymer and hybrid vesicles and observed efficient membrane and content mixing. We determined bending rigidity and pore edge tension as key parameters for fusion and described its plausible progression through cryo-EM snapshots. These findings demonstrate that dynamic membrane phenomena can be reconstituted in synthetic materials, thereby providing new tools for the assembly of synthetic protocells.

Solubilization of artificial mitochondrial membranes by amphiphilic copolymers of different charge

Certain amphiphilic copolymers form lipid-bilayer nanodiscs from artificial and natural membranes, thereby rendering incorporated membrane proteins optimal for structural analysis. Recent studies have shown that the amphiphilicity of a copolymer strongly determines its solubilization efficiency. This is especially true for highly negatively charged membranes, which experience pronounced Coulombic repulsion with polyanionic polymers. Here, we present a systematic study on the solubilization of artificial multicomponent lipid vesicles that mimic inner mitochondrial membranes, which harbor essential membrane-protein complexes. In particular, we compared the lipid-solubilization efficiencies of established anionic with less densely charged or zwitterionic and even cationic copolymers in low- and high-salt concentrations. The nanodiscs formed under these conditions were characterized by dynamic light scattering and negative-stain electron microscopy, pointing to a bimodal distribution of nanodisc diameters with a considerable fraction of nanodiscs engaging in side-by-side interactions through their polymer rims. Overall, our results show that some recent, zwitterionic copolymers are best suited to solubilize negatively charged membranes at high ionic strengths even at low polymer/lipid ratios.

Nanoscale Model System for the Human Myelin Sheath

Neurodegenerative disorders are among the most common diseases in modern society. However, the molecular bases of diseases such as multiple sclerosis or Charcot-Marie-Tooth disease remain far from being fully understood. Research in this field is limited by the complex nature of native myelin and by difficulties in obtaining good in vitro model systems of myelin. Here, we introduce an easy-to-use model system of the myelin sheath that can be used to study myelin proteins in a native-like yet well-controlled environment. To this end, we present myelin-mimicking nanodiscs prepared through one of the amphiphilic copolymers styrene/maleic acid (SMA), diisobutylene/maleic acid (DIBMA), and styrene/maleimide sulfobetaine (SMA-SB). These nanodiscs were tested for their lipid composition using chromatographic (HPLC) and mass spectrometric (MS) methods and, utilizing spin probes within the nanodisc, their comparability with liposomes was studied. In addition, their binding behavior with bovine myelin basic protein (MBP) was scrutinized to ensure that the nanodiscs represent a suitable model system of myelin. Our results suggest that both SMA and SMA-SB are able to solubilize the myelin-like (cytoplasmic) liposomes without preferences for specific lipid headgroups or fatty acyl chains. In nanodiscs of both SMA and SMA-SB (called SMA(-SB)-lipid particles, short SMALPs or SMA-SBLPs, respectively), the polymers restrict the lipids' motion in the hydrophobic center of the bilayer. The headgroups of the lipids, however, are sterically less hindered in nanodiscs when compared with liposomes. Myelin-like SMALPs are able to bind bovine MBP, which can stack the lipid bilayers like in native myelin, showing the usability of these simple, well-controlled systems in further studies of protein-lipid interactions of native myelin.

Coupling proteomics and metabolomics for the unsupervised identification of protein-metabolite interactions in Chaetomium thermophilum

Protein-metabolite interactions play an important role in the cell's metabolism and many methods have been developed to screen them in vitro. However, few methods can be applied at a large scale and not alter biological state. Here we describe a proteometabolomic approach, using chromatography to generate cell fractions which are then analyzed with mass spectrometry for both protein and metabolite identification. Integrating the proteomic and metabolomic analyses makes it possible to identify protein-bound metabolites. Applying the concept to the thermophilic fungus Chaetomium thermophilum, we predict 461 likely protein-metabolite interactions, most of them novel. As a proof of principle, we experimentally validate a predicted interaction between the ribosome and isopentenyl adenine.

BRCA1-BRCT Mutations Alter the Subcellular Localization of BRCA1 In Vitro

BACKGROUND/AIM

Numerous missense mutations have been determined in the BRCT domain of the BRCA1 gene, affecting localization and interaction of BRCA1 with other proteins.

MATERIALS AND METHODS

We examined whether the M1775K and V1809F mutations in the BRCT domain affect BRCA1 cellular localization. Cells were transfected with pEGFP-C3-BRCA1 and detected by fluorescence microscopy.

RESULTS

Following induction of DNA damage, cytoplasmic mislocalization was observed for both M1775K and V1809F mutants compared to EGFP-BRCA1wt and the less common variant M1652I. These results indicate that M1775K and V1809F mutations may change the function of the protein by affecting BRCA1 localization.

CONCLUSION

There is a correlation between subcellular localization of BRCA1 and diminished DNA repair observed in breast cancer cells, which may be explained by structural variations and altered binding properties of phosphopeptides.

Detecting Protein Communities in Native Cell Extracts by Machine Learning: A Structural Biologist's Perspective

Native cell extracts hold great promise for understanding the molecular structure of ordered biological systems at high resolution. This is because higher-order biomolecular interactions, dubbed as protein communities, may be retained in their (near-)native state, in contrast to extensively purifying or artificially overexpressing the proteins of interest. The distinct machine-learning approaches are applied to discover protein-protein interactions within cell extracts, reconstruct dedicated biological networks, and report on protein community members from various organisms. Their validation is also important, e.g., by the cross-linking mass spectrometry or cell biology methods. In addition, the cell extracts are amenable to structural analysis by cryo-electron microscopy (cryo-EM), but due to their inherent complexity, sorting structural signatures of protein communities derived by cryo-EM comprises a formidable task. The application of image-processing workflows inspired by machine-learning techniques would provide improvements in distinguishing structural signatures, correlating proteomic and network data to structural signatures and subsequently reconstructed cryo-EM maps, and, ultimately, characterizing unidentified protein communities at high resolution. In this review article, we summarize recent literature in detecting protein communities from native cell extracts and identify the remaining challenges and opportunities. We argue that the progress in, and the integration of, machine learning, cryo-EM, and complementary structural proteomics approaches would provide the basis for a multi-scale molecular description of protein communities within native cell extracts.

Integrative structure of a 10-megadalton eukaryotic pyruvate dehydrogenase complex from native cell extracts

The pyruvate dehydrogenase complex (PDHc) is a giant enzymatic assembly involved in pyruvate oxidation. PDHc components have been characterized in isolation, but the complex's quaternary structure has remained elusive due to sheer size, heterogeneity, and plasticity. Here, we identify fully assembled Chaetomium thermophilum α-keto acid dehydrogenase complexes in native cell extracts and characterize their domain arrangements utilizing mass spectrometry, activity assays, crosslinking, electron microscopy (EM), and computational modeling. We report the cryo-EM structure of the PDHc core and observe unique features of the previously unknown native state. The asymmetric reconstruction of the 10-MDa PDHc resolves spatial proximity of its components, agrees with stoichiometric data (60 E2p:12 E3BP:∼20 E1p: ≤ 12 E3), and proposes a minimum reaction path among component enzymes. The PDHc shows the presence of a dynamic pyruvate oxidation compartment, organized by core and peripheral protein species. Our data provide a framework for further understanding PDHc and α-keto acid dehydrogenase complex structure and function.

Algorithmic robustness to preferred orientations in single particle analysis by CryoEM

The presence of preferred orientations in single particle analysis (SPA) by cryo-Electron Microscopy (cryoEM) is currently one of the hurdles preventing many structural analyses from yielding high-resolution structures. Although the existence of preferred orientations is mostly related to the grid preparation, in this technical note, we show that some image processing algorithms used for angular assignment and three-dimensional (3D) reconstruction are more robust than others to these detrimental conditions. We exemplify this argument with three different data sets in which the presence of preferred orientations hindered achieving a 3D reconstruction without artifacts or, even worse, a 3D reconstruction could never be achieved.

Insights on cross-species transmission of SARS-CoV-2 from structural modeling

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the ongoing global pandemic that has infected more than 31 million people in more than 180 countries worldwide. Like other coronaviruses, SARS-CoV-2 is thought to have been transmitted to humans from wild animals. Given the scale and widespread geographical distribution of the current pandemic and confirmed cases of cross-species transmission, the question of the extent to which this transmission is possible emerges, as well as what molecular features distinguish susceptible from non-susceptible animal species. Here, we investigated the structural properties of several ACE2 orthologs bound to the SARS-CoV-2 spike protein. We found that species known not to be susceptible to SARS-CoV-2 infection have non-conservative mutations in several ACE2 amino acid residues that disrupt key polar and charged contacts with the viral spike protein. Our models also allow us to predict affinity-enhancing mutations that could be used to design ACE2 variants for therapeutic purposes. Finally, our study provides a blueprint for modeling viral-host protein interactions and highlights several important considerations when designing these computational studies and analyzing their results.

Structural models of human ACE2 variants with SARS-CoV-2 Spike protein for structure-based drug design

Emergence of coronaviruses poses a threat to global health and economy. The current outbreak of SARS-CoV-2 has infected more than 28,000,000 people and killed more than 915,000. To date, there is no treatment for coronavirus infections, making the development of therapies to prevent future epidemics of paramount importance. To this end, we collected information regarding naturally-occurring variants of the Angiotensin-converting enzyme 2 (ACE2), an epithelial receptor that both SARS-CoV and SARS-CoV-2 use to enter the host cells. We built 242 structural models of variants of human ACE2 bound to the receptor binding domain (RBD) of the SARS-CoV-2 surface spike glycoprotein (S protein) and refined their interfaces with HADDOCK. Our dataset includes 140 variants of human ACE2 representing missense mutations found in genome-wide studies, 39 mutants with reported effects on the recognition of the RBD, and 63 predictions after computational alanine scanning mutagenesis of ACE2-RBD interface residues. This dataset will help accelerate the design of therapeutics against SARS-CoV-2, as well as contribute to prevention of possible future coronaviruses outbreaks.

Measurement(s)	Molecular Genetic Variation
Technology Type(s)	digital curation
Factor Type(s)	ACE2 variants
Sample Characteristic - Organism	Homo sapiens

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12902498

Unstructured regions of large enzymatic complexes control the availability of metabolites with signaling functions

Metabolites produced via traditional biochemical processes affect intracellular communication, inflammation, and malignancy. Unexpectedly, acetyl-CoA, α-ketoglutarate and palmitic acid, which are chemical species of reactions catalyzed by highly abundant, gigantic enzymatic complexes, dubbed as "metabolons", have broad "nonmetabolic" signaling functions. Conserved unstructured regions within metabolons determine the yield of these metabolites. Unstructured regions tether functional protein domains, act as spatial constraints to confine constituent enzyme communication, and, in the case of acetyl-CoA production, tend to be regulated by intricate phosphorylation patterns. This review presents the multifaceted roles of these three significant metabolites and describes how their perturbation leads to altered or transformed cellular function. Their dedicated enzymatic systems are then introduced, namely, the pyruvate dehydrogenase (PDH) and oxoglutarate dehydrogenase (OGDH) complexes, and the fatty acid synthase (FAS), with a particular focus on their structural characterization and the localization of unstructured regions. Finally, upstream metabolite regulation, in which spatial occupancy of unstructured regions within dedicated metabolons may affect metabolite availability and subsequently alter cell functions, is discussed. Video abstract.

Structural analysis of 70S ribosomes by cross-linking/mass spectrometry reveals conformational plasticity

The ribosome is not only a highly complex molecular machine that translates the genetic information into proteins, but also an exceptional specimen for testing and optimizing cross-linking/mass spectrometry (XL-MS) workflows. Due to its high abundance, ribosomal proteins are frequently identified in proteome-wide XL-MS studies of cells or cell extracts. Here, we performed in-depth cross-linking of the E. coli ribosome using the amine-reactive cross-linker disuccinimidyl diacetic urea (DSAU). We analyzed 143 E. coli ribosomal structures, mapping a total of 10,771 intramolecular distances for 126 cross-link-pairs and 3,405 intermolecular distances for 97 protein pairs. Remarkably, 44% of intermolecular cross-links covered regions that have not been resolved in any high-resolution E. coli ribosome structure and point to a plasticity of cross-linked regions. We systematically characterized all cross-links and discovered flexible regions, conformational changes, and stoichiometric variations in bound ribosomal proteins, and ultimately remodeled 2,057 residues (15,794 atoms) in total. Our working model explains more than 95% of all cross-links, resulting in an optimized E. coli ribosome structure based on the cross-linking data obtained. Our study might serve as benchmark for conducting biochemical experiments on newly modeled protein regions, guided by XL-MS. Data are available via ProteomeXchange with identifier PXD018935.

Constructing artificial respiratory chain in polymer compartments: Insights into the interplay between bo3 oxidase and the membrane

Cytochrome bo3 ubiquinol oxidase is a transmembrane protein, which oxidizes ubiquinone and reduces oxygen, while pumping protons. Apart from its combination with F1Fo-ATPase to assemble a minimal ATP regeneration module, the utility of the proton pump can be extended to other applications in the context of synthetic cells such as transport, signaling, and control of enzymatic reactions. In parallel, polymers have been speculated to be phospholipid mimics with respect to their ability to self-assemble in compartments with increased stability. However, their usability as interfaces for complex membrane proteins has remained questionable. In the present work, we optimized a fusion/electroformation approach to reconstitute bo3 oxidase in giant unilamellar vesicles made of PDMS-g-PEO and/or phosphatidylcholine (PC). This enabled optical access, while microfluidic trapping allowed for online analysis of individual vesicles. The tight polymer membranes and the inward oriented enzyme caused 1 pH unit difference in 30 min, with an initial rate of 0.35 pH·min-1 To understand the interplay in these composite systems, we studied the relevant mechanical and rheological membrane properties. Remarkably, the proton permeability of polymer/lipid hybrids decreased after protein insertion, while the latter also led to a 20% increase of the polymer diffusion coefficient in polymersomes. In addition, PDMS-g-PEO increased the activity lifetime and the resistance to free radicals. These advantageous properties may open diverse applications, ranging from cell-free biotechnology to biomedicine. Furthermore, the presented study serves as a comprehensive road map for studying the interactions between membrane proteins and synthetic membranes, which will be fundamental for the successful engineering of such hybrid systems.

Flexibility of intrinsically disordered degrons in AUX/IAA proteins reinforces auxin co-receptor assemblies

Cullin RING-type E3 ubiquitin ligases SCF^TIR1/AFB1-5 and their AUX/IAA targets perceive the phytohormone auxin. The F-box protein TIR1 binds a surface-exposed degron in AUX/IAAs promoting their ubiquitylation and rapid auxin-regulated proteasomal degradation. Here, by adopting biochemical, structural proteomics and in vivo approaches we unveil how flexibility in AUX/IAAs and regions in TIR1 affect their conformational ensemble allowing surface accessibility of degrons. We resolve TIR1·auxin·IAA7 and TIR1·auxin·IAA12 complex topology, and show that flexible intrinsically disordered regions (IDRs) in the degron’s vicinity, cooperatively position AUX/IAAs on TIR1. We identify essential residues at the TIR1 N- and C-termini, which provide non-native interaction interfaces with IDRs and the folded PB1 domain of AUX/IAAs. We thereby establish a role for IDRs in modulating auxin receptor assemblies. By securing AUX/IAAs on two opposite surfaces of TIR1, IDR diversity supports locally tailored positioning for targeted ubiquitylation, and might provide conformational flexibility for a multiplicity of functional states.

Auxin-mediated recruitment of AUX/IAAs by the F-box protein TIR1 prompts rapid AUX/IAA ubiquitylation and degradation. By resolving auxin receptor topology, the authors show that intrinsically disordered regions near the degrons of two Aux/IAA proteins reinforce complex assembly and position Aux/IAAs for ubiquitylation.

2.7 Å cryo-EM structure of vitrified M. musculus H-chain apoferritin from a compact 200 keV cryo-microscope

Here we present the structure of mouse H-chain apoferritin at 2.7 Å (FSC = 0.143) solved by single particle cryogenic electron microscopy (cryo-EM) using a 200 kV device, the Thermo Fisher Glacios®. This is a compact, two-lens illumination system with a constant power objective lens, without any energy filters or aberration correctors, often thought of as a "screening cryo-microscope". Coulomb potential maps reveal clear densities for main chain carbonyl oxygens, residue side chains (including alternative conformations) and bound solvent molecules. We used a quasi-crystallographic reciprocal space approach to fit model coordinates to the experimental cryo-EM map. We argue that the advantages offered by (a) the high electronic and mechanical stability of the microscope, (b) the high emission stability and low beam energy spread of the high brightness Field Emission Gun (X-FEG), (c) direct electron detection technology and (d) particle-based Contrast Transfer Function (CTF) refinement have contributed to achieving high resolution. Overall, we show that basic electron optical settings for automated cryo-electron microscopy imaging can be used to determine structures approaching atomic resolution.

Integrative biology of native cell extracts: a new era for structural characterization of life processes

Advances in electron microscopy have provided unprecedented access to the structural characterization of large, flexible and heterogeneous complexes. Until recently, cryo-electron microscopy (cryo-EM) has been applied to understand molecular organization in either highly purified, isolated biomolecules or in situ. An emerging field is developing, bridging the gap between the two approaches, and focuses on studying molecular organization in native cell extracts. This field has demonstrated its potential by resolving the structure of fungal fatty acid synthase (FAS) at 4.7 Å [Fourier shell correlation (FSC) = 0.143]; FAS was not only less than 50% enriched, but also retained higher-order binders, previously unknown. Although controversial in the sense that the lysis step might introduce artifacts, cell extracts preserve aspects of cellular function. In addition, cell extracts are accessible, besides cryo-EM, to modern proteomic methods, chemical cross-linking, network biology and biophysical modeling. We expect that automation in imaging cell extracts, along with the integration of molecular/cell biology approaches, will provide remarkable achievements in the study of closer-to-life biomolecular states of pronounced biotechnological and medical importance. Such steps will, eventually, bring us a step closer to the biophysical description of cellular processes in an integrative, holistic approach.

Capturing protein communities by structural proteomics in a thermophilic eukaryote

The arrangement of proteins into complexes is a key organizational principle for many cellular functions. Although the topology of many complexes has been systematically analyzed in isolation, their molecular sociology in situ remains elusive. Here, we show that crude cellular extracts of a eukaryotic thermophile, Chaetomium thermophilum, retain basic principles of cellular organization. Using a structural proteomics approach, we simultaneously characterized the abundance, interactions, and structure of a third of the C. thermophilum proteome within these extracts. We identified 27 distinct protein communities that include 108 interconnected complexes, which dynamically associate with each other and functionally benefit from being in close proximity in the cell. Furthermore, we investigated the structure of fatty acid synthase within these extracts by cryoEM and this revealed multiple, flexible states of the enzyme in adaptation to its association with other complexes, thus exemplifying the need for in situ studies. As the components of the captured protein communities are known—at both the protein and complex levels—this study constitutes another step forward toward a molecular understanding of subcellular organization.

Subnanometre-resolution structure of the doublet microtubule reveals new classes of microtubule-associated proteins

Cilia are ubiquitous, hair-like appendages found in eukaryotic cells that carry out functions of cell motility and sensory reception. Cilia contain an intriguing cytoskeletal structure, termed the axoneme that consists of nine doublet microtubules radially interlinked and longitudinally organized in multiple specific repeat units. Little is known, however, about how the axoneme allows cilia to be both actively bendable and sturdy or how it is assembled. To answer these questions, we used cryo-electron microscopy to structurally analyse several of the repeating units of the doublet at sub-nanometre resolution. This structural detail enables us to unambiguously assign α- and β-tubulins in the doublet microtubule lattice. Our study demonstrates the existence of an inner sheath composed of different kinds of microtubule inner proteins inside the doublet that likely stabilizes the structure and facilitates the specific building of the B-tubule.

Cilia are hair-like appendages involved in cell motility and sensory reception. Here, the authors report a high resolution cryo-EM structure of the microtubule doublet from motile cilia and identify microtubule inner proteins (MIPs) bound to the inner surface of the doublet that appear to stabilize its structure.

Sense and simplicity in HADDOCK scoring: Lessons from CASP-CAPRI round 1

Our information-driven docking approach HADDOCK is a consistent top predictor and scorer since the start of its participation in the CAPRI community-wide experiment. This sustained performance is due, in part, to its ability to integrate experimental data and/or bioinformatics information into the modelling process, and also to the overall robustness of the scoring function used to assess and rank the predictions. In the CASP-CAPRI Round 1 scoring experiment we successfully selected acceptable/medium quality models for 18/14 of the 25 targets - a top-ranking performance among all scorers. Considering that for only 20 targets acceptable models were generated by the community, our effective success rate reaches as high as 90% (18/20). This was achieved using the standard HADDOCK scoring function, which, thirteen years after its original publication, still consists of a simple linear combination of intermolecular van der Waals and Coulomb electrostatics energies and an empirically derived desolvation energy term. Despite its simplicity, this scoring function makes sense from a physico-chemical perspective, encoding key aspects of biomolecular recognition. In addition to its success in the scoring experiment, the HADDOCK server takes the first place in the server prediction category, with 16 successful predictions. Much like our scoring protocol, because of the limited time per target, the predictions relied mainly on either an ab initio center-of-mass and symmetry restrained protocol, or on a template-based approach whenever applicable. These results underline the success of our simple but sensible prediction and scoring scheme. Proteins 2017; 85:417-423. © 2016 Wiley Periodicals, Inc.

Prediction of homoprotein and heteroprotein complexes by protein docking and template-based modeling: A CASP-CAPRI experiment

We present the results for CAPRI Round 30, the first joint CASP-CAPRI experiment, which brought together experts from the protein structure prediction and protein-protein docking communities. The Round comprised 25 targets from amongst those submitted for the CASP11 prediction experiment of 2014. The targets included mostly homodimers, a few homotetramers, and two heterodimers, and comprised protein chains that could readily be modeled using templates from the Protein Data Bank. On average 24 CAPRI groups and 7 CASP groups submitted docking predictions for each target, and 12 CAPRI groups per target participated in the CAPRI scoring experiment. In total more than 9500 models were assessed against the 3D structures of the corresponding target complexes. Results show that the prediction of homodimer assemblies by homology modeling techniques and docking calculations is quite successful for targets featuring large enough subunit interfaces to represent stable associations. Targets with ambiguous or inaccurate oligomeric state assignments, often featuring crystal contact-sized interfaces, represented a confounding factor. For those, a much poorer prediction performance was achieved, while nonetheless often providing helpful clues on the correct oligomeric state of the protein. The prediction performance was very poor for genuine tetrameric targets, where the inaccuracy of the homology-built subunit models and the smaller pair-wise interfaces severely limited the ability to derive the correct assembly mode. Our analysis also shows that docking procedures tend to perform better than standard homology modeling techniques and that highly accurate models of the protein components are not always required to identify their association modes with acceptable accuracy. Proteins 2016; 84(Suppl 1):323-348. © 2016 Wiley Periodicals, Inc.

dMM-PBSA: A New HADDOCK Scoring Function for Protein-Peptide Docking

Molecular-docking programs coupled with suitable scoring functions are now established and very useful tools enabling computational chemists to rapidly screen large chemical databases and thereby to identify promising candidate compounds for further experimental processing. In a broader scenario, predicting binding affinity is one of the most critical and challenging components of computer-aided structure-based drug design. The development of a molecular docking scoring function which in principle could combine both features, namely ranking putative poses and predicting complex affinity, would be of paramount importance. Here, we systematically investigated the performance of the MM-PBSA approach, using two different Poisson-Boltzmann solvers (APBS and DelPhi), in the currently rising field of protein-peptide interactions (PPIs), identifying the correct binding conformations of 19 different protein-peptide complexes and predicting their binding free energies. First, we scored the decoy structures from HADDOCK calculation via the MM-PBSA approach in order to assess the capability of retrieving near-native poses in the best-scoring clusters and of evaluating the corresponding free energies of binding. MM-PBSA behaves well in finding the poses corresponding to the lowest binding free energy, however the built-in HADDOCK score shows a better performance. In order to improve the MM-PBSA-based scoring function, we dampened the MM-PBSA solvation and coulombic terms by 0.2, as proposed in the HADDOCK score and LIE approaches. The new dampened MM-PBSA (dMM-PBSA) outperforms the original MM-PBSA and ranks the decoys structures as the HADDOCK score does. Second, we found a good correlation between the dMM-PBSA and HADDOCK scores for the near-native clusters of each system and the experimental binding energies, respectively. Therefore, we propose a new scoring function, dMM-PBSA, to be used together with the built-in HADDOCK score in the context of protein-peptide docking simulations.

PRODIGY: a web server for predicting the binding affinity of protein-protein complexes

Gaining insights into the structural determinants of protein-protein interactions holds the key for a deeper understanding of biological functions, diseases and development of therapeutics. An important aspect of this is the ability to accurately predict the binding strength for a given protein-protein complex. Here we present PROtein binDIng enerGY prediction (PRODIGY), a web server to predict the binding affinity of protein-protein complexes from their 3D structure. The PRODIGY server implements our simple but highly effective predictive model based on intermolecular contacts and properties derived from non-interface surface.

AVAILABILITY AND IMPLEMENTATION

PRODIGY is freely available at: http://milou.science.uu.nl/services/PRODIGY CONTACT: a.m.j.j.bonvin@uu.nl, a.vangone@uu.nl.

New potential peptide therapeutics perturbing CK1δ/α-tubulin interaction

Members of the CK1 family are highly conserved serine/threonine specific kinases being expressed in all eukaryotes. They are involved in many cellular processes and therefore tightly regulated. A central mechanism to modulate CK1 activity is via interaction with cellular proteins. CK1δ interacts with α-/β-tubulin and is involved in the regulation of microtubule dynamics. Therefore, it is important to identify the structural elements responsible for the interaction between these proteins. Using a peptide library covering the human CK1δ amino acid sequence in SPR and ELISA analyses, we identified peptide 39 (P39), encompassing aa361-aa375 of CK1δ, as a prominent binding partner of α-tubulin. P39 decreases α-tubulin phosphorylation by CK1δ and reduces the thermodynamic stability of α-tubulin in fluorescence thermal shift assays. Furthermore, P39 induces an inhibition of mitotic progression and a disruption of cells entering mitosis in CV-1 cells. Taken together our data provide valuable information regarding the interaction of CK1δ and α-tubulin and a novel approach for the development of pharmacological tools to inhibit proliferation of cancer cells.

The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes

The prediction of the quaternary structure of biomolecular macromolecules is of paramount importance for fundamental understanding of cellular processes and drug design. In the era of integrative structural biology, one way of increasing the accuracy of modeling methods used to predict the structure of biomolecular complexes is to include as much experimental or predictive information as possible in the process. This has been at the core of our information-driven docking approach HADDOCK. We present here the updated version 2.2 of the HADDOCK portal, which offers new features such as support for mixed molecule types, additional experimental restraints and improved protocols, all of this in a user-friendly interface. With well over 6000 registered users and 108,000 jobs served, an increasing fraction of which on grid resources, we hope that this timely upgrade will help the community to solve important biological questions and further advance the field. The HADDOCK2.2 Web server is freely accessible to non-profit users at http://haddock.science.uu.nl/services/HADDOCK2.2.