Epidermal growth factor (EGF) regulates cell proliferation and differentiation by binding to the EGF receptor (EGFR) extra-cellular domains. Human EGF is a small, single-chain protein comprising three distinct loops (A, B, and C), which are connected by three disulfide bridges (Cys6-Cys20, Cys14-Cys31, and Cys33-Cys42). These disulfide bridges are essential for structural stability and biological activity. EGF was extensively studied by disulfide scrambling, an experimental technique for the conformational entrapment of intermediate states, which allows us to study the folding pathway of proteins containing disulfide bonds. The experimental results showed that there is a major 2-disulfide intermediate (denoted EGF-II) and that the native disulfide bonding pattern is less prevalent in one of the mutants. In this article, we investigated for the first time the solution conformations of wild-type EGF, EGF-II, and the mutant S9C through extensive molecular dynamics (MD) simulations in water using both the standard MD technique and a recently developed amplified-collective-motion (ACM) sampling method. Compared to standard MD simulations, we achieved a much more enhanced sampling by the ACM simulations, and the structures were sufficiently relaxed to estimate configurational entropies. The simulation results suggest a predominantly entropic folding pathway governed by the disorder of three functional loop regions. Although EGF-II exhibits two native disulfide bonds (Cys14-Cys31 and Cys33-Cys42), its large configurational entropy inhibits a direct transition to the native structure in the folding process. When Ser9 is mutated into Cys, a non-native disulfide bridge Cys9-Cys20 is slightly more favorable than the native Cys6-Cys20 because a less constrained N-terminus affords larger entropy. Isomers that are functionally less active also exhibit a more localized dynamics of the functional loop regions, which may suggest a possible mechanism for the modulation of EGF activity.