PocketGraph : graph representation of binding site volumes

The representation of small molecules as molecular graphs is a common technique in various fields of cheminformatics. This approach employs abstract descriptions of topology and properties for rapid analyses and comparison. Receptor-based methods in contrast mostly depend on more complex representations impeding simplified analysis and limiting the possibilities of property assignment. In this study we demonstrate that ligand-based methods can be applied to receptor-derived binding site analysis. We introduce the new method PocketGraph that translates representations of binding site volumes into linear graphs and enables the application of graph-based methods to the world of protein pockets. The method uses the PocketPicker algorithm for characterization of binding site volumes and employs a Growing Neural Gas procedure to derive graph representations of pocket topologies. Self-organizing map (SOM) projections revealed a limited number of pocket topologies. We argue that there is only a small set of pocket shapes realized in the known ligand-receptor complexes

What can we learn from molecular dynamics simulations for GPCR drug design?

Recent years have seen a tremendous progress in the elucidation of experimental structural information for G-protein coupled receptors (GPCRs). Although for the vast majority of pharmaceutically relevant GPCRs structural information is still accessible only by homology models the steadily increasing amount of structural information fosters the application of structure-based drug design tools for this important class of drug targets. In this article we focus on the application of molecular dynamics (MD) simulations in GPCR drug discovery programs. Typical application scenarios of MD simulations and their scope and limitations will be described on the basis of two selected case studies, namely the binding of small molecule antagonists to the human CC chemokine receptor 3 (CCR3) and a detailed investigation of the interplay between receptor dynamics and solvation for the binding of small molecules to the human muscarinic acetylcholine receptor 3 (hM3R)

Charge Recombination and Protein Dynamics in Bacterial Photosynthetic Reaction Centers Entrapped in a Sol-Gel Matrix

Many proteins can be immobilized in silica hydrogel matrices without compromising their function, making this a suitable technique for biosensor applications. Immobilization will in general affect protein structure and dynamics. To study these effects, we have measured the P(+)Q(A)(−) charge recombination kinetics after laser excitation of Q(B)-depleted wild-type photosynthetic reaction centers from Rhodobacter sphaeroides in a tetramethoxysilane (TMOS) sol-gel matrix and, for comparison, also in cryosolvent. The nonexponential electron transfer kinetics observed between 10 and 300 K were analyzed quantitatively using the spin boson model for the intrinsic temperature dependence of the electron transfer and an adiabatic change of the energy gap and electronic coupling caused by protein motions in response to the altered charge distributions. The analysis reveals similarities and differences in the TMOS-matrix and bulk-solvent samples. In both preparations, electron transfer is coupled to the same spectrum of low frequency phonons. As in bulk solvent, charge-solvating protein motions are present in the TMOS matrix. Large-scale conformational changes are arrested in the hydrogel, as evident from the nonexponential kinetics even at room temperature. The altered dynamics is likely responsible for the observed changes in the electronic coupling matrix element

Ligand binding and protein dynamics in neuroglobin

Neuroglobin (Ngb) is a recently discovered protein in vertebrate brain tissue that belongs to the globin family of proteins. It has been implicated in the neuronal response to hypoxia or ischemia, although its physiological role has been hitherto unknown. Ngb is hexacoordinate in the ferrous deoxy form under physiological conditions. To bind exogenous ligands like O(2) and CO, the His E7 endogenous ligand is displaced from the sixth coordination. By using infrared spectroscopy and nanosecond time-resolved visible spectroscopy, we have investigated the ligand-binding reaction over a wide temperature range (3–353 K). Multiple, intrinsically heterogeneous distal heme pocket conformations exist in NgbCO. Photolysis at cryogenic temperatures creates a five-coordinate deoxy species with very low geminate-rebinding barriers. The photodissociated CO is observed to migrate within the distal heme pocket even at 20 K. Flash photolysis near physiological temperature (275–353 K) exhibits four sequential kinetic features: (i) geminate rebinding (t < 1 μs); (ii) extremely fast bimolecular exogenous ligand binding (10 μs < t < 1 ms) with a nontrivial temperature dependence; (iii) endogenous ligand binding (100 μs < t < 10 ms), which can be studied by using flash photolysis on deoxy Ngb; and (iv) displacement of the endogenous by the exogenous ligand (10 ms < t < 10 ks). All four processes are markedly nonexponential, suggesting that Ngb fluctuates among different conformations on surprisingly long time scales

Probing Electric Fields in Protein Cavities by Using the Vibrational Stark Effect of Carbon Monoxide

To determine the magnitude and direction of the internal electric field in the Xe4 cavity of myoglobin mutant L29W-S108L, we have studied the vibrational Stark effect of carbon monoxide (CO) using infrared spectroscopy at cryogenic temperatures. CO was photodissociated from the heme iron and deposited selectively in Xe4. Its infrared spectrum exhibits Stark splitting into two bands associated with CO in opposite orientations. Two different photoproduct states can be distinguished, C′ and C″, with markedly different properties. For C′, characteristic temperature-dependent changes of the area, shift, and width were analyzed, based on a dynamic model in which the CO performs fast librations within a double-well model potential. For the barrier between the wells, a height of ∼1.8 kJ/mol was obtained, in which the CO performs oscillations at an angular frequency of ∼25 cm(−1). The magnitude of the electric field in the C′ conformation was determined as 11.1 MV/cm; it is tilted by an angle of 29° to the symmetry axis of the potential. Above 140 K, a protein relaxation leads to a significantly altered photoproduct, C″, with a smaller Stark splitting and a more confining potential (barrier >4 kJ/mol) governing the CO librations

Temperature dependence of the heat diffusivity of proteins

In a combined experimental–theoretical study, we investigated the transport of vibrational energy from the surrounding solvent into the interior of a heme protein, the sperm whale myoglobin double mutant L29W-S108L, and its dependence on temperature from 20 to 70 K. The hindered libration of a CO molecule that is not covalently bound to any part of the protein but is trapped in one of its binding pockets (the Xe4 pocket) was used as the local thermometer. Energy was deposited into the solvent by IR excitation. Experimentally, the energy transfer rate increased from (30 ps)−1 at 20 K to (8 ps)−1 at 70 K. This temperature trend is opposite to what is expected, assuming that the mechanism of heat transport is similar to that in glasses. In order to elucidate the mechanism and its temperature dependence, nonequilibrium molecular dynamics (MD) simulations were performed, which, however, predicted an essentially temperature-independent rate of vibrational energy flow. We tentatively conclude that the MD potentials overestimate the coupling between the protein and the CO molecule, which appears to be the rate-limiting step in the real system at low temperatures. Assuming that this coupling is anharmonic in nature, the observed temperature trend can readily be explained