Prediction of cross-recognition of peptide-HLA A2 by Melan-A-specific cytotoxic T lymphocytes using three-dimensional quantitative structure-activity relationships.

The cross-recognition of peptides by cytotoxic T lymphocytes is a key element in immunology and in particular in peptide based immunotherapy. Here we develop three-dimensional (3D) quantitative structure-activity relationships (QSARs) to predict cross-recognition by Melan-A-specific cytotoxic T lymphocytes of peptides bound to HLA A*0201 (hereafter referred to as HLA A2). First, we predict the structure of a set of self- and pathogen-derived peptides bound to HLA A2 using a previously developed ab initio structure prediction approach [Fagerberg et al., J. Mol. Biol., 521-46 (2006)]. Second, shape and electrostatic energy calculations are performed on a 3D grid to produce similarity matrices which are combined with a genetic neural network method [So et al., J. Med. Chem., 4347-59 (1997)] to generate 3D-QSAR models. The models are extensively validated using several different approaches. During the model generation, the leave-one-out cross-validated correlation coefficient (q (2)) is used as the fitness criterion and all obtained models are evaluated based on their q (2) values. Moreover, the best model obtained for a partitioned data set is evaluated by its correlation coefficient (r = 0.92 for the external test set). The physical relevance of all models is tested using a functional dependence analysis and the robustness of the models obtained for the entire data set is confirmed using y-randomization. Finally, the validated models are tested for their utility in the setting of rational peptide design: their ability to discriminate between peptides that only contain side chain substitutions in a single secondary anchor position is evaluated. In addition, the predicted cross-recognition of the mono-substituted peptides is confirmed experimentally in chromium-release assays. These results underline the utility of 3D-QSARs in peptide mimetic design and suggest that the properties of the unbound epitope are sufficient to capture most of the information to determine the cross-recognition.

Comparison of indices proposed as criteria for assigning skin notation

OBJECTIVES: Skin notations are used as a hazard identification tool to flag chemicals associated with a potential risk related to transdermal penetration. The transparency and rigorousness of the skin notation assignment process have recently been questioned. We compared different approaches proposed as criteria for these notations as a starting point for improving and systematizing current practice. METHODS: In this study, skin notations, dermal acute lethal dose 50 in mammals (LD(50)s) and two dermal risk indices derived from previously published work were compared using the lists of Swiss maximum allowable concentrations (MACs) and threshold limit values (TLVs) from the American Conference of Governmental Industrial Hygienists (ACGIH). The indices were both based on quantitative structure-activity relationship (QSAR) estimation of transdermal fluxes. One index compared the cumulative dose received through skin given specific exposure surface and duration to that received through lungs following inhalation 8 h at the MAC or TLV. The other index estimated the blood level increase caused by adding skin exposure to the inhalation route at kinetic steady state. Dermal-to-other route ratios of LD(50) were calculated as secondary indices of dermal penetrability. RESULTS: The working data set included 364 substances. Depending on the subdataset, agreement between the Swiss and ACGIH skin notations varied between 82 and 87%. Chemicals with a skin notation were more likely to have higher dermal risk indices and lower dermal LD(50) than chemicals without a notation (probabilities between 60 and 70%). The risk indices, based on cumulative dose and kinetic steady state, respectively, appeared proportional up to a constant independent of chemical-specific properties. They agreed well with dermal LD(50)s (Spearman correlation coefficients -0.42 to -0.43). Dermal-to-other routes LD(50) ratios were moderately associated with QSAR-based transdermal fluxes (Spearman correlation coefficients -0.2 to -0.3). CONCLUSIONS: The plausible but variable relationship between current skin notations and the different approaches tested confirm the need to improve current skin notations. QSAR-based risk indices and dermal toxicity data might be successfully integrated in a systematic alternative to current skin notations for detecting chemicals associated with potential dermal risk in the workplace. [Authors]

To Monty Kier, a friendly tribute.

