Benchmarking of TASSER in the ab initio limit

A significant number of protein sequences in a given proteome have no obvious evolutionarily related protein in the database of solved protein structures, the PDB. Under these conditions, ab initio or template-free modeling methods are the sole means of predicting protein structure. To assess its expected performance on proteomes, the TASSER structure prediction algorithm is benchmarked in the ab initio limit on a representative set of 1129 nonhomologous sequences ranging from 40 to 200 residues that cover the PDB at 30% sequence identity and which adopt alpha, alpha + beta, and beta secondary structures. For sequences in the 40-100 (100-200) residue range, as assessed by their root mean square deviation from native, RMSD, the best of the top five ranked models of TASSER has a global fold that is significantly close to the native structure for 25% (16%) of the sequences, and with a correct identification of the structure of the protein core for 59% (36%). In the absence of a native structure, the structural similarity among the top five ranked models is a moderately reliable predictor of folding accuracy. If we classify the sequences according to their secondary structure content, then 64% (36%) of alpha, 43% (24%) of alpha + beta, and 20% (12%) of beta sequences in the 40-100 (100-200) residue range have a significant TM-score (TM-score > or = 0.4). TASSER performs best on helical proteins because there are less secondary structural elements to arrange in a helical protein than in a beta protein of equal length, since the average length of a helix is longer than that of a strand. In addition, helical proteins have shorter loops and dangling tails. If we exclude these flexible fragments, then TASSER has similar accuracy for sequences containing the same number of secondary structural elements, irrespective of whether they are helices and/or strands. Thus, it is the effective configurational entropy of the protein that dictates the average likelihood of correctly arranging the secondary structure elements.     

CONFOLD: Residue‐residue contact‐guided ab initio protein folding

Predicted protein residue–residue contacts can be used to build three‐dimensional models and consequently to predict protein folds from scratch. A considerable amount of effort is currently being spent to improve contact prediction accuracy, whereas few methods are available to construct protein tertiary structures from predicted contacts. Here, we present an ab initio protein folding method to build three‐dimensional models using predicted contacts and secondary structures. Our method first translates contacts and secondary structures into distance, dihedral angle, and hydrogen bond restraints according to a set of new conversion rules, and then provides these restraints as input for a distance geometry algorithm to build tertiary structure models. The initially reconstructed models are used to regenerate a set of physically realistic contact restraints and detect secondary structure patterns, which are then used to reconstruct final structural models. This unique two‐stage modeling approach of integrating contacts and secondary structures improves the quality and accuracy of structural models and in particular generates better β‐sheets than other algorithms. We validate our method on two standard benchmark datasets using true contacts and secondary structures. Our method improves TM‐score of reconstructed protein models by 45% and 42% over the existing method on the two datasets, respectively. On the dataset for benchmarking reconstructions methods with predicted contacts and secondary structures, the average TM‐score of best models reconstructed by our method is 0.59, 5.5% higher than the existing method. The CONFOLD web server is available at http://protein.rnet.missouri.edu/confold/ . Proteins 2015; 83:1436–1449. © 2015 Wiley Periodicals, Inc.     

Toward optimal fragment generations for ab initio protein structure assembly

Fragment assembly using structural motifs excised from other solved proteins has shown to be an efficient method for ab initio protein‐structure prediction. However, how to construct accurate fragments, how to derive optimal restraints from fragments, and what the best fragment length is are the basic issues yet to be systematically examined. In this work, we developed a gapless‐threading method to generate position‐specific structure fragments. Distance profiles and torsion angle pairs are then derived from the fragments by statistical consistency analysis, which achieved comparable accuracy with the machine‐learning‐based methods although the fragments were taken from unrelated proteins. When measured by both accuracies of the derived distance profiles and torsion angle pairs, we come to a consistent conclusion that the optimal fragment length for structural assembly is around 10, and at least 100 fragments at each location are needed to achieve optimal structure assembly. The distant profiles and torsion angle pairs as derived by the fragments have been successfully used in QUARK for ab initio protein structure assembly and are provided by the QUARK online server at http://zhanglab.ccmb. med.umich.edu/QUARK/ . Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Use of residual dipolar couplings as restraints in ab initio protein structure prediction

