In this thesis several important proteins are investigated from a structural perspective. Some of the proteins are disease related while other have important but not completely characterised functions. The techniques used are general as demonstrated by applications on metabolic proteins (CYP21, CYP11B1, IAPP, ADH3), regulatory proteins (p53, GDNF) and a transporter protein (ANTR1). When the protein CYP21 (steroid 21-hydroxylase) is deficient it causes CAH (congenital adrenal hyperplasia). For this protein, there are about 60 known mutations with characterised clinical phenotypes. Using manual structural analysis we managed to explain the severity of all but one of the mutations. By observing the properties of these mutations we could perform good predictions on, at the time, not classified mutations. For the cancer suppressor protein p53, there are over thousand mutations with known activity. To be able to analyse such a large number of mutations we developed an automated method for evaluation of the mutation effect called PREDMUT. In this method we include twelve different prediction parameters including two energy parameters calculated using an energy minimization procedure. The method manages to differentiate severe mutations from non-severe mutations with 77% accuracy on all possible single base substitutions and with 88% on mutations found in breast cancer patients. The automated prediction was further applied to CYP11B1 (steroid 11-beta-hydroxylase), which in a similar way as CYP21 causes CAH when deficient. A generalized method applicable to any kind of globular protein was developed. The method was subsequently evaluated on nine additional proteins for which mutants were known with annotated disease phenotypes. This prediction achieved 84% accuracy on CYP11B1 and 81% accuracy in total on the evaluation proteins while leaving 8% as unclassified. By increasing the number of unclassified mutations the accuracy of the remaining mutations could be increased on the evaluation proteins and substantially increase the classification quality as measured by the Matthews correlation coefficient. Servers with predictions for all possible single based substitutions are provided for p53, CYP21 and CYP11B1. The amyloid formation of IAPP (islet amyloid polypeptide) is strongly connected to diabetes and has been studied using both molecular dynamics and Monte Carlo energy minimization. The effects of mutations on the amount and speed of amyloid formation were investigated using three approaches. Applying a consensus of the three methods on a number of interesting mutations, 94% of the mutations could be correctly classified as amyloid forming or not, evaluated with in vitro measurements. In the brain there are many proteins whose functions and interactions are largely unknown. GDNF (glial cell line-derived neurotrophic factor) and NCAM (neural cell adhesion molecule) are two such neuron connected proteins that are known to interact. The form of interaction was studied using protein--protein docking where a docking interface was found mediated by four oppositely charged residues in respective protein. This interface was subsequently confirmed by mutagenesis experiments. The NCAM dimer interface upon binding to the GDNF dimer was also mapped as well as an additional interacting protein, GFRα1, which was successfully added to the protein complex without any clashes. A large and well studied protein family is the alcohol dehydrogenase family, ADH. A class of this family is ADH3 (alcohol dehydrogenase class III) that has several known substrates and inhibitors. By using virtual screening we tried to characterize new ligands. As some ligands were already known we could incorporate this knowledge when the compound docking simulations were scored and thereby find two new substrates and two new inhibitors which were subsequently successfully tested in vitro. ANTR1 (anion transporter 1) is a membrane bound transporter important in the photosynthesis in plants. To be able to study the amino acid residues involved in inorganic phosphate transportation a homology model of the protein was created. Important residues were then mapped onto the structure using conservation analysis and we were in this way able to propose roles of amino acid residues involved in the transportation of inorganic phosphate. Key residues were subsequently mutated in vitro and a transportation process could be postulated. To conclude, we have used several molecular modelling techniques to find functional clues, interaction sites and new ligands. Furthermore, we have investigated the effect of muations on the function and structure of a multitude of disease related proteins.