In 2008 we published the first set of guidelines for standardizing research in autophagy. Since then, research on this topic has continued to accelerate, and many new scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Accordingly, it is important to update these guidelines for monitoring autophagy in different organisms. Various reviews have described the range of assays that have been used for this purpose. Nevertheless, there continues to be confusion regarding acceptable methods to measure autophagy, especially in multicellular eukaryotes. A key point that needs to be emphasized is that there is a difference between measurements that monitor the numbers or volume of autophagic elements (e.g., autophagosomes or autolysosomes) at any stage of the autophagic process vs. those that measure flux through the autophagy pathway (i.e., the complete process); thus, a block in macroautophagy that results in autophagosome accumulation needs to be differentiated from stimuli that result in increased autophagic activity, defined as increased autophagy induction coupled with increased delivery to, and degradation within, lysosomes (in most higher eukaryotes and some protists such as Dictyostelium) or the vacuole (in plants and fungi). In other words, it is especially important that investigators new to the field understand that the appearance of more autophagosomes does not necessarily equate with more autophagy. In fact, in many cases, autophagosomes accumulate because of a block in trafficking to lysosomes without a concomitant change in autophagosome biogenesis, whereas an increase in autolysosomes may reflect a reduction in degradative activity. Here, we present a set of guidelines for the selection and interpretation of methods for use by investigators who aim to examine macroautophagy and related processes, as well as for reviewers who need to provide realistic and reasonable critiques of papers that are focused on these processes. These guidelines are not meant to be a formulaic set of rules, because the appropriate assays depend in part on the question being asked and the system being used. In addition, we emphasize that no individual assay is guaranteed to be the most appropriate one in every situation, and we strongly recommend the use of multiple assays to monitor autophagy. In these guidelines, we consider these various methods of assessing autophagy and what information can, or cannot, be obtained from them. Finally, by discussing the merits and limits of particular autophagy assays, we hope to encourage technical innovation in the field.

The activation process of the alpha1B-adrenergic receptor: potential role of protonation and hydrophobicity of a highly conserved aspartate.

In this study, a quantitative approach was used to investigate the role of D142, which belongs to the highly conserved E/DRY sequence, in the activation process of the alpha1B-adrenergic receptor (alpha1B-AR). Experimental and computer-simulated mutagenesis were performed by substituting all possible natural amino acids at the D142 site. The resulting congeneric set of proteins together with the finding that all the receptor mutants show various levels of constitutive (agonist-independent) activity enabled us to quantitatively analyze the relationships between structural/dynamic features and the extent of constitutive activity. Our results suggest that the hydrophobic/hydrophilic character of D142, which could be regulated by protonation/deprotonation of this residue, is an important modulator of the transition between the inactive (R) and active (R*) state of the alpha1B-AR. Our study represents an example of quantitative structure-activity relationship analysis of the activation process of a G protein-coupled receptor.

Rational design of 4-aryl-1,2,3-triazoles for indoleamine 2,3-dioxygenase 1 inhibition.

Indoleamine 2,3-dioxygenase 1 (IDO1) is an important therapeutic target for the treatment of diseases such as cancer that involve pathological immune escape. Starting from the scaffold of our previously discovered IDO1 inhibitor 4-phenyl-1,2,3-triazole, we used computational structure-based methods to design more potent ligands. This approach yielded highly efficient low molecular weight inhibitors, the most active being of nanomolar potency both in an enzymatic and in a cellular assay, while showing no cellular toxicity and a high selectivity for IDO1 over tryptophan 2,3-dioxygenase (TDO). A quantitative structure-activity relationship based on the electrostatic ligand-protein interactions in the docked binding modes and on the quantum chemically derived charges of the triazole ring demonstrated a good explanatory power for the observed activities.

Chemodiversity and molecular plasticity: recognition processes as explored by property spaces.

In the last few years, a need to account for molecular flexibility in drug-design methodologies has emerged, even if the dynamic behavior of molecular properties is seldom made explicit. For a flexible molecule, it is indeed possible to compute different values for a given conformation-dependent property and the ensemble of such values defines a property space that can be used to describe its molecular variability; a most representative case is the lipophilicity space. In this review, a number of applications of lipophilicity space and other property spaces are presented, showing that this concept can be fruitfully exploited: to investigate the constraints exerted by media of different levels of structural organization, to examine processes of molecular recognition and binding at an atomic level, to derive informative descriptors to be included in quantitative structure--activity relationships and to analyze protein simulations extracting the relevant information. Much molecular information is neglected in the descriptors used by medicinal chemists, while the concept of property space can fill this gap by accounting for the often-disregarded dynamic behavior of both small ligands and biomacromolecules. Property space also introduces some innovative concepts such as molecular sensitivity and plasticity, which appear best suited to explore the ability of a molecule to adapt itself to the environment variously modulating its property and conformational profiles. Globally, such concepts can enhance our understanding of biological phenomena providing fruitful descriptors in drug-design and pharmaceutical sciences.

Structure-function relationships of the alpha1b-adrenergic receptor.

