Structure of stable protein folding intermediates by equilibrium phi-analysis: the apoflavodoxin thermal intermediate.

Protein intermediates in equilibrium with native states may play important roles in protein dynamics but, in cases, can initiate harmful aggregation events. Investigating equilibrium protein intermediates is thus important for understanding protein behaviour (useful or pernicious) but it is hampered by difficulties in gathering structural information. We show here that the phi-analysis techniques developed to investigate transition states of protein folding can be extended to determine low-resolution three-dimensional structures of protein equilibrium intermediates. The analysis proposed is based solely on equilibrium data and is illustrated by determination of the structure of the apoflavodoxin thermal unfolding intermediate. In this conformation, a large part of the protein remains close to natively folded, but a 40 residue region is clearly unfolded. This structure is fully consistent with the NMR data gathered on an apoflavodoxin mutant designed specifically to stabilise the intermediate. The structure shows that the folded region of the intermediate is much larger than the proton slow-exchange core at 25 degrees C. It also reveals that the unfolded region is made of elements whose packing surface is more polar than average. In addition, it constitutes a useful guide to rationally stabilise the native state relative to the intermediate state, a far from trivial task.

Effect of pyrazinamidase activity on pyrazinamide resistance in Mycobacterium tuberculosis.

Resistance of Mycobacterium tuberculosis to pyrazinamide is associated with mutations in the pncA gene, which codes for pyrazinamidase. The association between the enzymatic activity of mutated pyrazinamidases and the level of pyrazinamide resistance remains poorly understood. Twelve M. tuberculosis clinical isolates resistant to pyrazinamide were selected based on Wayne activity and localization of pyrazinamidase mutation. Recombinant pyrazinamidases were expressed and tested for their kinetic parameters (activity, k(cat), K(m), and efficiency). Pyrazinamide resistance level was measured by Bactec-460TB and 7H9 culture. The linear correlation between the resistance level and the kinetic parameters of the corresponding mutated pyrazinamidase was calculated. The enzymatic activity and efficiency of the mutated pyrazinamidases varied with the site of mutation and ranged widely from low to high levels close to the corresponding of the wild type enzyme. The level of resistance was significantly associated with pyrazinamidase activity and efficiency, but only 27.3% of its statistical variability was explained. Although pyrazinamidase mutations are indeed associated with resistance, the loss of pyrazinamidase activity and efficiency as assessed in the recombinant mutated enzymes is not sufficient to explain a high variability of the level of pyrazinamide resistance, suggesting that complementary mechanisms for pyrazinamide resistance in M. tuberculosis with mutations in pncA are more important than currently thought.

A double-deletion method to quantifying incremental binding energies in proteins from experiment: example of a destabilizing hydrogen bonding pair.

The contribution of a specific hydrogen bond in apoflavodoxin to protein stability is investigated by combining theory, experiment and simulation. Although hydrogen bonds are major determinants of protein structure and function, their contribution to protein stability is still unclear and widely debated. The best method so far devised to estimate the contribution of side-chain interactions to protein stability is double mutant cycle analysis, but the interaction energies so derived are not identical to incremental binding energies (the energies quantifying net contributions of two interacting groups to protein stability). Here we introduce double-deletion analysis of 'isolated' residue pairs as a means to precisely quantify incremental binding. The method is exemplified by studying a surface-exposed hydrogen bond in a model protein (Asp96/Asn128 in apoflavodoxin). Combined substitution of these residues by alanines slightly destabilizes the protein due to a decrease in hydrophobic surface burial. Subtraction of this effect, however, clearly indicates that the hydrogen-bonded groups in fact destabilize the native conformation. In addition, molecular dynamics simulations and classic double mutant cycle analysis explain quantitatively that, due to frustration, the hydrogen bond must form in the native structure because when the two groups get approximated upon folding their binding becomes favorable. We would like to remark that 1), this is the first time the contribution of a specific hydrogen bond to protein stability has been measured by experiment; and 2), more hydrogen bonds need to be analyzed to draw general conclusions on protein hydrogen bond energetics. To that end, the double-deletion method should be of help.

Design of ligand binding to an engineered protein cavity using virtual screening and thermal up-shift evaluation.

