Efficient calculation of two‐dimensional adiabatic and free energy maps: Application to the isomerization of the C13C14 and C15N16 bonds in the retinal of bacteriorhodopsin

Accurate calculation of potential energy and free‐energy profiles along reaction coordinates of biological processes such as enzymatic reactions or conformational changes is fundamental to the obtention of theoretical insight into protein function. We describe here the practical implementation of the Automatic Map Refinement Procedure (AMRP) and two‐dimensional Weighted Histogram Analysis Method (WHAM) for efficient computation of adiabatic potential energy and free‐energy maps, respectively. Methods for efficiently sampling configuration space with high‐energy barriers and for removing hysteresis in the case of periodic reaction coordinates are presented. The application of these techniques to the isomerization of the C13C14 and C15N16 bonds in the retinal of bacteriorhodopsin is described. In dark‐adapted bacteriorhodopsin (bR), the retinal moiety exists in two conformers, all‐trans and (13,15)cis, with the latter making ≃67% of the population. This experimental free energy difference is reproduced here to within kBT. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1644–1658, 1999

[1]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[2]  H. Berendsen,et al.  ALGORITHMS FOR MACROMOLECULAR DYNAMICS AND CONSTRAINT DYNAMICS , 1977 .

[3]  W. Stoeckenius,et al.  Retinal isomer ratio in dark-adapted purple membrane and bacteriorhodopsin monomers. , 1989, Biochemistry.

[4]  Hiroshi Akima,et al.  A Method of Bivariate Interpolation and Smooth Surface Fitting for Irregularly Distributed Data Points , 1978, TOMS.

[5]  G. Torrie,et al.  Monte Carlo free energy estimates using non-Boltzmann sampling: Application to the sub-critical Lennard-Jones fluid , 1974 .

[6]  C. Brooks,et al.  Constant-temperature free energy surfaces for physical and chemical processes , 1993 .

[7]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[8]  Klaus Schulten,et al.  Quantum Chemistry: Molecular Dynamics Study of the Dark-Adaptation Process in Bacteriorhodopsin , 1996 .

[9]  M. El-Sayed,et al.  RETINAL ISOMER COMPOSITION IN SOME BACTERIORHODOPSIN MUTANTS UNDER LIGHT AND DARK ADAPTATION CONDITIONS , 1995 .

[10]  Jeremy C. Smith,et al.  Thermodynamic stability of water molecules in the bacteriorhodopsin proton channel: a molecular dynamics free energy perturbation study. , 1996, Biophysical journal.

[11]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[12]  Charles L. Brooks,et al.  Promotion of helix formation in peptides dissolved in alcohol and water-alcohol mixtures , 1993 .

[13]  Jeremy C. Smith,et al.  Ab initio quantum chemical analysis of Schiff base-water interactions in bacteriorhodopsin , 1993 .

[14]  A. Schulte,et al.  Equilibrium Composition of Retinal Isomers in Dark-Adapted Bacteriorhodopsin and Effect of High Pressure Probed by Near-Infrared Raman Spectroscopy , 1995 .

[15]  B. Roux.  The calculation of the potential of mean force using computer simulations , 1995 .

[16]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[17]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[18]  Jeremy C. Smith,et al.  Simulation analysis of the retinal conformational equilibrium in dark-adapted bacteriorhodopsin. , 1999, Biophysical journal.

[19]  S. O. Smith,et al.  Dark-adapted bacteriorhodopsin contains 13-cis, 15-syn and all-trans, 15-anti retinal Schiff bases. , 1984, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Jeremy C. Smith,et al.  Quantum Chemical and Free Energy Simulation Analysis of Retinal Conformational Energetics , 1997, J. Chem. Inf. Comput. Sci..

[21]  D. Oesterhelt,et al.  Rhodopsin-like protein from the purple membrane of Halobacterium halobium. , 1971, Nature: New biology.

[22]  Jeremy C. Smith,et al.  Functional interactions in bacteriorhodopsin: a theoretical analysis of retinal hydrogen bonding with water. , 1995, Biophysical journal.

[23]  R Henderson,et al.  Electron-crystallographic refinement of the structure of bacteriorhodopsin. , 1996, Journal of molecular biology.

[24]  E. Pebay-Peyroula,et al.  X-ray structure of bacteriorhodopsin at 2.5 angstroms from microcrystals grown in lipidic cubic phases. , 1997, Science.