Automated tertiary structure prediction with accurate local model quality assessment using the IntFOLD-TS method

The IntFOLD-TS method was developed according to the guiding principle that the model quality assessment would be the most critical stage for our template based modelling pipeline. Thus, the IntFOLD-TS method firstly generates numerous alternative models, using in-house versions of several different sequence-structure alignment methods, which are then ranked in terms of global quality using our top performing quality assessment method – ModFOLDclust2. In addition to the predicted global quality scores, the predictions of local errors are also provided in the resulting coordinate files, using scores that represent the predicted deviation of each residue in the model from the equivalent residue in the native structure. The IntFOLD-TS method was found to generate high quality 3D models for many of the CASP9 targets, whilst also providing highly accurate predictions of their per-residue errors. This important information may help to make the 3D models that are produced by the IntFOLD-TS method more useful for guiding future experimental work

gamma-zein secondary structure in solution by circular dichroism

The proline-rich N-terminal domain of gamma-zein has been reported in relevant process, which include its ability to cross the cell membranes. Evidences indicate that synthetic hexapeptide (PPPVHL), naturally found in N-terminal portion of gamma-zein, can adopt the polyproline II (PPII) conformation in aqueous solution. The secondary structure of gamma-zein in maize protein bodies had been analyzed by solid state Fourier transform infrared and nuclear magnetic resonance spectroscopies. However, it was not possible to measure PPII content in physiological environment since the beta-sheet and PPII signals overlap in both solid state techniques. Here, the secondary structure of gamma-zein has been analyzed by circular dichroism in SDS aqueous solution with and without ditiothreitol (DTT), and in 60% of 2-propanol and water with DTT The results show that gamma-zein has high helical content in all solutions. The PPII conformation was present at about 7% only in water/DTT solution. (c) 2007 Wiley Periodicals, Inc.

Structure and function in rhodopsin: Packing of the helices in the transmembrane domain and folding to a tertiary structure in the intradiscal domain are coupled*

A previous study of the retinitis pigmentosa mutation L125R and two designed mutations at this site, L125A and L125F, showed that these mutations cause partial or total misfolding of the opsins expressed in COS cells from the corresponding mutant opsin genes. We now report on expression and characterization of the opsins from the following retinitis pigmentosa mutants in the transmembrane domain of rhodopsin that correspond to six of the seven helices: G51A and G51V (helix A), G89D (helix B), A164V (helix D), H211P (helix E), P267L and P267R (helix F), and T297R (helix G). All the mutations caused partial misfolding of the opsins as observed by the UV/visible absorption characteristics and by separation of the expressed opsins into fractions that bound 11-cis-retinal to form the corresponding mutant rhodopsins and those that did not bind 11-cis-retinal. Further, all the mutant rhodopsins prepared from the above mutants, except for G51A, showed strikingly abnormal bleaching behavior with abnormal metarhodopsin II photointermediates. The results show that retinitis pigmentosa mutations in every one of the transmembrane helices can cause misfolding of the opsin. Therefore, on the basis of these and previous results, we conclude that defects in the packing of the transmembrane helices resulting from these mutations are relayed to the intradiscal domain, where they cause misfolding of the opsin by inducing the formation of a disulfide bond other than the native Cys-110—Cys-187 disulfide bond. Thus, there is coupling between packing of the helices in the transmembrane domain and folding to a tertiary structure in the intradiscal domain.

Conformational changes in tertiary structure near the ligand binding site of an integrin I domain

