The HADDOCK2.4 web server for integrative modeling of biomolecular complexes.

Interactions between macromolecules, such as proteins and nucleic acids, are essential for cellular functions. Experimental methods can fail to provide all the information required to fully model biomolecular complexes at atomic resolution, particularly for large and heterogeneous assemblies. Integrative computational approaches have, therefore, gained popularity, complementing traditional experimental methods in structural biology. Here, we introduce HADDOCK2.4, an integrative modeling platform, and its updated web interface ( https://wenmr.science.uu.nl/haddock2.4 ). The platform seamlessly integrates diverse experimental and theoretical data to generate high-quality models of macromolecular complexes. The user-friendly web server offers automated parameter settings, access to distributed computing resources, and pre- and post-processing steps that enhance the user experience. To present the web server's various interfaces and features, we demonstrate two different applications: (i) we predict the structure of an antibody-antigen complex by using NMR data for the antigen and knowledge of the hypervariable loops for the antibody, and (ii) we perform coarse-grained modeling of PRC1 with a nucleosome particle guided by mutagenesis and functional data. The described protocols require some basic familiarity with molecular modeling and the Linux command shell. This new version of our widely used HADDOCK web server allows structural biologists and non-experts to explore intricate macromolecular assemblies encompassing various molecule types.

ARCTIC-3D: automatic retrieval and clustering of interfaces in complexes from 3D structural information.

The formation of a stable complex between proteins lies at the core of a wide variety of biological processes and has been the focus of countless experiments. The huge amount of information contained in the protein structural interactome in the Protein Data Bank can now be used to characterise and classify the existing biological interfaces. We here introduce ARCTIC-3D, a fast and user-friendly data mining and clustering software to retrieve data and rationalise the interface information associated with the protein input data. We demonstrate its use by various examples ranging from showing the increased interaction complexity of eukaryotic proteins, 20% of which on average have more than 3 different interfaces compared to only 10% for prokaryotes, to associating different functions to different interfaces. In the context of modelling biomolecular assemblies, we introduce the concept of "recognition entropy", related to the number of possible interfaces of the components of a protein-protein complex, which we demonstrate to correlate with the modelling difficulty in classical docking approaches. The identified interface clusters can also be used to generate various combinations of interface-specific restraints for integrative modelling. The ARCTIC-3D software is freely available at github.com/haddocking/arctic3d and can be accessed as a web-service at wenmr.science.uu.nl/arctic3d.

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.

An antibiotic from an uncultured bacterium binds to an immutable target.

Antimicrobial resistance is a leading mortality factor worldwide. Here, we report the discovery of clovibactin, an antibiotic isolated from uncultured soil bacteria. Clovibactin efficiently kills drug-resistant Gram-positive bacterial pathogens without detectable resistance. Using biochemical assays, solid-state nuclear magnetic resonance, and atomic force microscopy, we dissect its mode of action. Clovibactin blocks cell wall synthesis by targeting pyrophosphate of multiple essential peptidoglycan precursors (C55PP, lipid II, and lipid IIIWTA). Clovibactin uses an unusual hydrophobic interface to tightly wrap around pyrophosphate but bypasses the variable structural elements of precursors, accounting for the lack of resistance. Selective and efficient target binding is achieved by the sequestration of precursors into supramolecular fibrils that only form on bacterial membranes that contain lipid-anchored pyrophosphate groups. This potent antibiotic holds the promise of enabling the design of improved therapeutics that kill bacterial pathogens without resistance development.

A new antibiotic from an uncultured bacterium binds to an immutable target.

Antimicrobial resistance is a leading mortality factor worldwide. Here we report the discovery of clovibactin, a new antibiotic, isolated from uncultured soil bacteria. Clovibactin efficiently kills drug-resistant bacterial pathogens without detectable resistance. Using biochemical assays, solid-state NMR, and atomic force microscopy, we dissect its mode of action. Clovibactin blocks cell wall synthesis by targeting pyrophosphate of multiple essential peptidoglycan precursors (C 55 PP, Lipid II, Lipid WTA ). Clovibactin uses an unusual hydrophobic interface to tightly wrap around pyrophosphate, but bypasses the variable structural elements of precursors, accounting for the lack of resistance. Selective and efficient target binding is achieved by the irreversible sequestration of precursors into supramolecular fibrils that only form on bacterial membranes that contain lipid-anchored pyrophosphate groups. Uncultured bacteria offer a rich reservoir of antibiotics with new mechanisms of action that could replenish the antimicrobial discovery pipeline.

Structural Biology in the Clouds: The WeNMR-EOSC Ecosystem.

Structural biology aims at characterizing the structural and dynamic properties of biological macromolecules at atomic details. Gaining insight into three dimensional structures of biomolecules and their interactions is critical for understanding the vast majority of cellular processes, with direct applications in health and food sciences. Since 2010, the WeNMR project (www.wenmr.eu) has implemented numerous web-based services to facilitate the use of advanced computational tools by researchers in the field, using the high throughput computing infrastructure provided by EGI. These services have been further developed in subsequent initiatives under H2020 projects and are now operating as Thematic Services in the European Open Science Cloud portal (www.eosc-portal.eu), sending >12 millions of jobs and using around 4,000 CPU-years per year. Here we review 10 years of successful e-infrastructure solutions serving a large worldwide community of over 23,000 users to date, providing them with user-friendly, web-based solutions that run complex workflows in structural biology. The current set of active WeNMR portals are described, together with the complex backend machinery that allows distributed computing resources to be harvested efficiently.

Prediction of protein assemblies, the next frontier: The CASP14-CAPRI experiment.

We present the results for CAPRI Round 50, the fourth joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of twelve targets, including six dimers, three trimers, and three higher-order oligomers. Four of these were easy targets, for which good structural templates were available either for the full assembly, or for the main interfaces (of the higher-order oligomers). Eight were difficult targets for which only distantly related templates were found for the individual subunits. Twenty-five CAPRI groups including eight automatic servers submitted ~1250 models per target. Twenty groups including six servers participated in the CAPRI scoring challenge submitted ~190 models per target. The accuracy of the predicted models was evaluated using the classical CAPRI criteria. The prediction performance was measured by a weighted scoring scheme that takes into account the number of models of acceptable quality or higher submitted by each group as part of their five top-ranking models. Compared to the previous CASP-CAPRI challenge, top performing groups submitted such models for a larger fraction (70-75%) of the targets in this Round, but fewer of these models were of high accuracy. Scorer groups achieved stronger performance with more groups submitting correct models for 70-80% of the targets or achieving high accuracy predictions. Servers performed less well in general, except for the MDOCKPP and LZERD servers, who performed on par with human groups. In addition to these results, major advances in methodology are discussed, providing an informative overview of where the prediction of protein assemblies currently stands.

Information-driven modeling of biomolecular complexes.

Proteins play crucial roles in every cellular process by interacting with each other, nucleic acids, metabolites, and other molecules. The resulting assemblies can be very large and intricate and pose challenges to experimental methods. In the current era of integrative modeling, it is often only by a combination of various experimental techniques and computations that three-dimensional models of those molecular machines can be obtained. Among the various computational approaches available, molecular docking is often the method of choice when it comes to predicting three-dimensional structures of complexes. Docking can generate particularly accurate models when taking into account the available information on the complex of interest. We review here the use of experimental and bioinformatics data in protein-protein docking, describing recent software developments and highlighting applications for the modeling of antibody-antigen complexes and membrane protein complexes, and the use of evolutionary and shape information.

Nucleosome binding peptide presents laudable biophysical and in vivo effects.

Chromatin state is highly dependent on the nucleosome binding proteins. Herein, we used a multipronged approach employing biophysical and in vivo experiments to characterize the effects of Nucleosome Binding Peptides (NBPeps) on nucleosome and cell activity. We performed a series of structure-based calculations on the nucleosome surface interaction with GMIP1 (a novel NBPep generated in silico), and HMGN2 (nucleosome binding motif of HMGN2), which contains sites that bind DNA and the acid patch, and also LANA and H4pep (nucleosome binding motif of H4 histone tail) that only bind to the acidic patch. Biochemical assays shows that H4pep, but not HMGN2, GMIP1 and LANA, is highly specific for targeting the nucleosome, with important effects on the final nucleosome structure and robust in vivo effects. These findings suggest that NBPeps might have important therapeutic implications and relevance as tools for chromatin investigation.

An overview of data-driven HADDOCK strategies in CAPRI rounds 38-45.

Our information-driven docking approach HADDOCK has demonstrated a sustained performance since the start of its participation to CAPRI. This is due, in part, to its ability to integrate data into the modeling process, and to the robustness of its scoring function. We participated in CAPRI both as server and manual predictors. In CAPRI rounds 38-45, we have used various strategies depending on the available information. These ranged from imposing restraints to a few residues identified from literature as being important for the interaction, to binding pockets identified from homologous complexes or template-based refinement/CA-CA restraint-guided docking from identified templates. When relevant, symmetry restraints were used to limit the conformational sampling. We also tested for a large decamer target a new implementation of the MARTINI coarse-grained force field in HADDOCK. Overall, we obtained acceptable or better predictions for 13 and 11 server and manual submissions, respectively, out of the 22 interfaces. Our server performance (acceptable or higher-quality models when considering the top 10) was better (59%) than the manual (50%) one, in which we typically experiment with various combinations of protocols and data sources. Again, our simple scoring function based on a linear combination of intermolecular van der Waals and electrostatic energies and an empirical desolvation term demonstrated a good performance in the scoring experiment with a 63% success rate across all 22 interfaces. An analysis of model quality indicates that, while we are consistently performing well in generating acceptable models, there is room for improvement for generating/identifying higher quality models.

MARTINI-Based Protein-DNA Coarse-Grained HADDOCKing.

Modeling biomolecular assemblies is an important field in computational structural biology. The inherent complexity of their energy landscape and the computational cost associated with modeling large and complex assemblies are major drawbacks for integrative modeling approaches. The so-called coarse-graining approaches, which reduce the degrees of freedom of the system by grouping several atoms into larger "pseudo-atoms," have been shown to alleviate some of those limitations, facilitating the identification of the global energy minima assumed to correspond to the native state of the complex, while making the calculations more efficient. Here, we describe and assess the implementation of the MARTINI force field for DNA into HADDOCK, our integrative modeling platform. We combine it with our previous implementation for protein-protein coarse-grained docking, enabling coarse-grained modeling of protein-nucleic acid complexes. The system is modeled using MARTINI topologies and interaction parameters during the rigid body docking and semi-flexible refinement stages of HADDOCK, and the resulting models are then converted back to atomistic resolution by an atom-to-bead distance restraints-guided protocol. We first demonstrate the performance of this protocol using 44 complexes from the protein-DNA docking benchmark, which shows an overall ~6-fold speed increase and maintains similar accuracy as compared to standard atomistic calculations. As a proof of concept, we then model the interaction between the PRC1 and the nucleosome (a former CAPRI target in round 31), using the same information available at the time the target was offered, and compare all-atom and coarse-grained models.

Less Is More: Coarse-Grained Integrative Modeling of Large Biomolecular Assemblies with HADDOCK.

Predicting the 3D structure of protein interactions remains a challenge in the field of computational structural biology. This is in part due to difficulties in sampling the complex energy landscape of multiple interacting flexible polypeptide chains. Coarse-graining approaches, which reduce the number of degrees of freedom of the system, help address this limitation by smoothing the energy landscape, allowing an easier identification of the global energy minimum. They also accelerate the calculations, allowing for modeling larger assemblies. Here, we present the implementation of the MARTINI coarse-grained force field for proteins into HADDOCK, our integrative modeling platform. Docking and refinement are performed at the coarse-grained level, and the resulting models are then converted back to atomistic resolution through a distance restraints-guided morphing procedure. Our protocol, tested on the largest complexes of the protein docking benchmark 5, shows an overall ∼7-fold speed increase compared to standard all-atom calculations, while maintaining a similar accuracy and yielding substantially more near-native solutions. To showcase the potential of our method, we performed simultaneous 7 body docking to model the 1:6 KaiC-KaiB complex, integrating mutagenesis and hydrogen/deuterium exchange data from mass spectrometry with symmetry restraints, and validated the resulting models against a recently published cryo-EM structure.

Glutaminase Affects the Transcriptional Activity of Peroxisome Proliferator-Activated Receptor γ (PPARγ) via Direct Interaction.

Phosphate-activated glutaminases catalyze the deamidation of glutamine to glutamate and play key roles in several physiological and pathological processes. In humans, GLS encodes two multidomain splicing isoforms: KGA and GAC. In both isoforms, the canonical glutaminase domain is flanked by an N-terminal region that is folded into an EF-hand-like four-helix bundle. However, the splicing event replaces a well-structured three-repeat ankyrin domain in KGA with a shorter, unordered C-terminal stretch in GAC. The multidomain architecture, which contains putative protein-protein binding motifs, has led to speculation that glutaminases are involved in cellular processes other than glutamine metabolism; in fact, some proteins have been identified as binding partners of KGA and the isoforms of its paralogue gene, GLS2. Here, a yeast two-hybrid assay identified nuclear receptor peroxisome proliferator-activated receptor γ (PPARγ) as a new binding partner of the glutaminase. We show that KGA and GAC directly bind PPARγ with a low-micromolar dissociation constant; the interaction involves the N-terminal and catalytic domains of glutaminases as well as the ligand-binding domain of the nuclear receptor. The interaction occurs within the nucleus, and by sequestering PPARγ from its responsive element DR1, the glutaminases decreased nuclear receptor activity as assessed by a luciferase reporter assay. Altogether, our findings reveal an unexpected glutaminase-binding partner and, for the first time, directly link mitochondrial glutaminases to an unanticipated role in gene regulation.

The mechanism by which a distinguishing arabinofuranosidase can cope with internal di-substitutions in arabinoxylans.

BACKGROUND: Arabinoxylan is an abundant polysaccharide in industrially relevant biomasses such as sugarcane, corn stover and grasses. However, the arabinofuranosyl di-substitutions that decorate the xylan backbone are recalcitrant to most known arabinofuranosidases (Abfs). RESULTS: In this work, we identified a novel GH51 Abf (XacAbf51) that forms trimers in solution and can cope efficiently with both mono- and di-substitutions at terminal or internal xylopyranosyl units of arabinoxylan. Using mass spectrometry, the kinetic parameters of the hydrolysis of 33-α-l-arabinofuranosyl-xylotetraose and 23,33-di-α-l-arabinofuranosyl-xylotetraose by XacAbf51 were determined, demonstrating the capacity of this enzyme to cleave arabinofuranosyl linkages of internal mono- and di-substituted xylopyranosyl units. Complementation studies of fungal enzyme cocktails with XacAbf51 revealed an increase of up to 20% in the release of reducing sugars from pretreated sugarcane bagasse, showing the biotechnological potential of a generalist GH51 in biomass saccharification. To elucidate the structural basis for the recognition of internal di-substitutions, the crystal structure of XacAbf51 was determined unveiling the existence of a pocket strategically arranged near to the - 1 subsite that can accommodate a second arabinofuranosyl decoration, a feature not described for any other GH51 Abf structurally characterized so far. CONCLUSIONS: In summary, this study reports the first kinetic characterization of internal di-substitution release by a GH51 Abf, provides the structural basis for this activity and reveals a promising candidate for industrial processes involving plant cell wall depolymerization.

Structural basis of exo-ß-mannanase activity in the GH2 family.

The classical microbial strategy for depolymerization of ß-mannan polysaccharides involves the synergistic action of at least two enzymes, endo-1,4-ß-mannanases and ß-mannosidases. In this work, we describe the first exo-ß-mannanase from the GH2 family, isolated from Xanthomonas axonopodis pv. citri (XacMan2A), which can efficiently hydrolyze both manno-oligosaccharides and ß-mannan into mannose. It represents a valuable process simplification in the microbial carbon uptake that could be of potential industrial interest. Biochemical assays revealed a progressive increase in the hydrolysis rates from mannobiose to mannohexaose, which distinguishes XacMan2A from the known GH2 ß-mannosidases. Crystallographic analysis indicates that the active-site topology of XacMan2A underwent profound structural changes at the positive-subsite region, by the removal of the physical barrier canonically observed in GH2 ß-mannosidases, generating a more open and accessible active site with additional productive positive subsites. Besides that, XacMan2A contains two residue substitutions in relation to typical GH2 ß-mannosidases, Gly439 and Gly556, which alter the active site volume and are essential to its mode of action. Interestingly, the only other mechanistically characterized mannose-releasing exo-ß-mannanase so far is from the GH5 family, and its mode of action was attributed to the emergence of a blocking loop at the negative-subsite region of a cleft-like active site, whereas in XacMan2A, the same activity can be explained by the removal of steric barriers at the positive-subsite region in an originally pocket-like active site. Therefore, the GH2 exo-ß-mannanase represents a distinct molecular route to this rare activity, expanding our knowledge about functional convergence mechanisms in carbohydrate-active enzymes.

Thioredoxin-1 Negatively Modulates ADAM17 Activity Through Direct Binding and Indirect Reductive Activity.

AIMS: A disintegrin and metalloprotease 17 (ADAM17) modulates signaling events by releasing surface protein ectodomains such as TNFa and the EGFR-ligands. We have previously characterized cytoplasmic thioredoxin-1 (Trx-1) as a partner of ADAM17 cytoplasmic domain. Still, the mechanism of ADAM17 regulation by Trx-1 is unknown, and it has become of paramount importance to assess the degree of influence that Trx-1 has on metalloproteinase ADAM17. RESULTS: Combining discovery and targeted proteomic approaches, we uncovered that Trx-1 negatively regulates ADAM17 by direct and indirect effect. We performed cell-based assays with synthetic peptides and site-directed mutagenesis, and we demonstrated that the interaction interface of Trx-1 and ADAM17 is important for the negative regulation of ADAM17 activity. However, both Trx-1K72A and catalytic site mutant Trx-1C32/35S rescued ADAM17 activity, although the interaction with Trx-1C32/35S was unaffected, suggesting an indirect effect of Trx-1. We confirmed that the Trx-1C32/35S mutant showed diminished reductive capacity, explaining this indirect effect on increasing ADAM17 activity through oxidant levels. Interestingly, Trx-1K72A mutant showed similar oxidant levels to Trx-1C32/35S, even though its catalytic site was preserved. We further demonstrated that the general reactive oxygen species inhibitor, Nacetylcysteine (NAC), maintained the regulation of ADAM17 dependent of Trx-1 reductase activity levels; whereas the electron transport chain modulator, rotenone, abolished Trx-1 effect on ADAM17 activity. INNOVATION: We show for the first time that the mechanism of ADAM17 regulation, Trx-1 dependent, can be by direct interaction and indirect effect, bringing new insights into the cross-talk between isomerases and mammalian metalloproteinases. CONCLUSION: This unexpected Trx-1K72A behavior was due to more dimer formation and, consequently, the reduction of its Trx-1 reductase activity, evaluated through dimer verification, by gel filtration and mass spectrometry analysis. Antioxid. Redox Signal. 29, 717-734.

Pyrrole-indolinone SU11652 targets the nucleoside diphosphate kinase from Leishmania parasites.

Nucleoside diphosphate kinases (NDKs) are key enzymes in the purine-salvage pathway of trypanosomatids and have been associated with the maintenance of host-cell integrity for the benefit of the parasite, being potential targets for rational drug discovery and design. The NDK from Leishmania major (LmNDK) and mutants were expressed and purified to homogeneity. Thermal shift assays were employed to identify potential inhibitors for LmNDK. Calorimetric experiments, site-directed mutagenesis and molecular docking analysis were performed to validate the interaction and to evaluate the structural basis of ligand recognition. Furthermore, the anti-leishmanial activity of the newly identified and validated compound was tested in vitro against different Leishmania species. The molecule SU11652, a Sunitinib analog, was identified as a potential inhibitor for LmNDK and structural studies indicated that this molecule binds to the active site of LmNDK in a similar conformation to nucleotides, mimicking natural substrates. Isothermal titration calorimetry experiments combined with site-directed mutagenesis revealed that the residues H50 and H117, considered essential for catalysis, play an important role in ligand binding. In vitro cell studies showed that SU11652 had similar efficacy to Amphotericin b against some Leishmania species. Together, our results indicate the pyrrole-indolinone SU11652 as a promising scaffold for the rational design of new drugs targeting the enzyme NDK from Leishmania parasites.

Calcium and magnesium ions modulate the oligomeric state and function of mitochondrial 2-Cys peroxiredoxins in Leishmania parasites.

Leishmania parasites have evolved a number of strategies to cope with the harsh environmental changes during mammalian infection. One of these mechanisms involves the functional gain that allows mitochondrial 2-Cys peroxiredoxins to act as molecular chaperones when forming decamers. This function is critical for parasite infectivity in mammals, and its activation has been considered to be controlled exclusively by the enzyme redox state under physiological conditions. Herein, we have revealed that magnesium and calcium ions play a major role in modulating the ability of these enzymes to act as molecular chaperones, surpassing the redox effect. These ions are directly involved in mitochondrial metabolism and participate in a novel mechanism to stabilize the decameric form of 2-Cys peroxiredoxins in Leishmania mitochondria. Moreover, we have demonstrated that a constitutively dimeric Prx1m mutant impairs the survival of Leishmania under heat stress, supporting the central role of the chaperone function of Prx1m for Leishmania parasites during the transition from insect to mammalian hosts.

Negative regulation of bacterial killing and inflammation by two novel CD16 ligands.

Sepsis, a leading cause of death worldwide, involves exacerbated proinflammatory responses and inefficient bacterial clearance. Phagocytic cells play a crucial part in the prevention of sepsis by clearing bacteria through host innate receptors. Here, we used a phage display library to identify two peptides in Escherichia coli that interact with host innate receptors. One of these peptides, encoded by the wzxE gene of E. coli K-12, was involved in the transbilayer movement of a trisaccharide-lipid intermediate in the assembly of enterobacterial common antigen. Peptide-receptor interactions induced CD16-mediated inhibitory immunoreceptor tyrosine-based activating motif signaling, blocking the production of ROS and bacterial killing. This CD16-mediated inhibitory signaling was abrogated in a WzxE(-/-) mutant of E. coli K-12, restoring the production of ROS and bacterial killing. Taken together, the two novel CD16 ligands identified negatively regulate bacterial killing and inflammation. Our findings may contribute toward the development of new immunotherapies for E. coli-mediated infectious diseases and inflammation.