The specific interactions between base pairs and amino acids were studied by the multicanonical Monte Carlo method. We sampled numerous interaction configurations and side‐chain conformations of the amino acid by the multicanonical algorithm, and calculated the free energies of the interactions between an amino acid at given Cα positions and a fixed base pair. The contour maps of free energy derived from this calculation represent the preferred Cα position of the amino acid around the base, and these maps of various combinations of bases and amino acids can be used to quantify the specificity of intrinsic base–amino acid interactions. Similarly, enthalpy and entropy maps will provide further details of the specific interactions. We have also calculated the free‐energy map of the orientations of the CαCβ bond vector, which indicates the preferential orientation of the amino acid against the base. We compared the results obtained by the multicanonical method with those of the exhaustive sampling and canonical Monte Carlo methods. The free‐energy map of the base–amino acid interaction obtained by the multicanonical simulation method was nearly identical to the accurate result derived from the exhaustive sampling method. This indicates that a single multicanonical Monte Carlo simulation can produce an accurate free‐energy map. Multicanonical Monte Carlo sampling produced free‐energy maps that were more accurate than those produced by canonical Monte Carlo sampling. Thus, the multicanonical Monte Carlo method can serve as a powerful tool for estimating the free‐energy landscape of base–amino acid interactions and for elucidating the mechanism by which amino acids of proteins recognize particular DNA base pairs. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 954–962, 2000
[1]
H. Kono,et al.
Structure‐based prediction of DNA target sites by regulatory proteins
,
1999,
Proteins.
[2]
Yuko Okamoto,et al.
Prediction of peptide conformation by multicanonical algorithm: New approach to the multiple‐minima problem
,
1993,
J. Comput. Chem..
[3]
Y. Okamoto,et al.
Finite-size scaling of helix–coil transitions in poly-alanine studied by multicanonical simulations
,
1998
.
[4]
N. Metropolis,et al.
Equation of State Calculations by Fast Computing Machines
,
1953,
Resonance.
[5]
P. Kollman,et al.
A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules
,
1995
.
[6]
R. Sauer,et al.
Transcription factors: structural families and principles of DNA recognition.
,
1992,
Annual review of biochemistry.
[7]
M. Michael Gromiha,et al.
Free-Energy Maps of Base−Amino Acid Interactions for DNA−Protein Recognition
,
1999
.
[8]
Berg,et al.
Multicanonical ensemble: A new approach to simulate first-order phase transitions.
,
1992,
Physical review letters.
[9]
Yuko Okamoto,et al.
The Generalized-Ensemble Approach for Protein Folding Simulations
,
1999
.
[10]
B. Berg,et al.
Multicanonical algorithms for first order phase transitions
,
1991
.
[11]
Yuko Okamoto,et al.
Comparative Study of Multicanonical and Simulated Annealing Algorithms in the Protein Folding Problem
,
1994
.