The alpha1b-adrenergic receptor (AR) is a member of the large superfamily of seven transmembrane domain (TMD) G protein-coupled receptors (GPCR). Combining site-directed mutagenesis of the alpha1b-AR with computational simulations of receptor dynamics, we have explored the conformational changes underlying the process of receptor activation, i.e. the transition between the inactive and active states. Our findings suggest that the structural constraint stabilizing the alpha1b-AR in the inactive form is a network of H-bonding interactions amongst conserved residues forming a polar pocket and R143 of the DRY sequence at the end of TMDIII. We have recently reported that point mutations of D142, of the DRY sequence and of A293 in the distal portion of the third intracellular loop resulted in ligand-independent (constitutive) activation of the alpha1b-AR. These constitutively activating mutations could induce perturbations resulting in the shift of R143 out of the polar pocket. The main role of R143 may be to mediate receptor activation by triggering the exposure of several basic amino acids of the intracellular loops towards the G protein. Our investigation has been extended also to the biochemical events involved in the desensitization process of alpha1b-AR. Our results indicate that immediately following agonist-induced activation, the alpha1b-AR can undergo rapid agonist-induced phosphorylation and desensitization. Different members of the G protein coupled receptor kinase family can play a role in agonist-induced regulation of the alpha1b-AR. In addition, constitutively active alpha1b-AR mutants display different phosphorylation and internalization features. The future goal is to further elucidate the molecular mechanism underlying the complex equilibrium between activation and inactivation of the alpha1b-AR and its regulation by pharmacological substances. These findings can help to elucidate the mechanism of action of various agents displaying properties of agonists or inverse agonists at the adrenergic system.

Catabolic ornithine carbamoyltransferase of Pseudomonas aeruginosa. Importance of the N-terminal region for dodecameric structure and homotropic carbamoylphosphate cooperativity.

Pseudomonas aeruginosa has an anabolic (ArgF) and a catabolic (ArcB) ornithine carbamoyltransferase (OTCase). Despite extensive sequence similarities, these enzymes function unidirectionally in vivo. In the dodecameric catabolic OTCase, homotropic cooperativity for carbamoylphosphate strongly depresses the anabolic reaction; the residue Glu1O5 and the C-terminus are known to be essential for this cooperativity. When Glu1O5 and nine C-terminal amino acids of the catabolic OTCase were introduced, by in vitro genetic manipulation, into the closely related, trimeric, anabolic (ArgF) OTCase of Escherichia coli, the enzyme displayed Michaelis-Menten kinetics and no cooperativity was observed. This indicates that additional amino acid residues are required to produce homotropic cooperativity and a dodecameric assembly. To localize these residues, we constructed several hybrid enzymes by fusing, in vivo or in vitro, the E. coli argF gene to the P. aeruginosa arcB gene. A hybrid enzyme consisting of 101 N-terminal ArgF amino acids fused to 233 C-terminal ArcB residues and the reciprocal ArcB-ArgF hybrid were both trimers with little or no cooperativity. Replacing the seven N-terminal residues of the ArcB enzyme by the corresponding six residues of E. coli ArgF enzyme produced a dodecameric enzyme which showed a reduced affinity for carbamoylphosphate and an increase in homotropic cooperativity. Thus, the N-terminal amino acids of catabolic OTCase are important for interaction with carbamoylphosphate, but do not alone determine dodecameric assembly. Hybrid enzymes consisting of either 26 or 42 N-terminal ArgF amino acids and the corresponding C-terminal ArcB residues were both trimeric, yet they retained some homotropic cooperativity. Within the N-terminal ArcB region, a replacement of motif 28-33 by the corresponding ArgF segment destabilized the dodecameric structure and the enzyme existed in trimeric and dodecameric states, indicating that this region is important for dodecameric assembly. These findings were interpreted in the light of the three-dimensional structure of catabolic OTCase, which allows predictions about trimer-trimer interactions. Dodecameric assembly appears to require at least three regions: the N- and C-termini (which are close to each other in a monomer), residues 28-33 and residues 147-154. Dodecameric structure correlates with high carbamoylphosphate cooperativity and thermal stability, but some trimeric hybrid enzymes retain cooperativity, and the dodecameric Glu1O5-->Ala mutant gives hyperbolic carbamoylphosphate saturation, indicating that dodecameric structure is neither necessary nor sufficient to ensure cooperativity.

Structure-function analyses point to a polynucleotide-accommodating groove essential for APOBEC3A restriction activities.

Members of the human APOBEC3 family of editing enzymes can inhibit various mobile genetic elements. APOBEC3A (A3A) can block the retrotransposon LINE-1 and the parvovirus adeno-associated virus type 2 (AAV-2) but does not inhibit retroviruses. In contrast, APOBEC3G (A3G) can block retroviruses but has only limited effects on AAV-2 or LINE-1. What dictates this differential target specificity remains largely undefined. Here, we modeled the structure of A3A based on its homology with the C-terminal domain of A3G and further compared the sequence of human A3A to those of 11 nonhuman primate orthologues. We then used these data to perform a mutational analysis of A3A, examining its ability to restrict LINE-1, AAV-2, and foreign plasmid DNA and to edit a single-stranded DNA substrate. The results revealed an essential functional role for the predicted single-stranded DNA-docking groove located around the A3A catalytic site. Within this region, amino acid differences between A3A and A3G are predicted to affect the shape of the polynucleotide-binding groove. Correspondingly, transferring some of these A3A residues to A3G endows the latter protein with the ability to block LINE-1 and AAV-2. These results suggest that the target specificity of APOBEC3 family members is partly defined by structural features influencing their interaction with polynucleotide substrates.

