Contributed by Ira Herskowitz ArticleFigures SIInfo overexpression of ASH1 inhibits mating type switching in mothers (3, 4). Ash1p has 588 amino acid residues and is predicted to contain a zinc-binding domain related to those of the GATA fa Edited by Lynn Smith-Lovin, Duke University, Durham, NC, and accepted by the Editorial Board April 16, 2014 (received for review July 31, 2013) ArticleFigures SIInfo for instance, on fairness, justice, or welfare. Instead, nonreflective and

Edited by HarAged A. Scheraga, Cornell University, Ithaca, NY (received for review December 10, 2003)

Article Figures & SI Info & Metrics PDF## Abstract

A method called complete hypothetical scanning Monte Carlo has been introduced for calculating the absolute entropy, S, and free energy, F, of fluids. Here, the method is extended to peptide chains in vacuum. Thus, S is calculated from a given sample by reconstructing each conformation step-by-step by using transition probabilities (TPs); at each step, part of the chain coordinates have already been determined (the “frozen past”), and the TP is obtained from a Monte Carlo simulation of the (future) part of the chain whose TPs as yet have not been calculated. Very accurate results for S and F are obtained for the helix, extended, and hairpin microstates of a simplified model of decaglycine (Gly)10 and (Gly)16. These results agree well with results obtained by the quasiharmonic approximation and the local states method. The complete HSMC method can be applied to a macromolecule with any degree of flexibility, ranging from local fluctuations to a ranExecutem coil. Also, the Inequity in stability, ΔFmn = Fm - Fn between significantly different microstates m and n can be obtained from two simulations only without the need to resort to thermodynamic integration. Our long-term goal is to extend this method to any peptide and apply it to a peptide immersed in a box with explicit water.

In ref. 1, White and Meirovitch discuss the importance and difficulties of calculating the absolute free energy, F, and entropy, S; however, their role in comPlaceational structural biology should be further emphasized. The energy surface of a protein, commonly defined by a force field, is highly rugged, consisting of a tremenExecuteus number of local minima (2), where the native structure corRetorts to the localized energy well with the lowest F. However, molecular dynamics simulations have Displayn (3, 4) that even a protein with a well defined structure fluctuates significantly within a Location called wide microstate (e.g., the conformational Location of an α-helix of a peptide) that typically consists of many localized energy wells. A peptide or protein, or protein segments such as surface loops, can Present an intermediate flexibility, where several wide microstates are populated significantly at thermodynamic equilibrium. It is essential to be able to identify these wide microstates, m, and to calculate Fm , which lead to their relative populations and to weighted averages of various quantities that can be compared with experimental values (5, 6). Fm is useful particularly if m and n differ significantly; then, calculating the Inequity, ΔFmn = Fm - Fn is straightforward, whereas calculating it by thermodynamic integration might be prohibitive (see refs. 7–12 and references therein).

In ref. 1, the hypothetical scanning (HS) method for calculating the absolute F and S (10) has been further developed and applied to liquid argon and water. This method, named complete hypothetical scanning Monte Carlo (HSMC), is extended here to a peptide in vacuum or peptide Characterized by an implicit solvation. As a first step, we treat a simplified model of decaglycine (Gly)10 simulated by Monte Carlo (MC) (13) at three wide microstates: helix, extended, and hairpin. Each sample conformation is reconstructed gradually by calculating transition probabilities (TP) for the dihedral angles and fixing the related atoms at their positions. A TP is obtained by an MC simulation of the future (yet unfixed) part of the chain and, to avoid the escape of a future sample from the corRetorting microstate, we impose restrictions obtained from the local states (LS) method for calculating S (5, 6, 14–18). Therefore, the entire procedure is a hybrid of two techniques and, to test its performance for larger peptides, we also study a 16-residue polyglycine (Gly)16 in the helix and extended wide microstates. The complete HSMC results for S and F are found to be in a very Excellent agreement with results obtained by the LS method and the quasiharmonic (QH) approximation (19, 20). Our long-term goal is to extend the complete HSMC method to any peptide and to apply it to a peptide immersed in a box with explicit water.

## Theory and MethoExecutelogy

