How long is a piece of loop?

Challenges in Peptide Solubilization - Amyloids Case Study

Peptide science has been a rapidly growing research field because of the enormous potential application of these biocompatible and bioactive molecules. However, many factors limit the widespread use of peptides in medicine, and low solubility is among the most common problems that hamper drug development in the early stages of research. Solubility is a crucial, albeit poorly understood, feature that determines peptide behavior. Several different solubility predictors have been proposed, and many strategies and protocols have been reported to dissolve peptides, but none of them is a one-size-fits-all method for solubilization of even the same peptide. In this review, we look for the reasons behind the difficulties in dissolving peptides, analyze the factors influencing peptide aggregation, conduct a critical analysis of solubilization strategies and protocols available in the literature, and give some tips on how to deal with the so-called difficult sequences. We focus on amyloids, which are particularly difficult to dissolve and handle such as amyloid beta (Aβ), insulin, and phenol-soluble modulins (PSMs).

Influence of the interaction between p53 and ZNF568 on mitochondrial oxidative phosphorylation

The tumor suppressor p53 plays important roles in suppressing the development and progression of cancer by responding to various stress signals. In addition, p53 can regulate the metabolic pathways of cancer cells by regulating energy metabolism and oxidative phosphorylation. Here, we present a mechanism for the interaction between p53 and ZNF568. Initially, we used X-ray crystallography to determine the irregular loop structure of the ZNF568 KRAB domain; this loop plays an important role in the interaction between p53 and ZNF568. In addition, Cryo-EM was used to examine how the p53 DBD and ZNF568 KRAB domains bind together. The function of ZNF568 on p53-mediated mitochondrial respiration was confirmed by measuring glucose consumption and lactate production. These findings show that ZNF568 can reduce p53-mediated mitochondrial respiratory activity by binding to p53 and inhibiting the transcription of SCO2. SIGNIFICANCE: ZNF568 can directly bind to the p53 DBD and transcriptionally regulate the SCO2 gene. SCO2 transcriptional regulation by interaction between ZNF568 and p53 may regulate the balance between mitochondrial respiration and glycolysis.

Hydrodynamic Radii of Intrinsically Disordered Proteins: Fast Prediction by Minimum Dissipation Approximation and Experimental Validation

The diffusion coefficients of globular and fully unfolded proteins can be predicted with high accuracy solely from their mass or chain length. However, this approach fails for intrinsically disordered proteins (IDPs) containing structural domains. We propose a rapid predictive methodology for estimating the diffusion coefficients of IDPs. The methodology uses accelerated conformational sampling based on self-avoiding random walks and includes hydrodynamic interactions between coarse-grained protein subunits, modeled using the generalized Rotne-Prager-Yamakawa approximation. To estimate the hydrodynamic radius, we rely on the minimum dissipation approximation recently introduced by Cichocki et al. Using a large set of experimentally measured hydrodynamic radii of IDPs over a wide range of chain lengths and domain contributions, we demonstrate that our predictions are more accurate than the Kirkwood approximation and phenomenological approaches. Our technique may prove to be valuable in predicting the hydrodynamic properties of both fully unstructured and multidomain disordered proteins.

Conformation-selective rather than avidity-based binding to tumor associated antigen derived peptide-MHC enables targeting of WT1-pMHC low expressing cancer cells by anti-WT1-pMHC/CD3 T cell engagers

T cell engagers, a category of T cell-retargeting immunotherapy, are rapidly transforming clinical cancer care. However, the lack of tumor-specific targets poses a significant roadblock for broad adaptation of this therapeutic modality in many indications, often resulting in systemic on-target off-tumor toxicity. Though various tumor-derived intracellular mutations provide a massive pool of potential tumor-specific antigens, targeting them is extremely challenging, partly due to the low copy number of tumor associated antigen (TAA)-derived pMHC on tumor cell surface. Further, the interplay of binding geometry and format valency in relation to the capacity of a T cell engager to efficiently target low density cell-surface pMHC is not well understood. Using the Wilms' tumor 1 (WT1) oncoprotein as a proof-of-principle TAA, combined with an array of IgG-like T cell engager modalities that differ in their anti-TAA valency and binding geometry, we show that the ability to induce an immunological synapse formation, resulting in potent killing of WT1 positive cancer cell lines is primarily dependent on the distinct geometrical conformations between the Fab arms of anti-WT1-HLA-A*02:01 and anti-CD3. The augmented avidity conferred by the binding of two anti-WT1-HLA-A*02:01 Fab arms has only minimal influence on cell killing potency. These findings demonstrate the need for careful examination of key design parameters for the development of next-generation T cell engagers targeting low density TAA-pMHCs on tumor cells.

Sequence-independent activity of a predicted long disordered segment of the human papillomavirus type 16 L2 capsid protein during virus entry

