Assessment of the bound conformation of bombesin to the BB1 and BB2 receptors

Bombesin is an endogenous peptide involved in a wide spectrum of physiological activities ranging from satiety, control of circadian rhythm and thermoregulation in the central nervous system, to stimulation of gastrointestinal hormone release, activation of macrophages and effects on development in peripheral tissues. Actions of the peptide are mediated through the two high affinity G-protein coupled receptors BB1R and BB2R. Under pathophysiological conditions, these receptors are overexpressed in many different types of tumors, such as prostate cancer, breast cancer, small and non-small cell lung cancer and pancreatic cancer. This observation has been used for designing cell markers, but it has not been yet exploited for therapeutical purposes. Despite the enormous biological interest of the peptide, little is known about the stereochemical features that contribute to their activity. On the one hand, mutagenesis studies identified a few receptor residues important for high bombesin affinity and on the other, a few studies focused on the relevance of diverse residues of the peptide for receptor activation. Models of the peptide bound to BB1R and BB2R can be helpful to improve our understanding of the stereochemical features granting bombesin activity. Accordingly, the present study describes the computational process followed to construct such models by means of Steered Molecular Dynamics, using models of the peptide and its receptors. Present results provide new insights into the structure-activity relationships of bombesin and its receptors, as well as render an explanation for the differential binding affinity observed towards BB1R and BB2R. Finally, these models can be further exploited to help for designing novel small molecule peptidomimetics with improved pharmacokinetics profile.

A Step Forward toward Selective Activation/Inhibition of Bak, a Pro-Apoptotic Member of the Bcl-2 Protein Family: Discovery of New Prospective Allosteric Sites Using Molecular Dynamics

Bak is a pro-apoptotic protein and a member of the Bcl-2 family that plays a key role in apoptosis, a programmed cell death mechanism of multicellular organisms. Its activation under death stimuli triggers the permeabilization of the mitochondrial outer membrane that represents a point of no return in the apoptotic pathway. This process is deregulated in many tumors where Bak is inactivated, whereas in other cases like in neurodegeneration, it exhibits an excessive response leading to disorders such as the Alzheimer disease. Members of the Bcl-2 family share a common 3D structure, exhibiting an extremely similar orthosteric binding site, a place where both pro and antiapoptotic proteins bind. This similarity raises a selectivity issue that hampers the identification of new drugs, capable of altering Bak activation in a selective manner. An alternative activation site triggered by antibodies has been recently identified, opening the opportunity to undertake new drug discovery studies. Despite this recent identification, an exhaustive study to identify cryptic pockets as prospective allosteric sites has not been yet performed. Thus, the present study aims to characterize novel hotspots in the Bak structure. For this purpose, we have carried out extensive molecular dynamics simulations using three different Bak systems including Bak in its apo form, Bak in complex with its endogen activator Bim and an intermediate form, set up by removing Bim from the previous complex. The results reported in the present work shed some light on future docking studies on Bak through the identification of new prospective allosteric sites, not previously described in this protein.

Effect of Sodium Valproate on the Conformational Stability of the Visual G Protein-Coupled Receptor Rhodopsin

Rhodopsin is the G protein-coupled receptor of rod photoreceptor cells that mediates vertebrate vision at low light intensities. Mutations in rhodopsin cause inherited retinal degenerative diseases such as retinitis pigmentosa. Several therapeutic strategies have attempted to address and counteract the deleterious effect of rhodopsin mutations on the conformation and function of this photoreceptor protein, but none has been successful in efficiently preventing retinal degeneration in humans. These approaches include, among others, the use of small molecules, known as pharmacological chaperones, that bind to the receptor stabilizing its proper folded conformation. Valproic acid, in its sodium valproate form, has been used as an anticonvulsant in epileptic patients and in the treatment of several psychiatric disorders. More recently, this compound has been tested as a potential therapeutic agent for the treatment of retinal degeneration associated with retinitis pigmentosa caused by rhodopsin mutations. We now report on the effect of sodium valproate on the conformational stability of heterologously expressed wild-type rhodopsin and a rhodopsin mutant, I307N, which has been shown to be an appropriate model for studying retinal degeneration in mice. We found no sign of enhanced stability for the dark inactive conformation of the I307N mutant. Furthermore, the photoactivated conformation of the mutant appears to be destabilized by sodium valproate as indicated by a faster decay of its active conformation. Therefore, our results support a destabilizing effect of sodium valproate on rhodopsin I307N mutant associated with retinal degeneration. These findings, at the molecular level, agree with recent clinical studies reporting negative effects of sodium valproate on the visual function of retinitis pigmentosa patients.

Molecular Determinants for the Activation/Inhibition of Bak Protein by BH3 Peptides

Apoptosis is a key cell death pathway in mammalian cells. Understanding this process and its regulation has been a subject of study in the last three decades. Members of the Bcl-2 family of proteins are involved in the regulation of apoptosis through mitochondrial poration with the subsequent initiation of apoptosis. Deregulation of proapoptotic proteins contributes to the progression of many tumor processes. Understanding how these pore-forming Bcl-2 proteins Bak and Bax are activated is key to find new anticancer treatments. As no drug capable of activating Bak has been disclosed yet, the study of the structural features of BH3 peptides-known as Bak activators-relevant for binding along with its binding energy decomposition analysis, becomes essential for designing novel small-molecule mimics of BH3. Interestingly, a BH3 Bim analogue-inactivating Bak has recently been discovered, opening a question on the molecular features that determine the functions of BH3 peptides. Therefore, the present work is aimed at understanding the way BH3 peptides activate or inactivate Bak in order to identify differential structural features that can be used in drug design. For this purpose, complexes of Bak with an activator and an inhibitor have been subjected to a molecular dynamics study. Structural differences were assessed by means of the fluctuations of the corresponding principal components. Moreover, the MMPB/GBSA approach was used to compute the binding free energy of the diverse complexes to identify those residues of the BH3 peptide that exhibit the larger contributions to complex formation. The results obtained in this work show differences between activators and inhibitors, both in structural and energetic terms, which can be used in the design of new molecules that can activate or inactivate proapoptotic Bak.