Understanding ligand binding, selectivity and functions on the G protein-coupled receptors: A molecular modeling approach

The assessment of target protein molecular structure provides a distinct advantage in the rational drug design process. The increasing number of available G protein-coupled receptor crystal structures has enabled utilization of a varied number of computational approaches for understanding the ligand-receptor interactions, ligand selectivity and even receptor response upon ligand binding. The following dissertation examines the results from three different projects with varied objectives – i) structural modeling of human C-C chemokine receptor type 5 (CCR5) and assessment of the ligand binding pocket of the receptor, ii) assessment of the selectivity profile of naltrexone derivatives on the three opioid receptors (μ-opioid, κ-opioid, δ-opioid) with an aim towards designing selective μ-opioid receptor antagonists, and iii) structural modeling of the ‘active’ state conformation of the κ-opioid receptor in response to agonist binding and determination of a plausible molecular mechanism involved in activation ‘switch’ of the κ-opioid receptor. In absence of a crystal-based molecular structure of CCR5, a homology model of the receptor was built and the ligand binding pocket was validated. On the basis of evaluation of the ligand-receptor interactions on the validated binding pocket, structural and chemical modifications to anibamine, a natural plant product, were proposed to enhance its receptor binding. The selectivity of naltrexone (a universal antagonist) was assessed with respect to the three opioid receptors by employing ligand docking studies and the ‘message-address’ concept. Multiple address sites were identified on the opioid receptors and structural modifications were proposed for the naltrexone derivatives for their enhanced selectivity. In the third project, structural modeling of the active state conformation of the κ-opioid receptor covalently bound to a salvinorin A derivative (agonist) was attempted via molecular dynamics simulations. Although the obtained molecular model lacked the signature ‘agonist-like’ conformations, the result provides a template for such studies in the future

A combined experimental and computational approach to investigate emergent network dynamics based on large-scale neuronal recordings

Sviluppo di un approccio integrato computazionale-sperimentale per lo studio di reti neuronali mediante registrazioni elettrofisiologich

The role of dynamic hydrogen bond networks in protonation coupled dynamics of retinal proteins

