Introduction

The combination of improved experimental capability, great advances in computer performance, and the development of new algorithms from computer science have led to quite sophisticated methods for the study of certain biomolecules, in particular of folded protein structures. One such technique, called ‘threading’, picks out small pieces of the primary structure of a protein whose structure is unknown and examines extensive databases of known protein structures to find similar pieces of primary structure. One then guesses that this piece will have the same folded structure as that in the known structure. Since pieces do not all fit together perfectly, an effective force field is used to ‘optimize’ the resultant structure, and Monte Carlo methods have already begun to play a role in this approach. (There are substantial similarities to ‘homology modeling’ approaches to the same, or similar, problems.) Of course, the certainty that the structure is correct comes primarily from comparison with experimental structure determination of crystallized proteins. One limitation is thus that not all proteins can be crystallized, and, even if they can, there is no assurance that the structure will be the same in vivo. Threading algorithms have, in some cases, been extraordinarily successful, but since they do not make use of the interactions between atoms it would be useful to complement this approach by atomistic simulations. (For an introductory overview of protein structure prediction, see Wooley and Ye (2007).) Biological molecules are extremely large and complex; moreover, they are usually surrounded by a large number of water molecules. Thus, realistic simulations that include water explicitly and take into account polarization effects are inordinately difficult. There have also been many attempts to handle this task by means of molecular dynamics simulations, but the necessity of performing very long runs of very large systems makes it extremely difficult (if not impossible) to reach equilibrium. We are thus in the quite early stages in the study of biological molecules via Monte Carlo methods, but, as we shall see, results are quite promising for the future. The field is developing rapidly, and in some ways we are passing through the same kind of maturation period as Monte Carlo enthusiasts in physics did 30–40 years ago in which simulations were, at first, not taken too seriously by experimentalists. This will surely change.