Geometric and energy parameters in lysine-retinal chromophores

The parameters used in the computer program ECEPP (Empirical Conformational Energy Program for Peptides) have been expanded to cover some key elements in retinal-containing proteins. These elements are 'all-trans retinal lysine with unprotonated imine', 'all-trans retinal lysine with protonated imine', '13-cis retinal lysine with unprotonated imine' and '13-cis retinal lysine with protonated imine' respectively. The geometric parameters of these four new 'amino acid residues' were derived by optimizing their molecular structures with the AM1 Hamiltonian included in MOPAC (Molecular Orbital PACkage), and their partial atomic charges were determined with a CNDO/2 (Complete Neglect of Differential Overlap) calculation. The parameters for nonbonded interactions and torsional potentials were obtained from the existing ECEPP parameters through a logical extension. The augmented ECEPP system thus obtained can be employed to investigate the conformation of bacteriorhodopsin and its proton-pumping mechanism from an energetic point of view. The computer modeling study on bacteriorhodopsin and other seven-helix membrane proteins, e.g. serotonin receptor and dopamine receptor, is under way in the Upjohn Laboratories.