NMR residual dipolar couplings (RDCs), in the form of the projection angles between the respective internuclear bond vectors, are used as structural restraints in the ab initio structure prediction of a test set of six proteins. The restraints are applied using a recently developed SICHO (SIde-CHain-Only) lattice protein model that employs a replica exchange Monte Carlo (MC) algorithm to search conformational space. Using a small number of RDC restraints, the quality of the predicted structures is improved as reflected by lower RMSD/dRMSD (root mean square deviation/distance root mean square deviation) values from the corresponding native structures and by the higher correlation of the most cooperative mode of motion of each predicted structure with that of the native structure. The latter, in particular, has possible implications for the structure-based functional analysis of predicted structures.     

Distance geometry generates native-like folds for small helical proteins using the consensus distances of predicted protein structures.

For successful ab initio protein structure prediction, a method is needed to identify native-like structures from a set containing both native and non-native protein-like conformations. In this regard, the use of distance geometry has shown promise when accurate inter-residue distances are available. We describe a method by which distance geometry restraints are culled from sets of 500 protein-like conformations for four small helical proteins generated by the method of Simons et al. (1997). A consensus-based approach was applied in which every inter-Calpha distance was measured, and the most frequently occurring distances were used as input restraints for distance geometry. For each protein, a structure with lower coordinate root-mean-square (RMS) error than the mean of the original set was constructed; in three cases the topology of the fold resembled that of the native protein. When the fold sets were filtered for the best scoring conformations with respect to an all-atom knowledge-based scoring function, the remaining subset of 50 structures yielded restraints of higher accuracy. A second round of distance geometry using these restraints resulted in an average coordinate RMS error of 4.38 A.     

Ab initio protein modelling reveals novel human MIT domains

Database searches can fail to detect all truly homologous sequences, particularly when dealing with short, highly sequence diverse protein families. Here, using microtubule interacting and transport (MIT) domains as an example, we have applied an approach of profile-profile matching followed by ab initio structure modelling to the detection of true homologues in the borderline significant zone of database searches. Novel MIT domains were confidently identified in USP54, containing an apparently inactive ubiquitin carboxyl-terminal hydrolase domain, a katanin-like ATPase KATNAL1, and an uncharacterized protein containing a VPS9 domain. As a proof of principle, we have confirmed the novel MIT annotation for USP54 by in vitro profiling of binding to CHMP proteins.

Structured summary

USP8 binds:CHMPs 1A 1B 2A 2B 4CUSP54 binds:CHMPs 1B 2A 2B 4C 6     

Contact order and ab initio protein structure prediction

Although much of the motivation for experimental studies of protein folding is to obtain insights for improving protein structure prediction, there has been relatively little connection between experimental protein folding studies and computational structural prediction work in recent years. In the present study, we show that the relationship between protein folding rates and the contact order (CO) of the native structure has implications for ab initio protein structure prediction. Rosetta ab initio folding simulations produce a dearth of high CO structures and an excess of low CO structures, as expected if the computer simulations mimic to some extent the actual folding process. Consistent with this, the majority of failures in ab initio prediction in the CASP4 (critical assessment of structure prediction) experiment involved high CO structures likely to fold much more slowly than the lower CO structures for which reasonable predictions were made. This bias against high CO structures can be partially alleviated by performing large numbers of additional simulations, selecting out the higher CO structures, and eliminating the very low CO structures; this leads to a modest improvement in prediction quality. More significant improvements in predictions for proteins with complex topologies may be possible following significant increases in high-performance computing power, which will be required for thoroughly sampling high CO conformations (high CO proteins can take six orders of magnitude longer to fold than low CO proteins). Importantly for such a strategy, simulations performed for high CO structures converge much less strongly than those for low CO structures, and hence, lack of simulation convergence can indicate the need for improved sampling of high CO conformations. The parallels between Rosetta simulations and folding in vivo may extend to misfolding: The very low CO structures that accumulate in Rosetta simulations consist primarily of local up-down beta-sheets that may resemble precursors to amyloid formation.     

