Molecular mechanism of spontaneous pigment activation in retinal cones

Spontaneous current and voltage fluctuations (dark noise) in the photoreceptor cells of the retina limit the ability of the visual system to detect dim light. We recorded the dark current noise of individual salamander L cones. Previous work showed that the dark noise in these cells arises from thermal activation of the visual pigment. From the temperature dependence of the rate of occurrence of elementary noise events, we found an Arrhenius activation energy E(a) of 25 +/- 7 kcal/mol (mean +/- SD). This E(a) is similar to that reported for the thermal isomerization of 11-cis retinal in solution, suggesting that the cone pigment noise results from isomerization of the retinal chromophore. E(a) for the cone noise is similar to that previously reported for the "photon-like" noise of rods, but the preexponential factor is five orders of magnitude higher. To test the hypothesis that thermal isomerization can only occur in molecules whose Schiff base linkage is unprotonated, we changed the pH of the solution bathing the cone outer segment. This had little effect on the rate of occurrence of elementary noise events. The rate was also unchanged when the cone was exposed to Ringer solution made up from heavy water, whose solvent isotope effect should reduce the probability, that the Schiff base nitrogen is naked.     

Spectral tuning in salamander visual pigments studied with dihydroretinal chromophores.

In visual pigments, opsin proteins regulate the spectral absorption of a retinal chromophore by mechanisms that change the energy level of the excited electronic state relative to the ground state. We have studied these mechanisms by using photocurrent recording to measure the spectral sensitivities of individual red rods and red (long-wavelength-sensitive) and blue (short-wavelength-sensitive) cones of salamander before and after replacing the native 3-dehydro 11-cis retinal chromophore with retinal analogs: 11-cis retinal, 3-dehydro 9-cis retinal, 9-cis retinal, and 5,6-dihydro 9-cis retinal. The protonated Schiff's bases of analogs with unsaturated bonds in the ring had broader spectra than the same chromophores bound to opsins. Saturation of the bonds in the ring reduced the spectral bandwidths of the protonated Schiff's bases and the opsin-bound chromophores and made them similar to each other. This indicates that torsion of the ring produces spectral broadening and that torsion is limited by opsin. Saturating the 5,6 double bond in retinal reduced the perturbation of the chromophore by opsin in red and in blue cones but not in red rods. Thus an interaction between opsin and the chromophoric ring shifts the spectral maxima of the red and blue cone pigments, but not that of the red rod pigment.     

An ab initio study of valence isomerization in the HOCl-HClO system

Hartree-Fock level and post Hartree-Fock level molecular orbital calculations have been completed for HOCl and its valence isomer, HClO. Ground state geometries have been determined for each molecule. The energy change for the reaction HOCl→HClO is estimated to be 67±5 kcal/mol endothermic and the activation energy for the process is 74±5 kcal/ mol. The transition state for the reaction is identified and discussed. A vibrational analysis has been carried out for both HOCl and HClO. The calculated photoelectron spectrum is discussed for each of the molecules.     

Protein-bound water molecules in primate red- and green-sensitive visual pigments