The activity of proteins is thought to be invariably determined by their amino acid sequence or composition, but we show that a long segment of a viral protein can support infection independent of its sequence or composition. During virus entry, the papillomavirus L2 capsid protein protrudes through the endosome membrane into the cytoplasm to bind cellular factors such as retromer required for intracellular virus trafficking. Here, we show that an ~110 amino acid segment of L2 is predicted to be disordered and that large deletions in this segment abolish infectivity of HPV16 pseudoviruses by inhibiting cytoplasmic protrusion of L2, association with retromer, and proper virus trafficking. The activity of these mutants can be restored by insertion of protein segments with diverse sequences, compositions, and chemical properties, including scrambled amino acid sequences, a tandem array of a short sequence, and the intrinsically disordered region of an unrelated cellular protein. The infectivity of mutants with small in-frame deletions in this segment directly correlates with the size of the segment. These results indicate that the length of the disordered segment, not its sequence or composition, determines its activity during HPV16 pseudovirus infection. We propose that a minimal length of L2 is required for it to protrude far enough into the cytoplasm to bind cytoplasmic trafficking factors, but the sequence of this segment is largely irrelevant. Thus, protein segments can carry out complex biological functions such as Human papillomavirus pseudovirus infection in a sequence-independent manner. This finding has important implications for protein function and evolution.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

Whole genome sequence analysis of CPV-2 isolates from 1998 to 2020

Canine parvovirus-2 (CPV-2) is a virus with worldwide spread causing canine gastroenteritis. New strains of this virus have unique characteristics and are resistant to some vaccine strains. Therefore, understanding the root causes of resistance has proven to be of increasing concern to many scientists. This study collected 126 whole genome sequences of CPV-2 subtypes with specific collection dates from the NCBI data bank. The whole genome sequences of CPV-2 collected from different countries were analyzed to detect the new substitutions and update these mutations. The result indicated 12, 7, and 10 mutations in NS1, VP1, and VP2, in that respective order. Moreover, the A5G and Q370R mutations of VP2 are the most common changes in the recent isolates of the CPV-2C subtype, and the new N93K residue of VP2 is speculated to be the cause of vaccine failure. To summarize, the observed mutations, which are increasing over time, causes several changes in viral characteristic. A comprehensive understanding of these mutations can lead us to control potential future epidemics associated with this virus more efficiently.

Sequence independent activity of a predicted long disordered segment of the human papillomavirus L2 capsid protein during virus entry

The papillomavirus L2 capsid protein protrudes through the endosome membrane into the cytoplasm during virus entry to bind cellular factors required for intracellular virus trafficking. Cytoplasmic protrusion of HPV16 L2, virus trafficking, and infectivity are inhibited by large deletions in an ∼110 amino acid segment of L2 that is predicted to be disordered. The activity of these mutants can be restored by inserting protein segments with diverse compositions and chemical properties into this region, including scrambled sequences, a tandem array of a short sequence, and the intrinsically disordered region of a cellular protein. The infectivity of mutants with small in-frame insertions and deletions in this segment directly correlates with the size of the segment. These results indicate that the length of the disordered segment, not its sequence or its composition, determines its activity during virus entry. Sequence independent but length dependent activity has important implications for protein function and evolution.

Vaccine design based on 16 epitopes of SARS-CoV-2 spike protein

The global outbreak of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) urgently requires an effective vaccine for prevention. In this study, 66 epitopes containing pentapeptides of SARS‐CoV‐2 spike protein in the IEDB database were compared with the amino acid sequence of SARS‐CoV‐2 spike protein, and 66 potentially immune‐related peptides of SARS‐CoV‐2 spike protein were obtained. Based on the single‐nucleotide polymorphisms analysis of spike protein of 1218 SARS‐CoV‐2 isolates, 52 easily mutated sites were identified and used for vaccine epitope screening. The best vaccine candidate epitopes in the 66 peptides of SARS‐CoV‐2 spike protein were screened out through mutation and immunoinformatics analysis. The best candidate epitopes were connected by different linkers in silico to obtain vaccine candidate sequences. The results showed that 16 epitopes were relatively conservative, immunological, nontoxic, and nonallergenic, could induce the secretion of cytokines, and were more likely to be exposed on the surface of the spike protein. They were both B‐ and T‐cell epitopes, and could recognize a certain number of HLA molecules and had high coverage rates in different populations. Moreover, epitopes 897‐913 were predicted to have possible cross‐immunoprotection for SARS‐CoV and SARS‐CoV‐2. The results of vaccine candidate sequences screening suggested that sequences (without linker, with linker GGGSGGG, EAAAK, GPGPG, and KK, respectively) were the best. The proteins translated by these sequences were relatively stable, with a high antigenic index and good biological activity. Our study provided vaccine candidate epitopes and sequences for the research of the SARS‐CoV‐2 vaccine.