For efficient ligand binding, integrins must be activated. Specifically, a conformational change has been proposed in a ligand binding domain present within some integrins, the inserted (I) domain [Lee, J., Bankston, L., Arnaout, M. & Liddington, R. C. (1995) Structure (London) 3, 1333–1340]. This proposal remains controversial, however, despite extensive crystal structure studies on the I domain [Lee, J., Bankston, L., Arnaout, M. & Liddington, R. C. (1995) Structure (London) 3, 1333–1340; Liddington, R. & Bankston, L. (1998) Structure (London) 6, 937–938; Qu, A. & Leahy, D. J. (1996) Structure (London) 4, 931–942; and Baldwin, E. T., Sarver, R. W., Bryant, G. L., Jr., Curry, K. A., Fairbanks, M. B., Finzel, B. C., Garlick, R. L., Heinrikson, R. L., Horton, N. C. & Kelly, L. L. (1998) Structure (London) 6, 923–935]. By defining the residues present in the epitope of a mAb against the human Mac-1 integrin (αMβ2, CD11b/CD18) that binds only the active receptor, we provide biochemical evidence that the I domain itself undergoes a conformational change with activation. This mAb, CBRM1/5, binds the I domain very close to the ligand binding site in a region that is widely exposed regardless of activation as judged by reactivity with other antibodies. The conformation of the epitope differs in two crystal forms of the I domain, previously suggested to represent active and inactive receptor. Our data suggests that conformational differences in the I domain are physiologically relevant and not merely a consequence of different crystal lattice interactions. We also demonstrate that the transition between the two conformational states depends on species-specific residues at the bottom of the I domain, which are proposed to be in an interface with another integrin domain, and that this transition correlates with functional activity.

Tertiary structure of an amyloid immunoglobulin light chain protein: a proposed model for amyloid fibril formation.

An immunoglobulin light chain protein was isolated from the urine of an individual (BRE) with systemic amyloidosis. Complete amino acid sequence of the variable region of the light chain (VL) protein established it as a kappa I, which when compared with other kappa I amyloid associated proteins had unique residues, including Ile-34, Leu-40, and Tyr-71. To study the tertiary structure, BRE VL was expressed in Escherichia coli by using a PCR product amplified from the patient BRE's bone marrow DNA. The PCR product was ligated into pCZ11, a thermal-inducible replication vector. Recombinant BRE VL was isolated, purified to homogeneity, and crystallized by using ammonium sulfate as the precipitant. Two crystal forms were obtained. In crystal form I the BRE VL kappa domain crystallizes as a dimer with unit cell constants isomorphous to previously published kappa protein structures. Comparison with a nonamyloid VL kappa domain from patient REI, identified significant differences in position of residues in the hypervariable segments plus variations in framework region (FR) segments 40-46 (FR2) and 66-67 (FR3). In addition, positional differences can be seen along the two types of local diads, corresponding to the monomer-monomer and dimer-dimer interfaces. From the packing diagram, a model for the amyloid light chain (AL) fibril is proposed based on a pseudohexagonal spiral structure with a rise of approximately the width of two dimers per 360 degree turn. This spiral structure could be consistent with the dimensions of amyloid fibrils as determined by electron microscopy.

Phocid seal leptin: tertiary structure and hydrophobic receptor binding site preservation during distinct leptin gene evolution

The cytokine hormone leptin is a key signalling molecule in many pathways that control physiological functions. Although leptin demonstrates structural conservation in mammals, there is evidence of positive selection in primates, lagomorphs and chiropterans. We previously reported that the leptin genes of the grey and harbour seals (phocids) have significantly diverged from other mammals. Therefore we further investigated the diversification of leptin in phocids, other marine mammals and terrestrial taxa by sequencing the leptin genes of representative species. Phylogenetic reconstruction revealed that leptin diversification was pronounced within the phocid seals with a high dN/dS ratio of 2.8, indicating positive selection. We found significant evidence of positive selection along the branch leading to the phocids, within the phocid clade, but not over the dataset as a whole. Structural predictions indicate that the individual residues under selection are away from the leptin receptor (LEPR) binding site. Predictions of the surface electrostatic potential indicate that phocid seal leptin is notably different to other mammalian leptins, including the otariids. Cloning the grey seal leptin binding domain of LEPR confirmed that this was structurally conserved. These data, viewed in toto, support a hypothesis that phocid leptin divergence is unlikely to have arisen by random mutation. Based upon these phylogenetic and structural assessments, and considering the comparative physiology and varying life histories among species, we postulate that the unique phocid diving behaviour has produced this selection pressure. The Phocidae includes some of the deepest diving species, yet have the least modified lung structure to cope with pressure and volume changes experienced at depth. Therefore, greater surfactant production is required to facilitate rapid lung re-inflation upon surfacing, while maintaining patent airways. We suggest that this additional surfactant requirement is met by the leptin pulmonary surfactant production pathway which normally appears only to function in the mammalian foetus.