Protein-bound water molecules play crucial roles in the structure and function of proteins. The functional role of water molecules has been discussed for rhodopsin, the light sensor for twilight vision, on the basis of X-ray crystallography, Fourier transform infrared (FTIR) spectroscopy, and a radiolytic labeling method, but nothing is known about the protein-bound waters in our color visual pigments. Here we apply low-temperature FTIR spectroscopy to monkey red (MR)- and green (MG)-sensitive color pigments at 77 K and successfully identify water vibrations using D(2)O and D(2)(18)O in the whole midinfrared region. The observed water vibrations are 6-8 for MR and MG, indicating that several water molecules are present near the retinal chromophore and change their hydrogen bonds upon retinal photoisomerization. In this sense, color visual pigments possess protein-bound water molecules essentially similar to those of rhodopsin. The absence of strongly hydrogen-bonded water molecules (O-D stretch at <2400 cm(-1)) is common between rhodopsin and color pigments, which greatly contrasts with the case of proton-pumping microbial rhodopsins. On the other hand, two important differences are observed in water signal between rhodopsin and color pigments. First, the water vibrations are identical between the 11-cis and 9-cis forms of rhodopsin, but different vibrational bands are observed at >2550 cm(-1) for both MR and MG. Second, strongly hydrogen-bonded water molecules (2303 cm(-1) for MR and 2308 cm(-1) for MG) are observed for the all-trans form after retinal photoisomerization, which is not the case for rhodopsin. These specific features of MR and MG can be explained by the presence of water molecules in the Cl(-)-biding site, which are located near positions C11 and C9 of the retinal chromophore. The averaged frequencies of the observed water O-D stretching vibrations for MR and MG are lower as the λ(max) is red-shifted, suggesting that water molecules are involved in the color tuning of our vision.     

Complete assignment of the hydrogen out-of-plane wagging vibrations of bathorhodopsin: chromophore structure and energy storage in the primary photoproduct of vision

Resonance Raman vibrational spectra of the retinal chromophore in bathorhodopsin have been obtained after regenerating bovine visual pigments with an extensive series of 13C- and deuterium-labeled retinals. A low-temperature spinning cell technique was used to produce high-quality bathorhodopsin spectra exhibiting resolved hydrogen out-of-plane wagging vibrations at 838, 850, 858, 875, and 921 cm-1. The isotopic shifts and a normal coordinate analysis permit the assignment of these lines to the HC7 = C8H Bg, C14H, C12H, C10H, and C11H hydrogen out-of-plane wagging modes, respectively. The coupling constant between the C11H and C12H wags as well as the C12H wag force constant are unusually low compared to those of retinal model compounds. This quantitatively confirms the lack of coupling between the C11H and C12H wags and the low C12H wag vibrational frequency noted earlier by Eyring et al. [(1982) Biochemistry 21, 384]. The force constants for the C10H and C14H wags are also significantly below the values observed in model compounds. We suggest that the perturbed hydrogen out-of-plane wagging and C-C stretching force constants for the C10-C11 = C12-C13 region of the chromophore in bathorhodopsin result from electrostatic interactions with a charged protein residue. This interaction may also contribute to the 33 kcal/mol energy storage in bathorhodopsin.     

Large-amplitude bending motions in phenylalanine transfer RNA

Conformational energy calculations on yeast tRNA^Phe reveal several likely modes of intramolecular bending, including both hingelike motions (rotations about a discrete point) and distributed flexibility (deformations that bend a double-helical segment along a smooth curve). By combining these modes of motion, the molecule can be bent from the L-shaped crystallographic structure to two extremes. It can be straightened into a nearly linear conformation at an energy cost of about 50 kcal/mol, and it can be doubled over to a conformation where the anticodon and the amino acid acceptor terminus are separated by about 40 Å at an energy cost of less than 100 kcal/mol. A bending range of over 100° can be covered for 50 kcal/mol, and we estimate that this value could be cut in half with a minimization algorithm that produced optimum stereochemistry. These energies are comparable to those that would be associated with changes in solvation due to changes in surface area as the molecule bends, indicating that there are no major steric barriers to tRNA flexibility and that variations in solvent conditions and interactions with other molecules may produce large changes in the overall conformation of tRNA.     

Effect of the Vibrational/Rotational Energy Partitioning on the Energy Transfer in Atom–Triatomic Molecule Collisions

