Quantum Chemical and Free Energy Simulation Analysis of Retinal Conformational Energetics

Ab initio quantum chemical and free energy molecular dynamics calculations are performed to examine energy differences between the all-trans and (13,15)cis conformers of retinal, i.e., those populated in the dark-adapted state of bacteriorhodopsin. The quantum chemical results are used to derive an empirical force field that is used to calculate an adiabatic potential energy map for rotation about the bonds concerned. The same potential function is used in restrained molecular dynamics free energy calculations, with and without umbrella sampling. The simulation model yields a free energy of the (13,15)cis isomer ∼2.1 kcal/mol higher than that of the all-trans species.

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

[2]  S. Seltzer.  MNDO barrier heights for catalyzed bicycle-pedal, hula-twist, and ordinary cis-trans isomerizations of protonated retinal Schiff base , 1987 .

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

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

[5]  K. Schulten,et al.  Molecular dynamics study of the proton pump cycle of bacteriorhodopsin. , 1993, Biochemistry.

[6]  K Schulten,et al.  Molecular dynamics study of bacteriorhodopsin and artificial pigments. , 1994, Biochemistry.

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

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

[9]  Richard Henderson,et al.  A model for the structure of bacteriorhodopsin based on high resolution electron cryomicroscopy , 1990 .

[10]  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.

[11]  Arieh Warshel,et al.  Bicycle-pedal model for the first step in the vision process , 1976, Nature.

[12]  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 .

[13]  R. Henderson,et al.  Model for the structure of bacteriorhodopsin based on high-resolution electron cryo-microscopy. , 1990, Journal of molecular biology.

[14]  K. Schulten,et al.  Molecular dynamics study of the 13-cis form (bR548) of bacteriorhodopsin and its photocycle. , 1995, Biophysical journal.

[15]  S. Seltzer.  MNDO Barrier Heights for Catalyzed Bicycle-Pedal, Hula-Twist, and Ordinary cis-trans Isomerizations of Protonated Retinal Schiff Base. , 1987 .

[16]  M. Rees,et al.  The nucleus of Centaurus A , 1976, Nature.

[17]  A. Schulte,et al.  High-pressure near-infrared Raman spectroscopy of bacteriorhodopsin light to dark adaptation. , 1995, Biophysical journal.

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

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

[20]  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.

[21]  Molecular dynamics study of the early intermediates in the bacteriorhodopsin photocycle , 1995 .

[22]  Structure and dynamics of bacteriorhodopsin , 1993, FEBS letters.

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

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