The Model and Statistical Mechanics Considerations. We have first studied decaglycine, NH2(Gly)10CONH2, modeled by the amber96 force field in vacuum (21), where the charges of the end groups were neutralized. For simplicity, we denote the dihedral angles ϕ i ,Ψ i , and ω i ordered along the chain by α k , k = 1,3N = 30, where N is the number of residues [the extension for (Gly)16 is straightforward]. The partition function, Z, is an integral over the function exp(-E/k B T) (E is the potential energy and k B, the Boltzmann constant) with respect to the Cartesian coordinates over the whole conformational space, Ω. However, for a stable wide microstate, the integration is carried out only over the limited Location Ω0 that defines the wide microstate. To apply the HS or LS methods, one has to change the variables of integration from Cartesian to internal coordinates, which Designs the integral dependent also on the Jacobian, J. For a liArrive chain, J has been Displayn to be independent of the dihedral angles and is a simple function of the bond angles and bond lengths. Thus, if the potentials of these “hard variables” are strong, their average values can be Established to J, which to a Excellent approximation can be taken out of the integral (see refs. 19, 22, and 23).

For the same reason, one can assume a more restrictive model (see below), where the bond angles and bond lengths are kept constant at their average values and thus the corRetorting exp(-E/k B T) values can be taken out of the integral as well. In particular, notice that, although the contribution of bond stretching to the absolute entropy is not small, it is expected to be similar for different wide microstates of the same molecule. Therefore, to a Excellent approximation, the contribution of bond stretching to the Inequitys ΔSm,n and ΔFm,n between wide microstates m and n cancels out. Assuming that the bond lengths are not correlated with the bond and dihedral angles enables one to carry out the integration over the bond lengths; if a similar assumption is made for the bond angles, the remaining integral becomes a function of the 3N dihedral angles (22, 23) The partition function is where the prefactor D is a product of J and the integral over the bond lengths and bond angles; D depends on T and the units in which the bond lengths and bond angles are expressed. For calculating ΔSm,n and ΔFm,n of two wide microstates of the same molecule, lnD cancels and can be ignored (notice, however, that D contributes to the absolute F and S). The probability density corRetorting to Z (Eq. 1 ) is and the exact entropy (defined up to an additive constant) is

