Combining Structural Aggregation Propensity and Stability Predictions To Redesign Protein Solubility.

The aggregation propensity of each particular protein seems to be shaped by evolution according to its natural abundance in the cell. The production and downstream processing of recombinant polypeptides implies attaining concentrations that are orders of magnitude above their natural levels, often resulting in their aggregation; a phenomenon that precludes the marketing of many globular proteins for biomedical or biotechnological applications. Therefore, there is a huge interest in methods aimed to increase the proteins solubility above their natural limits. Here, we demonstrate that an updated version of our AGGRESCAN 3D structural aggregation predictor, that now takes into account protein stability, allows for designing mutations at specific positions in the structure that improve the solubility of proteins without compromising their conformation. Using this approach, we have designed a highly soluble variant of the green fluorescent protein and a human single-domain VH antibody displaying significantly reduced aggregation propensity. Overall, our data indicate that the solubility of unrelated proteins can be easily tuned by in silico-designed nondestabilizing amino acid changes at their surfaces.

Ensemble-based evaluation for protein structure models.

MOTIVATION: Comparing protein tertiary structures is a fundamental procedure in structural biology and protein bioinformatics. Structure comparison is important particularly for evaluating computational protein structure models. Most of the model structure evaluation methods perform rigid body superimposition of a structure model to its crystal structure and measure the difference of the corresponding residue or atom positions between them. However, these methods neglect intrinsic flexibility of proteins by treating the native structure as a rigid molecule. Because different parts of proteins have different levels of flexibility, for example, exposed loop regions are usually more flexible than the core region of a protein structure, disagreement of a model to the native needs to be evaluated differently depending on the flexibility of residues in a protein. RESULTS: We propose a score named FlexScore for comparing protein structures that consider flexibility of each residue in the native state of proteins. Flexibility information may be extracted from experiments such as NMR or molecular dynamics simulation. FlexScore considers an ensemble of conformations of a protein described as a multivariate Gaussian distribution of atomic displacements and compares a query computational model with the ensemble. We compare FlexScore with other commonly used structure similarity scores over various examples. FlexScore agrees with experts' intuitive assessment of computational models and provides information of practical usefulness of models. AVAILABILITY AND IMPLEMENTATION: https://bitbucket.org/mjamroz/flexscore CONTACT: dkihara@purdue.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

CABS-dock web server for the flexible docking of peptides to proteins without prior knowledge of the binding site.

Protein-peptide interactions play a key role in cell functions. Their structural characterization, though challenging, is important for the discovery of new drugs. The CABS-dock web server provides an interface for modeling protein-peptide interactions using a highly efficient protocol for the flexible docking of peptides to proteins. While other docking algorithms require pre-defined localization of the binding site, CABS-dock does not require such knowledge. Given a protein receptor structure and a peptide sequence (and starting from random conformations and positions of the peptide), CABS-dock performs simulation search for the binding site allowing for full flexibility of the peptide and small fluctuations of the receptor backbone. This protocol was extensively tested over the largest dataset of non-redundant protein-peptide interactions available to date (including bound and unbound docking cases). For over 80% of bound and unbound dataset cases, we obtained models with high or medium accuracy (sufficient for practical applications). Additionally, as optional features, CABS-dock can exclude user-selected binding modes from docking search or to increase the level of flexibility for chosen receptor fragments. CABS-dock is freely available as a web server at http://biocomp.chem.uw.edu.pl/CABSdock.

AGGRESCAN3D (A3D): server for prediction of aggregation properties of protein structures.

Protein aggregation underlies an increasing number of disorders and constitutes a major bottleneck in the development of therapeutic proteins. Our present understanding on the molecular determinants of protein aggregation has crystalized in a series of predictive algorithms to identify aggregation-prone sites. A majority of these methods rely only on sequence. Therefore, they find difficulties to predict the aggregation properties of folded globular proteins, where aggregation-prone sites are often not contiguous in sequence or buried inside the native structure. The AGGRESCAN3D (A3D) server overcomes these limitations by taking into account the protein structure and the experimental aggregation propensity scale from the well-established AGGRESCAN method. Using the A3D server, the identified aggregation-prone residues can be virtually mutated to design variants with increased solubility, or to test the impact of pathogenic mutations. Additionally, A3D server enables to take into account the dynamic fluctuations of protein structure in solution, which may influence aggregation propensity. This is possible in A3D Dynamic Mode that exploits the CABS-flex approach for the fast simulations of flexibility of globular proteins. The A3D server can be accessed at http://biocomp.chem.uw.edu.pl/A3D/.