ASIC1a oligomerization structure-function of its extracellular vestibule and functional characterization of αS243Pβɣ ENaC

Quatre cristaux du canal ASIC1a ont été publiés et soutiennent une stoechiométrie trimérique. Cependant, ces données contredisant de précédentes analyses fonctionnelles effectuées sur des canaux de la même famille, notre intérêt fut porté sur l'oligomérisation d'ASIC1a. Dans ce sens, un nouvel essai couplant la méthode d'analyse par substitution de cystéines (SCAM) avec l'utilisation de réactifs sulfhydryls bifonctionnels (crosslinkers) a été mis en place. Le but étant de stabiliser, puis sélectionner les canaux fonctionnels, pour ensuite les séparer selon leur taille par SDS-PAGE. Grâce à cette technique, nous avons démontré que le complexe stabilisé a une taille coïncidant avec une organisation tétramérique. En plus de son oligomérisation, le chemin emprunté par les ions pour traverser le canal n'est pas clairement défini dans ces structures. De ce fait, utilisant une approche électrophysiologique, nous avons étudié le lien entre la structure et la fonction du vestibule extracellulaire d'ASIC1a. Dans ce but, nous nous sommes intéressés l'accessibilité de cystéines spécifiques localisées dans ce vestibule pour des réactifs méthanethiosulfonates (MTS). Ainsi, nous avons pu corréler les cinétiques de modification de ces cystéines par les MTS avec les effets sur le courant sodique, et donc avoir des informations supplémentaires sur la voie empruntée par les ions. De plus, la simulation informatique de liaison de ces réactifs illustre le remplissage total de ce vestibule. Fonctionnellement, cette interaction ne perturbe pas le passage de ions, c'est pourquoi il nous apparaît probable que le vestibule présente une taille plus large que celle illustrée par les cristaux. Dans un deuxième temps, notre intérêt fut porté sur ENaC. Ce canal est composé des trois sous-unités (a, ß et y) et est exprimé dans divers épithéliums, dont les tubules des reins. Il participe à l'homéostasie sodique et est essentiellement régulé par voie hormonale via l'aldostérone et la Vasopressine, mais également par des sérines protéases ou le Na+. Nous avons étudié la répercussion fonctionnelle de la mutation aS243P, découverte chez un nouveau-né prématuré atteint de pseudohypoaldostéronisme de type 1. Cette maladie autosomale récessive se caractérise, généralement, par une hyponatrémie liée à d'importantes pertes de sel dans les urines, une hyperkaliémie, ainsi qu'un niveau élevé d'aldostérone. Tout d'abord aucune des expériences biochimiques et électrophysiologiques n'a pu démontrer un défaut d'expression ou une forte diminution de l'activité soutenant les données cliniques. Cependant, en challengeant aS243PßyENaC avec une forte concentration de Na+ externe, une hypersensibilité de canal fut observée. En effet, ni les phénomènes régulateurs de « feedback inhibition » ou de « Na+ self-inhibition » n'étaient semblables au canal sauvage. De ce fait, ils apparaissaient exacerbés en présence de la mutation, amenant ainsi à une diminution de la réabsorption de Na+. Ceci corrobore entièrement l'hyponatrémie diagnostiquée. Le rein d'un prématuré étant immature, la quantité de Na+ atteignant la partie distale du néphron est plus élevée, du fait que les autres mécanismes de réabsorption en amont ne sont probablement pas encore en place. Cette hypothèse est renforcée par l'existence d'un frère présentant la même mutation, mais qui, né à terme, ne présentait aucun signe d'hyponatrémie. - The main topic of my thesis is the structure-function relationship of the ENaC/Deg family of ion channels, namely the Acid-Sensing Ion Channel ASIC1a and the Epithelial Na Channel ENaC. The primary part of this research is dedicated to the structure of ASIC1a. Four channel crystals have been published, which support a trimeric stoichiometry, although these data contradict previous functional experiments on other ENaC/Deg members. We are therefore interested in ASIC1a oligomerization and have set up a new assay combining the Substituted- Cysteine Accessibility Method (SCAM) with Afunctional sulfhydryl reagents (crosslinkers) allowing its study. The aim was to first stabilize the channels, then select those that are functional and then resolve them according to their size on SDS-PAGE. We demonstrated that the stabilized complex has a molecular weight corresponding to a tetrameric stoichiometry. In addition to our interest in the oligomerization of the ENaC/Deg family of ion channels, we also wanted to investigate the thus far undefined way of permeation for these channels. Therefore, taking the advantage of a more electrophysiological approach, we studied the accessibility of specific cysteines for methanethiosulfonate reagents (MTS) and were able to correlate the MTS association kinetics on cysteine residues with Na+ currents. These results have given us an insight into ion permeation and our functional evidence indicates that the extracellular is larger than that depicted by the crystal structures. As a side project, we focused on ENaC, which is made up of three subunits (a, ß and y) and is expressed in various epithelia, especially in the distal nephron of the kidneys. It plays a role in Na+ homeostasis and is essentially regulated by hormones via aldosterone and vasopressin, but also by serine proteases or Na+. We have studied the functional impact of the aS243P mutation, discovered in a premature baby suffering from pseudohypoaldosteronism of type 1. This autosomal recessive disease is characterized by hyponatremia, hyperkalemia and high aldosterone levels. Firstly, neither biochemical nor electrophysiological experiments indicated an expression defect or a strong decrease in activity. However, challenging aS243PßyENaC with increased external Na+ concentration showed channel hypersensitivity. Indeed, both the "feedback inhibition" and the "Na+ self-inhibition" regulatory mechanisms are impaired, leading to a decrease in Na+ reabsorption, entirely supports the diagnosis. The kidneys in preterm infants are immature and Na+ levels reaching the distal nephron are higher than normally observed. We hypothesize that the upstream reabsorption machinery is unlikely to be sufficiently matured and this assumption is supported by an asymptomatic sibling carrying the same mutation, but born at term. - La cellule, unité fonctionnelle du corps humain, est délimitée par une membrane plasmique servant de barrière biologique entre les milieux intra et extracellulaires. Une communication entre cellules est indispensable pour un fonctionnement adéquat. Sa survie dépend, entre autres, du maintien de la teneur en ions dans chacun des milieux qui doivent pouvoir être réabsorbés, ou sécrétés, selon les besoins. Les protéines insérées dans la membrane forment un canal et sont un moyen de communication permettant spécifiquement à des ions tel que le sodium (Na+) de traverser. Le Na+ se trouve dans la plupart des aliments et le sel, et est spécifiquement réabsorbé au niveau des reins grâce au canal sodique épithélial ENaC. Cette réabsorption se fait de l'urine primaire vers l'intérieur de la cellule, puis est transporté vers le sang. Pour maintenir un équilibre, une régulation de ce canal est nécessaire. En effet, des dysfonctionnements impliquant la régulation ou l'activité d'ENaC lui-même sont à l'origine de maladies telles que la mucoviscidose, l'hypertension ou encore, le pseudohypoaldostéronisme (PHA). Cette maladie est caractérisée, notamment, par d'importantes pertes de sel dans les urines. Des pédiatres ont diagnostiqué un PHA chez un nouveau-né, ce dernier présentant une modification du canal ENaC, nous avons recréé cette protéine afin d'étudier l'impact de ce changement sur son activité. Nous avons démontré que la régulation d'ENaC était effectivement perturbée, conduisant ainsi à une forte réduction de la réabsorption sodique. Afin de développer des molécules capables de moduler l'activité de protéines. Il est nécessaire d'en connaître la structure. Celle du canal sodique sensible à l'acidification ASIC1, un canal cousin d'ENaC, est connue. Ces données structurales contredisant cependant les analyses fonctionnelles, nous nous sommes penchés une nouvelle fois sur ASIC1. Une protéine est une macromolécule biologique composée d'une chaîne d'acides aminés (aa). De l'enchaînement d'aa à la protéine fonctionnelle, quatre niveaux de structuration existent. Chaque aa donne une indication quant au repliement et plus particulièrement la cystéine. Arborant un groupe sulfhydryle (SH) capable de former une liaison spécifique et stable avec un autre SH, celle-ci est souvent impliquée dans la structure tridimensionnelle de la protéine. Ce type de liaison intervient également dans la stabilisation de la structure quaternaire, qui est l'association de plusieurs protéines identiques (homomère), ou pas (hétéromère). Dans cette partie, nous avons remplacé des aa par des cystéines à des endroits spécifiques. Le but était de stabiliser plusieurs homomères d'ASICl ensemble avec des réactifs créant des ponts entre deux SH. Ainsi, nous avons pu déterminer le nombre de protéines ASIC1 participant à la formation d'un canal fonctionnel. Nos résultats corroborent les données fonctionnelles soutenant un canal tétramérique. Nous avons également étudié l'accessibilité de ces nouvelles cystéines afin d'obtenir des informations supplémentaires sur la structure du chemin emprunté par le Na+ à travers ASIC1 et plus particulièrement du vestibule extracellulaire.

