Search

The HADDOCK team participated in CAPRI rounds 47-55 as server, manual predictor, and scorers. Throughout these CAPRI rounds, we used a plethora of computational strategies to predict the structure of protein complexes. Of the 10 targets comprising 24 interfaces, we achieved acceptable or better models for 3 targets in the human category and 1 in the server category. Our performance in the scoring challenge was slightly better, with our simple scoring protocol being the only one capable of identifying an acceptable model for Target 234. This result highlights the robustness of the simple, fully physics-based HADDOCK scoring function, especially when applied to highly flexible antibody-antigen complexes. Inspired by the significant advances in machine learning for structural biology and the dramatic improvement in our success rates after the public release of Alphafold2, we identify the integration of classical approaches like HADDOCK with AI-driven structure prediction methods as a key strategy for improving the accuracy of model generation and scoring., (© 2024 The Author(s). PROTEINS: Structure, Function, and Bioinformatics published by Wiley Periodicals LLC.)

In the recent years, therapeutic use of antibodies has seen a huge growth, "due to their inherent proprieties and technological advances in the methods used to study and characterize them. Effective design and engineering of antibodies for therapeutic purposes are heavily dependent on knowledge of the structural principles that regulate antibody-antigen interactions. Several experimental techniques such as X-ray crystallography, cryo-electron microscopy, NMR, or mutagenesis analysis can be applied, but these are usually expensive and time-consuming. Therefore computational approaches like molecular docking may offer a valuable alternative for the characterization of antibody-antigen complexes.Here we describe a protocol for the prediction of the 3D structure of antibody-antigen complexes using the integrative modelling platform HADDOCK. The protocol consists of (1) the identification of the antibody residues belonging to the hypervariable loops which are known to be crucial for the binding and can be used to guide the docking and (2) the detailed steps to perform docking with the HADDOCK 2.4 webserver following different strategies depending on the availability of information about epitope residues., (© 2023. Springer Science+Business Media, LLC, part of Springer Nature.)

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., (© 2021 Wiley Periodicals LLC.)

Molecular docking excels at creating a plethora of potential models of protein-protein complexes. To correctly distinguish the favorable, native-like models from the remaining ones remains, however, a challenge. We assessed here if a protocol based on molecular dynamics (MD) simulations would allow distinguishing native from non-native models to complement scoring functions used in docking. To this end, the first models for 25 protein-protein complexes were generated using HADDOCK. Next, MD simulations complemented with machine learning were used to discriminate between native and non-native complexes based on a combination of metrics reporting on the stability of the initial models. Native models showed higher stability in almost all measured properties, including the key ones used for scoring in the Critical Assessment of PRedicted Interaction (CAPRI) competition, namely the positional root mean square deviations and fraction of native contacts from the initial docked model. A random forest classifier was trained, reaching a 0.85 accuracy in correctly distinguishing native from non-native complexes. Reasonably modest simulation lengths of the order of 50-100 ns are sufficient to reach this accuracy, which makes this approach applicable in practice.

Background: Intravenous drug use (IVDU) represents the major factor of HCV transmission, but the treatment uptake among people who inject drugs (PWID) remains low owing to a false presumption of low efficacy. The aim of our study was to assess treatment efficacy in PWID and factors determining adherence to therapy., Methods: A total of 278 consecutive patients starting DAA (direct-acting antivirals) therapy were included, divided into two groups: individuals with a history of IVDU, PWID group (N = 101) and the control group (N = 177) without a history of IVDU., Results: Sustained virological response 12 weeks after the end of therapy (SVR12) was achieved by 99/101 (98%) and 172/177 (98%) patients in the PWID and control group, respectively; in PWID group, two patients were lost to follow-up, and in the control group, four patients relapsed and one was lost to follow-up. PWID patients postponed appointments significantly more often, 29 (28.7%) in PWID versus 7 (4%) in the control group, p = 0.001. Thirteen of 101 (12.9%) and six of 177 (3.4%) patients in the PWID and in the control group, respectively, missed at least one visit (p < 0.01). However, postponing visits led to a lack of medication in only one PWID. In the PWID group, older age (p < 0.05; OR 1.07, 95% CI 1.00-1.20) and stable housing (p < 0.01; OR 9.70, 95% CI 2.10-56.20) were factors positively contributing to adherence. Contrarily, a stable job was a factor negatively influencing adherence (p < 0.05; OR 0.24, 95% CI 0.06-0.81). In the control group, none of the analyzed social and demographic factors had an impact on adherence to therapy., Conclusions: In PWID, treatment efficacy was excellent and was comparable with SVR of the control group. Stable housing and older age contributed to a better adherence to therapy.