Depression is a major psychological disorder that affects a person's mental and physical abilities. The National Institute of Mental Health (NIMH) classified it as a serious medical illness. It causes huge economic, as well as financial impact on the people, and it is also becoming a major public health issue. Antidepressant drugs are prescribed to mitigate the suffering caused by this disorder. Different generations of antidepressants have been developed with dissimilar mechanisms of action. According to the Center for Disease Control, the usage of antidepressants has skyrocketed by 400 percent increase over 2005- 2008 survey period. This dramatic rise in usage indicates that these are the most prescribed drugs in the US. Even with the FDA mandated "black box" warning of increased suicidal thoughts upon use of selected antidepressants, these drugs are still being used at a higher rate. <br>All classes of antidepressants are plagued by side effects with mainly sexual dysfunction common among them. To avoid the adverse effects, an emphasis is to discover novel structural drug scaffolds that can be further developed as a new generation of antidepressants. The importance of this research is to discover structurally novel antidepressants by performing in silico virtual screening (VS) of chemical databases using the serotonin transporter (SERT). In the absence of a SERT crystal structure, a homology model was developed. The homology model was utilized to develop the first structure-based pharmacophore for the extracellular facing secondary ligand binding pocket. The pharmacophore captured the necessary drug-SERT interaction pattern for SERT inhibitory action. This pharmacophore was employed as one of the filters for VS of candidate ligands. The ten compounds identified were purchased and tested pharmacologically. Out of the ten hits, three structurally novel ligands were identified as lead compounds. Two of these compounds exhibited selectivity towards SERT; the remaining lead compound was selective towards the dopamine transporter and displayed cocaine inhibition. The two SERT selective compounds will provide new opportunities in the development of novel therapeutics to treat depression. <br>For dopamine transporter (DAT), the study was based on recently developed structurally diverse photo probes. In an effort to better understand the binding profile similarities among these different scaffolds, the photo probes were docked into DAT. The finger print analysis of the interaction pattern of docked poses was performed to identify the inhibitor-binding sites. <br>For mitogen activated protein kinase 5 (MEK5), given the lack of structural information, a homology model of MEK5 was developed to guide the rational design of inhibitors. Docking of known MEK5 inhibitors into the homology model was performed to understand the inhibitory interaction profile. Several series of analogues were designed utilizing the generated interaction profile. / Bayer School of Natural and Environmental Sciences / Chemistry and Biochemistry / PhD / Dissertation

The D1-class receptors (D1R, D5R) each possess distinct signaling characteristics; however, pharmacological selectivity between them remains elusive. The third intracellular loops (IL3) of D1R and D5R harbour divergent residues that may contribute to their individual signalling phenotypes. Here we probe the function of central region of IL3 of D1R and D5R using deletion mutagenesis. Radioligand binding and whole cell cAMP assays suggest that the N-terminal and C-terminal moieties of the central IL3 oppositely contribute to the constitutive and agonist-dependant activity of D1-Class receptors. Whereas the N-terminal deletions ablated constitutive activity and decreased DA-induced activation, C-terminal deletions induced robust increases. These data, interpreted in concert with structural predictions generated from homology modeling implicate the central IL3 as playing an important role in the activation and subtype-specific characteristics of the D1-class receptors. This study may serve as a basis for the development of novel drugs targeting the central IL3 region.

<p>Rational conversion of biologically active peptides to nonpeptide compounds with retained activity is an appealing approach in drug development. One important objective of the work presented in this thesis was to use computational modeling to aid in such a conversion of the peptide angiotensin II (Ang II, Asp-Arg-Val-Tyr-Ile-His-Pro-Phe). An equally important objective was to gain an understanding of the requirements for ligand binding to the Ang II receptors, with a focus on interactions with the AT₂ receptor.</p><p>The bioactive conformation of a peptide can provide important guidance in peptidomimetic design. By designing and introducing well-defined secondary structure mimetics into Ang II the bioactive conformation can be addressed. In this work, both γ- and β-turn mimetic scaffolds have been designed and characterized for incorporation into Ang II. Using conformational analysis and the pharmacophore recognition method DISCO, a model was derived of the binding mode of the pseudopeptide Ang II analogues. This model indicated that the positioning of the Arg side chain was important for AT₂ receptor binding, which was also supported when the structure–activity relationship of Ang II was investigated by performing a glycine scan.</p><p>To further examine ligand binding, a 3D model of the AT₂ receptor was constructed employing homology modeling. Using this receptor model in a docking study of the ligands, binding modes were identified that were in agreement with data from point-mutation studies of the AT₂ receptor.</p><p>By investigating truncated Ang II analogues, small pseudopeptides were developed that were structurally similar to nonpeptide AT₂ receptor ligands. For further guidance in ligand design of nonpeptide compounds, three-dimensional quantitative structure–activity relationship models for AT₁ and AT₂ receptor affinity as well as selectivity were derived. </p>