The effect of the initial partitioning of the molecular energy between vibrational and rotational modes of a triatomic molecule on the collisional energy transfer is studied for a model atom–triatomic molecule system. We considered the collisions of thermal bath Ar atoms with SO2 molecules, and used the trajectory calculations for determining the energy transfer for three different samplings of initial conditions of the molecule. The first sampling method generated the microcanonical distribution over all states, entering into the vibrational and rotational manifolds, while two others produced distributions with relatively lower values of the rotational energies. It is shown that both the average energy transfer per collision and the mechanism of the energy exchange are significantly affected by the vibrational/rotational energy partitioning before the collisions. Relative decrease in the rotational energy results in the decrease of the averaged energy transfer and progressively emphasizes the role of active rotation as the gateway for translation-vibration energy exchange.     

Wavelength regulation in iodopsin, a cone pigment.

The opsin shift, the difference in wavenumber between the absorption peak of a visual pigment and the protonated Schiff base of the chromophore, represents the influence of the opsin binding site on the chromophore. The opsin shift for the chicken cone pigment iodopsin is much larger than that for rhodopsin. To understand the origin of this opsin shift and the mechanism of wavelength regulation in iodopsin, a series of synthetic 9-cis and 11-cis dehydro- and dihydro-retinals was used to regenerate iodopsin-based pigments. The opsin shifts of these pigments are quite similar to those found in bacteriorhodopsin-based artificial pigments. On the basis of these studies, a tentative model of wavelength regulation in iodopsin is proposed.     

Rapid release of retinal from a cone visual pigment following photoactivation

As part of the visual cycle, the retinal chromophore in both rod and cone visual pigments undergoes reversible Schiff base hydrolysis and dissociation following photobleaching. We characterized light-activated release of retinal from a short-wavelength-sensitive cone pigment (VCOP) in 0.1% dodecyl maltoside using fluorescence spectroscopy. The half-time (t(1/2)) of release of retinal from VCOP was 7.1 s, 250-fold faster than that of rhodopsin. VCOP exhibited pH-dependent release kinetics, with the t(1/2) decreasing from 23 to 4 s with the pH decreasing from 4.1 to 8, respectively. However, the Arrhenius activation energy (E(a)) for VCOP derived from kinetic measurements between 4 and 20 °C was 17.4 kcal/mol, similar to the value of 18.5 kcal/mol for rhodopsin. There was a small kinetic isotope (D(2)O) effect in VCOP, but this effect was smaller than that observed in rhodopsin. Mutation of the primary Schiff base counterion (VCOP(D108A)) produced a pigment with an unprotonated chromophore (λ(max) = 360 nm) and dramatically slowed (t(1/2) ~ 6.8 min) light-dependent retinal release. Using homology modeling, a VCOP mutant with two substitutions (S85D and D108A) was designed to move the counterion one α-helical turn into the transmembrane region from the native position. This double mutant had a UV-visible absorption spectrum consistent with a protonated Schiff base (λ(max) = 420 nm). Moreover, the VCOP(S85D/D108A) mutant had retinal release kinetics (t(1/2) = 7 s) and an E(a) (18 kcal/mol) similar to those of the native pigment exhibiting no pH dependence. By contrast, the single mutant VCOP(S85D) had an ~3-fold decreased retinal release rate compared to that of the native pigment. Photoactivated VCOP(D108A) had kinetics comparable to those of a rhodopsin counterion mutant, Rho(E113Q), both requiring hydroxylamine to fully release retinal. These results demonstrate that the primary counterion of cone visual pigments is necessary for efficient Schiff base hydrolysis. We discuss how the large differences in retinal release rates between rod and cone visual pigments arise, not from inherent differences in the rate of Schiff base hydrolysis but rather from differences in the properties of noncovalent binding of the retinal chromophore to the protein.     

Quantifying the Entropy of Binding for Water Molecules in Protein Cavities by Computing Correlations