Hemoglobin enhances the biological activity of synthetic and natural bacterial (endotoxic) virulence factors: a general principle.

Although hemoglobin (Hb) is mainly present in the cytoplasm of erythrocytes (red blood cells), lower concentrations of pure, cell-free Hb are released permanently into the circulation due to an inherent intravascular hemolytic disruption of erythrocytes. Previously it was shown that the interaction of Hb with bacterial endotoxins (lipopolysaccharides, LPS) results in a significant increase of the biological activity of LPS. There is clear evidence that the enhancement of the biological activity of LPS by Hb is connected with a disaggregation of LPS. From these findings one questions whether the property to enhance the biological activity of endotoxin, in most cases proven by the ability to increase the cytokine (tumor-necrosis-factor-alpha, interleukins) production in human mononuclear cells, is restricted to bacterial endotoxin or is a more general principle in nature. To elucidate this question, we investigated the interaction of various synthetic and natural virulence (pathogenicity) factors with hemoglobin of human or sheep origin. In addition to enterobacterial R-type LPS a synthetic bacterial lipopeptide and synthetic phospholipid-like structures mimicking the lipid A portion of LPS were analysed. Furthermore, we also tested endotoxically inactive LPS and lipid A compounds such as those from Chlamydia trachomatis. We found that the observations made for endotoxically active form of LPS can be generalized for the other synthetic and natural virulence factors: In every case, the cytokine-production induced by them is increased by the addition of Hb. This biological property of Hb is connected with its physical property to convert the aggregate structures of the virulence factors into one with cubic symmetry, accompanied with a considerable reduction of the size and number of the original aggregates.