Rational conversion of biologically active peptides to nonpeptide compounds with retained activity is an appealing approach in drug development. One important objective of the work presented in this thesis was to use computational modeling to aid in such a conversion of the peptide angiotensin II (Ang II, Asp-Arg-Val-Tyr-Ile-His-Pro-Phe). An equally important objective was to gain an understanding of the requirements for ligand binding to the Ang II receptors, with a focus on interactions with the AT2 receptor. The bioactive conformation of a peptide can provide important guidance in peptidomimetic design. By designing and introducing well-defined secondary structure mimetics into Ang II the bioactive conformation can be addressed. In this work, both γ- and β-turn mimetic scaffolds have been designed and characterized for incorporation into Ang II. Using conformational analysis and the pharmacophore recognition method DISCO, a model was derived of the binding mode of the pseudopeptide Ang II analogues. This model indicated that the positioning of the Arg side chain was important for AT2 receptor binding, which was also supported when the structure–activity relationship of Ang II was investigated by performing a glycine scan. To further examine ligand binding, a 3D model of the AT2 receptor was constructed employing homology modeling. Using this receptor model in a docking study of the ligands, binding modes were identified that were in agreement with data from point-mutation studies of the AT2 receptor. By investigating truncated Ang II analogues, small pseudopeptides were developed that were structurally similar to nonpeptide AT2 receptor ligands. For further guidance in ligand design of nonpeptide compounds, three-dimensional quantitative structure–activity relationship models for AT1 and AT2 receptor affinity as well as selectivity were derived.

The D1-class receptors (D1R, D5R) each possess distinct signaling characteristics; however, pharmacological selectivity between them remains elusive. The third intracellular loops (IL3) of D1R and D5R harbour divergent residues that may contribute to their individual signalling phenotypes. Here we probe the function of central region of IL3 of D1R and D5R using deletion mutagenesis. Radioligand binding and whole cell cAMP assays suggest that the N-terminal and C-terminal moieties of the central IL3 oppositely contribute to the constitutive and agonist-dependant activity of D1-Class receptors. Whereas the N-terminal deletions ablated constitutive activity and decreased DA-induced activation, C-terminal deletions induced robust increases. These data, interpreted in concert with structural predictions generated from homology modeling implicate the central IL3 as playing an important role in the activation and subtype-specific characteristics of the D1-class receptors. This study may serve as a basis for the development of novel drugs targeting the central IL3 region.

The D1-class receptors (D1R, D5R) each possess distinct signaling characteristics; however, pharmacological selectivity between them remains elusive. The third intracellular loops (IL3) of D1R and D5R harbour divergent residues that may contribute to their individual signalling phenotypes. Here we probe the function of central region of IL3 of D1R and D5R using deletion mutagenesis. Radioligand binding and whole cell cAMP assays suggest that the N-terminal and C-terminal moieties of the central IL3 oppositely contribute to the constitutive and agonist-dependant activity of D1-Class receptors. Whereas the N-terminal deletions ablated constitutive activity and decreased DA-induced activation, C-terminal deletions induced robust increases. These data, interpreted in concert with structural predictions generated from homology modeling implicate the central IL3 as playing an important role in the activation and subtype-specific characteristics of the D1-class receptors. This study may serve as a basis for the development of novel drugs targeting the central IL3 region.