Protein structural analysis demonstrates that water molecules are commonly found in the internal cavities of proteins. Analysis of experimental data on the entropies of inorganic crystals suggests that the entropic cost of transferring such a water molecule to a protein cavity will not typically be greater than 7.0 cal/mol/K per water molecule, corresponding to a contribution of approximately +2.0 kcal/mol to the free energy. In this study, we employ the statistical mechanical method of inhomogeneous fluid solvation theory to quantify the enthalpic and entropic contributions of individual water molecules in 19 protein cavities across five different proteins. We utilize information theory to develop a rigorous estimate of the total two-particle entropy, yielding a complete framework to calculate hydration free energies. We show that predictions from inhomogeneous fluid solvation theory are in excellent agreement with predictions from free energy perturbation (FEP) and that these predictions are consistent with experimental estimates. However, the results suggest that water molecules in protein cavities containing charged residues may be subject to entropy changes that contribute more than +2.0 kcal/mol to the free energy. In all cases, these unfavorable entropy changes are predicted to be dominated by highly favorable enthalpy changes. These findings are relevant to the study of bridging water molecules at protein-protein interfaces as well as in complexes with cognate ligands and small-molecule inhibitors.     

Ground state isomerization of a model green fluorescent protein chromophore

The relationship between ground state cis-trans isomerization and protonation state is explored for a model green fluorescent protein chromophore, 4-hydroxybenzylidene-1,2-dimethylimidazolinone (HBDI). We find that the protonation state has only a modest effect on the free energy differences between cis and trans isomers and on the activation energies for isomerization. Specifically, the experimental free energy differences are 3.3, 8.8, and 9.6 kJ/mol for cationic, neutral, and anionic forms of HBDI, respectively, and the activation energies are 48.9, 54.8, and 54.8 kJ/mol for cationic, neutral, and anionic forms, respectively. Furthermore, these activation energies are much smaller than might be expected based on comparison with similar systems. These results suggest that there may be a sub-population of the chromophore, which is nearly equally accessible to all three protonation states, through which thermal isomerization may proceed.     

all-trans-retinoids and dihydroretinoids as probes of the role of chromophore structure in rhodopsin activation

The absorption of a photon of light by rhodopsin results in the cis to trans isomerization of the 11-cis-retinal Schiff base chromophore. In the studies reported here, an attempt is made to determine the mechanism of the energization of rhodopsin as it relates to the chemistry of the isomerization process and the geometrical state of the chromophore. Studies were performed with vitamin A analogues to probe this mechanism. Both 11-cis-7,8-dihydroretinal and 9-cis-7,8-dihydroretinal form bleachable pigments when combined with opsin. Photolysis of these pigments in the presence of G-protein results in the activation of the latter as revealed by its GTPase activity. Phosphodiesterase is also activated when it is included in the incubation. Therefore, the possibility that rhodopsin is energized by mechanisms involving photochemically induced charge transfer from the protonated Schiff base to the beta-ionone ring can be discarded. Further studies were conducted with all-trans-vitamin A derivatives to determine if these compounds can form the GTPase-activating state R*, a situation that is possible, in principle, by microscopic reversibility. Neither all-trans-retinal nor its oxime, when incubated with bovine opsin in the dark, caused activation of the GTPase, requiring at least a 5 kcal/mol energy gap between them. Furthermore, stoichiometric adducts of all-trans-retinoids and opsin were also unable to mediate activation of the GTPase. Since both all-trans-15,16-dihydroretinylopsin and all-trans-retinoylopsin possess an all-trans-retinoid permanently adducted to opsin, it can be concluded that the all-trans-retinoid chromophore-opsin linkage may be necessary but not sufficient to achieve activation of the visual pigment.     

Molecular modelling studies unveil potential binding sites on human serum albumin for selected experimental and in silico COVID-19 drug candidate molecules

