Analysis of a large SARS-CoV-2 (Alpha) outbreak in a Catalan prison using conventional and genomic epidemiology.

Enforcing strict protocols that prevent transmission of airborne infections in prisons is challenging. We examine a large SARS-CoV-2 outbreak in a Catalan penitentiary center from February-March 2021, prior to vaccination deployment. The aim was to describe the evolution of the outbreak using classical and genomic epidemiology and the containment strategy applied. The outbreak was initially detected in one module but spread to four, infecting 7 staff members and 140 incarcerated individuals, 6 of whom were hospitalized (4.4%). Genomic analysis confirmed a single origin (B.1.1.7). Contact tracing identified transmission vectors between modules and prevented further viral spread. In future similar scenarios, the control strategy described here may help limiting transmission of airborne infections in correctional settings.

The multivalency of the glucocorticoid receptor ligand-binding domain explains its manifold physiological activities.

The glucocorticoid receptor (GR) is a ubiquitously expressed transcription factor that controls metabolic and homeostatic processes essential for life. Although numerous crystal structures of the GR ligand-binding domain (GR-LBD) have been reported, the functional oligomeric state of the full-length receptor, which is essential for its transcriptional activity, remains disputed. Here we present five new crystal structures of agonist-bound GR-LBD, along with a thorough analysis of previous structural work. We identify four distinct homodimerization interfaces on the GR-LBD surface, which can associate into 20 topologically different homodimers. Biologically relevant homodimers were identified by studying a battery of GR point mutants including crosslinking assays in solution, quantitative fluorescence microscopy in living cells, and transcriptomic analyses. Our results highlight the relevance of non-canonical dimerization modes for GR, especially of contacts made by loop L1-3 residues such as Tyr545. Our work illustrates the unique flexibility of GR's LBD and suggests different dimeric conformations within cells. In addition, we unveil pathophysiologically relevant quaternary assemblies of the receptor with important implications for glucocorticoid action and drug design.

Discriminating physiological from non-physiological interfaces in structures of protein complexes: A community-wide study.

Reliably scoring and ranking candidate models of protein complexes and assigning their oligomeric state from the structure of the crystal lattice represent outstanding challenges. A community-wide effort was launched to tackle these challenges. The latest resources on protein complexes and interfaces were exploited to derive a benchmark dataset consisting of 1677 homodimer protein crystal structures, including a balanced mix of physiological and non-physiological complexes. The non-physiological complexes in the benchmark were selected to bury a similar or larger interface area than their physiological counterparts, making it more difficult for scoring functions to differentiate between them. Next, 252 functions for scoring protein-protein interfaces previously developed by 13 groups were collected and evaluated for their ability to discriminate between physiological and non-physiological complexes. A simple consensus score generated using the best performing score of each of the 13 groups, and a cross-validated Random Forest (RF) classifier were created. Both approaches showed excellent performance, with an area under the Receiver Operating Characteristic (ROC) curve of 0.93 and 0.94, respectively, outperforming individual scores developed by different groups. Additionally, AlphaFold2 engines recalled the physiological dimers with significantly higher accuracy than the non-physiological set, lending support to the reliability of our benchmark dataset annotations. Optimizing the combined power of interface scoring functions and evaluating it on challenging benchmark datasets appears to be a promising strategy.

Impact of AlphaFold on structure prediction of protein complexes: The CASP15-CAPRI experiment.

We present the results for CAPRI Round 54, the 5th joint CASP-CAPRI protein assembly prediction challenge. The Round offered 37 targets, including 14 homodimers, 3 homo-trimers, 13 heterodimers including 3 antibody-antigen complexes, and 7 large assemblies. On average ~70 CASP and CAPRI predictor groups, including more than 20 automatics servers, submitted models for each target. A total of 21 941 models submitted by these groups and by 15 CAPRI scorer groups were evaluated using the CAPRI model quality measures and the DockQ score consolidating these measures. The prediction performance was quantified by a weighted score based on the number of models of acceptable quality or higher submitted by each group among their five best models. Results show substantial progress achieved across a significant fraction of the 60+ participating groups. High-quality models were produced for about 40% of the targets compared to 8% two years earlier. This remarkable improvement is due to the wide use of the AlphaFold2 and AlphaFold2-Multimer software and the confidence metrics they provide. Notably, expanded sampling of candidate solutions by manipulating these deep learning inference engines, enriching multiple sequence alignments, or integration of advanced modeling tools, enabled top performing groups to exceed the performance of a standard AlphaFold2-Multimer version used as a yard stick. This notwithstanding, performance remained poor for complexes with antibodies and nanobodies, where evolutionary relationships between the binding partners are lacking, and for complexes featuring conformational flexibility, clearly indicating that the prediction of protein complexes remains a challenging problem.

