The systematic design of functional peptides has technological and therapeutic applications. However, there is a need for pattern-based search engines that help locate desired functional motifs in primary sequences regardless of their evolutionary conservation. Existing databases such as The Protein Secondary Structure database (PSS) no longer serves the community, while the Dictionary of Protein Secondary Structure (DSSP) annotates the secondary structures when tertiary structures of proteins are provided. Here, we extract 1.7 million helices from the PDB and compile them into a database (Therapeutic Peptide Design database; TP-DB) that allows queries of compounded patterns to facilitate the identification of sequence motifs of helical structures. We show how TP-DB helps us identify a known purification-tag-specific antibody that can be repurposed into a diagnostic kit for Helicobacter pylori. We also show how the database can be used to design a new antimicrobial peptide that shows better Candida albicans clearance and lower hemolysis than its template homologs. Finally, we demonstrate how TP-DB can suggest point mutations in helical peptide blockers to prevent a targeted tumorigenic protein-protein interaction. TP-DB is made available at http://dyn.life.nthu.edu.tw/design/.

Plasmodium falciparum malaria, which degrades haemoglobin through falcipain-2 (FP2), is a serious disease killing 445 thousand people annually. Since the P. falciparum’s survival in humans depends on its ability to degrade human’s haemoglobin, stoppage or hindrance of FP2 has antimalarial effects. Therefore, we studied the atomic details of how E64 approaches, binds to, and inhibits FP2. We found that E64 (1) gradually approaches FP2 by first interacting with FP2’s D170 and Q171 or N81, N77, and K76; (2) binds FP2 tightly (ΔG_binding = −12.2 ± 1.1 kJ/mol); and (3) persistently blocks access to FP2’s catalytic residues regardless of whether or not E64 has already been able to form a covalent bond with FP2’s C42. Furthermore, the results suggest that S41, D234, D170, N38, N173, and L172 (which are located in or near the FP2’s catalytic site’s binding pocket) contribute the most towards the favourable binding of E64 to FP2. Their in silico mutations adversely affect E64-FP2 binding affinity with D234L/A, N173L/A, W43F/A, D234L/A, H174F/A, and N38L/A having the most significant adverse effects on E64-FP2 binding and interactions. The findings presented in this article, which has antimalarial implications, suggest that hydrogen bonding and electrostatic interactions play important roles in E64-FP2 binding, and that a potential FP2-blocking E64-based/E64-like antimalarial drug should be capable of being both hydrogen-bond donor and acceptor, and/or have the ability to favourably interact with polar amino acids (such as S41, S149, N38, N173, N77, Q171) and with charged amino acids (such as D234, D170, H174) of FP2. The abilities to favourably interact with ASN, ASP, and SER appears to be important characteristics that such potential drug should have.

Primary dystonia's prolonged muscle contractions and the associated abnormal postures and twisting movements remain incurable. Genetic mutation/deletion of GAG from TorsonA's gene resulting in ΔE303 (which weakens the binding between TorsinA and its activator, such as LULL1) primarily cause this neurodegenerative disorder. We studied TorsinA-LULL1 (or TorsinAΔE303-LULL1) bindings and interactions. For the first time, we show the atomic details of TorsinA-LULL1 dynamic interactions and TorsinAΔE303-LULL1 dynamic interactions and their binding affinities. Our results show extensive effects of ΔE303 on TorsinAΔE303-LULL1 interactions, and suggest that the differences between TorsinA-LULL1 interactions and TorsinAΔE303-LULL1 interactions are non-subtle. ΔE303 significantly weakens TorsinAΔE303-LULL1's binding affinity. We present pieces of evidence proving that the effects of ΔE303 (on the differences between TorsinA-LULL1 interactions and TorsinAΔE303-LULL1 interactions) are more pronounced than previously suggested, and that the nanobody used for achieving the X-ray crystallization in the previous study attenuated the differences between TorsinA-LULL1 and TorsinAΔE303-LULL1 interactions. Our accounts of the dynamic interactions between “TorsinA and LULL1” and between “TorsinAΔE303 and LULL1” and the detailed effects of ΔE303 on TorsinA-/TorsinAΔE303-LULL1 build on previous findings and offer new insights for a better understanding of the molecular basis of Primary Dystonia. Our results have long-term potentials of guiding the development of medications for the disease.

Secondary structures (SS) of proteins are of great importance to structural, molecular, and computational biology and chemistry. Accurate and reliable method for automatic SS assignment when only coarse-grained (CG) information is available is needed. RaFoSA, a novel, accurate, and reliable method for automatic SS assignment based on coordinates of alpha carbon (CAC) atoms alone is presented here. Results from RaFoSa have been rigorously compared to those from Dictionary of Protein SS (DSSP, the acclaimed gold-standard for automatic SS assignment) and STRIDE. Requiring only CAC, RaFoSA achieves an agreement of 96% (and 94%) with DSSP (and STRIDE) that require all-atom and hydrogen-bonding information. No known automatic SS assignment method based on CG system has ever achieved such agreement with DSSP and STRIDE. Furthermore, RaFoSA has been applied to a real-life problem and its possible use for ranking proteins in their order of SS-based stability is shown in this paper. Overall, RaFoSA’s abilities to accurately and reliably assign SS to CG or all-atom protein systems make this work important. Furthermore, it must be emphasized that SS assignment by RaFoSA is different from (and is more rigorous than) SS prediction from amino acids sequence. Indeed, SS assignment by RaFoSA can differentiate between frames from molecular dynamics simulations trajectories, while existing methods for SS prediction from amino acid sequence cannot. Source codes and a webserver implementation of RaFoSA are available at http://bioinformatics.center/RaFoSA.

Receptor-binding and subsequent signal-activation of interleukin-1 beta (IL-1β) are essential to immune and proinflammatory responses. We mutated 12 residues to identify sites important for biological activity and/or receptor binding. Four of these mutants with mutations in loop 9 (T117A, E118K, E118A, E118R) displayed significantly reduced biological activity. Neither T117A nor E118K mutants substantially affected receptor binding, whereas both mutants lack the IL-1β signaling in vitro but can antagonize wild-type (WT) IL-1β. Crystal structures of T117A, E118A, and E118K revealed that the secondary structure or surface charge of loop 9 is dramatically altered compared with that of wild-type chicken IL-1β. Molecular dynamics simulations of IL-1β bound to its receptor (IL-1RI) and receptor accessory protein (IL-1RAcP) revealed that loop 9 lies in a pocket that is formed at the IL-1RI/IL-1RAcP interface. This pocket is also observed in the human ternary structure. The conformations of above mutants in loop 9 may disrupt structural packing and therefore the stability in a chicken IL-1β/IL-1RI/IL-1RAcP signaling complex. We identify the hot spots in IL-1β that are essential to immune responses and elucidate a mechanism by which IL-1β activity can be inhibited. These findings should aid in the development of new therapeutics that neutralize IL-1 activity.