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.

Information-Driven Antibody-Antigen Modelling with HADDOCK.

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.

Native or Non-Native Protein-Protein Docking Models? Molecular Dynamics to the Rescue.

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.

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.

Therapy of chronic hepatitis C in people who inject drugs: focus on adherence.

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.

Redox Cofactor Rotates during Its Stepwise Decarboxylation: Molecular Mechanism of Conversion of Coproheme to Heme b.

Coproheme decarboxylase (ChdC) catalyzes the last step in the heme biosynthesis pathway of monoderm bacteria with coproheme acting both as redox cofactor and substrate. Hydrogen peroxide mediates the stepwise decarboxylation of propionates 2 and 4 of coproheme. Here we present the crystal structures of coproheme-loaded ChdC from Listeria monocytogenes (LmChdC) and the three-propionate intermediate, for which the propionate at position 2 (p2) has been converted to a vinyl group and is rotated by 90° compared to the coproheme complex structure. Single, double, and triple mutants of LmChdC, in which H-bonding interactions to propionates 2, 4, 6, and 7 were eliminated, allowed us to obtain the assignment of the coproheme propionates by resonance Raman spectroscopy and to follow the H2O2-mediated conversion of coproheme to heme b. Substitution of H2O2 by chlorite allowed us to monitor compound I formation in the inactive Y147H variant which lacks the catalytically essential Y147. This residue was demonstrated to be oxidized during turnover by using the spin-trap 2-methyl-2-nitrosopropane. Based on these findings and the data derived from molecular dynamics simulations of cofactor structures in distinct poses, we propose a reaction mechanism for the stepwise decarboxylation of coproheme that includes a 90° rotation of the intermediate three-propionate redox cofactor.

Free-Energy Calculations for Bioisosteric Modifications of A3 Adenosine Receptor Antagonists.

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 A3 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.

Binding Modes and Metabolism of Caffeine.

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.

Saturation Mutagenesis by Efficient Free-Energy Calculation.

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.

Free energy calculations on the stability of the 14-3-3ζ protein.

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.

Correction: Native Phytoremediation Potential of Urtica dioica for Removal of PCBs and Heavy Metals Can Be Improved by Genetic Manipulations Using Constitutive CaMV 35S Promoter.

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

Serogroup and Clonal Characterization of Czech Invasive Neisseria meningitidis Strains Isolated from 1971 to 2015.

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.

Native Phytoremediation Potential of Urtica dioica for Removal of PCBs and Heavy Metals Can Be Improved by Genetic Manipulations Using Constitutive CaMV 35S Promoter.

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.

Chemistry and Molecular Dynamics Simulations of Heme b-HemQ and Coproheme-HemQ.

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.

'Click chemistry' synthesis of 1-(α-D-mannopyranosyl)-1,2,3-triazoles for inhibition of α-mannosidases.

Three new triazole conjugates derived from d-mannose were synthesized and assayed in in vitro assays to investigate their ability to inhibit α-mannosidase enzymes from the glycoside hydrolase (GH) families 38 and 47. The triazole conjugates were more selective for a GH47 α-mannosidase (Aspergillus saitoi α1,2-mannosidase), showing inhibition at the micromolar level (IC50 values of 50-250 µM), and less potent towards GH38 mannosidases (IC50 values in the range of 0.5-6 mM towards jack bean α-mannosidase or Drosophila melanogaster lysosomal and Golgi α-mannosidases). The highest selectivity ratio [IC50(GH38)/IC50(GH47)] of 100 was exhibited by the phenyltriazole conjugate. To understand structure-activity properties of synthesized compounds, 3-D complexes of inhibitors with α-mannosidases were built using molecular docking calculations.