Hydrogen bonds (H-bonds) are an essential interaction in membrane proteins. Embedded in complex hydrated lipid bilayers, intramolecular interactions through the means of hydrogen bonding networks are often crucial for the function of the protein. Internal water molecules that occupy stable sites inside the protein, or water molecules that visit transiently from the bulk, can play an important role in shaping local conformational dynamics forming complex networks that bridge regions of the protein via water-mediated hydrogen bonds that can function as wires for the transferring of protons as a part of the protein’s function. For example, the membrane-embedded channelrhodopsins which are found in archaea are proteins that couple light induced isomerization of a retinal chromophore with proton transfer reactions and passive flow of cations through their pore. I contributed to the development of a new algorithm package that features a unique approach to H-bond analyses. I performed analyses of long Molecular Dynamics (MD) trajectories of channelrhodopsin variants embedded in hydrated lipid membranes and large data sets of static structures, to detect and dissect dynamic hydrogen-bond networks. The photocycle of channelrhodopsins begins with absorption and isomerization of the retinal from an all-trans state to a 13-cis state and followed by the deprotonation of the Schiff base. Thus, the retinal is found in the epicenter of the analyses. Through the use of 2-dimensional graphs of the protein H-bond networks I identified protein groups potentially important for the proton transfer activity. Local dynamics are highly affected by point mutations of amino acids important for function. The interior of channelrhodopsin C1C2 hosts extensive networks of protein and H-bonded-water molecules, and a never reported before, network that can bridge transiently the two retinal chromophores in channelrhodopsin dimers. In a recently identified inward proton pump, AntR, I applied centrality measures on MD trajectories of the homology model I generated, to assess the communication of the amino acid residues within the networks. I detected a frequently sampled long water chain that connects the retinal with a candidate proton acceptor, as well as a conserved serine in the vicinity of the retinal chromophore plays a significant role in the connectivity and communication of the H-bond networks upon isomerization. A similar water bridge is sampled in independent simulations of ChR2, where a participant for the proton donor group connects to the 13-cis,15-anti retinal. Proton transfer reactions often take place through certain amino acids, forming patterns. I analyzed H-bond patterns or motifs in large hand-curated datasets of static structures of α-transmembrane helix proteins, organized according to the superfamilies they belong, their function and an alternative classification method. The presence of motifs in TM proteins is tightly related to their families/superfamilies of the host protein and their position along the membrane normal.Wasserstoffbrücken (H-Brücke) sind eine wesentliche Wechselwirkung in Membranproteinen. Eingebettet in komplexe hydratisierte Lipiddoppelschichten sind intramolekulare Wechselwirkungen über Wasserstoffbrückenbindungsnetzwerke oft entscheidend für die Funktion des Proteins. Interne Wassermoleküle, die stabile Stellen im Inneren des Proteins besetzen, oder Wassermoleküle, die vorübergehend aus der Masse zu Besuch kommen, können eine wichtige Rolle bei der Gestaltung der lokalen Konformationsdynamik spielen, indem sie komplexe Netzwerke bilden, die Regionen des Proteins über wasservermittelte Wasserstoffbrückenbindungen überbrücken, die als Drähte für den Transfer von Protonen als Teil der Proteinfunktion funktionieren können. Die in Archaeen vorkommenden, in die Membran eingebetteten Kanalrhodopsine sind beispielsweise Proteine, die die lichtinduzierte Isomerisierung eines Retinachromophors mit Protonentransferreaktionen und dem passiven Fluss von Kationen durch ihre Pore verbinden. Ich habe an der Entwicklung eines neuen Algorithmuspakets mitgewirkt, das einen einzigartigen Ansatz für H-Bindungsanalysen bietet. Ich habe lange Molekulardynamik-Trajektorien von Kanalrhodopsine-Varianten, die in hydratisierte Lipidmembranen eingebettet sind, sowie große Datensätze statischer Strukturen analysiert, um dynamische Wasserstoffbrücken-bindungsnetzwerke zu erkennen und zu zerlegen. Der Photozyklus der Kanalrhodopsine beginnt mit der Absorption und Isomerisierung des Retinals von einem all-trans-Zustand zu einem 13-cis-Zustand, gefolgt von der Deprotonierung der Schiff-Base. Somit steht das Retinal im Mittelpunkt der Analysen. Durch die Verwendung von 2-dimensionalen Graphen der Protein- H-Brückenetzwerke identifizierte ich Proteingruppen, die für die Protonentransferaktivität wichtig sein könnten. Die lokale Dynamik wird durch Punktmutationen der für die Funktion wichtigen Aminosäuren stark beeinflusst. Das Innere von Kanalrhodopsine C1C2 beherbergt ausgedehnte Netzwerke von Protein- und H-Brücke-Wassermolekülen und ein bisher unbekanntes Netzwerk, das die beiden retinalen Chromophore in Kanalrhodopsine-Dimeren vorübergehend überbrücken kann. In einer kürzlich identifizierten Protonenpumpe, AntR, wendete ich Zentralitätsmaße auf MD-Trajektorien des von mir erstellten Homologiemodells an, um die Kommunikation der Aminosäurereste innerhalb der Netzwerke zu bewerten. Ich fand, dass eine häufig gesampelte lange Wasserkette, die das Retinal mit einem Protonenakzeptor verbindet, sowie ein konserviertes Serin in der Nähe des Retinal-Chromophors eine wichtige Rolle bei der Konnektivität und Kommunikation der H-Brückesnetzwerke bei der Isomerisierung spielt. Eine ähnliche Wasserbrücke ist in unabhängigen Simulationen von Kanalrhodopsine-2 zu finden, wo ein Teilnehmer für die Protonendonorgruppe mit dem 13-cis,15-anti-Retinal verbunden ist. Protonenübertragungsreaktionen finden oft über bestimmte Aminosäuren statt und bilden Muster. Ich analysierte H-Brückemuster oder -motive in großen, von Hand kuratierten Datensätzen statischer Strukturen von α-Transmembranhelix-Proteinen, geordnet nach den Superfamilien, zu denen sie gehören, ihrer Funktion und einer alternativen Klassifizierungsmethode. Das Vorhandensein von Motiven in TM-Proteinen steht in engem Zusammenhang mit ihren Familien/Superfamilien des Wirtsproteins und ihrer Position entlang der Membrannormale

The Structure and Function of Photosystem I and Photosystem I – Hydrogenase Protein Fusions: An Experimental and Computational Study

Photosystem I (PSI) is a membrane protein involved in the photosynthetic cycle of plants, algae, and cyanobacteria that is of specific interest due to its ability to harness solar energy to generate reducing power. This work seeks to form an in vitro hybrid protein fusion between the membrane integral PSI protein and the membrane-bound hydrogenase (MBH) enzyme, in an effort to improve electron transport between these two proteins. Small-angle neutron scattering (SANS) was used to characterize the detergent-solubilized solution structure of trimeric PSI from the cyanobacterium Thermosynechococcus elongatus, which showed that the detergent interacts primarily with the hydrophobic periphery of PSI. The SANS results were used as a guide to constructing a model of trimeric PSI embedded in a detergent belt. Subsequent all-atom molecular dynamics (MD) simulations of the PSI-detergent complex suggested that the detergent environment could negatively impact the long-term stability of PSI, but is not likely to affect PSI activity or hinder its ligation to the MBH. Having verified that the solution structure of the PSI-detergent complex will not affect formation of PSI-MBH fusions, the membrane-bound [NiFe]-hydrogenase of Ralstonia eutropha was genetically engineered to express a Gly3 [Gly-Gly-Gly] tag on the N-terminus of the small subunit to allow for site-specific ligation to the psaE subunit of PSI. H2 [hydrogen] uptake activity results show a complete loss of activity in the mutant R. eutropha strain, possibly due to mutations introduced during previous genetic engineering work. In parallel, MD simulations of the PSI-MBH fusion protein indicate this ligation strategy is not optimal for electron transport between these proteins. This MD approach can be used to evaluate other PSI-MBH fusion strategies, possibly targeting other stromal subunits of PSI. Finally, MD simulations of previously studied PSI-[FeFe]-hydrogenase fusions were conducted, revealing significant distortion of the protein structure that could limit their long-term stability

