A new approach incorporating flexible docking simulations and NMR data is presented for calculating the bound conformation of a ligand that interacts weakly with an enzyme. This approach consists of sampling directly the conformation of a flexible ligand inside a receptor active site containing surrounding flexible loops. To make this sampling efficient, a ligand-growing procedure has been adopted. Optimization of the ECEPP/3-plus-NOE constraint function is carried out by using a collective variable Monte Carlo minimization technique. Numerous energy minimizations are made possible for such a large system by using a Bezier splines energy grid technique. This new flexible docking approach was applied to determine the structure of a fibrinogen Aalpha-like peptide (7DFLAEGGGVRGPRV20) bound to an active site mutant of thrombin [thrombin(S195A)]. Structure calculations of the bound ligand, using 2D-transferred NOESY distance constraints in the DIANA program, showed that the N-terminal portion of the peptide (D7-R16) involves a chain reversal, whereas the C-terminal portion (G17-V20) adopts a fold that exists in several different orientations. In addition, the ECEPP/3 flexible docking package was used to assess the conformational variability of the ligand and surrounding 60D-insertion loop of thrombin. Amino acid residues (17-20) of the peptide interact with a region of the enzyme that exhibits broad specificity, with a preferred direction between the 60D-insertion loop and Pro37 of thrombin.