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.

Mutations in the photoreceptor protein rhodopsin are known as one of the leading causes of retinal degeneration in humans. Two rhodopsin mutations, Y102H and I307N, obtained in chemically mutagenized mice, are currently the subject of increased interest as relevant models for studying the process of retinal degeneration in humans. Here, we report on the biochemical and functional characterization of the structural and functional alterations of these two rhodopsin mutants and we compare them with the G90V mutant previously analyzed, as a basis for a better understanding of in vivo studies. This mechanistic knowledge is fundamental to use it for developing novel therapeutic approaches for the treatment of inherited retinal degeneration in retinitis pigmentosa. We find that Y102H and I307N mutations affect the inactive–active equilibrium of the receptor. In this regard, the mutations reduce the stability of the inactive conformation but increase the stability of the active conformation. Furthermore, the initial rate of the functional activation of transducin, by the I307N mutant is reduced, but its kinetic profile shows an unusual increase with time suggesting a profound effect on the signal transduction process. This latter effect can be associated with a change in the flexibility of helix 7 and an indirect effect of the mutation on helix 8 and the C-terminal tail of rhodopsin, whose potential role in the functional activation of the receptor has been usually underestimated. In the case of the Y102H mutant, the observed changes can be associated with conformational alterations affecting the folding of the rhodopsin intradiscal domain, and its presumed involvement in the retinal binding process by the receptor.

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.

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.