Molecular dynamics of rhodopsin and free opsin: computer simulation

Quantum Mechanical and Molecular Mechanics Modeling of Membrane-Embedded Rhodopsins

Computational chemistry provides versatile methods for studying the properties and functioning of biological systems at different levels of precision and at different time scales. The aim of this article is to review the computational methodologies that are applicable to rhodopsins as archetypes for photoactive membrane proteins that are of great importance both in nature and in modern technologies. For each class of computational techniques, from methods that use quantum mechanics for simulating rhodopsin photophysics to less-accurate coarse-grained methodologies used for long-scale protein dynamics, we consider possible applications and the main directions for improvement.

Femtosecond spectroscopic study of photochromic reactions of bacteriorhodopsin and visual rhodopsin

Photochromic ultrafast reactions of bacteriorhodopsin (H. salinarum) and bovine rhodopsin were conducted with a femtosecond two-pump probe pulse setup with the time resolution of 20-25fs. The dynamics of the forward and reverse photochemical reactions for both retinal-containing proteins was compared. It is demonstrated that when retinal-containing proteins are excited by femtosecond pulses, dynamics pattern of the vibrational coherent wave packets in the course of the reaction is different for bacteriorhodopsin and visual rhodopsin. As shown in these studies, the low-frequencies that form a wave packets experimentally observed in the dynamics of primary products formation as a result of retinal photoisomerization have different intensities and are clearer for bovine rhodopsin. Photo-reversible reactions for both retinal proteins were performed from the stage of the relatively stable photointermediates that appear within 3-5ps after the light pulse impact. It is demonstrated that the efficiency of the reverse phototransition K-form→bacteriorhodopsin is almost five-fold higher than that of the Batho-intermediate→visual rhodopsin phototransition. The results obtained indicate that in the course of evolution the intramolecular mechanism of the chromophore-protein interaction in visual rhodopsin becomes more perfect and specific. The decrease in the probability of the reverse chromophore photoisomerization (all-trans→11-cis retinal) in primary photo-induced rhodopsin products causes an increase in the efficiency of the photoreception process.

Model of Abnormal Chromophore-Protein Interaction for Е181К Rhodopsin Mutation: Computer Molecular Dynamics Study

The interaction of the 11-cis-retinal chromophore with the surrounding amino acid residues in the chromophore center of the rhodopsin protein has been investigated for the Е181К mutant form using molecular dynamics simulation. A comparative analysis of the arrangement of the amino acid residues in the chromophore center has been performed for both wild (native) and mutant rhodopsins. It is shown that for the Е181К mutant rhodopsin there is no proper binding of 11-cis-retinal with the surrounding amino acid residues. The distortion of the conformation states in the mutant rhodopsin molecule takes place in both the chromophore center and cytoplasmic domain. Our simulations suggest that a stable covalent linkage of 11-cis-retinal with the protein part (viz. opsin) of the rhodopsin molecule will not form. This, on the other hand, implies that the protein's active site in the cytoplasmic domain, which is responsible for the G-protein binding (so-called transducin), may not be completely blocked.Based on our molecular simulation data, we discuss the possible correlation between retinitis pigmentosa pathogenesis and the structural and functional properties of the rhodopsin protein.

Conserved water-mediated hydrogen bond network between TM-I, -II, -VI, and -VII in 7TM receptor activation

Five highly conserved polar residues connected by a number of structural water molecules together with two rotamer micro-switches, TrpVI:13 and TyrVII:20, constitute an extended hydrogen bond network between the intracellular segments of TM-I, -II, -VI, and -VII of 7TM receptors. Molecular dynamics simulations showed that, although the fewer water molecules in rhodopsin were relatively movable, the hydrogen bond network of the beta2-adrenergic receptor was fully loaded with water molecules that were surprisingly immobilized between the two rotamer switches, both apparently being in their closed conformation. Manipulations of the rotamer state of TyrVII:20 and TrpVI:13 demonstrated that these residues served as gates for the water molecules at the intracellular and extracellular ends of the hydrogen bond network, respectively. TrpVI:13 at the bottom of the main ligand-binding pocket was shown to apparently function as a catching trap for water molecules. Mutational analysis of the beta2-adrenergic receptor demonstrated that the highly conserved polar residues of the hydrogen bond network were all important for receptor signaling but served different functions, some dampening constitutive activity (AsnI:18, AspII:10, and AsnVII:13), whereas others (AsnVII:12 and AsnVII:16) located one helical turn apart and sharing a water molecule were shown to be essential for agonist-induced signaling. It is concluded that the conserved water hydrogen bond network of 7TM receptors constitutes an extended allosteric interface between the transmembrane segments being of crucial importance for receptor signaling and that part of the function of the rotamer micro-switches, TyrVII:20 and TrpVI:13, is to gate or trap the water molecules.