A new class of isothiocyanate-based irreversible inhibitors of macrophage migration inhibitory factor.

Macrophage migration inhibitory factor (MIF) is a homotrimeric multifunctional proinflammatory cytokine that has been implicated in the pathogenesis of several inflammatory and autoimmune diseases. Current therapeutic strategies for targeting MIF focus on developing inhibitors of its tautomerase activity or modulating its biological activities using anti-MIF neutralizing antibodies. Herein we report a new class of isothiocyanate (ITC)-based irreversible inhibitors of MIF. Modification by benzyl isothiocyanate (BITC) and related analogues occurred at the N-terminal catalytic proline residue without any effect on the oligomerization state of MIF. Different alkyl and arylalkyl ITCs modified MIF with nearly the same efficiency as BITC. To elucidate the mechanism of action, we performed detailed biochemical, biophysical, and structural studies to determine the effect of BITC and its analogues on the conformational state, quaternary structure, catalytic activity, receptor binding, and biological activity of MIF. Light scattering, analytical ultracentrifugation, and NMR studies on unmodified and ITC-modified MIF demonstrated that modification of Pro1 alters the tertiary, but not the secondary or quaternary, structure of the trimer without affecting its thermodynamic stability. BITC induced drastic effects on the tertiary structure of MIF, in particular residues that cluster around Pro1 and constitute the tautomerase active site. These changes in tertiary structure and the loss of catalytic activity translated into a reduction in MIF receptor binding activity, MIF-mediated glucocorticoid overriding, and MIF-induced Akt phosphorylation. Together, these findings highlight the role of tertiary structure in modulating the biochemical and biological activities of MIF and present new opportunities for modulating MIF biological activities in vivo.

The C-terminal Region Of The Human P23 Chaperone Modulates Its Structure And Function.

The p23 protein is a chaperone widely involved in protein homeostasis, well known as an Hsp90 co-chaperone since it also controls the Hsp90 chaperone cycle. Human p23 includes a β-sheet domain, responsible for interacting with Hsp90; and a charged C-terminal region whose function is not clear, but seems to be natively unfolded. p23 can undergo caspase-dependent proteolytic cleavage to form p19 (p231-142), which is involved in apoptosis, while p23 has anti-apoptotic activity. To better elucidate the function of the human p23 C-terminal region, we studied comparatively the full-length human p23 and three C-terminal truncation mutants: p23₁₋₁₁₇; p23₁₋₁₃₁ and p23₁₋₁₄₂. Our data indicate that p23 and p19 have distinct characteristics, whereas the other two truncations behave similarly, with some differences to p23 and p19. We found that part of the C-terminal region can fold in an α-helix conformation and slightly contributes to p23 thermal-stability, suggesting that the C-terminal interacts with the β-sheet domain. As a whole, our results suggest that the C-terminal region of p23 is critical for its structure-function relationship. A mechanism where the human p23 C-terminal region behaves as an activation/inhibition module for different p23 activities is proposed.

Solution structure of the sodium channel antagonist conotoxin GS: A new molecular caliper for probing sodium channel geometry