KnotProt: a database of proteins with knots and slipknots.

The protein topology database KnotProt, http://knotprot.cent.uw.edu.pl/, collects information about protein structures with open polypeptide chains forming knots or slipknots. The knotting complexity of the cataloged proteins is presented in the form of a matrix diagram that shows users the knot type of the entire polypeptide chain and of each of its subchains. The pattern visible in the matrix gives the knotting fingerprint of a given protein and permits users to determine, for example, the minimal length of the knotted regions (knot's core size) or the depth of a knot, i.e. how many amino acids can be removed from either end of the cataloged protein structure before converting it from a knot to a different type of knot. In addition, the database presents extensive information about the biological functions, families and fold types of proteins with non-trivial knotting. As an additional feature, the KnotProt database enables users to submit protein or polymer chains and generate their knotting fingerprints.

Coarse-Grained Simulations of Membrane Insertion and Folding of Small Helical Proteins Using the CABS Model.

The CABS coarse-grained model is a well-established tool for modeling globular proteins (predicting their structure, dynamics, and interactions). Here we introduce an extension of the CABS representation and force field (CABS-membrane) to the modeling of the effect of the biological membrane environment on the structure of membrane proteins. We validate the CABS-membrane model in folding simulations of 10 short helical membrane proteins not using any knowledge about their structure. The simulations start from random protein conformations placed outside the membrane environment and allow for full flexibility of the modeled proteins during their spontaneous insertion into the membrane. In the resulting trajectories, we have found models close to the experimental membrane structures. We also attempted to select the correctly folded models using simple filtering followed by structural clustering combined with reconstruction to the all-atom representation and all-atom scoring. The CABS-membrane model is a promising approach for further development toward modeling of large protein-membrane systems.

Institutional analysis of health promotion for older people in Europe - concept and research tool.

BACKGROUND: European societies are ageing rapidly and thus health promotion for older people (HP4OP) is becoming an increasingly relevant issue. Crucial here is not only the clinical aspect of health promotion but also its organisational and institutional dimension. The latter has been relatively neglected in research on HP4OP. This issue is addressed in this study, constituting a part of the EU project ProHealth65+, engaging ten member countries. This paper is based on two intertwining research goals: (1) exploring which institutions/organisations are performing HP4OP activities in selected European countries (including sectors involved, performed roles of these institutions, organisation of those activities); (2) developing an institutional approach to HP4OP. Thus, the paper provides a description of the analytical tools for further research in this area. METHODS: The mentioned aims were addressed through the mutual use of two complementary methods: (a) a literature review of scientific and grey literature; and (b) questionnaire survey with selected expert respondents from 10 European countries. The expert respondents, selected by the project's collaborating partners, were asked to fill in a custom designed questionnaire concerning HP4OP institutional aspects. RESULTS: The literature review provided an overview of the organisational arrangements in different HP4OP initiatives. It also enabled the development of functional institutional definitions of health promotion, health promotion activities and interventions, as well as an institutional definition adequate to the health promotion context. The distinctions between sectors were also clarified. The elaborated questionnaires provided in-depth information on countries specifically indicating the key sectors involved in HP4OP in those selected countries. These are: health care, regional/local authorities, NGO's/voluntary institutions. The questionnaire and literature review both resulted in the indication of a significant level of cross-sectorial cooperation in HP4OP. CONCLUSIONS: The inclusion of the institutional analysis within the study of HP4OP provides a valuable opportunity to analyse, in a systematic way, good practices in this respect, also in terms of institutional arrangements. A failure to address this aspect in policymaking might potentially cause organisational failure even in evidence-based programmes. This paper frames the perception of this problem.

CABS-flex predictions of protein flexibility compared with NMR ensembles.

MOTIVATION: Identification of flexible regions of protein structures is important for understanding of their biological functions. Recently, we have developed a fast approach for predicting protein structure fluctuations from a single protein model: the CABS-flex. CABS-flex was shown to be an efficient alternative to conventional all-atom molecular dynamics (MD). In this work, we evaluate CABS-flex and MD predictions by comparison with protein structural variations within NMR ensembles. RESULTS: Based on a benchmark set of 140 proteins, we show that the relative fluctuations of protein residues obtained from CABS-flex are well correlated to those of NMR ensembles. On average, this correlation is stronger than that between MD and NMR ensembles. In conclusion, CABS-flex is useful and complementary to MD in predicting protein regions that undergo conformational changes as well as the extent of such changes. AVAILABILITY AND IMPLEMENTATION: The CABS-flex is freely available to all users at http://biocomp.chem.uw.edu.pl/CABSflex. CONTACT: sekmi@chem.uw.edu.pl SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

