Folding studies on proteases by the conventional hydrogen exchange experiments are severely hampered because of interference from the autolytic reaction in the interpretation of the exchange data. We report here NMR identification of the hierarchy of early conformational transitions (folding propensities) in HIV-1 protease by systematic monitoring of the changes in the state of the protein as it is subjected to different degrees of denaturation by guanidine hydrochloride. Secondary chemical shifts, H-N-H-alpha coupling constants, H-1-N-15 nuclear Overhauser effects, and N-15 transverse relaxation parameters have been used to report on the residual structural propensities, motional restrictions, conformational transitions, etc., and the data suggest that even under the strongest denaturing conditions (6 M guanidine) hydrophobic clusters as well as different native and non-native secondary structural elements are transiently formed. These constitute the folding nuclei, which include residues spanning the active site, the hinge region, and the dimerization domain. Interestingly, the proline residues influence the structural propensities, and the small amino acids, Gly and Ala, enhance the flexibility of the protein. On reducing the denaturing conditions, partially folded forms appear. The residues showing high folding propensities are contiguous along the sequence at many locations or are in close proximity on the native protein structure, suggesting a certain degree of local cooperativity in the conformational transitions. The dimerization domain, the flaps, and their hinges seem to exhibit the highest folding propensities. The data suggest that even the early folding events may involve many states near the surface of the folding funnel.