A comparison of the predicted and X-ray structures of angiogenin. Implications for further studies of model building of homologous proteins

Exploring the RNase A scaffold to combine catalytic and antimicrobial activities. Structural characterization of RNase 3/1 chimeras

Design of novel antibiotics to fight antimicrobial resistance is one of the first global health priorities. Novel protein-based strategies come out as alternative therapies. Based on the structure-function knowledge of the RNase A superfamily we have engineered a chimera that combines RNase 1 highest catalytic activity with RNase 3 unique antipathogen properties. A first construct (RNase 3/1-v1) was successfully designed with a catalytic activity 40-fold higher than RNase 3, but alas in detriment of its anti-pathogenic activity. Next, two new versions of the original chimeric protein were created showing improvement in the antimicrobial activity. Both second generation versions (RNases 3/1-v2 and -v3) incorporated a loop characteristic of RNase 3 (L7), associated to antimicrobial activity. Last, removal of an RNase 1 flexible loop (L1) in the third version enhanced its antimicrobial properties and catalytic efficiency. Here we solved the 3D structures of the three chimeras at atomic resolution by X-ray crystallography. Structural analysis outlined the key functional regions. Prediction by molecular docking of the protein chimera in complex with dinucleotides highlighted the contribution of the C-terminal region to shape the substrate binding cavity and determine the base selectivity and catalytic efficiency. Nonetheless, the structures that incorporated the key features related to RNase 3 antimicrobial activity retained the overall RNase 1 active site conformation together with the essential structural elements for binding to the human ribonuclease inhibitor (RNHI), ensuring non-cytotoxicity. Results will guide us in the design of the best RNase pharmacophore for anti-infective therapies.

dissectHMMER: a HMMER-based score dissection framework that statistically evaluates fold-critical sequence segments for domain fold similarity

Background

Annotation transfer for function and structure within the sequence homology concept essentially requires protein sequence similarity for the secondary structural blocks forming the fold of a protein. A simplistic similarity approach in the case of non-globular segments (coiled coils, low complexity regions, transmembrane regions, long loops, etc.) is not justified and a pertinent source for mistaken homologies. The latter is either due to positional sequence conservation as a result of a very simple, physically induced pattern or integral sequence properties that are critical for function. Furthermore, against the backdrop that the number of well-studied proteins continues to grow at a slow rate, it necessitates for a search methodology to dive deeper into the sequence similarity space to connect the unknown sequences to the well-studied ones, albeit more distant, for biological function postulations.

Results

Based on our previous work of dissecting the hidden markov model (HMMER) based similarity score into fold-critical and the non-globular contributions to improve homology inference, we propose a framework-dissectHMMER, that identifies more fold-related domain hits from standard HMMER searches. Subsequent statistical stratification of the fold-related hits into cohorts of functionally-related domains allows for the function postulation of the query sequence. Briefly, the technical problems as to how to recognize non-globular parts in the domain model, resolve contradictory HMMER2/HMMER3 results and evaluate fold-related domain hits for homology, are addressed in this work. The framework is benchmarked against a set of SCOP-to-Pfam domain models. Despite being a sequence-to-profile method, dissectHMMER performs favorably against a profile-to-profile based method-HHsuite/HHsearch. Examples of function annotation using dissectHMMER, including the function discovery of an uncharacterized membrane protein Q9K8K1_BACHD (WP_010899149.1) as a lactose/H+ symporter, are presented. Finally, dissectHMMER webserver is made publicly available at http://dissecthmmer.bii.a-star.edu.sg.

Conclusions

The proposed framework-dissectHMMER, is faithful to the original inception of the sequence homology concept while improving upon the existing HMMER search tool through the rescue of statistically evaluated false-negative yet fold-related domain hits to the query sequence. Overall, this translates into an opportunity for any novel protein sequence to be functionally characterized.

Reviewers

This article was reviewed by Masanori Arita, Shamil Sunyaev and L. Aravind.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0068-3) contains supplementary material, which is available to authorized users.

On the necessity of dissecting sequence similarity scores into segment-specific contributions for inferring protein homology, function prediction and annotation

Background

Protein sequence similarities to any types of non-globular segments (coiled coils, low complexity regions, transmembrane regions, long loops, etc. where either positional sequence conservation is the result of a very simple, physically induced pattern or rather integral sequence properties are critical) are pertinent sources for mistaken homologies. Regretfully, these considerations regularly escape attention in large-scale annotation studies since, often, there is no substitute to manual handling of these cases. Quantitative criteria are required to suppress events of function annotation transfer as a result of false homology assignments.

Results

The sequence homology concept is based on the similarity comparison between the structural elements, the basic building blocks for conferring the overall fold of a protein. We propose to dissect the total similarity score into fold-critical and other, remaining contributions and suggest that, for a valid homology statement, the fold-relevant score contribution should at least be significant on its own. As part of the article, we provide the DissectHMMER software program for dissecting HMMER2/3 scores into segment-specific contributions. We show that DissectHMMER reproduces HMMER2/3 scores with sufficient accuracy and that it is useful in automated decisions about homology for instructive sequence examples. To generalize the dissection concept for cases without 3D structural information, we find that a dissection based on alignment quality is an appropriate surrogate. The approach was applied to a large-scale study of SMART and PFAM domains in the space of seed sequences and in the space of UniProt/SwissProt.

Conclusions

Sequence similarity core dissection with regard to fold-critical and other contributions systematically suppresses false hits and, additionally, recovers previously obscured homology relationships such as the one between aquaporins and formate/nitrite transporters that, so far, was only supported by structure comparison.

