The flap region in aspartic proteases is a unique structural feature to this class of enzymes, and found to have a profound impact on protein overall structure, function, and dynamics. Understanding the structure and dynamic behavior of the flap regions is crucial in the design of selective inhibitors against aspartic proteases. Cathepsin-D, an aspartic protease enzyme, has been implicated in a long list of degenerative diseases as well as breast cancer progression. Presented herein, for the first time, is a comprehensive description of the conformational flap dynamics of cathepsin-D using a comparative 50 ns "multiple" molecular dynamics simulations. Diverse collective metrics were proposed to accurately define flap dynamics. These are distance d1 between the flap tips residues (Gly79 and Met301); dihedral angle ϕ; in addition to TriCα angles Gly79-Asp33-Asp223, θ1 , and Gly79-Asp223-Met301, θ2 . The maximum distance attained throughout the simulation was 17.42 and 11.47 Å for apo and bound cathepsin-D, respectively, while the minimum distance observed was 8.75 and 6.32 Å for apo and bound cathepsin-D, respectively. The movement of the flap as well as the twist of the active pocket can properly be explained by measuring the angle, θ1 , between Gly79-Asp33-Met301 and correlating it with the distance Cα of the flap tip residues. The asymmetrical opening of the binding cavity was best described by the large shift of -6.26° to +20.94° in the dihedral angle, ϕ, corresponding to the full opening of the flap at a range of 31-33 ns. A wide-range of post-dynamic analyses was also applied in this report to supplement our findings. We believe that this report would augment current efforts in designing potent structure-based inhibitors against cathepsin-D in the treatment of breast cancer and other degenerative diseases. J. Cell. Biochem. 117: 2643-2657, 2016. © 2016 Wiley Periodicals, Inc.