The venom of cone snails (predatory marine molluscs of the genus Conus) has yielded a rich source of novel neuroactive peptides or “conotoxins”. Conotoxins are bioactive peptides found in the venom duct of Conus spp. Like other neuropeptides, conotoxins are expressed as propeptides that undergo posttranslational proteolytic processing. Peptides derived from propeptides are typically cleaved at a pair of dibasic residues (Lys-Arg, Arg-Arg, Lys-Lys or Arg-Lys) by proteases found in secretory vesicles. However, many precursor peptides contain multiple sets of basic residues, suggesting that highly substrate specific or differentially expressed proteases can determine processing outcomes. As many of the substrate-specific proteases remain unidentified, predicting new bioactive peptides from cDNA sequences is presently difficult, if not impossible. In order to understand more about the substrate specificity of conotoxin substrate-specific proteases a characterisation study of one such endoprotease isolated from the venom duct of Conus textile was undertaken. The C. textile mollusc was chosen as a good source from which to isolate the endoprotease for two reasons; firstly, these cone shells are found in great abundance on the Great Barrier Reef (Queensland, Australia) and are readily obtainable and secondly, a number of conotoxin precursors and their cleavage products have been previously identified in the venom duct. In order to purify the endoprotease an activity-guided fractionation protocol that included a para-nitroanilide (p-NA) substrate assay was developed. The p-NA substrate mimicked the cleavage site of the conotoxin TxVIA, a member of the C. textile O-superfamily of toxins. The protocol included a number of chromatographic techniques including ion exchange, size-exclusion and reverse-phased HPLC and resulted in isolation of an active protease, termed Tex31, to >95% purity. The purification of microgram quantities of Tex31 made it possible to characterise the proteolytic nature of Tex31 and to further characterise the O-superfamily conopeptide propeptide cleavage site specificity. Specificity experiments showed Tex31 requires a minimum of four residues including a leucine in the P4 position (LNKR↓) for efficient substrate processing. The complete sequence of Tex31 was determined from cDNA. A BLAST search revealed Tex31 to have high amino acid sequence similarity to the CAP (abbreviated from CRISP (Cysteine-rich secretory protein), Antigen 5 and PR-1 (pathogenesis-related protein)) superfamily and most closely related to the CRISP family of mammalian and venom proteins that, like Tex31, have a cysteine-rich C-terminal domain. The CAP superfamily is widely distributed in the animal, plant and fungal kingdoms, and is implicated in processes as diverse as human brain tumour growth and plant pathogenesis. This is the first report of a biological role for the N-terminal domain of CAP proteins. A homology model of Tex31 constructed from two PR-1 proteins, Antigen 5 and P14a, revealed the highly conserved and likely catalytic residues, His78, Ser99 and Glu115. These three amino acids fall within a structurally conserved N-terminal domain found in all CAP proteins. It is possible that other CAP proteins are also substrate-specific proteases. With no homology to any known proteases, Tex31 may belong to a new class of protease. The sequence alignment of five Tex31-like proteins cloned from C. marmoreus, C. litteratus, C. arentus, C. planboris, and C. omaria show very high sequence similarity to Tex31 (~80%), but only one weakly conserved serine residue was identified when the conserved residues of the new Tex31-like protein sequences were aligned with members of the CAP superfamily. Future work to identify members of catalytic diad or triad, e.g. by site-directed mutagenesis, will rely on the expression of active recombinant Tex31. In this study neither Escherichia coli nor Pichia pastoris expression systems yielded active recombinant Tex31 protein, possibly due to the number of cysteine residues hindering the expression of correctly folded active Tex31. This study has shown Tex31 to be highly sequence specific in its cleavage site and it is likely that this high substrate specificity has confounded previous attempts to identify the proteolytic nature of other CAP proteins. With the proteolytic nature of one member of the CAP protein family confirmed, it is hoped this important discovery may lead the way to discovering the role of other CAP family members.