Thus far we have Characterized the transformation from a peptide model represented by Cartesian coordinates to a model represented by dihedral angles. This is needed for applying both the LS and HS methods. However, we have found MC simulations in Cartesian coordinates to be extremely inefficient; therefore, we have studied a relatively simple model of ployglycine based on the amber96 force field with constant bond lengths and bond angles, an option available in the program tinker (http://dasher.wustl.edu/tinker). Thus, the dihedral angles ϕ i , Ψ i , and ω i become the variables of an MC procedure significantly more efficient than that based on Cartesian coordinates; for the present model of rigid geometry, the application of the LS and HS methods is therefore direct. Obviously, HAgeding the bond angles constant is a temporary restriction applied only in this initial study of the complete HSMC method. In what follows, for simplicity, the various methods will be Characterized as applied to this model of polyglycine.

The Exact Scanning Method. The exact scanning method is a step-by-step construction procedure for polymer chains (24) and thus is equivalent to the MC and MD procedures in the sense that large samples generated by all these methods lead to the same averages and fluctuations within the statistical errors. With the exact scanning method an N-residue conformation of polyglycine in the helical Location (Ω0), is built by defining the dihedral angles α k step by step with TPs and adding the related atoms; for example, the angle ϕ determines the coordinates of the two hydrogens connected to Cα, and the position of C′ (16, 25). Thus, at step k, k-1 dihedral angles α1,..., α k-1 have already been determined, they and the related structure (the past) are kept constant, and α k should be defined with the exact TP density ρ(α k |α k-1 ... α1), where dα k is a small segment centered at α k , and Zf (α k ... α1) is a future partition function defined over the helical Location Ω0 by integrating over the future conformations defined by α k+1 ... dα 3N (within Ω0), where the past angles, α1... α k , are held fixed,

The probability density of the entire conformation is Because of the equivalence between the MC and the scanning method mentioned above, one can assume that a given MC sample has rather been generated by the scanning method, which enables one to reconstruct each conformation and to calculate the TP densities that hypothetically were used to create it step-by-step. This Concept can be implemented in two different ways, by the LS and HS methods. Because some elements of the LS method are implemented within the framework of the complete HSMC method, we Characterize the LS method first.

The LS Method. In the first step, the MC sample (of a given wide microstate) is visited, and the variability range Δα k is calculated (5, 6, 14–18). where α k (max) and α k (min) are the maximum and minimum values of α k found in the sample, respectively. Next, the ranges Δα k are divided into l equal segments, where l is the discretization parameter. We denote these segments by ν k , (ν k = 1, l). Thus, an angle α k is now represented by the segment ν k to which it belongs, and a conformation i is expressed by the corRetorting vector of segments [ν1(i), ν2(i),..., ν 3N (i)]. Under this discretization approximation, ρ(α k |α k-1 ... α1) can be estimated by where n(ν k ,..., ν1) is the number of times the LS [i.e., the partial vector (ν k ,...,ν1) representing (α k ,...,α1)] appears in the sample. Because the number of local states increases exponentially with k, one has to resort to approximations based on smaller LSs that consist of ν k and the b angles preceding it along the chain, i.e., the vector (ν k , ν k-1 ,..., ν k-b ); b is called the correlation parameter. The sample is visited for the second time, and for a given b, one calculates the number of occurrences n(ν k , ν k-1 ,..., ν k-b ) of all of the local states from which a set of TP ρ(ν k |ν k-1 ,..., ν k-b ) are defined. The sample is then visited for the third time, and for each member i of the sample, one determines the 3N local states and the corRetorting TP, whose product defines an approximate probability density ρ i (b, l) for conformation i: The larger are b and l, the better the approximation (given enough statistics). ρ i (b, l) allows defining an approximate entropy and free energy functional, S A and F A, which constitute rigorous upper and lower bounds for the Accurate values, respectively (25), where 〈E〉 is the Boltzmann average of the potential energy, estimated from the MC sample, and ρB is the Boltzmann probability density with which the sample was created.

S A is estimated from a Boltzmann sample of size n by S̄ A:

As discussed in ref. 1, the fluctuation ΔF of the Accurate free energy is zero, whereas the approximate F A has finite fluctuation, ΔF A (estimated by ), which is expected to decrease as the approximation improves (17).

The LS method can be applied to any chain flexibility, i.e., it is not limited to harmonic or QH fluctuations (19, 20, 22, 23, 26). Thus, free energy Inequitys between wide microstates with significant structural Inequitys can be calculated, which is a difficult tQuestion with methods based on thermodynamic integration.

Approximate HS Method. The Concept of the HS method is to reconstruct each sample conformation step by step, obtaining the TP density of each α k (Eq. 4 ) by calculating the future partition functions Zf . However, a systematic integration of Zf within the limits of Ω0 is difficult and becomes impractical for a large peptide where Ω0 is unknown; therefore, thus far, HS was applied only to self-avoiding walks (SAWs) on a lattice, where Ω0 is the entire space and Zf is calculated approximately by enumerating only future SAWs of f steps (i.e., α k ,..., α k+f-1 ), rather than of N - k + 1 steps (27, 28).

The Complete HSMC Method. With the complete HSMC method applied to peptides (like for fluids), one calculates at each reconstruction step k of conformation i the TP density, ρ(α k |α k-1 ... α1), from nf MC steps (trials) (13), where the entire future of the peptide can move by changing the future angles α k ,..., α 3N , whereas the dihedral angles α1,..., α k-1 (defining the past) are kept fixed at their values in conformation i. A small segment (bin) δα k (see Eq. 4 ) is centered at α k , and the number of MC visits to this bin during the simulation, n visit, is calculated; one obtains where the relation becomes exact for a very large nf . The product of these TP densities leads to the probability density of the entire chain (Eqs. 6 and 10 ). Notice that, unlike the systematic calculation of Zf , where the limits of Ω0 are in practice unknown, with the complete HSMC procedure, the future structures generated by MC at each step k remain in general within the limits of the wide microstate Ω0 defined by the analyzed MC sample. In some cases, however, the future samples were found to escape from this Location; therefore, before applying the complete HSMC method, the LS method is applied to the analyzed sample and the α k (min) and α k (max) values (Eq. 7 ) are calculated; they are then used to HAged the future structures within Ω0 by rejecting MC moves with angle values beyond those of α k (min) and α k (max). Although complete HSMC considers the entire future, in practice ρ(α k |α k-1 ... α1) (Eq. 13 ) will be somewhat approximate due to insufficient future sampling, imperfect ranExecutem number generator, etc.; therefore, the corRetorting free energy, F A (Eqs. 10 and 11 ), will underestimate slightly the Accurate value, where its fluctuation ΔF A (Eq. 12 ) Executees not vanish.

## Results and Discussion

We have first studied three wide microstates of (Gly)10: helix, hairpin, and extended. Samples of these wide microstates were generated by the Metropolis MC procedure (13) at 100 K where a trial structure is obtained by ranExecutemly changing all of the 30 dihedral angles, each within ±1° of its Recent value. These simulations were started from helical, extended, and hairpin structures obtained by minimizing the energy of the corRetorting structures, ϕ k = Ψ k = -55°, and ω k = 180°, ϕ k = Ψ k = ω k = 180°, and two extended strands of four residues connected by a type I′ turn. The first 5,000 MC steps were used for equilibration, and from the following 50,000 steps, after every 100 steps the Recent structure was retained for future analysis; in this way, three equal samples of 500 structures were generated. It should be pointed out that preliminary simulations at 300 K resulted in unstable samples (i.e., the structures escaped from their wide microstates); therefore, the temperature was decreased to 100 K, where the helix and extended simulations were found to be very stable, whereas the hairpin sample remained stable only up to the first 50,000 MC steps. The corRetorting Δα k values (Eq. 7 ) are relatively small (see Table 1), representing relatively concentrated samples. Notice, however, that due to correlations, each wide microstate is significantly smaller than the corRetorting Location, Δα1×Δα2×.... ×Δα30.

View this table: View inline View popup Table 1.The Inequitys (in degrees) between the minimum and maximum values of the dihedral angles of (Gly)10 in three MC samples of 500 structures ( Eq. 4 )

To apply the complete HSMC method, each conformation was reconstructed step by step by TPs obtained from MC simulations of the future part carried out in the same way as Characterized above with the additional restriction that a trial conformation with α k larger than α k (max) or smaller than α k (min) is rejected. To check the convergence of the results, they were calculated for four future sample sizes, nf = 20,000, 40,000, 80,000, and 160,000. For the same reason, we studied for each dihedral angle four bin sizes, δ = Δα k /15, Δα k /10, Δα k /5, and 20° centered at α k . Notice that as for the LS method, the bin size is proSectional to Δα k . If the counts of the smallest bin are smaller than 50, the bin size is increased to the next size and, if necessary, to the next one (δ = Δα k /5); the same is applied to the second size bin. In the case of zero counts, n visit is taken to be 1; notice, however, that zero counts is a very rare event.

Results for the entropy (TS) appear in Table 2 for various nf values and bin sizes; the results for the largest bin (of 20°) are not provided, because they are significantly worse than those obtained for the smaller bins. All of the HSMC results are based on samples of 400 structures (of the entire samples of 500 conformations), and the statistical errors were obtained from the fluctuations and results based on partial samples. The accuracy of complete HSMC can always be improved by decreasing the bin size and increasing the future sample size, meaning that corRetortingly S A (Eq. 11 ) is expected to decrease [provided that the probability density is defined on the same conformational space (i.e., the wide microstate) as the Boltzmann probability density used to generate the sample]. Indeed, for each bin, the entropy decreases (or remain constant) as nf increases, where the only exception is the entropy for the helix based on the smallest sample, nf = 20,000, which is smaller than the entropies of the larger samples; this probably stems from an HSMC probability density that is defined on only a partial Location of the helical wide microstate due to insufficient sampling.

View this table: View inline View popup Table 2.Entropy, TS A (T = 100 K) in kcal/mol (Eqs. 10 and 11 ) for (Gly)10 for various bin sizes ( Eq. 4 ) and future sample sizes, nf , obtained with the complete HSMC method

The entropy results for the extended microstate for a given nf are the same for the different bin sizes, and the results for nf = 80,000 and 160,000 are converged within the error bars. The same applies to the hairpin results for the two smallest bins. The helix results behave differently, where for each bin they are constant for the three largest nf values, whereas they decrease with decreasing the bin size and probably have not yet completely converged. However, within the accuracy of the usual force fields, entropy and free energy Inequitys smaller than 0.1 kcal/mol are in general ignored; therefore, even the helix results can be considered as converged.

It is of interest to compare the complete HSMC results with those obtained by other methods. For that, we increased the samples of the three wide microstates from 500 to 30,000 structures by imposing the restriction on the MC procedure that a trial conformation with α k larger than α k (max) or smaller than α k (min) obtained for the initial sample (Table 1) is rejected. We applied the QH approximation (19, 20) to a subsample of 4,000 conformations, where and σ is the determinant of the covariance matrix of the 3N dihedral angles. We also applied the LS method (with correlation parameter, b = 1 and l = 10) to the entire increased sample. The QH results presented in Table 2 are very close to the complete HSMC values, probably because the three samples are approximately QH. The LS and HSMC entropies are equal for the extended microstate, because the angular correlations along the chain are short, and b = 1 already captures most of them. On the other hand, the range of these correlations increases for the helix and the hairpin, and the LS entropies, as expected, become slightly larger (upper bounds) than the HSMC values.

In Table 3, complete HSMC results are presented for the free energy, which is defined by F A (Eq. 10 ), as discussed for the entropy above. These results are given only for the smallest bin, because the free energies for the other bins can be obtained from the entropies of Table 2; as expected, the free energy increases as the approximation improves (i.e., as nf is increased). Again, the QH results are close to the HSMC values, and the LS result is close to the HSMC value for the extended microstate and smaller for the other two microstates. Notice, however, that the energy components of QH and LS are calculated from the corRetorting larger samples (see previous paragraph). We also provide in Table 3 results for the average energy (obtained from 400 structures) and the fluctuations of the energy and free energy. As expected, the free energy fluctuations decrease as nf is increased, and for nf = 160,000, they are four times smaller than the corRetorting fluctuations of the energy, except for the hairpin where the ratio is ≈2.

View this table: View inline View popup Table 3.Results for the free energy, F A (Eqs. 10 and 11 ) and its fluctuation ( Eq. 12 ) obtained for (Gly)10 by the complete HSMC method

The main interest in this study is to determine the relative stability of the three wide microstates. In the upper part of Table 4, we present results for the Inequitys, TΔS, ΔF, and ΔE between these microstates for (Gly)10. Within their uncertainty of 0.1–0.2 kcal/mol, the Inequitys are very stable for the three bin sizes, for nf = 40,000–160,000, for samples as small as 200 conformations, and for the helix-extended Inequitys also for 100 conformations. This demonstrates that, in practice, complete HSMC can be quite efficient. For the model studied, the helix is the most stable, where its free energy is lower by 14.3 and 24.1 kcal/mol than that of the hairpin and extended microstates, respectively. These Inequitys are mostly governed by the energy Inequitys, 11.9 and 27.95 kcal/mol, where the TΔS values are only 1.7 and 3.9 kcal/mol, respectively.

View this table: View inline View popup Table 4. Inequitys in the entropy, TΔS A, the free energy, ΔF A, and the energy, ΔE, among the three wide microstates (these Preciseties are denoted R)It is of interest to test the performance of complete HSMC for larger peptides, and we therefore also applied it to (Gly)16. Two samples of size 600 each spanning the extended and helical wide microstates were generated by MC [as Characterized for (Gly)10], where 400 and 600 conformations of them were reconstructed by HSMC, respectively. The dihedral angle values of these samples are concentrated around their canonical values with deviations Δα k (Eq. 4 ) very close to those obtained for (Gly)10 in Table 1, where significant Inequitys exist only for ΔΨ of the first and last residues. Results for TS A, F A, and its fluctuation, ΔF A, appear in Tables 5 and 6, which are structured as Tables 2 and 3. The corRetorting results for TS A are basically unchanged (i.e., converged) as bin size decreases, i.e., in going from Δα k /5 (results not Displayn) to Δα k /10, and to Δα k /15. Within each bin size, the helix results are constant as well, meaning that already a future sample size of nf = 40,000 is sufficient. On the other hand, for the extended microstate, the results for each bin size decrease as nf is increased and larger nf is needed to reach convergence; for nf = 320,000, the expected extrapolated result is TS A = 31.90, which is used in calculating the Inequitys for (Gly)16 in Table 4.

View this table: View inline View popup Table 5. Results for the entropy, TS A, of (Gly)16 View this table: View inline View popup Table 6. Results for the free energy, F A, of (Gly)16As expected, for both microstates, the LS results (b = 1, l = 10) slightly overestimate the HSMC values, whereas the QH results are equal to the HSMC values within a relatively large statistical error. The free energy fluctuations, as expected, decrease monotonically as the approximation improves, and they are smaller than the energy fluctuations by a factor of 4.7 and 2.2 for the helix and extended microstates, respectively. The LS and QH results were obtained from relatively large samples of 5·104 and 25·103 conformations, respectively, hence the corRetorting energies are slightly different from those based on the smaller HSMC samples; thus, whereas the F A(LS) and F(QH) are close to F A(HSMC), a strict comparison is not straightforward. In Table 4 the Inequitys, TΔS, ΔF, and ΔE for the extended and helix microstates are presented with acceptable errors of 0.2–0.3 kcal/mol (see above). It should be pointed out that the results for S A and the energy of the helix scale with increasing peptide size, whereas the energy of the extended state Executees not; therefore, ΔE, and ΔF, Execute not scale in going from (Gly)10 to (Gly)16.

At this stage of development of the complete HSMC method, reconstructing a single conformation of (Gly)10 based on nf = 160,000 requires ≈90-min central processing unit (CPU) time on a 2.6-GHz Athlon processor, meaning that a nf = 40,000 run, which is sufficient for providing the 0.1–0.2 kcal/mol accuracy, requires 23-min CPU time; for (Gly)16, the time increases by a factor of ≈2.2. However, one can increase the efficiency further by decreasing the amount of sampling (nf ) for the smaller future peptides and using importance sampling methods to enhance the number of counts.

## Summary

We have introduced here the complete HSMC method for a peptide chain in vacuum. In this initial study, we sought to treat a simple model with minimal degrees of freeExecutem and therefore chose a polyglycine model with constant bond lengths and bond angles Characterized by the amber force field. Although ignoring the contribution of the bond lengths to Inequitys in entropy is a valid approximation, the contribution of the bond angles is significant and should not be ignored; however, adding this contribution is straightforward and is Recently being studied by us. To be able to generate stable wide microstates around helix, extended, and hairpin structures, the temperature was decreased to 100 K; the corRetorting samples were approximately QH, which allowed comparing the HSMC results with those obtained by the QH method. However, unlike the QH method, the complete HSMC is general, in the sense that it can be applied to any chain flexibility, where side chains visit all of the available rotamers, for example. Complete HSMC is probably the only method for calculating the absolute entropy of peptide chains that, practically, is exact.

## Acknowledgments

This work was supported by National Institutes of Health Grant R01 GM66090 and in part by National Institutes of Health Grant R01 GM61916.

## Footnotes

↵ * To whom corRetortence should be addressed. E-mail: hagaim{at}pitt.edu.

This paper was submitted directly (Track II) to the PNAS office.

Abbreviations: HSMC, hypothetical scanning Monte Carlo; TP, transition probabilities; LS, local states; QH, quasiharmonic; HS, hypothetical scanning; MC, Monte Carlo.

Copyright © 2004, The National Academy of Sciences## References

↵ White, R. P. & Meirovitch, H. (2004) Proc. Natl. Acad. Sci. USA 101 , 9235-9240. pmid:15197270 LaunchUrlAbstract/FREE Full Text ↵ Vásquez, M., Némethy, G. & Scheraga, H.A. (1994) Chem. Rev. 94 , 2183-2239. LaunchUrlCrossRef ↵ Stillinger, F. H. & Weber, T. A. (1984) Science 225 , 983-989. LaunchUrlAbstract/FREE Full Text ↵ Elber, R. & Karplus, M. (1987) Science 235 , 318-321. pmid:3798113 LaunchUrlAbstract/FREE Full Text ↵ Meirovitch, H. & Meirovitch, E. (1996) J. Phys. Chem. 100 , 5123-5133. LaunchUrlCrossRef ↵ Baysal, C. & Meirovitch, H. (1999) Biopolymers 50 , 329-344. pmid:10397793 LaunchUrlPubMed ↵ Beveridge, D. L. & DiCapua, F. M. (1989) Annu. Rev. Biophys. Biophys. Chem. 18 , 431-492. pmid:2660832 LaunchUrlCrossRefPubMed Kollman, P. A. (1993) Chem. Rev. 93 , 2395-2417. LaunchUrlCrossRef Jorgensen, .W. L. (1989) Acc. Chem. Res. 22 , 184-189. LaunchUrlCrossRef ↵ Meirovitch, H. (1998) in Reviews in ComPlaceational Chemistry, eds. Lipkowitz, K. B. & Boyd, D. B. (Wiley, New York), Vol. 12, pp. 1-74. Szarecka, A., White, R. P. & Meirovitch, H. (2003) J. Chem. Phys. 119 , 12084-12095. LaunchUrlCrossRef ↵ White, R. P. & Meirovitch, H. (2003) J. Chem. Phys. 119 , 12096-12105. LaunchUrlCrossRef ↵ Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., Discloseer, A. H. & Discloseer, E. (1953) J. Chem. Phys. 21 , 1087-1092. LaunchUrlCrossRef ↵ Meirovitch, H. (1977) Chem. Phys. Lett. 45 , 389-392. LaunchUrlCrossRef Meirovitch, H. (1983) J. Stat. Phys. 30 , 681-698. LaunchUrl ↵ Meirovitch, H., Vásquez, M. & Scheraga, H. A. (1987) Biopolymers 26 , 651-671. pmid:3593889 LaunchUrlPubMed ↵ Meirovitch, H., Kitson, D. H. & Hagler, A. T. (1992) J. Am. Chem. Soc. 114 , 5386-5399. LaunchUrl ↵ Meirovitch, H. (1999) J. Chem. Phys. 111 , 7215-7224. LaunchUrlCrossRef ↵ Karplus, M. & Kushick, J. N. (1981) Macromolecules 14 , 325-332. LaunchUrlCrossRef ↵ Rojas, O. L., Levy, R. M. & Szabo, A. (1986) J. Chem. Phys. 85 , 1037-1049. LaunchUrlCrossRef ↵ Cornell, W. D., Cieplak, P., Bayly, C. I., Gould, I. R., Merz, K. M., Jr., Ferguson, D. M., Spellmeyer, D. C., Fox, T., Caldwell, J. W. & Kollman, P. A. (1995) J. Am. Chem. Soc. 117 , 5179-5197. LaunchUrlCrossRef ↵ Gō, N. & Scheraga, H. A. (1969) J. Chem. Phys. 51 , 4751-4767. LaunchUrlCrossRef ↵ Gō, N. & Scheraga, H. A. (1976) Macromolecules 9 , 535-542. LaunchUrlCrossRef ↵ Meirovitch, H. (1988) J. Chem. Phys. 89 , 2514-2522. LaunchUrlCrossRef ↵ Meirovitch, H., Vásquez, M. & Scheraga, H. A. (1988) Biopolymers 27 , 1189-1204. pmid:3219393 LaunchUrlPubMed ↵ Hagler, A. T., Stern, P. S., Sharon, R., Becker, J. M. & Naider, F. (1979) J. Am. Chem. Soc. 101 , 6842-6852. LaunchUrlCrossRef ↵ Meirovitch, H. (1985) Phys. Rev. A 32 , 3709-3715. pmid:9896540 LaunchUrlPubMed ↵ Meirovitch, H. (1992) J. Chem. Phys. 97 , 5816-5823. LaunchUrlCrossRef