Human serum albumin (HSA) is the most prevalent protein in the blood plasma which binds an array of exogenous compounds. Drug binding to HSA is an important consideration when developing new therapeutic molecules, and it also aids in understanding the underlying mechanisms that govern their pharmacological effects. This study aims to investigate the molecular binding of coronavirus disease 2019 (COVID-19) therapeutic candidate molecules to HSA and to identify their putative binding sites. Binding energies and interacting residues were used to evaluate the molecular interaction. Four drug candidate molecules (β-D-N4-hydroxycytidine, Chloroquine, Disulfiram, and Carmofur) demonstrate weak binding to HSA, with binding energies ranging from −5 to −6.7 kcal/mol. Ivermectin, Hydroxychloroquine, Remdesivir, Arbidol, and other twenty drug molecules with binding energies ranging from −6.9 to −9.5 kcal/mol demonstrated moderate binding to HSA. The strong HSA binding drug candidates consist of fourteen molecules (Saquinavir, Ritonavir, Dihydroergotamine, Daclatasvir, Paritaprevir etc.) with binding energies ranging from −9.7 to −12.1 kcal/mol. All these molecules bind to different HSA subdomains (IA, IB, IIA, IIB, IIIA, and IIIB) through molecular forces such as hydrogen bonds and hydrophobic interactions. Various pharmacokinetic properties (gastrointestinal absorption, blood-brain barrier permeation, P-glycoprotein substrate, and cytochrome P450 inhibitor) of each molecule were determined using SwissADME program. Further, the stability of the HSA-ligand complexes was analyzed through 100 ns molecular dynamics simulations considering various geometric properties. The binding free energy between free HSA and compounds were calculated using Molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) and molecular mechanics generalized Born surface area (MM/GBSA) approach. The findings of this study might be useful in understanding the mechanism of COVID-19 drug candidates binding to serum albumin protein, as well as their pharmacodynamics and pharmacokinetics.Keyword: Human serum albumin, HAS, Serum protein, COVID-19, Molecular docking, Molecular dynamics simulation, Pharmacokinetics, Pharmacodynamics     

B3LYP/6-311++G** study of monohydrates of alpha- and beta-D-glucopyranose: hydrogen bonding, stress energies, and effect of hydration on internal coordinates

Twenty-six monohydrates of alpha- and beta-D-glucopyranose were studied using gradient methods at the B3LYP/6-311++G** level of theory. Geometry optimization was carried out with the water molecules at different configurations around the glucose molecule. A new nomenclature for hydrated carbohydrates was developed to describe the water configurations. Zero-point vibrational energy, enthalpy, entropy, and relative free energy were obtained using the harmonic approximation. Hydrogen-bond energies for the monohydrates range from approximately -5 to -12 kcal/mol, and the average relative free energy is approximately 5 kcal/mol. The 1-hydroxy position is the most energetically favored site for hydration, and the region between the two and three positions is the next-most favored site. A water molecule approaching alpha-D-glucose between the 1- and 2-hydroxy positions pulls the 2-hydroxyl hydrogen atom away from the 1-hydroxy oxygen atom, thus increasing the hydrogen-bond length and also increasing the alpha-D-glucose energy. The increase in energy that occurs with a similar interaction on the beta-anomer is much less effective since the hydrogen bond is much longer. Using the calculated free energies of all 26 configurations, the anomer population (alpha/beta) increases in the beta-anomer population relative to the in vacuo case by approximately 10% at the expense of the alpha-anomer, giving an (alpha/beta) ratio of approximately 50/50. This result arises from entropy contributions favoring the beta-anomer more than the alpha-anomer. From analysis of donor and acceptor hydrogen-bond lengths, excellent correlation is found between the DFT calculated distances and those taken from carbohydrate structures in the Cambridge Crystallographic Data Bank.     

Comparative study on the chromophore binding sites of rod and red-sensitive cone visual pigments by use of synthetic retinal isomers and analogues