Background: The venoms of Conus snails contain small, disulfide-rich inhibitors of voltage-dependent sodium channels. Conotoxin GS is a 34-residue polypeptide isolated from Conus geographus that interacts with the extracellular entrance of skeletal muscle sodium channels to prevent sodium ion conduction. Although conotoxin GS binds competitively with mu conotoxin GIIIA to the sodium channel surface, the two toxin types have little sequence identity with one another, and conotoxin GS has a four-loop structural framework rather than the characteristic three-loop mu-conotoxin framework. The structural study of conotoxin GS will form the basis for establishing a structure-activity relationship and understanding its interaction with the pore region of sodium channels. Results: The three-dimensional structure of conotoxin GS was determined using two-dimensional NMR spectroscopy. The protein exhibits a compact fold incorporating a beta hairpin and several turns. An unusual feature of conotoxin GS is the exceptionally high proportion (100%) of cis-imide bond geometry for the three proline or hydroxyproline residues. The structure of conotoxin GS bears little resemblance to the three-loop mu conotoxins, consistent with the low sequence identity between the two toxin types and their different structural framework. However, the tertiary structure and cystine-knot motif formed by the three disulfide bonds is similar to that present in several other polypeptide ion channel inhibitors. Conclusions: This is the first three-dimensional structure of a 'four-loop' sodium channel inhibitor, and it represents a valuable new structural probe for the pore region of voltage-dependent sodium channels. The distribution of amino acid sidechains in the structure creates several polar and charged patches, and comparison with the mu conotoxins provides a basis for determining the binding surface of the conotoxin GS polypeptide.

The solution structure of C1-T1, a two-domain proteinase inhibitor derived from a circular precursor protein from Nicotiana alata

A two-domain portion of the proteinase inhibitor precursor from Nicotiana alata (NaProPI) has been expressed and its structure determined by NMR spectroscopy. NaProPI contains six almost identical 53 amino acid repeats that fold into six highly similar domains; however, the sequence repeats do nut coincide with the structural domains. Five of the structural domains comprise the C-terminal portion of one repeat and the N-terminal portion of the next. The sixth domain contains the C-terminal portion of the sixth repeat and the N-terminal portion of the first repeat. Disulphide bonds link these C and N-terminal fragments to generate the clasped-bracelet fold of NaProPI. The three-dimensional structure of NaProPI is not known, but it is conceivable that adjacent domains in NaProPI interact to generate the circular bracelet with the N and C termini in close enough proximity to facilitate formation of the disulphide bonds that form the clasp The expressed protein, examined in the current study, comprises residues 25-135 of NaProPI and encompasses the first two contiguous structural domains, namely the chymotrypsin inhibitor C1 and the trypsin inhibitor T1, joined by a five-residue linker, and is referred to as C1-T1. The tertiary structure of each domain in C1-T1 is identical to that found in the isolated inhibitors. However, no nuclear Overhauser effect contacts are observed between the two domains and the five-residue linker adopts an extended conformation. The absence of interactions between the domains indicates that adjacent domains do not specifically interact to drive the circularisation of NaProPI. These results are in agreement with recent data which describe similar PI precursors from other members of the Solanaceae having two, three, or four repeats. The lack of strong interdomain association is likely to be important for the function of individual inhibitors by ensuring that there is no masking of reactive sites upon release from the precursor. (C) 2001 Academic Press.

Discovery of catalases in members of the Chlamydiales order.