Structural basis for the dominant or recessive character of GLIALCAM mutations found in leukodystrophies.

Megalencephalic leukoencephalopathy with subcortical cysts (MLC) is a type of leukodystrophy characterized by white matter edema, and it is caused mainly by recessive mutations in MLC1 and GLIALCAM genes. These variants are called MLC1 and MLC2A with both types of patients sharing the same clinical phenotype. In addition, dominant mutations in GLIALCAM have also been identified in a subtype of MLC patients with a remitting phenotype. This variant has been named MLC2B. GLIALCAM encodes for an adhesion protein containing two immunoglobulin (Ig) domains and it is needed for MLC1 targeting to astrocyte-astrocyte junctions. Most mutations identified in GLIALCAM abolish GlialCAM targeting to junctions. However, it is unclear why some mutations behave as recessive or dominant. Here, we used a combination of biochemistry methods with a new developed anti-GlialCAM nanobody, double-mutants and cysteine cross-links experiments, together with computer docking, to create a structural model of GlialCAM homo-interactions. Using this model, we suggest that dominant mutations affect different GlialCAM-GlialCAM interacting surfaces in the first Ig domain, which can occur between GlialCAM molecules present in the same cell (cis) or present in neighbouring cells (trans). Our results provide a framework that can be used to understand the molecular basis of pathogenesis of all identified GLIALCAM mutations.

UEP: an open-source and fast classifier for predicting the impact of mutations in protein-protein complexes.

MOTIVATION: Single protein residue mutations may reshape the binding affinity of protein-protein interactions. Therefore, predicting its effects is of great interest in biotechnology and biomedicine. Unfortunately, the availability of experimental data on binding affinity changes upon mutation is limited, which hampers the development of new and more precise algorithms. Here, we propose UEP, a classifier for predicting beneficial and detrimental mutations in protein-protein complexes trained on interactome data. RESULTS: Regardless of the simplicity of the UEP algorithm, which is based on a simple three-body contact potential derived from interactome data, we report competitive results with the gold standard methods in this field with the advantage of being faster in terms of computational time. Moreover, we propose a consensus selection procedure by involving the combination of three predictors that showed higher classification accuracy in our benchmark: UEP, pyDock and EvoEF1/FoldX. Overall, we demonstrate that the analysis of interactome data allows predicting the impact of protein-protein mutations using UEP, a fast and reliable open-source code. AVAILABILITY AND IMPLEMENTATION: UEP algorithm can be found at: https://github.com/pepamengual/UEP. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

New Insights into the Evolution of the Electron Transfer from Cytochrome f to Photosystem I in the Green and Red Branches of Photosynthetic Eukaryotes.

In cyanobacteria and most green algae of the eukaryotic green lineage, the copper-protein plastocyanin (Pc) alternatively replaces the heme-protein cytochrome c6 (Cc6) as the soluble electron carrier from cytochrome f (Cf) to photosystem I (PSI). The functional and structural equivalence of 'green' Pc and Cc6 has been well established, representing an example of convergent evolution of two unrelated proteins. However, plants only produce Pc, despite having evolved from green algae. On the other hand, Cc6 is the only soluble donor available in most species of the red lineage of photosynthetic organisms, which includes, among others, red algae and diatoms. Interestingly, Pc genes have been identified in oceanic diatoms, probably acquired by horizontal gene transfer from green algae. However, the mechanisms that regulate the expression of a functional Pc in diatoms are still unclear. In the green eukaryotic lineage, the transfer of electrons from Cf to PSI has been characterized in depth. The conclusion is that in the green lineage, this process involves strong electrostatic interactions between partners, which ensure a high affinity and an efficient electron transfer (ET) at the cost of limiting the turnover of the process. In the red lineage, recent kinetic and structural modeling data suggest a different strategy, based on weaker electrostatic interactions between partners, with lower affinity and less efficient ET, but favoring instead the protein exchange and the turnover of the process. Finally, in diatoms the interaction of the acquired green-type Pc with both Cf and PSI may not yet be optimized.