Adenosine receptors are a family of G protein-coupled receptors with increased attention as drug targets on different indications. We investigate the thermodynamics of ligand binding to the A 3 adenosine receptor subtype, focusing on a recently reported series of diarylacetamidopyridine inhibitors via molecular dynamics simulations. With a combined approach of thermodynamic integration and one-step perturbation, we characterize the impact of the charge distribution in a central heteroaromatic ring on the binding affinity prediction. Standard charge distributions according to the GROMOS force field yield values in good agreement with the experimental data and previous free energy calculations. Subsequently, we examine the thermodynamics of inhibitor binding in terms of the energetic and entropic contributions. The highest entropy penalties are found for inhibitors with methoxy substituents in meta position of the aryl groups. This bulky group restricts rotation of aromatic rings attached to the pyrimidine core which leads to two distinct poses of the ligand. Our predictions support the previously proposed binding pose for the o-methoxy ligand, yielding in this case a very good correlation with the experimentally measured affinities with deviations below 4 kJ/mol.

A correct estimate of ligand binding modes and a ratio of their occupancies is crucial for calculations of binding free energies. The newly developed method BLUES combines molecular dynamics with nonequilibrium candidate Monte Carlo. Nonequilibrium candidate Monte Carlo generates a plethora of possible binding modes and molecular dynamics enables the system to relax. We used BLUES to investigate binding modes of caffeine in the active site of its metabolizing enzyme Cytochrome P450 1A2 with the aim of elucidating metabolite-formation profiles at different concentrations. Because the activation energies of all sites of metabolism do not show a clear preference for one metabolite over the others, the orientations in the active site must play a key role. In simulations with caffeine located in a spacious pocket above the I-helix, it points N3 and N1 to the heme iron, whereas in simulations where caffeine is in close proximity to the heme N7 and C8 are preferably oriented toward the heme iron. We propose a mechanism where at low caffeine concentrations caffeine binds to the upper part of the active site, leading to formation of the main metabolite paraxanthine. On the other hand, at high concentrations two molecules are located in the active site, forcing one molecule into close proximity to the heme and yielding metabolites theophylline and trimethyluretic acid. Our results offer an explanation of previously published experimental results.

Mutations of cysteine are often introduced to e.g. avoid formation of non-physiological inter-molecular disulfide bridges in in-vitro experiments, or to maintain specificity in labeling experiments. Alanine or serine is typically preferred, which usually do not alter the overall protein stability, when the original cysteine was surface exposed. However, selecting the optimal mutation for cysteines in the hydrophobic core of the protein is more challenging. In this work, the stability of selected Cys mutants of 14-3-3ζ was predicted by free-energy calculations and the obtained data were compared with experimentally determined stabilities. Both the computational predictions as well as the experimental validation point at a significant destabilization of mutants C94A and C94S. This destabilization could be attributed to the formation of hydrophobic cavities and a polar solvation of a hydrophilic side chain. A L12E, M78K double mutant was further studied in terms of its reduced dimerization propensity. In contrast to naïve expectations, this double mutant did not lead to the formation of strong salt bridges, which was rationalized in terms of a preferred solvation of the ionic species. Again, experiments agreed with the calculations by confirming the monomerization of the double mutants. Overall, the simulation data is in good agreement with experiments and offers additional insight into the stability and dimerization of this important family of regulatory proteins., (Copyright © 2017 The Author(s). Published by Elsevier B.V. All rights reserved.)

Single-point mutations in proteins can greatly influence protein stability, binding affinity, protein function or its expression per se. Here, we present accurate and efficient predictions of the free energy of mutation of amino acids. We divided the complete mutational free energy into an uncharging step, which we approximate by a third-power fitting (TPF) approach, and an annihilation step, which we approximate using the one-step perturbation (OSP) method. As a diverse set of test systems, we computed the solvation free energy of all amino acid side chain analogues and obtained an excellent agreement with thermodynamic integration (TI) data. Moreover, we calculated mutational free energies in model tripeptides and established an efficient protocol involving a single reference state. Again, the approximate methods agreed excellently with the TI references, with a root-mean-square error of only 3.6 kJ/mol over 17 mutations. Our combined TPF+OSP approach does show not only a very good agreement but also a 2-fold higher efficiency than full blown TI calculations.