Alpha-helical cationic antimicrobial peptides: relationships of structure and function

Antimicrobial peptides (AMPs), with their extraordinary properties, such as broad-spectrum activity, rapid action and difficult development of resistance, have become promising molecules as new antibiotics. Despite their various mechanisms of action, the interaction of AMPs with the bacterial cell membrane is the key step for their mode of action. Moreover, it is generally accepted that the membrane is the primary target of most AMPs, and the interaction between AMPs and eukaryotic cell membranes (causing toxicity to host cells) limits their clinical application. Therefore, researchers are engaged in reforming or de novo designing AMPs as a ‘single-edged sword’ that contains high antimicrobial activity yet low cytotoxicity against eukaryotic cells. To improve the antimicrobial activity of AMPs, the relationship between the structure and function of AMPs has been rigorously pursued. In this review, we focus on the current knowledge of α-helical cationic antimicrobial peptides, one of the most common types of AMPs in nature.     

家蝇抗菌肽的研究与应用

近年来,对家蝇Muscadomestica体内外抗菌活性物质的研究受到人们广泛的关注,其中研究较多的是抗菌肽。家蝇抗菌肽对细菌、真菌、肿瘤癌细胞和病毒等具有生物活性,可通过针刺、带菌针刺、超声波、放射性射线、高频电磁场和生理盐水等方法诱导增量产生,能耐极端的温度和pH值的溶液,在高浓度盐溶液中也很稳定。文章还对家蝇抗菌肽的结构特点、作用机理与分子生物学研究做了概述,同时对家蝇抗菌肽的应用前景进行了讨论。     

Design and characterization of novel hybrid antimicrobial peptides based on cecropin A, LL-37 and magainin II

Antimicrobial peptides (AMPs) are a naturally occurring component of the innate immune response of many organisms and can have activity against both Gram-negative and Gram-positive bacterial species. In order to optimize and improve the direct antimicrobial effect of AMPs against a broad spectrum of bacterial species, novel synthetic hybrids were rationally designed from cecropin A, LL-37 and magainin II. AMPs were selected based on their α-helical secondary structure and fragments of these were analyzed and combined in silico to determine which hybrid peptides would form the best amphipathic cationic α-helices. Four hybrid peptides were synthesized (CaLL, CaMA, LLaMA and MALL) and evaluated for direct antimicrobial activity against a range of bacterial species (Bacillus anthracis, Burkholderia cepacia, Francisella tularensis LVS and Yersinia pseudotuberculosis) alongside the original 'parent' AMPs. The hybrid peptides showed greater antimicrobial effects than the parent AMPs (in one case a parent is completely ineffective while a hybrid based on it removes all traces of bacteria by 3h), although they also demonstrated higher hemolytic properties. Modifications were then carried out to the most toxic hybrid AMP (CaLL) to further improve the therapeutic index. Modifications made to the hybrid lowered hemolytic activity and also lowered antimicrobial activity by various degrees. Overall, this work highlights the potential for rational design and synthesis of improved AMPs that have the capability to be used therapeutically for treatment of bacterial infections.     

Structure-function studies of chemokine-derived carboxy-terminal antimicrobial peptides