A comparative study on the chromophore (retinal) binding sites of the opsin (R-photopsin) from chicken red-sensitive cone visual pigment (iodopsin) and that scotopsin) from bovine rod pigment (rhodopsin) was made by the aid of geometric isomers of retinal (all-trans, 13-cis, 11-cis, 9-cis, and 7-cis) and retinal analogues including fluorinated (14-F, 12-F, 10-F, and 8-F) and methylated (12-methyl) 11-cis-retinals. The stereoselectivity of R-photopsin for the retinal isomers and analogues was almost identical with that of scotopsin, indicating that the shapes of the chromophore binding sites of both opsins are similar, although the former appears to be somewhat more restricted than the latter. The rates of pigment formation from R-photopsin were considerably greater than those from scotopsin. In addition, all the iodopsin isomers and analogues were more susceptible to hydroxylamine than were the rhodopsin ones. These observations suggest that the retinal binding site of iodopsin is located near the protein surface. On the basis of the spectral properties of fluorinated analogues, a polar group in the chromophore binding site of iodopsin as well as rhodopsin was estimated to be located near the hydrogen atom at the C10 position of the retinylidene chromophore. A large difference in wavelength between the absorption maxima of iodopsin and rhodopsin was significantly reduced in the 9-cis and 7-cis pigments. On the assumption that the retinylidene chromophore is anchored rigidly at the alpha-carbon of the lysine residue and loosely at the cyclohexenyl ring, each of the two isomers would have the Schiff-base nitrogen at a position altered from that of the 11-cis pigments.(ABSTRACT TRUNCATED AT 250 WORDS)     

Raman microscope and quantum yield studies on the primary photochemistry of A2-visual pigments.

The 77-K resonance Raman vibrational spectrum of intact goldfish rod photoreceptors containing 3,4-dehydro (A2) retinal is dominated by scattering from the 9-cis component of the steady state at all excitation wavelengths. Intact goldfish photoreceptors were regenerated with an A1-retinal chromophore to determine whether this behavior is caused by the protein or the chromophore. The resulting Raman spectrum was typical of an A1-pigment exhibiting significant scattering from all three components of the steady state: rhodopsin, bathorhodopsin, and isorhodopsin. Furthermore, regeneration of bovine opsin with A2-retinal produces a characteristic "A2-Raman spectrum" that is dominated by scattering from the 9-cis pigment. We conclude that the differences between the Raman spectra of the A1-and A2-pigments are caused by some intrinsic difference in the photochemical properties of the retinal chromophores. To quantitate these observations, the 77-K adsorption spectra and the photochemical quantum yields (phi) of the native A2-goldfish and the regenerated A2-bovine pigments were measured. In the goldfish A2-pigment, the value of phi 4 (9-cis----trans) is 0.05; phi 3 (trans----9-cis) is 0.10; and phi 2 (trans----11-cis) is 0.35. By contrast, in the bovine A1-pigment, these quantum yields are 0.10, 0.053, and 0.50, respectively. The reduced value of phi 4 and the increased value of phi 3 in the goldfish pigment confirms that the 9-cis isomer is photochemically more stable in A2-pigments.     

Deoxymyoglobin studied by the conformational normal mode analysis. II. The conformational change upon oxygenation

The conformational change taking place in myoglobin concomitantly with the observed geometrical change at the heme-His(F8) linkage upon oxygenation is studied by normal mode analysis, which is based on the quadratic approximation of the conformational energy function. The heme-globin interaction energy increases for this change by 8.114 kcal/mol when both the heme group and the globin molecule are held rigid. When they are permitted flexibility, the interaction energy relaxes by 7.038 kcal/mol, and the difference (1.076 kcal/mol) is distributed as strain energy within the molecule. This increase is the work necessary for the heme group to move against the force exerted by the globin. If the force is assumed to be invariable during this move, the work is small, 0.276 kcal/mol, meaning that the force is strongly variable. Furthermore, this means that the heme group is located near the equilibrium point of the potential energy of the heme-globin interaction. The change in the localized strain energy stored in the force field at the linkage between the heme and the imidazole of HisF8 is estimated to be of the same order of magnitude as the distributed energy. The largest atomic displacements are observed at the region from the F helix to the beginning of the G helix, and secondary large displacements occur at several regions, i.e, the A helix, from the C helix to the CD corner, the E helix, and the C-terminal side of the H helix. All of these regions have strong dynamic interactions with the heme group, either directly or indirectly. Their secondary structures show complex deformations. In other parts, relatively rigid segments undergo rotational and/or bending changes in a way consistent with the large changes described above and close atomic packing within the molecule. The calculated conformational change is decomposed to vibrational normal modes of deoxymyoglobin. The magnitude of the conformational change, measured by the mass-weighted mean-square atomic displacement, is accounted for up to 92.0% by the 151 normal modes with frequencies lower than 40 cm-1. In descending order of contribution, the first six modes, each of which has a frequency lower than 12 cm-1, account for up to 57.4%. This means that the functionally important conformational change can well be expressed in terms of a relatively small number of collective low frequency normal modes.     