Catalase is an important virulence factor for survival in macrophages and other phagocytic cells. In Chlamydiaceae, no catalase had been described so far. With the sequencing and annotation of the full genomes of Chlamydia-related bacteria, the presence of different catalase-encoding genes has been documented. However, their distribution in the Chlamydiales order and the functionality of these catalases remain unknown. Phylogeny of chlamydial catalases was inferred using MrBayes, maximum likelihood, and maximum parsimony algorithms, allowing the description of three clade 3 and two clade 2 catalases. Only monofunctional catalases were found (no catalase-peroxidase or Mn-catalase). All presented a conserved catalytic domain and tertiary structure. Enzymatic activity of cloned chlamydial catalases was assessed by measuring hydrogen peroxide degradation. The catalases are enzymatically active with different efficiencies. The catalase of Parachlamydia acanthamoebae is the least efficient of all (its catalytic activity was 2 logs lower than that of Pseudomonas aeruginosa). Based on the phylogenetic analysis, we hypothesize that an ancestral class 2 catalase probably was present in the common ancestor of all current Chlamydiales but was retained only in Criblamydia sequanensis and Neochlamydia hartmannellae. The catalases of class 3, present in Estrella lausannensis and Parachlamydia acanthamoebae, probably were acquired by lateral gene transfer from Rhizobiales, whereas for Waddlia chondrophila they likely originated from Legionellales or Actinomycetales. The acquisition of catalases on several occasions in the Chlamydiales suggests the importance of this enzyme for the bacteria in their host environment.

Caractérisation de l’ubiquitin-fold modifier (UFM1) dans un modèle C. elegans

L’ubiquitin-fold modifier (UFM1) fait partie de la classe 1 de la famille de protéine ubiquitin-like (Ubl). UFM1 et Ub ont très peu d’homologie de séquence, mais partagent des similarités remarquables au niveau de leur structure tertiaire. Tout comme l’Ub et la majorité des autres Ubls, UFM1 se lie de façon covalente à ses substrats par l’intermédiaire d’une cascade enzymatique. Il est de plus en plus fréquemment rapporté que les protéines Ubls sont impliquées dans des maladies humaines. Le gène Ufm1 est surexprimé chez des souris de type MCP développant une ischémie myocardique et dans les îlots de Langerhans de patients atteints du diabète de type 2. UFM1 et ses enzymes spécifiques, UBA5, UFL1 et UFC1, sont conservés chez les métazoaires et les plantes suggérant un rôle important pour les organismes multicellulaires. Le Caenorhabditis elegans est le modèle animal le plus simple utilisé en biologie. Sa morphologie, ses phénotypes visibles et ses lignées cellulaires ont été décrits de façon détaillée. De plus, son cycle de vie court permet de rapidement observer les effets de certains gènes sur la longévité. Ce modèle nous permet de facilement manipuler l’expression du gène Ufm1 et de mieux connaître ses fonctions. En diminuant l’expression du gène ufm-1 chez le C.elegans, par la technique de l’ARN interférence par alimentation, nous n’avons observé aucun problème morphologique grave. Les vers ressemblaient aux vers sauvages et possédaient un nombre de progéniture normal. Cependant, les vers sauvage exposés à l’ARNi d’ufm-1 vivent significativement moins longtemps que les contrôles et ce, de façon indépendante de la voie de signalisation de l’insuline/IGF. Chez le C. elegans la longévité et la résistance au stress cellulaire sont intimement liées. Nous n’avons remarqué aucun effet d’ufm-1 sur le stress thermal, osmotique ou oxydatif, mais il est requis pour la protection contre le stress protéotoxique. Il est également nécessaire au maintien de l’intégrité neuronale au cours du vieillissement des animaux. L’ensemble de nos données nous renseigne sur les fonctions putatives du gène Ufm1.

Quantification de la relation séquence-activité de l’ARN par prédiction de structure tridimensionnelle