Enhancing fragment-based protein structure prediction by customising fragment cardinality according to local secondary structure

BACKGROUND

Whenever suitable template structures are not available, usage of fragment-based protein structure prediction becomes the only practical alternative as pure ab initio techniques require massive computational resources even for very small proteins. However, inaccuracy of their energy functions and their stochastic nature imposes generation of a large number of decoys to explore adequately the solution space, limiting their usage to small proteins. Taking advantage of the uneven complexity of the sequence-structure relationship of short fragments, we adjusted the fragment insertion process by customising the number of available fragment templates according to the expected complexity of the predicted local secondary structure. Whereas the number of fragments is kept to its default value for coil regions, important and dramatic reductions are proposed for beta sheet and alpha helical regions, respectively.

RESULTS

The evaluation of our fragment selection approach was conducted using an enhanced version of the popular Rosetta fragment-based protein structure prediction tool. It was modified so that the number of fragment candidates used in Rosetta could be adjusted based on the local secondary structure. Compared to Rosetta's standard predictions, our strategy delivered improved first models, + 24% and + 6% in terms of GDT, when using 2000 and 20,000 decoys, respectively, while reducing significantly the number of fragment candidates. Furthermore, our enhanced version of Rosetta is able to deliver with 2000 decoys a performance equivalent to that produced by standard Rosetta while using 20,000 decoys. We hypothesise that, as the fragment insertion process focuses on the most challenging regions, such as coils, fewer decoys are needed to explore satisfactorily conformation spaces.

CONCLUSIONS

Taking advantage of the high accuracy of sequence-based secondary structure predictions, we showed the value of that information to customise the number of candidates used during the fragment insertion process of fragment-based protein structure prediction. Experimentations conducted using standard Rosetta showed that, when using the recommended number of decoys, i.e. 20,000, our strategy produces better results. Alternatively, similar results can be achieved using only 2000 decoys. Consequently, we recommend the adoption of this strategy to either improve significantly model quality or reduce processing times by a factor 10.

Simple Selection Procedure to Distinguish between Static and Flexible Loops

Loops are the most variable and unorganized elements of the secondary structure of proteins. Their ability to shift their shape can play a role in the binding of small ligands, enzymatic catalysis, or protein–protein interactions. Due to the loop flexibility, the positions of their residues in solved structures show the largest B-factors, or in a worst-case scenario can be unknown. Based on the loops’ movements’ timeline, they can be divided into slow (static) and fast (flexible). Although most of the loops that are missing in experimental structures belong to the flexible loops group, the computational tools for loop reconstruction use a set of static loop conformations to predict the missing part of the structure and evaluate the model. We believe that these two loop types can adopt different conformations and that using scoring functions appropriate for static loops is not sufficient for flexible loops. We showed that common model evaluation methods, are insufficient in the case of flexible solvent-exposed loops. Instead, we recommend using the potential energy to evaluate such loop models. We provide a novel model selection method based on a set of geometrical parameters to distinguish between flexible and static loops without the use of molecular dynamics simulations. We have also pointed out the importance of water network and interactions with the solvent for the flexible loop modeling.

Predicting loop conformational ensembles

Motivation

Protein function is often facilitated by the existence of multiple stable conformations. Structure prediction algorithms need to be able to model these different conformations accurately and produce an ensemble of structures that represent a target's conformational diversity rather than just a single state. Here, we investigate whether current loop prediction algorithms are capable of this. We use the algorithms to predict the structures of loops with multiple experimentally determined conformations, and the structures of loops with only one conformation, and assess their ability to generate and select decoys that are close to any, or all, of the observed structures.

Results

We find that while loops with only one known conformation are predicted well, conformationally diverse loops are modelled poorly, and in most cases the predictions returned by the methods do not resemble any of the known conformers. Our results contradict the often-held assumption that multiple native conformations will be present in the decoy set, making the production of accurate conformational ensembles impossible, and hence indicating that current methodologies are not well suited to prediction of conformationally diverse, often functionally important protein regions.

Contact

marks@stats.ox.ac.uk.

Supplementary information

Supplementary data are available at Bioinformatics online.

GradDock: rapid simulation and tailored ranking functions for peptide-MHC Class I docking

Motivation

The identification of T-cell epitopes has many profound translational applications in the areas of transplantation, disease diagnosis, vaccine/therapeutic protein development and personalized immunotherapy. While data-driven methods have been widely used for the prediction of peptide binders with notable successes, the structural modeling of peptide binding to MHC molecules is crucial for understanding the underlying molecular mechanism of the immunological processes.

Results