pyDockEneRes: per-residue decomposition of protein-protein docking energy.

MOTIVATION: Protein-protein interactions are key to understand biological processes at the molecular level. As a complement to experimental characterization of protein interactions, computational docking methods have become useful tools for the structural and energetics modeling of protein-protein complexes. A key aspect of such algorithms is the use of scoring functions to evaluate the generated docking poses and try to identify the best models. When the scoring functions are based on energetic considerations, they can help not only to provide a reliable structural model for the complex, but also to describe energetic aspects of the interaction. This is the case of the scoring function used in pyDock, a combination of electrostatics, desolvation and van der Waals energy terms. Its correlation with experimental binding affinity values of protein-protein complexes was explored in the past, but the per-residue decomposition of the docking energy was never systematically analyzed. RESULTS: Here, we present pyDockEneRes (pyDock Energy per-Residue), a web server that provides pyDock docking energy partitioned at the residue level, giving a much more detailed description of the docking energy landscape. Additionally, pyDockEneRes computes the contribution to the docking energy of the side-chain atoms. This fast approach can be applied to characterize a complex structure in order to identify energetically relevant residues (hot-spots) and estimate binding affinity changes upon mutation to alanine. AVAILABILITY AND IMPLEMENTATION: The server does not require registration by the user and is freely accessible for academics at https://life.bsc.es/pid/pydockeneres. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Integrative modeling of protein-protein interactions with pyDock for the new docking challenges.

The seventh CAPRI edition imposed new challenges to the modeling of protein-protein complexes, such as multimeric oligomerization, protein-peptide, and protein-oligosaccharide interactions. Many of the proposed targets needed the efficient integration of rigid-body docking, template-based modeling, flexible optimization, multiparametric scoring, and experimental restraints. This was especially relevant for the multimolecular assemblies proposed in the CASP12-CAPRI37 and CASP13-CAPRI46 joint rounds, which were described and evaluated elsewhere. Focusing on the purely CAPRI targets of this edition (rounds 38-45), we have participated in all 17 assessed targets (considering heteromeric and homomeric interfaces in T125 as two separate targets) both as predictors and as scorers, by using integrative modeling based on our docking and scoring approaches: pyDock, IRaPPA, and LightDock. In the protein-protein and protein-peptide targets, we have also participated with our webserver (pyDockWeb). On these 17 CAPRI targets, we submitted acceptable models (or better) within our top 10 models for 10 targets as predictors, 13 targets as scorers, and 4 targets as servers. In summary, our participation in this CAPRI edition confirmed the capabilities of pyDock for the scoring of docking models, increasingly used within the context of integrative modeling of protein interactions and multimeric assemblies.

SKEMPI 2.0: an updated benchmark of changes in protein-protein binding energy, kinetics and thermodynamics upon mutation.

Motivation: Understanding the relationship between the sequence, structure, binding energy, binding kinetics and binding thermodynamics of protein-protein interactions is crucial to understanding cellular signaling, the assembly and regulation of molecular complexes, the mechanisms through which mutations lead to disease, and protein engineering. Results: We present SKEMPI 2.0, a major update to our database of binding free energy changes upon mutation for structurally resolved protein-protein interactions. This version now contains manually curated binding data for 7085 mutations, an increase of 133%, including changes in kinetics for 1844 mutations, enthalpy and entropy changes for 443 mutations, and 440 mutations, which abolish detectable binding. Availability and implementation: The database is available as supplementary data and at https://life.bsc.es/pid/skempi2/. Supplementary information: Supplementary data are available at Bioinformatics online.

Blind prediction of homo- and hetero-protein complexes: The CASP13-CAPRI experiment.