Dans un premier temps, nous avons modélisé la structure d’une famille d’ARN avec une grammaire de graphes afin d’identifier les séquences qui en font partie. Plusieurs autres méthodes de modélisation ont été développées, telles que des grammaires stochastiques hors-contexte, des modèles de covariance, des profils de structures secondaires et des réseaux de contraintes. Ces méthodes de modélisation se basent sur la structure secondaire classique comparativement à nos grammaires de graphes qui se basent sur les motifs cycliques de nucléotides. Pour exemplifier notre modèle, nous avons utilisé la boucle E du ribosome qui contient le motif Sarcin-Ricin qui a été largement étudié depuis sa découverte par cristallographie aux rayons X au début des années 90. Nous avons construit une grammaire de graphes pour la structure du motif Sarcin-Ricin et avons dérivé toutes les séquences qui peuvent s’y replier. La pertinence biologique de ces séquences a été confirmée par une comparaison des séquences d’un alignement de plus de 800 séquences ribosomiques bactériennes. Cette comparaison a soulevée des alignements alternatifs pour quelques unes des séquences que nous avons supportés par des prédictions de structures secondaires et tertiaires. Les motifs cycliques de nucléotides ont été observés par les membres de notre laboratoire dans l'ARN dont la structure tertiaire a été résolue expérimentalement. Une étude des séquences et des structures tertiaires de chaque cycle composant la structure du Sarcin-Ricin a révélé que l'espace des séquences dépend grandement des interactions entre tous les nucléotides à proximité dans l’espace tridimensionnel, c’est-à-dire pas uniquement entre deux paires de bases adjacentes. Le nombre de séquences générées par la grammaire de graphes est plus petit que ceux des méthodes basées sur la structure secondaire classique. Cela suggère l’importance du contexte pour la relation entre la séquence et la structure, d’où l’utilisation d’une grammaire de graphes contextuelle plus expressive que les grammaires hors-contexte. Les grammaires de graphes que nous avons développées ne tiennent compte que de la structure tertiaire et négligent les interactions de groupes chimiques spécifiques avec des éléments extra-moléculaires, comme d’autres macromolécules ou ligands. Dans un deuxième temps et pour tenir compte de ces interactions, nous avons développé un modèle qui tient compte de la position des groupes chimiques à la surface des structures tertiaires. L’hypothèse étant que les groupes chimiques à des positions conservées dans des séquences prédéterminées actives, qui sont déplacés dans des séquences inactives pour une fonction précise, ont de plus grandes chances d’être impliqués dans des interactions avec des facteurs. En poursuivant avec l’exemple de la boucle E, nous avons cherché les groupes de cette boucle qui pourraient être impliqués dans des interactions avec des facteurs d'élongation. Une fois les groupes identifiés, on peut prédire par modélisation tridimensionnelle les séquences qui positionnent correctement ces groupes dans leurs structures tertiaires. Il existe quelques modèles pour adresser ce problème, telles que des descripteurs de molécules, des matrices d’adjacences de nucléotides et ceux basé sur la thermodynamique. Cependant, tous ces modèles utilisent une représentation trop simplifiée de la structure d’ARN, ce qui limite leur applicabilité. Nous avons appliqué notre modèle sur les structures tertiaires d’un ensemble de variants d’une séquence d’une instance du Sarcin-Ricin d’un ribosome bactérien. L’équipe de Wool à l’université de Chicago a déjà étudié cette instance expérimentalement en testant la viabilité de 12 variants. Ils ont déterminé 4 variants viables et 8 létaux. Nous avons utilisé cet ensemble de 12 séquences pour l’entraînement de notre modèle et nous avons déterminé un ensemble de propriétés essentielles à leur fonction biologique. Pour chaque variant de l’ensemble d’entraînement nous avons construit des modèles de structures tertiaires. Nous avons ensuite mesuré les charges partielles des atomes exposés sur la surface et encodé cette information dans des vecteurs. Nous avons utilisé l’analyse des composantes principales pour transformer les vecteurs en un ensemble de variables non corrélées, qu’on appelle les composantes principales. En utilisant la distance Euclidienne pondérée et l’algorithme du plus proche voisin, nous avons appliqué la technique du « Leave-One-Out Cross-Validation » pour choisir les meilleurs paramètres pour prédire l’activité d’une nouvelle séquence en la faisant correspondre à ces composantes principales. Finalement, nous avons confirmé le pouvoir prédictif du modèle à l’aide d’un nouvel ensemble de 8 variants dont la viabilité à été vérifiée expérimentalement dans notre laboratoire. En conclusion, les grammaires de graphes permettent de modéliser la relation entre la séquence et la structure d’un élément structural d’ARN, comme la boucle E contenant le motif Sarcin-Ricin du ribosome. Les applications vont de la correction à l’aide à l'alignement de séquences jusqu’au design de séquences ayant une structure prédéterminée. Nous avons également développé un modèle pour tenir compte des interactions spécifiques liées à une fonction biologique donnée, soit avec des facteurs environnants. Notre modèle est basé sur la conservation de l'exposition des groupes chimiques qui sont impliqués dans ces interactions. Ce modèle nous a permis de prédire l’activité biologique d’un ensemble de variants de la boucle E du ribosome qui se lie à des facteurs d'élongation.

