AbstractA method is proposed to perform computer simulations of protein dynamics in the long‐time regime. The method is based upon a Monte Carlo technique. The only molecular degrees of freedom considered are bond rotations. All other degrees of freedom including the amide plane torsions are kept rigid. These constraints approximately account for all interactions related to chemical bonding. An individual Monte Carlo step adopts the Go and Scheraga algorithm where local conformational changes in a small window of the protein backbone are performed. By using correlated rotations, the conformation of residues outside the window remains invariant. To test the reliability of the method, the nonbonded interactions are turned off in the present application. Exact statistical averages are compared with values obtained from data of computer simulation involving 2 × 106 scans of the window along the protein backbone. Time is related to the number of scans of the window along the protein backbone. End‐to‐end distance autocorrelation functions decay to 1/e of its initial value in about 103–104 scans of the window algorithm. Time decay follows a stretched exponential Kohlrausch decay law. © 1993 John Wiley & Sons, Inc.