We present the results for CAPRI Round 46, the third joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of 20 targets including 14 homo-oligomers and 6 heterocomplexes. Eight of the homo-oligomer targets and one heterodimer comprised proteins that could be readily modeled using templates from the Protein Data Bank, often available for the full assembly. The remaining 11 targets comprised 5 homodimers, 3 heterodimers, and two higher-order assemblies. These were more difficult to model, as their prediction mainly involved "ab-initio" docking of subunit models derived from distantly related templates. A total of ~30 CAPRI groups, including 9 automatic servers, submitted on average ~2000 models per target. About 17 groups participated in the CAPRI scoring rounds, offered for most targets, submitting ~170 models per target. The prediction performance, measured by the fraction of models of acceptable quality or higher submitted across all predictors groups, was very good to excellent for the nine easy targets. Poorer performance was achieved by predictors for the 11 difficult targets, with medium and high quality models submitted for only 3 of these targets. A similar performance "gap" was displayed by scorer groups, highlighting yet again the unmet challenge of modeling the conformational changes of the protein components that occur upon binding or that must be accounted for in template-based modeling. Our analysis also indicates that residues in binding interfaces were less well predicted in this set of targets than in previous Rounds, providing useful insights for directions of future improvements.

LightDock: a new multi-scale approach to protein-protein docking.

Motivation: Computational prediction of protein-protein complex structure by docking can provide structural and mechanistic insights for protein interactions of biomedical interest. However, current methods struggle with difficult cases, such as those involving flexible proteins, low-affinity complexes or transient interactions. A major challenge is how to efficiently sample the structural and energetic landscape of the association at different resolution levels, given that each scoring function is often highly coupled to a specific type of search method. Thus, new methodologies capable of accommodating multi-scale conformational flexibility and scoring are strongly needed. Results: We describe here a new multi-scale protein-protein docking methodology, LightDock, capable of accommodating conformational flexibility and a variety of scoring functions at different resolution levels. Implicit use of normal modes during the search and atomic/coarse-grained combined scoring functions yielded improved predictive results with respect to state-of-the-art rigid-body docking, especially in flexible cases. Availability and implementation: The source code of the software and installation instructions are available for download at https://life.bsc.es/pid/lightdock/. Contact: juanf@bsc.es. Supplementary information: Supplementary data are available at Bioinformatics online.

Structural and Computational Characterization of Disease-Related Mutations Involved in Protein-Protein Interfaces.

One of the known potential effects of disease-causing amino acid substitutions in proteins is to modulate protein-protein interactions (PPIs). To interpret such variants at the molecular level and to obtain useful information for prediction purposes, it is important to determine whether they are located at protein-protein interfaces, which are composed of two main regions, core and rim, with different evolutionary conservation and physicochemical properties. Here we have performed a structural, energetics and computational analysis of interactions between proteins hosting mutations related to diseases detected in newborn screening. Interface residues were classified as core or rim, showing that the core residues contribute the most to the binding free energy of the PPI. Disease-causing variants are more likely to occur at the interface core region rather than at the interface rim (p < 0.0001). In contrast, neutral variants are more often found at the interface rim or at the non-interacting surface rather than at the interface core region. We also found that arginine, tryptophan, and tyrosine are over-represented among mutated residues leading to disease. These results can enhance our understanding of disease at molecular level and thus contribute towards personalized medicine by helping clinicians to provide adequate diagnosis and treatments.

IRaPPA: information retrieval based integration of biophysical models for protein assembly selection.