Germin and germin-like proteins: evolution, structure, and function

Germin and germin-like proteins (GLPs) are encoded by a family of genes found in all plants. They are part of the cupin superfamily of biochemically diverse proteins, a superfamily that has a conserved tertiary structure, though with limited similarity in primary sequence. The subgroups of GLPs have different enzyme functions that include the two hydrogen peroxide-generating enzymes, oxalate oxidase (OxO) and superoxide dismutase. This review summarizes the sequence and structural details of GLPs and also discusses their evolutionary progression, particularly their amplification in gene number during the evolution of the land plants. In terms of function, the GLPs are known to be differentially expressed during specific periods of plant growth and development, a pattern of evolutionary subfunctionalization. They are also implicated in the response of plants to biotic (viruses, bacteria, mycorrhizae, fungi, insects, nematodes, and parasitic plants) and abiotic (salt, heat/cold, drought, nutrient, and metal) stress. Most detailed data come from studies of fungal pathogenesis in cereals. This involvement with the protection of plants from environmental stress of various types has led to numerous plant breeding studies that have found links between GLPs and QTLs for disease and stress resistance. In addition the OxO enzyme has considerable commercial significance, based principally on its use in the medical diagnosis of oxalate concentration in plasma and urine. Finally, this review provides information on the nutritional importance of these proteins in the human diet, as several members are known to be allergenic, a feature related to their thermal stability and evolutionary connection to the seed storage proteins, also members of the cupin superfamily.

Rapid model quality assessment for protein structure predictions using the comparison of multiple models without structural alignments

MOTIVATION: The accurate prediction of the quality of 3D models is a key component of successful protein tertiary structure prediction methods. Currently, clustering or consensus based Model Quality Assessment Programs (MQAPs) are the most accurate methods for predicting 3D model quality; however they are often CPU intensive as they carry out multiple structural alignments in order to compare numerous models. In this study, we describe ModFOLDclustQ - a novel MQAP that compares 3D models of proteins without the need for CPU intensive structural alignments by utilising the Q measure for model comparisons. The ModFOLDclustQ method is benchmarked against the top established methods in terms of both accuracy and speed. In addition, the ModFOLDclustQ scores are combined with those from our older ModFOLDclust method to form a new method, ModFOLDclust2, that aims to provide increased prediction accuracy with negligible computational overhead. RESULTS: The ModFOLDclustQ method is competitive with leading clustering based MQAPs for the prediction of global model quality, yet it is up to 150 times faster than the previous version of the ModFOLDclust method at comparing models of small proteins (<60 residues) and over 5 times faster at comparing models of large proteins (>800 residues). Furthermore, a significant improvement in accuracy can be gained over the previous clustering based MQAPs by combining the scores from ModFOLDclustQ and ModFOLDclust to form the new ModFOLDclust2 method, with little impact on the overall time taken for each prediction. AVAILABILITY: The ModFOLDclustQ and ModFOLDclust2 methods are available to download from: http://www.reading.ac.uk/bioinf/downloads/ CONTACT: l.j.mcguffin@reading.ac.uk.