Mapping the energetics of water-protein and water-ligand interactions with the "natural" HINT forcefield: predictive tools for characterizing the roles of water in biomolecules

The energetics and hydrogen bonding pattern of water molecules bound to proteins were mapped by analyzing structural data (resolution better than 2.3A) for sets of uncomplexed and ligand-complexed proteins. Water-protein and water-ligand interactions were evaluated using hydropatic interactions (HINT), a non-Newtonian forcefield based on experimentally determined logP(octanol/water) values. Potential water hydrogen bonding ability was assessed by a new Rank algorithm. The HINT-derived binding energies and Ranks for second shell water molecules were -0.04 kcal mol(-1) and 0.0, respectively, for first shell water molecules -0.38 kcal mol(-1) and 1.6, for active site water molecules -0.45 kcal mol(-1) and 2.3, for cavity water molecules -0.55 kcal mol(-1) and 3.3, and for buried water molecules -0.56 kcal mol(-1) and 4.4. For the last four classes, similar energies indicate that internal and external water molecules interact with protein almost equally, despite different degrees of hydrogen bonding. The binding energies and Ranks for water molecules bridging ligand-protein were -1.13 kcal mol(-1) and 4.5, respectively. This energetic contribution is shared equally between protein and ligand, whereas Rank favors the protein. Lastly, by comparing the uncomplexed and complexed forms of proteins, guidelines were developed for prediction of the roles played by active site water molecules in ligand binding. A water molecule with high Rank and HINT score is unlikely to make further interactions with the ligand and is largely irrelevant to the binding process, while a water molecule with moderate Rank and high HINT score is available for ligand interaction. Water molecule displaced for steric reasons were characterized by lower Rank and HINT score. These guidelines, tested by calculating HINT score and Rank for 50 water molecules bound in the active site of four uncomplexed proteins (for which the structures of the liganded forms were also available), correctly predicted the ultimate roles (in the complex) for 76% of water molecules. Some failures were likely due to ambiguities in the structural data.     

Membrane phospholipids as an energy source in the operation of the visual cycle

Biology depends on the coupling of the free energy of hydrolysis of phosphate esters, such as ATP, to drive processes which would otherwise be thermodynamically unfavorable. Carboxyl esters are like phosphate esters in their ability to hydrolyze with substantial negative free energies, enabling them to participate in group transfer processes as well. In particular, membrane phospholipids constitute an enormous store of potential energy that could be used to fuel energetically unfavorable processes. One such process involves the biosynthesis of 11-cis-retinal, the chromophore of rhodopsin, from all-trans-retinol (vitamin A). The difference in free energy between an all-trans retinoid and its corresponding 11-cis retinoid is approximately 4 kcal/mol. This energy is provided for in a minimally two-step process involving membrane phospholipids as the energy source. First, all-trans-retinol is esterified in the retinal pigment epithelium by lecithin retinol acyl transferase (LRAT) to produce an all-trans-retinyl ester. Second, this ester is transformed into 11-cis-retinol by an isomerohydrolase in a process that couples the negative free energy of hydrolysis of the acyl ester to the formation of the strained 11-cis-retinol.