CABS-flex: Server for fast simulation of protein structure fluctuations.

The CABS-flex server (http://biocomp.chem.uw.edu.pl/CABSflex) implements CABS-model-based protocol for the fast simulations of near-native dynamics of globular proteins. In this application, the CABS model was shown to be a computationally efficient alternative to all-atom molecular dynamics--a classical simulation approach. The simulation method has been validated on a large set of molecular dynamics simulation data. Using a single input (user-provided file in PDB format), the CABS-flex server outputs an ensemble of protein models (in all-atom PDB format) reflecting the flexibility of the input structure, together with the accompanying analysis (residue mean-square-fluctuation profile and others). The ensemble of predicted models can be used in structure-based studies of protein functions and interactions.

CABS-fold: Server for the de novo and consensus-based prediction of protein structure.

The CABS-fold web server provides tools for protein structure prediction from sequence only (de novo modeling) and also using alternative templates (consensus modeling). The web server is based on the CABS modeling procedures ranked in previous Critical Assessment of techniques for protein Structure Prediction competitions as one of the leading approaches for de novo and template-based modeling. Except for template data, fragmentary distance restraints can also be incorporated into the modeling process. The web server output is a coarse-grained trajectory of generated conformations, its Jmol representation and predicted models in all-atom resolution (together with accompanying analysis). CABS-fold can be freely accessed at http://biocomp.chem.uw.edu.pl/CABSfold.

Structure prediction of the second extracellular loop in G-protein-coupled receptors.

G-protein-coupled receptors (GPCRs) play key roles in living organisms. Therefore, it is important to determine their functional structures. The second extracellular loop (ECL2) is a functionally important region of GPCRs, which poses significant challenge for computational structure prediction methods. In this work, we evaluated CABS, a well-established protein modeling tool for predicting ECL2 structure in 13 GPCRs. The ECL2s (with between 13 and 34 residues) are predicted in an environment of other extracellular loops being fully flexible and the transmembrane domain fixed in its x-ray conformation. The modeling procedure used theoretical predictions of ECL2 secondary structure and experimental constraints on disulfide bridges. Our approach yielded ensembles of low-energy conformers and the most populated conformers that contained models close to the available x-ray structures. The level of similarity between the predicted models and x-ray structures is comparable to that of other state-of-the-art computational methods. Our results extend other studies by including newly crystallized GPCRs.

ClusCo: clustering and comparison of protein models.

BACKGROUND: The development, optimization and validation of protein modeling methods require efficient tools for structural comparison. Frequently, a large number of models need to be compared with the target native structure. The main reason for the development of Clusco software was to create a high-throughput tool for all-versus-all comparison, because calculating similarity matrix is the one of the bottlenecks in the protein modeling pipeline. RESULTS: Clusco is fast and easy-to-use software for high-throughput comparison of protein models with different similarity measures (cRMSD, dRMSD, GDT_TS, TM-Score, MaxSub, Contact Map Overlap) and clustering of the comparison results with standard methods: K-means Clustering or Hierarchical Agglomerative Clustering. CONCLUSIONS: The application was highly optimized and written in C/C++, including the code for parallel execution on CPU and GPU, which resulted in a significant speedup over similar clustering and scoring computation programs.

The 2021 plan for hospital care centralization in Poland - When politics overwhelms the policy process.

The health system in Poland is characterized by oversized hospital infrastructure, with simultaneous deficits in the ambulatory and long-term care sectors. The main challenges of the hospital sector involve i.a. weak stewardship and fragmented governance with a concurrent problem of persistent hospital debts as well as huge workforce deficits. The objective of this paper is to present the government's 2021 plan for hospital care centralization. The reform project aimed i.a. at improving hospital service coordination by implementing a professional and centralized system for hospital sector supervision and effective restructuration processes. The proposed regulation project focused on three major issues: (1) adjusting the existing hospital network towards better concentration of specialized services; (2) launching an independent central agency responsible for monitoring public hospital financial standing as well as initiating and/or managing hospitals restructuration processes; and (3) introducing a formal certification of hospital managers competencies. The reform plans were developed in a relatively short time frame, with a top-down approach and strongly pushed towards the adoption in 2022. Many of the health system stakeholders were strongly opposed to the project which, in connection with new challenges faced by the health system in 2022 (the economic crisis) led the reform suspension. At the same time, a new restructuration and debt relief programme for public hospitals was announced.

Structural features that predict real-value fluctuations of globular proteins.

It is crucial to consider dynamics for understanding the biological function of proteins. We used a large number of molecular dynamics (MD) trajectories of nonhomologous proteins as references and examined static structural features of proteins that are most relevant to fluctuations. We examined correlation of individual structural features with fluctuations and further investigated effective combinations of features for predicting the real value of residue fluctuations using the support vector regression (SVR). It was found that some structural features have higher correlation than crystallographic B-factors with fluctuations observed in MD trajectories. Moreover, SVR that uses combinations of static structural features showed accurate prediction of fluctuations with an average Pearson's correlation coefficient of 0.669 and a root mean square error of 1.04 Å. This correlation coefficient is higher than the one observed in predictions by the Gaussian network model (GNM). An advantage of the developed method over the GNMs is that the former predicts the real value of fluctuation. The results help improve our understanding of relationships between protein structure and fluctuation. Furthermore, the developed method provides a convienient practial way to predict fluctuations of proteins using easily computed static structural features of proteins.

Chirality measures of α-amino acids.

To measure molecular chirality, the molecule is treated as a finite set of points in the Euclidean R(3) space supplemented by k properties, p(1)((i)), p(2)((i)), ..., p(k)((i)) assigned to the ith atom, which constitute a point in the Property P(k) space. Chirality measures are described as the distance between a molecule and its mirror image minimized over all its arbitrary orientation-preserving isometries in the R(3) × P(k) Cartesian product space. Following this formalism, different chirality measures can be estimated by taking into consideration different sets of atomic properties. Here, for α-amino acid zwitterionic structures taken from the Cambridge Structural Database and for all 1684 neutral conformers of 19 biogenic α-amino acid molecules, except glycine and cystine, found at the B3LYP/6-31G** level, chirality measures have been calculated by a CHIMEA program written in this project. It is demonstrated that there is a significant correlation between the measures determined for the α-amino acid zwitterions in crystals and the neutral forms in the gas phase. Performance of the studied chirality measures with changes of the basis set and computation method was also checked. An exemplary quantitative structure­activity relationship (QSAR) application of the chirality measures was presented by an introductory model for the benchmark Cramer data set of steroidal ligands of the sex-hormone binding globulin.

On stability, chirality measures, and theoretical VCD spectra of the chiral C58X2 fullerenes (X = N, B).

The stability of all 23 C(58)N(2) and C(58)B(2) heterofullerenes in the singlet and triplet states was determined at the B3LYP/6-31G** level. In equilibrium mixture the achiral (1,4) C(58)N(2) isomer would be populated in ca. 95.8%, the chiral (1,16) one in ca. 3.3%, and the achiral (1,4) C(58)B(2) in 100%, whereas all triplet state isomers are less stable. Fourteen out of 23 C(58)X(2) are chiral. Four different chirality measures were calculated by our own CHIMEA program: pure geometrical, labeled, mass, and charge. Intercorrelations between the measures for all chiral compounds indicate that the pure geometrical chirality measure is unstable and should not be used in QSAR predictions of the other molecular properties, while the labeled and mass-weighted ones are promising QSAR descriptors. For each chiral C(58)N(2) molecule, some very strong VCD bands, of intensity comparable with that in the IR spectra, can serve in identification and characterization of the isomers.

On vibrational circular dichroism chirality transfer in electron donor-acceptor complexes: a prediction for the quinine···BF3 system.

Vibrational circular dichroism (VCD) chirality transfer occurs when an achiral molecule interacts with a chiral one and becomes VCD-active. Unlike for H-bonds, for organic electron donor-acceptor (EDA) complexes this phenomenon remains almost unknown. Here, the VCD chirality transfer from chiral quinine to achiral BF3 is studied at the B3LYP/aug-cc-pVDZ level. Accessibility of four quinine electron donor sites changes with conformation. Therefore, the quinine conformational landscape was explored and a considerable agreement between X-ray and the most stable conformer geometries was achieved. The BF3 complex through the aliphatic quinuclidine N atom is definitely dominating and is predicted to be easily recognizable in the VCD spectrum. Out of several VCD chirality transfer modes, the ν(s)(BF3) mode, the most intense in the entire VCD spectrum, satisfies the VCD mode robustness criterion and can be used for monitoring the chirality transfer phenomenon in quinine···BF3 system.

A Social Network Analysis of Tweets Related to Mandatory COVID-19 Vaccination in Poland.

Poland's efforts to combat COVID-19 were hindered by endemic vaccination hesitancy and the prevalence of opponents to pandemic restrictions. In this environment, the policy of a COVID-19 vaccination mandate faces strong resistance in the public debate. Exploring the discourse around this resistance could help uncover the motives and develop an understanding of vaccination hesitancy in Poland. This paper aims to conduct a social network analysis and content analysis of Twitter discussions around the intention of the Polish Ministry of Health to introduce mandatory vaccinations for COVID-19. Twitter was chosen as a platform to study because of the critical role it played during the global health crisis. Twitter data were retrieved from 26 July to 9 December 2021 through the API v2 for Academic Research, and analysed using NodeXL and Gephi. When conducting social network analysis, nodes were ranked by their betweenness centrality. Clustering analysis with the Clauset-Newman-Moore algorithm revealed two important groups of users: advocates and opponents of mandatory vaccination. The temporal trends of tweets, the most used hashtags, the sentiment expressed in the most popular tweets, and correlations with epidemiological data were also studied. The results reveal a substantial degree of polarisation, a high intensity of the discussion, and a high degree of involvement of Twitter users. Vaccination mandate advocates were consistently more numerous, but less engaged and less mobilised to "preach" their own stances. Vaccination mandate opponents were vocal and more mobilised to participate: either as original authors or as information diffusers. Our research leads to the conclusion that systematic monitoring of the public debate on vaccines is essential not only in counteracting misinformation, but also in crafting evidence-based as well as emotionally motivating narratives.

Modeling of loops in proteins: a multi-method approach.

BACKGROUND: Template-target sequence alignment and loop modeling are key components of protein comparative modeling. Short loops can be predicted with high accuracy using structural fragments from other, not necessairly homologous proteins, or by various minimization methods. For longer loops multiscale approaches employing coarse-grained de novo modeling techniques should be more effective. RESULTS: For a representative set of protein structures of various structural classes test predictions of loop regions have been performed using MODELLER, ROSETTA, and a CABS coarse-grained de novo modeling tool. Loops of various length, from 4 to 25 residues, were modeled assuming an ideal target-template alignment of the remaining portions of the protein. It has been shown that classical modeling with MODELLER is usually better for short loops, while coarse-grained de novo modeling is more effective for longer loops. Even very long missing fragments in protein structures could be effectively modeled. Resolution of such models is usually on the level 2-6 A, which could be sufficient for guiding protein engineering. Further improvement of modeling accuracy could be achieved by the combination of different methods. In particular, we used 10 top ranked models from sets of 500 models generated by MODELLER as multiple templates for CABS modeling. On average, the resulting molecular models were better than the models from individual methods. CONCLUSIONS: Accuracy of protein modeling, as demonstrated for the problem of loop modeling, could be improved by the combinations of different modeling techniques.

IR low-temperature matrix, X-ray and ab initio study on L-isoserine conformations.

The IR low-temperature Ar and Kr matrix spectra of l-isoserine were registered for the first time and interpreted by means of the anharmonic DFT frequencies calculated at the B3LYP/aug-cc-pVTZ and B3LYP/aug-cc-pVDZ levels. 54 l-isoserine conformers were predicted to be stable at the B3LYP/aug-cc-pVDZ level. Population of the 8 most stable conformers was based on the QCISD/aug-cc-pVDZ energies, corrected for thermal anharmonic factors obtained at the B3LYP/aug-cc-pVDZ level. We found several conformers to be present in the measured matrices and conformer 1 to be dominating. Presence of the conformer 2 is well confirmed by the nu(C=O) band at 1790 cm(-1) and two bands at 1380 and 1350 cm(-1). Presence of the conformer 4 is quite well confirmed by the nu(C-O) bands at 1120 and 1095 cm(-1). Slightly weaker arguments are found for the observation of conformers 6 and 3. Calculations on 54 neutral and 5 zwitterionic conformers in water at the IEF-PCM/B3LYP/aug-cc-pVDZ level suggest that one neutral and one zwitterionic conformer co-exist in the aqueous environment. The crystal structure of l-isoserine was solved by X-ray diffraction analysis. The compound crystallizes without solvent in the chiral P2(1)2(1)2 space group. The asymmetric unit contains a single molecule. The molecule is in its zwitterionic form with the CH(2)-NH(3) side chain in the gauche conformation with respect to the hydroxyl group and in the anti conformation with respect to the carboxylate group. The structure of l-isoserine is dominated by a set of intermolecular hydrogen bonds. The strongest one appears between the OH and COOH groups of two neighbouring molecules: the O...H contact is of 1.66(2) A, which is amongst the shortest H-bonds of this kind observed in amino acid crystal structures.