A method is described for calculating the reaction rate in globular proteins of activated processes such as ligand binding or enzymatic catalysis. The method is based on the determination of the probability that the system is in the transition state and of the magnitude of the reactive flux for transition-state systems. An "umbrella sampling" simulation procedure is outlined for evaluating the transition-state probability. The reactive flux is obtained from an approach described previously for calculating the dynamics of transition-state trajectories. An application to the rotational isomerization of an aromatic ring in the bovine pancreatic trypsin inhibitor is presented. The results demonstrate the feasibility of calculating rate constants for reactions in proteins and point to the importance of solvent effects for reactions that occur near the protein surface.