[This corrects the article DOI: 10.1371/journal.pone.0167927.].

Background: This study presents antigenic and genetic characteristics of Neisseria meningitidis strains recovered from invasive meningococcal disease (IMD) in the Czech Republic in 1971-2015., Material and Methods: A total of 1970 isolates from IMD, referred to the National Reference Laboratory for Meningococcal Infections in 1971-2015, were studied. All isolates were identified and characterized by conventional biochemical and serological tests. Most isolates (82.5%) were characterized by multilocus sequence typing method., Results: In the study period 1971-2015, the leading serogroup was B (52.4%), most often assigned to clonal complexes cc32, cc41/44, cc18, and cc269. A significant percentage of strains were of serogroup C (41.4%), with high clonal homogeneity due to hyperinvasive complex cc11, which played an important role in IMD in the Czech Republic in the mid-1990s. Serogroup Y isolates, mostly assigned to cc23, and isolates of clonally homogeneous serogroup W have also been recovered more often over the last years., Conclusion: The incidence of IMD and distribution of serogroups and clonal complexes of N. meningitidis in the Czech Republic varied over time, as can be seen from the long-term monitoring, including molecular surveillance data. Data from the conventional and molecular IMD surveillance are helpful in refining the antimeningococcal vaccination strategy in the Czech Republic., Competing Interests: The authors have declared that no competing interests exist.

Although stinging nettle (Urtica dioica) has been shown to reduce HM (heavy metal) content in soil, its wider phytoremediation potential has been neglected. Urtica dioica was cultivated in soils contaminated with HMs or polychlorinated biphenyls (PCBs). After four months, up to 33% of the less chlorinated biphenyls and 8% of HMs (Zn, Pb, Cd) had been removed. Bacteria were isolated from the plant tissue, with the endophytic bacteria Bacillus shackletonii and Streptomyces badius shown to have the most significant effect. These bacteria demonstrated not only benefits for plant growth, but also extreme tolerance to As, Zn and Pb. Despite these results, the native phytoremediation potential of nettles could be improved by biotechnologies. Transient expression was used to investigate the functionality of the most common constitutive promoter, CaMV 35S in Urtica dioica. This showed the expression of the CUP and bphC transgenes. Collectively, our findings suggest that remediation by stinging nettle could have a much wider range of applications than previously thought., Competing Interests: The authors have declared that no competing interests exist.

Recently, a novel pathway for heme b biosynthesis in Gram-positive bacteria has been proposed. The final poorly understood step is catalyzed by an enzyme called HemQ and includes two decarboxylation reactions leading from coproheme to heme b. Coproheme has been suggested to act as both substrate and redox active cofactor in this reaction. In the study presented here, we focus on HemQs from Listeria monocytogenes (LmHemQ) and Staphylococcus aureus (SaHemQ) recombinantly produced as apoproteins in Escherichia coli. We demonstrate the rapid and two-phase uptake of coproheme by both apo forms and the significant differences in thermal stability of the apo forms, coproheme-HemQ and heme b-HemQ. Reduction of ferric high-spin coproheme-HemQ to the ferrous form is shown to be enthalpically favored but entropically disfavored with standard reduction potentials of -205 ± 3 mV for LmHemQ and -207 ± 3 mV for SaHemQ versus the standard hydrogen electrode at pH 7.0. Redox thermodynamics suggests the presence of a pronounced H-bonding network and restricted solvent mobility in the heme cavity. Binding of cyanide to the sixth coproheme position is monophasic but relatively slow (∼1 × 10(4) M(-1) s(-1)). On the basis of the available structures of apo-HemQ and modeling of both loaded forms, molecular dynamics simulation allowed analysis of the interaction of coproheme and heme b with the protein as well as the role of the flexibility at the proximal heme cavity and the substrate access channel for coproheme binding and heme b release. Obtained data are discussed with respect to the proposed function of HemQ in monoderm bacteria.