From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions

Over the last few decades, conical intersections (CoIns) have grown from theoretical curiosities into common mechanistic features of photochemical reactions, whose function is to funnel electronically excited molecules back to their ground state in regions where the potential energy surfaces (PESs) of two electronic states become degenerate. Analogous to transition states in thermal chemistry, CoIns appear as transient structures providing a kinetic bottleneck along a reaction coordinate. However, such a bottleneck is not associated with the probability of crossing an energy barrier but rather with an excited state decay probability along a full "line" of transient structures connected by non-reactive modes, the intersection space (IS). This article will review our understanding of the factors controlling CoIn mediated ultrafast photochemical reactions, taking a physical organic chemist approach by discussing a number of case studies for small organic molecules and photoactive proteins. Such discussion will be carried out by first introducing the "standard" one-mode model based on Landau-Zener (LZ) theory to describe a reactive excited state decay event intercepting, locally, a single CoIn along a single direction, and then by providing a modern perspective based on the effects of the phase matching of multiple modes on the same local event, thus redefining and expanding the description of the excited state reaction coordinate. The direct proportionality between the slope (or velocity) along one mode and decay probability at a single CoIn is a widely applied fundamental principle that follows from the LZ model, yet it fails to provide a complete understanding of photochemical reactions whose local reaction coordinate changes along the IS. We show that in these situations, in particular by focussing on rhodopsin double bond photoisomerization, it is mandatory to consider additional molecular modes and their phase relationship approaching the IS, hence providing a key mechanistic principle of ultrafast photochemistry based on the phase matching of those modes. We anticipate that this qualitative mechanistic principle should be considered in the rational design of any ultrafast excited state process, impacting various fields of research ranging from photobiology to light-driven molecular devices

On the potential role of lateral connectivity in retinal anticipation

International audienceWe analyse the potential effects of lateral connectivity (amacrine cells and gap junctions) on motion anticipation in the retina. Our main result is that lateral connectivity can-under conditions analysed in the paper-trigger a wave of activity enhancing the anticipation mechanism provided by local gain control [8, 17]. We illustrate these predictions by two examples studied in the experimental literature: differential motion sensitive cells [1] and direction sensitive cells where direction sensitivity is inherited from asymmetry in gap junctions connectivity [73]. We finally present reconstructions of retinal responses to 2D visual inputs to assess the ability of our model to anticipate motion in the case of three different 2D stimuli

BIOPHYSICAL CHARACTERIZATION OF CHEMICALLY UNFOLDED STATES OF THE MEMBRANE PROTEIN RHODOPSIN

Membrane proteins function as important communication channels of the cell and its environment that aid in regulating the overall homeostasis of organisms. Understanding the pathways by which these proteins adopt their three-dimensional structures can provide us with key insights into their functions. Failure of a membrane protein to fold into its native structure can lead to disruption of their functions and cause diseases. Through an understanding of the folding mechanisms of membrane proteins it may be possible to identify avenues for the treatment of such diseases. Towards these goals, this thesis describes the biophysical characterization of denatured states of rhodopsin, a model system selected to study helical membrane protein folding. The first contribution of this thesis was to establish approaches that can be used to identify suitable conditions for studying membrane protein folding in vitro. This required screening different denaturing conditions to obtain maximum unfolding without causing aggregation of rhodopsin. 30% SDS and 3% SDS + 8 M urea were found to be the most suitable denaturing conditions. Next, structural features of largely unfolded states of rhodopsin under optimized denaturing conditions were systematically characterized focussing on three levels of structural resolution: global, local and site-specific. Global tertiary structure changes upon SDS denaturation were observed to correlate with SDS micellar structure changes and also hinted at formation of compact intermediate states. Local structural dynamics, probed by NMR spectroscopy, showed that the cytoplasmic domain is more flexible than extracellular and transmembrane domains taken together in spite of an overall increase in flexibility with denaturation. Mobility studies probing site-specific changes by EPR spectroscopy, showed that specific extracellular residues retain more rigidity than cytoplasmic residues in denatured states. These results indicate that the former domain is involved in more stable interactions forming a possible folding core like structure, the location of which correlates with that described by the long-range interaction model of folding. Finally, the importance of dynamics in understanding folding mechanisms of rhodopsin led us to contribute to the development of two novel methodologies: terahertz spectroscopy to detect global motions and 19F NMR using new monofluoro labels to quantify residue specific motions