We developed GradDock, a structure-based method for the rapid and accurate modeling of peptide binding to MHC Class I (pMHC-I). GradDock explicitly models diverse unbound peptides in vacuo and inserts them into the MHC-I groove through a steered gradient descent with a topological correction process. The simulation process yields diverse structural conformations including native-like peptides. We completely revised the Rosetta score terms and developed a new ranking function specifically for pMHC-I. Using the diverse peptides, a linear programming approach is applied to find the optimal weights for the individual Rosetta score terms. Our examination revealed that a refinement of the dihedral angles and a modification of the repulsion can dramatically improve the modeling quality. GradDock is five-times faster than a Rosetta-based docking approach for pMHC-I. We also demonstrate that the predictive capability of GradDock with the re-weighted Rosetta ranking function is consistently more accurate than the Rosetta-based method with the standard Rosetta score (approximately three-times better for a cross-docking set).

Availability and implementation

GradDock is freely available for academic purposes. The program and the ranking score weights for Rosetta are available at http://bel.kaist.ac.kr/research/GradDock.

Contact

hskim76@kaist.ac.kr.

Supplementary information

Supplementary data are available at Bioinformatics online.

The soluble loop BC region guides, but not dictates, the assembly of the transmembrane cytochrome b6

Studying folding and assembly of naturally occurring α-helical transmembrane proteins can inspire the design of membrane proteins with defined functions. Thus far, most studies have focused on the role of membrane-integrated protein regions. However, to fully understand folding pathways and stabilization of α–helical membrane proteins, it is vital to also include the role of soluble loops. We have analyzed the impact of interhelical loops on folding, assembly and stability of the heme-containing four-helix bundle transmembrane protein cytochrome b₆ that is involved in charge transfer across biomembranes. Cytochrome b₆ consists of two transmembrane helical hairpins that sandwich two heme molecules. Our analyses strongly suggest that the loop connecting the helical hairpins is not crucial for positioning the two protein “halves” for proper folding and assembly of the holo-protein. Furthermore, proteolytic removal of any of the remaining two loops, which connect the two transmembrane helices of a hairpin structure, appears to also not crucially effect folding and assembly. Overall, the transmembrane four-helix bundle appears to be mainly stabilized via interhelical interactions in the transmembrane regions, while the soluble loop regions guide assembly and stabilize the holo-protein. The results of this study might steer future strategies aiming at designing heme-binding four-helix bundle structures, involved in transmembrane charge transfer reactions.

An amino acid code for irregular and mixed protein packing

To advance our understanding of protein tertiary structure, the development of the knob-socket model is completed in an analysis of the packing in irregular coil and turn secondary structure packing as well as between mixed secondary structure. The knob-socket model simplifies packing based on repeated patterns of two motifs: a three-residue socket for packing within secondary (2°) structure and a four-residue knob-socket for tertiary (3°) packing. For coil and turn secondary structure, knob-sockets allow identification of a correlation between amino acid composition and tertiary arrangements in space. Coil contributes almost as much as α-helices to tertiary packing. In irregular sockets, Gly, Pro, Asp, and Ser are favored, while in irregular knobs, the preference order is Arg, Asp, Pro, Asn, Thr, Leu, and Gly. Cys, His,Met, and Trp are not favored in either. In mixed packing, the knob amino acid preferences are a function of the socket that they are packing into, whereas the amino acid composition of the sockets does not depend on the secondary structure of the knob. A unique motif of a coil knob with an XYZ β-sheet socket may potentially function to inhibit β-sheet extension. In addition, analysis of the preferred crossing angles for strands within a β-sheet and mixed α-helice/β-sheet identifies canonical packing patterns useful in protein design. Lastly, the knob-socket model abstracts the complexity of protein tertiary structure into an intuitive packing surface topology map.

Building a better fragment library for de novo protein structure prediction

Fragment-based approaches are the current standard for de novo protein structure prediction. These approaches rely on accurate and reliable fragment libraries to generate good structural models. In this work, we describe a novel method for structure fragment library generation and its application in fragment-based de novo protein structure prediction. The importance of correct testing procedures in assessing the quality of fragment libraries is demonstrated. In particular, the exclusion of homologs to the target from the libraries to correctly simulate a de novo protein structure prediction scenario, something which surprisingly is not always done. We demonstrate that fragments presenting different predominant predicted secondary structures should be treated differently during the fragment library generation step and that exhaustive and random search strategies should both be used. This information was used to develop a novel method, Flib. On a validation set of 41 structurally diverse proteins, Flib libraries presents both a higher precision and coverage than two of the state-of-the-art methods, NNMake and HHFrag. Flib also achieves better precision and coverage on the set of 275 protein domains used in the two previous experiments of the the Critical Assessment of Structure Prediction (CASP9 and CASP10). We compared Flib libraries against NNMake libraries in a structure prediction context. Of the 13 cases in which a correct answer was generated, Flib models were more accurate than NNMake models for 10. “Flib is available for download at: http://www.stats.ox.ac.uk/research/proteins/resources”.