En l’absence de structures tridimensionnelles expérimentales des cytochromes P450 CYP26A1, CYP26B1 et CYP26C1, la caractérisation de leur substrats et ligands s’est basée sur l’analyse des modèles structuraux obtenus par modélisation par homologie avec la structure expérimentale du cytochrome P450 CYP120. La justesse des modèles a été validée par l’amarrage de l’acide rétinoïque all-trans dans des configurations compatibles avec les métabolites attendus. L’amarrage d’agonistes et d’antagonistes des récepteurs nucléaires RARs prédirent l’acide tazaroténique (TA) et l’adapalène comme des substrats potentiels. Les expériences in vitro confirmèrent la métabolisation de ces 2 médicaments par les CYP26s. L’analyse de la cinétique de sulfoxidation du TA par CYP26A1 and CYP26B1 a permis d’établir le TA comme la référence contrôle de l’activité de ces enzymes. Puis, la comparaison des modèles des CYP26s avec la structure cristalline de CYP2C8 a permis d’identifier des similarités structurales de leurs inhibiteurs. Une corrélation entre l’inhibition de CYP26A1 et de CYP2C8 par des inhibiteurs connus de CYP2C8 a été démontrée après détermination de leurs IC50 pour CYP26A1 et CYP26B1 en utilisant le TA comme substrat de référence. La mesure de l’inhibition in vitro fut ensuite utilisée pour évaluer la possibilité que les CYP26s soient impliquées dans des interactions médicamenteuses observées pour certaines molécules. Cette thèse caractérise et appuie le rôle encore mal connu des CYP26s dans la métabolisation in vivo de certains xénobiotiques ainsi que l’effet potentiel de leur inhibition qui favoriserait la survenue d'effets indésirables. / Without crystal structures to study the CYP26 family of drug metabolizing enzymes, homology models were used to characterize CYP26A1, CYP26B1 and CYP26C1 and to identify substrates and inhibitors of the enzymes. Computational models of each isoform based on structural homology to CYP120 were validated by docking all-trans retinoic acid, an endogenous ligand of CYP26. Docking of retinoic acid receptor agonists and antagonists suggested that tazarotenic acid (TA) and adapalene may be metabolic substrates for CYP26, data which was confirmed using in vitro metabolite identification assays. Phenotyping experiments determined that CYP26s played a major role in the metabolism of these compounds in vitro. The kinetics of TA sulfoxidation by CYP26A1 and CYP26B1 were characterized and the compound was proposed as an in vitro probe of CYP26 activity in single enzyme expression systems. Structural characterization efforts identified similarities between the CYP26 homology models and the known crystal structure of CYP2C8, in agreement with previously published reports. Using TA as a probe, the IC50’s of known CYP2C8 inhibitors was measured against CYP26A1 and CYP26B1, with a statistically significant correlation observed between CYP26A1 and CYP2C8. Additional in vitro and computational experiments were used to characterize the inhibition mechanism for the most potent inhibitors. The observed in vitro inhibition was then used to predict the likelihood of CYP26 inhibition being involved in clinically relevant drug interactions. As a whole, the results presented support the role of the CYP26s in the metabolism of xenobiotic compounds as well as in potential in vivo drug interactions.

Foot-and-mouth disease (FMD) is a highly contagious and economically important disease that affects even-toed hoofed mammals. The FMD virus (FMDV) is the causative agent of FMD, of which there are seven clinically indistinguishable serotypes. Three serotypes, namely, South African Territories (SAT)1, SAT2 and SAT3 are endemic to southern Africa and are the most antigenically diverse among the FMDV serotypes. A negative consequence of this antigenic variation is that infection or vaccination with one virus may not provide immune protection from other strains or it may only confer partial protection. The identification of B-cell epitopes is therefore key to rationally designing cross-reactive vaccines that recognize the immunologically distinct serotypes present within the population. Computational epitope prediction methods that exploit the inherent physicochemical properties of epitopes in their algorithms have been proposed as a cost and time-effective alternative to the classical experimental methods. The aim of this project is to employ in silico epitope prediction programmes to predict B-cell epitopes on the capsids of the SAT serotypes. Sequence data for 18 immunologically distinct SAT1, SAT2 and SAT3 strains from across southern Africa were collated. Since, only one SAT1 virus has had its structure elucidated by X-ray crystallography (PDB ID: 2WZR), homology models of the 18 virus capsids were built computationally using Modeller v9.12. They were then subjected to energy minimizations using the AMBER force field. The quality of the models was evaluated and validated stereochemically and energetically using the PROMOTIF and ANOLEA servers respectively. The homology models were subsequently used as input to two different epitope prediction servers, namely Discotope1.0 and Ellipro. Only those epitopes predicted by both programmes were defined as epitopes. Both previously characterised and novel epitopes were predicted on the SAT strains. Some of the novel epitopes are located on the same loops as experimentally derived epitopes, while others are located on a putative novel antigenic site, which is located close to the five-fold axis of symmetry. A consensus set of 11 epitopes that are common on at least 15 out of 18 SAT strains was collated. In future work, the epitopes predicted in this study will be experimentally validated using mutagenesis studies. Those found to be true epitopes may be used in the rational design of broadly reactive SAT vaccines / Life and Consumer Sciences / M. Sc. (Life Sciences)