Transmembrane helix: simple or complex

Transmembrane helical segments (TMs) can be classified into two groups of so-called ‘simple’ and ‘complex’ TMs. Whereas the first group represents mere hydrophobic anchors with an overrepresentation of aliphatic hydrophobic residues that are likely attributed to convergent evolution in many cases, the complex ones embody ancestral information and tend to have structural and functional roles beyond just membrane immersion. Hence, the sequence homology concept is not applicable on simple TMs. In practice, these simple TMs can attract statistically significant but evolutionarily unrelated hits during similarity searches (whether through BLAST- or HMM-based approaches). This is especially problematic for membrane proteins that contain both globular segments and TMs. As such, we have developed the transmembrane helix: simple or complex (TMSOC) webserver for the identification of simple and complex TMs. By masking simple TM segments in seed sequences prior to sequence similarity searches, the false-discovery rate decreases without sacrificing sensitivity. Therefore, TMSOC is a novel and necessary sequence analytic tool for both the experimentalists and the computational biology community working on membrane proteins. It is freely accessible at http://tmsoc.bii.a-star.edu.sg or available for download.

Not all transmembrane helices are born equal: Towards the extension of the sequence homology concept to membrane proteins

BACKGROUND

Sequence homology considerations widely used to transfer functional annotation to uncharacterized protein sequences require special precautions in the case of non-globular sequence segments including membrane-spanning stretches composed of non-polar residues. Simple, quantitative criteria are desirable for identifying transmembrane helices (TMs) that must be included into or should be excluded from start sequence segments in similarity searches aimed at finding distant homologues.

RESULTS

We found that there are two types of TMs in membrane-associated proteins. On the one hand, there are so-called simple TMs with elevated hydrophobicity, low sequence complexity and extraordinary enrichment in long aliphatic residues. They merely serve as membrane-anchoring device. In contrast, so-called complex TMs have lower hydrophobicity, higher sequence complexity and some functional residues. These TMs have additional roles besides membrane anchoring such as intra-membrane complex formation, ligand binding or a catalytic role. Simple and complex TMs can occur both in single- and multi-membrane-spanning proteins essentially in any type of topology. Whereas simple TMs have the potential to confuse searches for sequence homologues and to generate unrelated hits with seemingly convincing statistical significance, complex TMs contain essential evolutionary information.

CONCLUSION

For extending the homology concept onto membrane proteins, we provide a necessary quantitative criterion to distinguish simple TMs (and a sufficient criterion for complex TMs) in query sequences prior to their usage in homology searches based on assessment of hydrophobicity and sequence complexity of the TM sequence segments.

Structure of the virus capsid protein VP1 of enterovirus 71 predicted by some homology modeling and molecular docking studies

The homology modeling technique has been used to construct the structure of enterovirus 71 (EV 71) capsid protein VP1. The protein is consisted of 297 amino acid residues and treated as the target. The amino acid sequence identity between the target protein and sequences of template proteins 1EAH, 1PIV, and 1D4M searched from NCBI protein BLAST and WorkBench protein tools were 38, 37, and 36%, respectively. Based on these template structures, the protein model was constructed by using the InsightII/Homology program. The protein model was briefly refined by energy minimization and molecular dynamics (MD) simulation steps. The protein model was validated using some web available servers such as ERRAT, PROCHECK, PROVE, and PROSA2003. However, an inconsistency between the docking scores and the measured activity was observed for a series of EV 71 VP1 inhibitors synthesized by Shia et al. (J Med Chem 2002, 45, 1644) and docked into the binding pocket of the protein model using the DOCK 4.0.2 program. The protein model with an EV 71 VP1 inhibitor docked and engulfed was then refined further by some MD simulation steps in the presence of water molecules. The docking scores obtained for these inhibitors after such a MD refinement were well correlated with the activities. The structure-activity relationships for the ligand-protein model system was also analyzed using the GRID-VOLSURF programs and the corresponding noncrossvalidated and crossvalidated (by leave-one-out) r2 and q2 were 0.99 and 0.61, respectively. The hydrophobic nature of the binding pocket of the protein model was also examined using the GRID21 program. The possibility of improving the potency of the current series of EV 71 VP1 inhibitors was discussed based on all the studies presented.

Evolution of physics-based methodology for exploring the conformational energy landscape of proteins

The evolution of our physics-based computational methods for determining protein conformation without the introduction of secondary-structure predictions, homology modeling, threading, or fragment coupling is described. Initial use of a hard-sphere potential captured much of the structural properties of polypeptide chains, and subsequent more refined force fields, together with efficient methods of global optimization provide indications that progress is being made toward an understanding of the interresidue interactions that underlie protein folding.

Force probe measurements of antibody-antigen interactions

The surface force apparatus has been used to quantify directly the forces that govern the interactions between proteins and ligands. In this work, we describe the measured interactions between the antigen fluorescein and the Fab' fragment of the monoclonal 4-4-20 anti-fluorescyl IgG antibody. Here we first describe the use of the surface force apparatus to demonstrate directly the impact of the charge composition in the region of the antibody binding site on the antibody interactions. Several approaches are described for immobilizing antigens, antibodies, and proteins in general for direct force measurements. The measured force profiles presented are accompanied by an extensive discussion of protocols used to analyze the force-distance curves and to interpret them in terms of the antibody structure. In addition to long-range electrostatic forces, we also consider short-range forces that can affect the strength of adhesion between the Fab' and immobilized fluorescein. The latter investigations demonstrate the influence of interfacial properties on the recognition of surface-bound antigens.