Proteins could be used to carry and deliver small compounds. As a tool for designing ligand binding sites in protein cores, a three-step virtual screening method is presented that has been optimised using existing data on T4 lysozyme complexes and tested in a newly engineered cavity in flavodoxin. The method can pinpoint, in large databases, ligands of specific protein cavities. In the first step, physico-chemical filters are used to screen the library and discard a majority of compounds. In the second step, a flexible, fast docking procedure is used to score and select a smaller number of compounds as potential binders. In the third step, a finer method is used to dock promising molecules of the hit list into the protein cavity, and an optimised free energy function allows discarding the few false positives by calculating the affinity of the modelled complexes. To demonstrate the portability of the method, several cavities have been designed and engineered in the flavodoxin from Anabaena PCC 7119, and the W66F/L44A double mutant has been selected as a suitable host protein. The NCI database has then been screened for potential binders, and the binding to the engineered cavity of five promising compounds and three tentative non-binders has been experimentally tested by thermal up-shift assays and spectroscopic titrations. The five tentative binders (some apolar and some polar), unlike the three tentative non-binders, are shown to bind to the host mutant and, importantly, not to bind to the wild type protein. The three-step virtual screening method developed can thus be used to identify ligands of buried protein cavities. We anticipate that the method could also be used, in a reverse manner, to identify natural or engineerable protein cavities for the hosting of ligands of interest.

The long and short flavodoxins: I. The role of the differentiating loop in apoflavodoxin structure and FMN binding.

Flavodoxins are well known one-domain alpha/beta electron-transfer proteins that, according to the presence or absence of a approximately 20-residue loop splitting the fifth beta-strand of the central beta-sheet, have been classified in two groups: long and short-chain flavodoxins, respectively. Although the flavodoxins have been extensively used as models to study electron transfer, ligand binding, protein stability and folding issues, the role of the loop has not been investigated. We have constructed two shortened versions of the long-chain Anabaena flavodoxin in which the split beta-strand has been spliced to remove the original loop. The two variants have been carefully analyzed using various spectroscopic and hydrodynamic criteria, and one of them is clearly well folded, indicating that the long loop is a peripheral element of the structure of long flavodoxins. However, the removal of the loop (which is not in contact with the cofactor in the native structure) markedly decreases the affinity of the apoflavodoxin-FMN complex. This seems related to the fact that, in long flavodoxins, the adjacent tyrosine-bearing FMN binding loop (which is longer and thus more flexible than in short flavodoxins) is stabilized in its competent conformation by interactions with the excised loop. The modest role played by the long loop of long flavodoxins in the structure of these proteins (and in its conformational stability, see Lopez-Llano, J., Maldonado, S., Jain, S., Lostao, A., Godoy-Ruiz, R., Sanchez-Ruiz, Cortijo, M., Fernandez-Recio, J., and Sancho, J. (2004) J. Biol. Chem. 279, 47184-47191) opens the possibility that its conservation in so many species is related to a functional role yet to be discovered. In this respect, we discuss the possibility that the long loop is involved in the recognition of some flavodoxin partners. In addition, we report on a structural feature of flavodoxins that could indicate that the short flavodoxins derive from the long ones.

The long and short flavodoxins: II. The role of the differentiating loop in apoflavodoxin stability and folding mechanism.

Flavodoxins are classified in two groups according to the presence or absence of a approximately 20-residue loop of unknown function. In the accompanying paper (36), we have shown that the differentiating loop from the long-chain Anabaena PCC 7119 flavodoxin is a peripheral structural element that can be removed without preventing the proper folding of the apoprotein. Here we investigate the role played by the loop in the stability and folding mechanism of flavodoxin by comparing the equilibrium and kinetic behavior of the full-length protein with that of loop-lacking, shortened variants. We show that, when the loop is removed, the three-state equilibrium thermal unfolding of apoflavodoxin becomes two-state. Thus, the loop is responsible for the complexity shown by long-chain apoflavodoxins toward thermal denaturation. As for the folding reaction, both shortened and wild type apoflavodoxins display three-state behavior but their folding mechanisms clearly differ. Whereas the full-length protein populates an essentially off-pathway transient intermediate, the additional state observed in the folding of the shortened variant analyzed seems to be simply an alternative native conformation. This finding suggests that the long loop may also be responsible for the accumulation of the kinetic intermediate observed in the full-length protein. Most revealing, however, is that the influence of the loop on the overall conformational stability of apoflavodoxin is quite low and the natively folded shortened variant Delta(120-139) is almost as stable as the wild type protein. The fact that the loop, which is not required for a proper folding of the polypeptide, does not even play a significant role in increasing the conformational stability of the protein supports our proposal (36) that the differentiating loop of long-chain flavodoxins may be related to a recognition function, rather than serving a structural purpose.