Structure-function characterization and optimization of a plant-derived antibacterial peptide.

Crushed seeds of the Moringa oleifera tree have been used traditionally as natural flocculants to clarify drinking water. We previously showed that one of the seed peptides mediates both the sedimentation of suspended particles such as bacterial cells and a direct bactericidal activity, raising the possibility that the two activities might be related. In this study, the conformational modeling of the peptide was coupled to a functional analysis of synthetic derivatives. This indicated that partly overlapping structural determinants mediate the sedimentation and antibacterial activities. Sedimentation requires a positively charged, glutamine-rich portion of the peptide that aggregates bacterial cells. The bactericidal activity was localized to a sequence prone to form a helix-loop-helix structural motif. Amino acid substitution showed that the bactericidal activity requires hydrophobic proline residues within the protruding loop. Vital dye staining indicated that treatment with peptides containing this motif results in bacterial membrane damage. Assembly of multiple copies of this structural motif into a branched peptide enhanced antibacterial activity, since low concentrations effectively kill bacteria such as Pseudomonas aeruginosa and Streptococcus pyogenes without displaying a toxic effect on human red blood cells. This study thus identifies a synthetic peptide with potent antibacterial activity against specific human pathogens. It also suggests partly distinct molecular mechanisms for each activity. Sedimentation may result from coupled flocculation and coagulation effects, while the bactericidal activity would require bacterial membrane destabilization by a hydrophobic loop.

BAFF, APRIL and their receptors: structure, function and signaling.

BAFF, APRIL and their receptors play important immunological roles, especially in the B cell arm of the immune system. A number of splice isoforms have been described for both ligands and receptors in this subfamily, some of which are conserved between mouse and human, while others are species-specific. Structural and mutational analyses have revealed key determinants of receptor-ligand specificity. BAFF-R has a strong selectivity for BAFF; BCMA has a higher affinity for APRIL than for BAFF, while TACI binds both ligands equally well. The molecular signaling events downstream of BAFF-R, BCMA and TACI are still incompletely characterized. Survival appears to be mediated by upregulation of Bcl-2 family members through NF-kappaB activation, degradation of the pro-apototic Bim protein, and control of subcellular localization of PCKdelta. Very little is known about other signaling events associated with receptor engagement by BAFF and APRIL that lead for example to B cell activation or to CD40L-independent Ig switch.

Cheminformatic studies of molecular properties and their influence on bioactivity in the context of drug discovery