Recent reports which show that several chemokines can act as direct microbicidal agents have drawn renewed attention to these chemotactic signalling proteins. Here we present a structure-function analysis of peptides derived from the human chemokines macrophage inflammatory protein-3α (MIP-3α/CCL20), interleukin-8 (IL-8), neutrophil activating protein-2 (NAP-2) and thrombocidin-1 (TC-1). These peptides encompass the C-terminal α-helices of these chemokines, which have been suggested to be important for the direct antimicrobial activities. Far-UV CD spectroscopy showed that the peptides are unstructured in aqueous solution and that a membrane mimetic solvent is required to induce a helical secondary structure. A co-solvent mixture was used to determine solution structures of the peptides by two-dimensional ¹H-NMR spectroscopy. The highly cationic peptide, MIP-3α_51-70, had the most pronounced antimicrobial activity and displayed an amphipathic structure. A shorter version of this peptide, MIP-3α_59-70, remained antimicrobial but its structure and mechanism of action were unlike that of the former peptide. The NAP-2 and TC-1 proteins differ in their sequences only by the deletion of two C-terminal residues in TC-1, but intact TC-1 is a very potent antimicrobial while NAP-2 is inactive. The corresponding C-terminal peptides, NAP-2_50-70 and TC-1_50-68, had very limited and no bactericidal activity, respectively. This suggests that other regions of TC-1 contribute to its bactericidal activity. Altogether, this work provides a rational structural basis for the biological activities of these peptides and proteins and highlights the importance of experimental characterization of peptide fragments as distinct entities because their activities and structural properties may differ substantially from their parent proteins.     

Encapsulin carrier proteins for enhanced expression of antimicrobial peptides

Antimicrobial peptides (AMPs) are regarded as attractive alternatives to conventional antibiotics, but their production in microbes remains challenging due to their inherent bactericidal nature. To address these limitations, we have developed a novel AMP fusion protein system based on an encapsulin nanocompartment protein and have demonstrated its utility in enhancing expression of HBCM2, an AMP with activity against Gram-negative bacteria. Here, HBCM2 was fused to the N-terminus of several Encapsulin monomer (Enc) variants engineered with multiple TEV protease recognition site insertions to facilitate proteolytic release of the fused HBCM2. Fusion of HBCM2 to the Enc variants, but not other common carrier proteins, enabled robust overexpression in Escherichia coli C43(DE3) cells. Interestingly, variants with a TEV site insertion following residue K71 in Enc exhibited the highest overexpression and HBCM2 release efficiencies compared to other variants but were deficient in cage formation. HBCM2 was purified from the highest expressing variant following TEV protease digestion and was found to be highly active in inhibiting E. coli growth (MIC = 5 μg/ml). Our study demonstrates the potential use of the Enc system to enhance expression of AMPs for biomanufacturing and therapeutic applications.     

Proteomic analysis reveals responsive proteins of Vibrio parahaemolyticus on exposure to cationic antimicrobial peptides

Aims: To investigate whether Vibrio parahaemolyticus can sense and directly respond to the presence of cationic antimicrobial peptides (AMPs). Methods and Results: We performed proteomic methodologies to investigate the responsive proteins of V. parahaemolyticus on exposure to AMP Q6. Differential expression patterns of the outer membrane, inner membrane and cytoplasmic proteins (OMPs, IMPs and CPs) from the bacteria with or without Q6 treatment were obtained using two‐dimensional gel electrophoresis (2‐DE). Three OMPs (maltoporin, flagellin and OmpV), two IMPs (ATP synthase F1, alpha subunit; and OmpV) and three CPs (pyruvate dehydrogenase subunit E1, glyceraldehyde‐3‐phosphate dehydrogenase and inositol‐5‐monophosphate dehydrogenase) were identified using LC‐ESI‐Q‐TOF MS/MS and Mascot program. Real‐time quantitative polymerase chain reaction was also performed to determine the mRNA expression level of the target genes. Conclusions: Our results suggested that V. parahaemolyticus may directly respond to AMPs through the upregulation of the efflux channel, increased yield of energy, effective repair of damaged membranes and downregulation of carbohydrate and nucleotide metabolism for energy preservation. Significance and Impact of the Study: AMP‐responsive proteins identified in the current study could serve as attractive targets for developing more effective antimicrobial agents.     

Ab initio quantum mechanical models of peptide helices and their vibrational spectra

