In this paper, we study nonlinear stability of the 3D plane Couette flow (y, 0,0) at high Reynolds number Re in a finite channel T x [-1, 1] x T. It is well known that the plane Couette flow is linearly stable for any Reynolds number. However, it could become nonlinearly unstable and transition to turbulence for small but finite perturbations at high Reynolds number. This is so-called Sommerfeld paradox. One resolution of this paradox is to study the transition threshold problem, which is concerned with how much disturbance will lead to the instability of the flow and the dependence of disturbance on the Reynolds number. This work shows that if the initial velocity v0 satisfies Ilv(0) - (y, 0, 0)IIH2 <= c(0)Re(-1) for some c(0) > 0 independent of Re, then the solution of the 3D Navier-Stokes equations is global in time and does not transit away from the Couette flow in the L infinity sense, and 1 rapidly converges to a streak solution for t >> Re-1/3 due to the mixing -enhanced dissipation effect. This result confirms the threshold result obtained by Chapman via an asymptotic analysis(JFM 2002). The most key ingredient of the proof is the resolvent estimates for the full linearized 3D Navier-Stokes system around the flow (V (y, z), 0, 0), where V (y, z) is a small perturbation(but independent of Re) of y.