Molecular shape has long been known to be an important property for the process of molecular recognition. Previous studies postulated the existence of a drug-like shape space that could be used to artificially bias the composition of screening libraries, with the aim to increase the chance of success in Hit Identification. In this work, it was analysed to which extend this assumption holds true. Normalized Principal Moments of Inertia Ratios (NPRs) have been used to describe the molecular shape of small molecules. It was investigated, whether active molecules of diverse targets are located in preferred subspaces of the NPR shape space. Results illustrated a significantly stronger clustering than could be expected by chance, with parts of the space unlikely to be occupied by active compounds. Furthermore, a strong enrichment of elongated, rather flat shapes could be observed, while globular compounds were highly underrepresented. This was confirmed for a wide range of small molecule datasets from different origins. Active compounds exhibited a high overlap in their shape distributions across different targets, making a purely shape­ based discrimination very difficult. An additional perspective was provided by comparing the shapes of protein binding pockets with those of their respective ligands. Although more globular than their ligands, it was observed that binding sites shapes exhibited a similarly skewed distribution in shape space: spherical shapes were highly underrepresented. This was different for unoccupied binding pockets of smaller size. These were on the contrary identified to possess a more globular shape. The relation between shape complementarity and exhibited bioactivity was analysed; a moderate correlation between bioactivity and parameters including pocket coverage, distance in shape space, and others could be identified, which reflects the importance of shape complementarity. However, this also suggests that other aspects are of relevance for molecular recognition. A subsequent analysis assessed if and how shape and volume information retrieved from pocket or respective reference ligands could be used as a pre-filter in a virtual screening approach. ln Lead Optimization compounds need to get optimized with respect to a variety of pararneters. Here, the availability of past success stories is very valuable, as they can guide medicinal chemists during their analogue synthesis plans. However, although of tremendous interest for the public domain, so far only large corporations had the ability to mine historical knowledge in their proprietary databases. With the aim to provide such information, the SwissBioisostere database was developed and released during this thesis. This database contains information on 21,293,355 performed substructural exchanges, corresponding to 5,586,462 unique replacements that have been measured in 35,039 assays against 1,948 molecular targets representing 30 target classes, and on their impact on bioactivity . A user-friendly interface was developed that provides facile access to these data and is accessible at http//www.swissbioisostere.ch. The ChEMBL database was used as primary data source of bioactivity information. Matched molecular pairs have been identified in the extracted and cleaned data. Success-based scores were developed and integrated into the database to allow re-ranking of proposed replacements by their past outcomes. It was analysed to which degree these scores correlate with chemical similarity of the underlying fragments. An unexpectedly weak relationship was detected and further investigated. Use cases of this database were envisioned, and functionalities implemented accordingly: replacement outcomes are aggregatable at the assay level, and it was shawn that an aggregation at the target or target class level could also be performed, but should be accompanied by a careful case-by-case assessment. It was furthermore observed that replacement success depends on the activity of the starting compound A within a matched molecular pair A-B. With increasing potency the probability to lose bioactivity through any substructural exchange was significantly higher than in low affine binders. A potential existence of a publication bias could be refuted. Furthermore, often performed medicinal chemistry strategies for structure-activity-relationship exploration were analysed using the acquired data. Finally, data originating from pharmaceutical companies were compared with those reported in the literature. It could be seen that industrial medicinal chemistry can access replacement information not available in the public domain. In contrast, a large amount of often-performed replacements within companies could also be identified in literature data. Preferences for particular replacements differed between these two sources. The value of combining different endpoints in an evaluation of molecular replacements was investigated. The performed studies highlighted furthermore that there seem to exist no universal substructural replacement that always retains bioactivity irrespective of the biological environment. A generalization of bioisosteric replacements seems therefore not possible. - La forme tridimensionnelle des molécules a depuis longtemps été reconnue comme une propriété importante pour le processus de reconnaissance moléculaire. Des études antérieures ont postulé que les médicaments occupent préférentiellement un sous-ensemble de l'espace des formes des molécules. Ce sous-ensemble pourrait être utilisé pour biaiser la composition de chimiothèques à cribler, dans le but d'augmenter les chances d'identifier des Hits. L'analyse et la validation de cette assertion fait l'objet de cette première partie. Les Ratios de Moments Principaux d'Inertie Normalisés (RPN) ont été utilisés pour décrire la forme tridimensionnelle de petites molécules de type médicament. Il a été étudié si les molécules actives sur des cibles différentes se co-localisaient dans des sous-espaces privilégiés de l'espace des formes. Les résultats montrent des regroupements de molécules incompatibles avec une répartition aléatoire, avec certaines parties de l'espace peu susceptibles d'être occupées par des composés actifs. Par ailleurs, un fort enrichissement en formes allongées et plutôt plates a pu être observé, tandis que les composés globulaires étaient fortement sous-représentés. Cela a été confirmé pour un large ensemble de compilations de molécules d'origines différentes. Les distributions de forme des molécules actives sur des cibles différentes se recoupent largement, rendant une discrimination fondée uniquement sur la forme très difficile. Une perspective supplémentaire a été ajoutée par la comparaison des formes des ligands avec celles de leurs sites de liaison (poches) dans leurs protéines respectives. Bien que plus globulaires que leurs ligands, il a été observé que les formes des poches présentent une distribution dans l'espace des formes avec le même type d'asymétrie que celle observée pour les ligands: les formes sphériques sont fortement sous­ représentées. Un résultat différent a été obtenu pour les poches de plus petite taille et cristallisées sans ligand: elles possédaient une forme plus globulaire. La relation entre complémentarité de forme et bioactivité a été également analysée; une corrélation modérée entre bioactivité et des paramètres tels que remplissage de poche, distance dans l'espace des formes, ainsi que d'autres, a pu être identifiée. Ceci reflète l'importance de la complémentarité des formes, mais aussi l'implication d'autres facteurs. Une analyse ultérieure a évalué si et comment la forme et le volume d'une poche ou de ses ligands de référence pouvaient être utilisés comme un pré-filtre dans une approche de criblage virtuel. Durant l'optimisation d'un Lead, de nombreux paramètres doivent être optimisés simultanément. Dans ce contexte, la disponibilité d'exemples d'optimisations réussies est précieuse, car ils peuvent orienter les chimistes médicinaux dans leurs plans de synthèse par analogie. Cependant, bien que d'un extrême intérêt pour les chercheurs dans le domaine public, seules les grandes sociétés pharmaceutiques avaient jusqu'à présent la capacité d'exploiter de telles connaissances au sein de leurs bases de données internes. Dans le but de remédier à cette limitation, la base de données SwissBioisostere a été élaborée et publiée dans le domaine public au cours de cette thèse. Cette base de données contient des informations sur 21 293 355 échanges sous-structuraux observés, correspondant à 5 586 462 remplacements uniques mesurés dans 35 039 tests contre 1948 cibles représentant 30 familles, ainsi que sur leur impact sur la bioactivité. Une interface a été développée pour permettre un accès facile à ces données, accessible à http:/ /www.swissbioisostere.ch. La base de données ChEMBL a été utilisée comme source de données de bioactivité. Une version modifiée de l'algorithme de Hussain et Rea a été implémentée pour identifier les Matched Molecular Pairs (MMP) dans les données préparées au préalable. Des scores de succès ont été développés et intégrés dans la base de données pour permettre un reclassement des remplacements proposés selon leurs résultats précédemment observés. La corrélation entre ces scores et la similarité chimique des fragments correspondants a été étudiée. Une corrélation plus faible qu'attendue a été détectée et analysée. Différents cas d'utilisation de cette base de données ont été envisagés, et les fonctionnalités correspondantes implémentées: l'agrégation des résultats de remplacement est effectuée au niveau de chaque test, et il a été montré qu'elle pourrait également être effectuée au niveau de la cible ou de la classe de cible, sous réserve d'une analyse au cas par cas. Il a en outre été constaté que le succès d'un remplacement dépend de l'activité du composé A au sein d'une paire A-B. Il a été montré que la probabilité de perdre la bioactivité à la suite d'un remplacement moléculaire quelconque est plus importante au sein des molécules les plus actives que chez les molécules de plus faible activité. L'existence potentielle d'un biais lié au processus de publication par articles a pu être réfutée. En outre, les stratégies fréquentes de chimie médicinale pour l'exploration des relations structure-activité ont été analysées à l'aide des données acquises. Enfin, les données provenant des compagnies pharmaceutiques ont été comparées à celles reportées dans la littérature. Il a pu être constaté que les chimistes médicinaux dans l'industrie peuvent accéder à des remplacements qui ne sont pas disponibles dans le domaine public. Par contre, un grand nombre de remplacements fréquemment observés dans les données de l'industrie ont également pu être identifiés dans les données de la littérature. Les préférences pour certains remplacements particuliers diffèrent entre ces deux sources. L'intérêt d'évaluer les remplacements moléculaires simultanément selon plusieurs paramètres (bioactivité et stabilité métabolique par ex.) a aussi été étudié. Les études réalisées ont souligné qu'il semble n'exister aucun remplacement sous-structural universel qui conserve toujours la bioactivité quel que soit le contexte biologique. Une généralisation des remplacements bioisostériques ne semble donc pas possible.

PPAR alpha structure-function relationships derived from species-specific differences in responsiveness to hypolipidemic agents.

The nuclear receptor PPAR alpha is a key regulatory transcription factor in lipid homeostasis, some liver detoxification processes and the control of inflammation. Recent findings suggest that many hypolipidemic drugs and anti-inflammatory agents can potentially act by binding to PPAR alpha and inducing its activity. Here, we identify some structure-function relationships in PPAR alpha, by using the species-specific responsiveness to the two hypolipidemic agents, Wy 14,643 and 5,8,11,14-eicosatetraynoic acid (ETYA). We first show that the species-specific differences are mediated primarily via the ligand binding domain of the receptor and that these two drugs are indeed ligands of PPAR alpha. By mutagenesis analyses we identify amino acid residues in the ligand binding domains of Xenopus, mouse and human PPAR alpha, that confer preferential responsiveness to ETYA and Wy 14,643. These findings will aid in the development of new synthetic PPAR alpha ligands as effective therapeutics for lipid-related diseases and inflammatory disorders.