Structural parameters for standard peptide helices (alpha, 3(10), 3(1) left-handed) were fully ab initio optimized for Ac-(L-Ala)(9)-NHMe and for Ac-(L-Pro)(9)-NHMe (poly-L-proline-PLP I and PLP II-forms), in order to better understand the relative stability and minimum energy geometries of these conformers and the dependence of the ir absorption and vibrational CD (VCD) spectra on detailed variation in these conformations. Only the 3(10)-helical Ala-based conformation was stable in vacuum for this decaamide structure, but both Pro-based conformers minimized successfully. Inclusion of solvent effects, by use of the conductor-like screening solvent model (COSMO), enabled ab initio optimizations [at the DFT/B3LYP/SV(P) level] without any constraints for the alpha- and 3(10)-helical Ala-based peptides as well as the two Pro-based peptides. The geometries obtained compare well with peptide chain torsion angles and hydrogen-bond distances found for these secondary structure types in x-ray structures of peptides and proteins. For the simulation of VCD spectra, force field and intensity response tensors were obtained ab initio for the complete Ala-based peptides in vacuum, but constrained to the COSMO optimized torsional angles, due to limitations of the solvent model. Resultant spectral patterns reproduce well many aspects of the experimental spectra and capture the differences observed for these various helical types.     

Prediction of protein tertiary structure using PROFESY, a novel method based on fragment assembly and conformational space annealing

A novel method for ab initio prediction of protein tertiary structures, PROFESY (PROFile Enumerating SYstem), is proposed. This method utilizes the secondary structure prediction information of a query sequence and the fragment assembly procedure based on global optimization. Fifteen-residue-long fragment libraries are constructed using the secondary structure prediction method PREDICT, and fragments in these libraries are assembled to generate full-length chains of a query protein. Tertiary structures of 50 to 100 conformations are obtained by minimizing an energy function for proteins, using the conformational space annealing method that enables one to sample diverse low-lying local minima of the energy. We apply PROFESY for benchmark tests to proteins with known structures to demonstrate its feasibility. In addition, we participated in CASP5 and applied PROFESY to four new-fold targets for blind prediction. The results are quite promising, despite the fact that PROFESY was in its early stages of development. In particular, PROFESY successfully provided us the best model-one structure for the target T0161.     

PackHelix: a tool for helix-sheet packing during protein structure prediction

The three‐dimensional structure of a protein is organized around the packing of its secondary structure elements. Predicting the topology and constructing the geometry of structural motifs involving α‐helices and/or β‐strands are therefore key steps for accurate prediction of protein structure. While many efforts have focused on how to pack helices and on how to sample exhaustively the topologies and geometries of multiple strands forming a β‐sheet in a protein, there has been little progress on generating native‐like packings of helices on sheets. We describe a method that can generate the packing of multiple helices on a given β‐sheet for αβα sandwich type protein folds. This method mines the results of a statistical analysis of the conformations of αβ₂ motifs in protein structures to provide input values for the geometric attributes of the packing of a helix on a sheet. It then proceeds with a geometric builder that generates multiple arrangements of the helices on the sheet of interest by sampling through these values and performing consistency checks that guarantee proper loop geometry between the helices and the strands, minimal number of collisions between the helices, and proper formation of a hydrophobic core. The method is implemented as a module of ProteinShop. Our results show that it produces structures that are within 4–6 Å RMSD of the native one, regardless of the number of helices that need to be packed, though this number may increase if the protein has several helices between two consecutive strands in the sequence that pack on the sheet formed by these two strands. Proteins 2011; Published 2011 Wiley‐Liss, Inc.     

The Implementation of the Iterative Diagonalization Scheme and ab initio Molecular Dynamics Simulation with the LAPW Method

Abstract

A new ab initio molecular dynamics method based on the full-potential linearized-augmented-plane-wave (LAPW) basis set has been implemented. The LAPW basis set has been successfully employed for systems containing localized electrons such as first row atoms and transition metals. In our implementation of the LAPW-MD scheme, iterative residual minimization algorithm is used to solve the electronic states problem. The atoms are moved according to forces derived from the Hellman–Feynman theorem and incomplete basis set correction terms. The performance of the program is further enhanced by parallelization. We will discuss technical details of the program implementation and present results obtained from this code to the equilibrium structures and vibrational properties of simple diatomic molecules.