An algorithm has been developed to estimate flexibility for potential hinge motion at specified residues, that is, the mutual movement of two domains by rotation around a set of main-chain dihedral angles with torsion angles of neighboring side chains as variables. Such conformational changes must occur without severe atomic collisions. Flexible hinges have been found that satisfy such criteria. Sequence flexibility charts were obtained by plotting the flexibility of each residue against the residue number. Such charts were calculated for 10 proteins (ovomucoid third domain, cytochrome c, lysozyme, hemoglobin /3-chain, a-chymotrypsin, elastase, carboxypeptidase A, dihydrofolate reductase, triosephosphate isomerase, and alcohol dehydrogenase) taken from the Protein Data Bank. The first step of unfolding is likely to occur at the hinge point with the largest flexibility. Following this idea, the polypeptide chain can be dissected into several folding units according to the sequence flexibility chart. When two domains are separated by conformational changes a t such a hinge, the sequence flexibility chart for each domain changes, and it is recalculated and used to indicate subsequent unfolding steps. In this process of iterative estimation of tlexibile hinges, some well-isolated hinges, or the border line between flexible and inflexible regions, were found to be directly at or close to the positions of splice junctions in the eukaryotic genes. Of a total of 45 splice junctions in the 10 proteins examined in this paper, 38 junctions can be identified as flexible hinges between folding units. We suggest that the iterative estimation of flexible hinges may define an array of possible folding/unfolding paths, and that the exon-intron arrangement in the gene may be closely correlated with the folding process of the protein.