MOTIVATION: In order to function, proteins frequently bind to one another and form 3D assemblies. Knowledge of the atomic details of these structures helps our understanding of how proteins work together, how mutations can lead to disease, and facilitates the designing of drugs which prevent or mimic the interaction. RESULTS: Atomic modeling of protein-protein interactions requires the selection of near-native structures from a set of docked poses based on their calculable properties. By considering this as an information retrieval problem, we have adapted methods developed for Internet search ranking and electoral voting into IRaPPA, a pipeline integrating biophysical properties. The approach enhances the identification of near-native structures when applied to four docking methods, resulting in a near-native appearing in the top 10 solutions for up to 50% of complexes benchmarked, and up to 70% in the top 100. AVAILABILITY AND IMPLEMENTATION: IRaPPA has been implemented in the SwarmDock server ( http://bmm.crick.ac.uk/∼SwarmDock/ ), pyDock server ( http://life.bsc.es/pid/pydockrescoring/ ) and ZDOCK server ( http://zdock.umassmed.edu/ ), with code available on request. CONTACT: moal@ebi.ac.uk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Structural and energy determinants in protein-RNA docking.

Deciphering the structural and energetic determinants of protein-RNA interactions harbors the potential to understand key cell processes at molecular level, such as gene expression and regulation. With this purpose, computational methods like docking aim to complement current biophysical and structural biology efforts. However, the few reported docking algorithms for protein-RNA interactions show limited predictive success rates, mainly due to incomplete sampling of the conformational space of both the protein and the RNA molecules, as well as to the difficulties of the scoring function in identifying the correct docking models. Here, we have tested the predictive value of a variety of knowledge-based and energetic scoring functions on a recently published protein-RNA docking benchmark and developed a scoring function able to efficiently discriminate docking decoys. We first performed docking calculations with the bound conformation, which allowed us to analyze the problem in optimal conditions. We found that geometry-based terms and electrostatics were the most important scoring terms, while binding propensities and desolvation were much less relevant for the scoring of protein-RNA models. This is in contrast with what we observed for protein-protein docking. The results also showed an interesting dependence of the predictive rates on the flexibility of the protein molecule, which arises from the observed higher positive charge of flexible interfaces and provides hints for future development of more efficient protein-RNA docking methods.

A systematic analysis of scoring functions in rigid-body protein docking: The delicate balance between the predictive rate improvement and the risk of overtraining.

Protein-protein interactions play fundamental roles in biological processes including signaling, metabolism, and trafficking. While the structure of a protein complex reveals crucial details about the interaction, it is often difficult to acquire this information experimentally. As the number of interactions discovered increases faster than they can be characterized, protein-protein docking calculations may be able to reduce this disparity by providing models of the interacting proteins. Rigid-body docking is a widely used docking approach, and is often capable of generating a pool of models within which a near-native structure can be found. These models need to be scored in order to select the acceptable ones from the set of poses. Recently, more than 100 scoring functions from the CCharPPI server were evaluated for this task using decoy structures generated with SwarmDock. Here, we extend this analysis to identify the predictive success rates of the scoring functions on decoys from three rigid-body docking programs, ZDOCK, FTDock, and SDOCK, allowing us to assess the transferability of the functions. We also apply set-theoretic measure to test whether the scoring functions are capable of identifying near-native poses within different subsets of the benchmark. This information can provide guides for the use of the most efficient scoring function for each docking method, as well as instruct future scoring functions development efforts. Proteins 2017; 85:1287-1297. © 2017 Wiley Periodicals, Inc.

pyDock scoring for the new modeling challenges in docking: Protein-peptide, homo-multimers, and domain-domain interactions.

The sixth CAPRI edition included new modeling challenges, such as the prediction of protein-peptide complexes, and the modeling of homo-oligomers and domain-domain interactions as part of the first joint CASP-CAPRI experiment. Other non-standard targets included the prediction of interfacial water positions and the modeling of the interactions between proteins and nucleic acids. We have participated in all proposed targets of this CAPRI edition both as predictors and as scorers, with new protocols to efficiently use our docking and scoring scheme pyDock in a large variety of scenarios. In addition, we have participated for the first time in the servers section, with our recently developed webserver, pyDockWeb. Excluding the CASP-CAPRI cases, we submitted acceptable models (or better) for 7 out of the 18 evaluated targets as predictors, 4 out of the 11 targets as scorers, and 6 out of the 18 targets as servers. The overall success rates were below those in past CAPRI editions. This shows the challenging nature of this last edition, with many difficult targets for which no participant submitted a single acceptable model. Interestingly, we submitted acceptable models for 83% of the evaluated protein-peptide targets. As for the 25 cases of the CASP-CAPRI experiment, in which we used a larger variety of modeling techniques (template-based, symmetry restraints, literature information, etc.), we submitted acceptable models for 56% of the targets. In summary, this CAPRI edition showed that pyDock scheme can be efficiently adapted to the increasing variety of problems that the protein interactions field is currently facing. Proteins 2017; 85:487-496. © 2016 Wiley Periodicals, Inc.

pyDockSAXS: protein-protein complex structure by SAXS and computational docking.

Structural characterization of protein-protein interactions at molecular level is essential to understand biological processes and identify new therapeutic opportunities. However, atomic resolution structural techniques cannot keep pace with current advances in interactomics. Low-resolution structural techniques, such as small-angle X-ray scattering (SAXS), can be applied at larger scale, but they miss atomic details. For efficient application to protein-protein complexes, low-resolution information can be combined with theoretical methods that provide energetic description and atomic details of the interactions. Here we present the pyDockSAXS web server (http://life.bsc.es/pid/pydocksaxs) that provides an automatic pipeline for modeling the structure of a protein-protein complex from SAXS data. The method uses FTDOCK to generate rigid-body docking models that are subsequently evaluated by a combination of pyDock energy-based scoring function and their capacity to describe SAXS data. The only required input files are structural models for the interacting partners and a SAXS curve. The server automatically provides a series of structural models for the complex, sorted by the pyDockSAXS scoring function. The user can also upload a previously computed set of docking poses, which opens the possibility to filter the docking solutions by potential interface residues or symmetry restraints. The server is freely available to all users without restriction.

Structural basis for the recruitment and activation of the Legionella phospholipase VipD by the host GTPase Rab5.

A challenge for microbial pathogens is to assure that their translocated effector proteins target only the correct host cell compartment during infection. The Legionella pneumophila effector vacuolar protein sorting inhibitor protein D (VipD) localizes to early endosomal membranes and alters their lipid and protein composition, thereby protecting the pathogen from endosomal fusion. This process requires the phospholipase A1 (PLA1) activity of VipD that is triggered specifically on VipD binding to the host cell GTPase Rab5, a key regulator of endosomes. Here, we present the crystal structure of VipD in complex with constitutively active Rab5 and reveal the molecular mechanism underlying PLA1 activation. An active site-obstructing loop that originates from the C-terminal domain of VipD is repositioned on Rab5 binding, thereby exposing the catalytic pocket within the N-terminal PLA1 domain. Substitution of amino acid residues located within the VipD-Rab5 interface prevented Rab5 binding and PLA1 activation and caused a failure of VipD mutant proteins to target to Rab5-enriched endosomal structures within cells. Experimental and computational analyses confirmed an extended VipD-binding interface on Rab5, explaining why this L. pneumophila effector can compete with cellular ligands for Rab5 binding. Together, our data explain how the catalytic activity of a microbial effector can be precisely linked to its subcellular localization.

Structural bases for the interaction and stabilization of the human amino acid transporter LAT2 with its ancillary protein 4F2hc.

Heteromeric amino acid transporters (HATs) are the unique example, known in all kingdoms of life, of solute transporters composed of two subunits linked by a conserved disulfide bridge. In metazoans, the heavy subunit is responsible for the trafficking of the heterodimer to the plasma membrane, and the light subunit is the transporter. HATs are involved in human pathologies such as amino acidurias, tumor growth and invasion, viral infection and cocaine addiction. However structural information about interactions between the heavy and light subunits of HATs is scarce. In this work, transmission electron microscopy and single-particle analysis of purified human 4F2hc/L-type amino acid transporter 2 (LAT2) heterodimers overexpressed in the yeast Pichia pastoris, together with docking analysis and crosslinking experiments, reveal that the extracellular domain of 4F2hc interacts with LAT2, almost completely covering the extracellular face of the transporter. 4F2hc increases the stability of the light subunit LAT2 in detergent-solubilized Pichia membranes, allowing functional reconstitution of the heterodimer into proteoliposomes. Moreover, the extracellular domain of 4F2hc suffices to stabilize solubilized LAT2. The interaction of 4F2hc with LAT2 gives insights into the structural bases for light subunit recognition and the stabilizing role of the ancillary protein in HATs.