Computational Methods in Biophysics and Medicinal Chemistry: Applications and Challenges

In this thesis I described the theory and application of several computational methods in solving medicinal chemistry and biophysical tasks. I pointed out to the valuable information which could be achieved by means of computer simulations and to the possibility to predict the outcome of traditional experiments. Nowadays, computer represents an invaluable tool for chemists. In particular, the main topics of my research consisted in the development of an automated docking protocol for the voltage-gated hERG potassium channel blockers, and the investigation of the catalytic mechanism of the human peptidyl-prolyl cis-trans isomerase Pin1.

Molekulare und supramolekulare Eigenschaften von neuen Donor-Akzeptor-Polymeren

Die organische Halbleitertechnologie befindet sich seit Jahrzehnten im Focus des Interesses, da sie eine kostengünstige und umweltverträgliche Alternative zu anorganischen Silizium-basierten Halbleitern darstellt. Die Möglichkeit der gezielten Funktionalisierung von definierten Strukturen durch synthetische Methoden, welche eine große Vielfalt an Materialien ermöglicht, steht dabei besonders im Vordergrund. Die Modifikation von physikalischen Eigenschaften ermöglicht dabei eine stark erleichterte Anpassung für den geplanten Anwendungsbereich. Im Rahmen der vorliegenden Dissertation wurden organische Halbleitermaterialien basierend auf Cyclopenta[2,1-b:3,4-b']dithiophen (CDT) dargestellt und hinsichtlich ihrer strukturellen und elektronischen Eigenschaften untersucht. In Kombination mit Benzo[c][1,2,5]thiadiazol (BTZ) und weiteren Akzeptoren wurden zunächst Donor(D)-Akzeptor(A)-Polymere synthetisiert und Struktur-Eigenschaft-Beziehungen aufgestellt. So konnte ein sehr hochmolekulares Polymer CDT-BTZ-Polymer (Mn = 36 kg mol-1, PDI = 2.6) erhalten werden, welches sich durch eine hohe lamellare Ordnung und eine gemessene Ladungsträgermobilität in FETs von über 5.0 cm2V-1s-1 bei Raumtemperatur auszeichnete; bei niedrigen Temperaturen (240 K) war letztgenannte 6.5 cm2V-1s-1. Aufgrund dieses hohen Ladungstransports und der Abwesenheit niedermolekularer Polymerketten innerhalb des Polymers konnte erstmals eine Messung eines HALL-Effektes bewerkstelligt werden. Dies war der erste Beweis eines Band-artigen Ladungstransportes an einem Polymerhalbleiter. Des Weiteren wurde durch synthetische Veränderung der Grundstruktur des Polymers zu längeren Alkylketten eine anisotrope Anordnung der Polymerketten erreicht und die Ladungsträgermobilität (6.5 cm2V-1s-1 bei Raumtemperatur) weiter gesteigert. Darauf aufbauend wurde der Einfluss von stereoisomeren Seitenketten an CDT-BTZ-Polymeren auf Packungsverhalten, Parametern (Sperrstrom, Einschaltstrom) in FETs und Löslichkeit in organischen Lösungsmitteln untersucht. Durch cis-trans-Isomerisierung der Seitenketten wurde hier eine neue Methode zur Optimierung des Packungsverhaltens von Polymeren in dünnen Filmen und Lösung gefunden. Zuletzt wurden D-π-A-Farbstoffen, welche CDT als Verbrückungseinheit (π) beinhalten, dargestellt. Durch Variation von D und A konnten Struktur-Eigenschaft-Beziehungen in der Anwendung in Solarzellen (Feststoffsolarzellen, Flüssigsolarzellen) gefunden werden. Die Untersuchungen der photoinduzierten Absorption und der Photolumisenzenzquantenausbeute lieferten dabei Erklärungen für physikalische Prozesse wie Ladungsinjektion- und rekombination.

Conformational analysis of a potent SSTR3-selective somatostatin analogue by NMR in water solution

The three-dimensional structure of a potent SSTR3-selective analogue of somatostatin, cyclo(3-14)H-Cys(3)-Phe(6)-Tyr(7)-D-Agl(8)(N(beta) Me, 2-naphthoyl)-Lys(9)-Thr(10)-Phe(11)-Cys(14)-OH (des-AA(1, 2, 4, 5, 12, 13)[Tyr(7), D-Agl(8)(N(beta) Me, 2-naphthoyl)]-SRIF) (peptide 1) has been determined by (1)H NMR in water and molecular dynamics (MD) simulations. The peptide exists in two conformational isomers differing mainly by the cis/trans isomerization of the side chain in residue 8. The structure of 1 is compared with the consensus structural motifs of other somatostatin analogues that bind predominantly to SSTR1, SSTR2/SSTR5 and SSTR4 receptors, and to the 3D structure of a non-selective SRIF analogue, cyclo(3-14)H-Cys(3)-Phe(6)-Tyr(7)-D-2Nal(8)-Lys(9)-Thr(10)-Phe(11)-Cys(14)-OH (des-AA(1, 2, 4, 5, 12, 13)[Tyr(7), D-2Nal(8)]-SRIF) (peptide 2). The structural determinant factors that could explain selectivity of peptide 1 for SSTR3 receptors are discussed.

The Aeromonas salmonicida subsp. salmonicida exoproteome: global analysis, moonlighting proteins and putative antigens for vaccination against furunculosis

BACKGROUND Aeromonas salmonicida subsp. salmonicida, the etiologic agent of furunculosis, is a major pathogen of fisheries worldwide. Despite the identification of several virulence factors the pathogenesis is still poorly understood. We have used high-throughput proteomics to display the differences between in vitro secretome of A. salmonicida wild-type (wt, hypervirulent, JF5054) and T3SS-deficient (isogenic ΔascV, extremely low-virulent, JF2747) strains in exponential (GP) and stationary (SP) phases of growth. RESULTS Among the different experimental conditions we obtained semi-quantitative values for a total of 2136 A. salmonicida proteins. Proteins of specific A. salmonicida species were proportionally less detected than proteins common to the Aeromonas genus or those shared with other Aeromonas species, suggesting that in vitro growth did not induce the expression of these genes. Four detected proteins which are unidentified in the genome of reference strains of A. salmonicida were homologous to components of the conjugative T4SS of A. hydrophila pRA1 plasmid. Polypeptides of three proteins which are specific to the 01-B526 strain were also discovered. In supernatants (SNs), the number of detected proteins was higher in SP (326 for wt vs 329 for mutant) than in GP (275 for wt vs 263 for mutant). In pellets, the number of identified proteins (a total of 1536) was approximately the same between GP and SP. Numerous highly conserved cytoplasmic proteins were present in A. salmonicida SNs (mainly EF-Tu, EF-G, EF-P, EF-Ts, TypA, AlaS, ribosomal proteins, HtpG, DnaK, peptidyl-prolyl cis-trans isomerases, GAPDH, Enolase, FbaA, TpiA, Pgk, TktA, AckA, AcnB, Mdh, AhpC, Tpx, SodB and PNPase), and several evidences support the theory that their extracellular localization was not the result of cell lysis. According to the Cluster of Orthologous Groups classification, 29% of excreted proteins in A. salmonicida SNs were currently poorly characterized. CONCLUSIONS In this part of our work we elucidated the whole in vitro exoproteome of hypervirulent A. salmonicida subsp. salmonicida and showed the secretion of several highly conserved cytoplasmic proteins with putative moonlighting functions and roles in virulence. All together, our results offer new information about the pathogenesis of furunculosis and point out potential candidates for vaccine development.

Prodigium & Elogium Perfidiae ac Ignaviae Strasburgensis : Olim Civitatis Imperialis; Nunc Municipii Gallici

A Rurali Solitario, C. F. a K. In Campis Silesiae Habitante, Celeri Calamo concinnatum; Nec non Perpetuae Recordationi, Orbi Germanico, Germanaeq[ue] Posteritati Publicatum

C1-C8 hydrocarbons in DSDP Leg 64 Holes sediments

Sediments from the Baja California Continental Margin Transect - Sites 474 and 476 - showed small amounts of C2-C8 hydrocarbons and functionalized compounds (alkenes) typical of organic-rich, Recent, cold (<30°C) marine sediments. In contrast, some samples from Sites 477, 478, 479, and Hole 481A in the Guaymas Basin, an active spreading center, showed the characteristics of thermally generated hydrocarbons. These include an increase (sometimes exponential) in amount and diversity of C2-C8 hydrocarbons and a decrease in alkenes in more thermally mature sediments. The results indicate that the injection of basaltic sills has minimal effect on C2-C8 hydrocarbon generation except in the immediate vicinity of the sill. The absence of light hydrocarbons close to the hottest sills suggests that the compounds distill away as they are formed in these areas of very active hydrothermal circulation. A sample of young sediment exposed to very high temperatures (>300°C) from deeper thermal sources at the hottest site, 477, showed a very limited hydrocarbon distribution, including primarily ethane, benzene, and toluene, together with smaller amounts of propane and butane.

C1-C7 volatile organic compounds in sediments at DSDP Legs 56-57 Holes

Volatile C1-C7 components in sediments were examined for Japan Trench DSDP Sites 438, 439, 435, 440, 434 and 436, proceeding from west to east. Levels of all components are lowest in the highly fractured sediments of Sites 440 and 434. A number of alkenes, furans, and sulfur compounds were detected in concentrations higher than noted in any other DSDP sediments examined to date. The types, amounts, and specificity of occurrence are similar to those for 1-meter gravity cores we have examined which bear a significant biological imprint. Site 436 shows high levels of saturated and aromatic hydrocarbons, as well as olefins, including traces of dimethycyclopentanes and the highest level of cyclohexene detected in any DSDP sediment we have examined to date. The results from Site 436 were unexpected, considering the low organic-carbon content, absence of biogenic methane, and evidence of an aerobic depositional environment at this site.

Shedding light on the dark and weakly fluorescent states of green fluorescent proteins

Recent experiments on various similar green fluorescent protein (GFP) mutants at the single-molecule level and in solution provide evidence of previously unknown short- and long-lived “dark” states and of related excited-state decay channels. Here, we present quantum chemical calculations on cis-trans photoisomerization paths of neutral, anionic, and zwitterionic GFP chromophores in their ground and first singlet excited states that explain the observed behaviors from a common perspective. The results suggest that favorable radiationless decay channels can exist for the different protonation states along these isomerizations, which apparently proceed via conical intersections. These channels are suggested to rationalize the observed dramatic reduction of fluorescence in solution. The observed single-molecule fast blinking is attributed to conversions between the fluorescent anionic and the dark zwitterionic forms whereas slow switching is attributed to conversions between the anionic and the neutral forms. The predicted nonadiabatic crossings are seen to rationalize the origins of a variety of experimental observations on a common basis and may have broad implications for photobiophysical mechanisms in GFP.

Juglone, an inhibitor of the peptidyl-prolyl isomerase Pin1, also directly blocks transcription

The C-terminal domain (CTD) of the large subunit of RNA polymerase II plays a role in transcription and RNA processing. Yeast ESS1, a peptidyl-prolyl cis/trans isomerase, is involved in RNA processing and can associate with the CTD. Using several types of assays we could not find any evidence of an effect of Pin1, the human homolog of ESS1, on transcription by RNA polymerase II in vitro or on the expression of a reporter gene in vivo. However, an inhibitor of Pin1, 5-hydroxy-1,4-naphthoquinone (juglone), blocked transcription by RNA polymerase II. Unlike N-ethylmaleimide, which inhibited all phases of transcription by RNA polymerase II, juglone disrupted the formation of functional preinitiation complexes by modifying sulfhydryl groups but did not have any significant effect on either initiation or elongation. Both RNA polymerases I and III, but not T7 RNA polymerase, were inhibited by juglone. The primary target of juglone has not been unambiguously identified, although a site on the polymerase itself is suggested by inhibition of RNA polymerase II during factor-independent transcription of single-stranded DNA. Because of its unique inhibitory properties juglone should prove useful in studying transcription in vitro.

Ligand-independent assembly of recombinant human CD1 by using oxidative refolding chromatography

CD1 is an MHC class I-like antigen-presenting molecule consisting of a heavy chain and β2-microglobulin light chain. The in vitro refolding of synthetic MHC class I molecules has always required the presence of ligand. We report here the use of a folding method using an immobilized chaperone fragment, a protein disulphide isomerase, and a peptidyl-prolyl cis-trans isomerase (oxidative refolding chromatography) for the fast and efficient assembly of ligand-free and ligand-associated CD1a and CD1b, starting with material synthesized in Escherichia coli. The results suggest that “empty” MHC class I-like molecules can assemble and remain stable at physiological temperatures in the absence of ligand. The use of oxidative refolding chromatography thus is extended to encompass complex multisubunit proteins and specifically to members of the extensive, functionally diverse and important immunoglobulin supergene family of proteins, including those for which a ligand has yet to be identified.

Escherichia coli trigger factor is a prolyl isomerase that associates with nascent polypeptide chains.

Correct folding of newly synthesized proteins is proposed to be assisted by molecular chaperones and folding catalysts. To identify cellular factors involved in the initial stages of this process we searched for proteins associated with nascent polypeptide chains. In an Escherichia coli transcription/translation system synthesizing beta-galactosidase we identified a 58-kDa protein which associated with translating ribosomes but dissociated from these ribosomes upon release of nascent beta-galactosidase. N-terminal sequencing identified it as trigger factor, previously implicated in protein secretion. Direct evidence for association of trigger factor with nascent polypeptide chains was obtained by crosslinking. In a wheat germ translation system complemented with E. coli lysates, epsilon-4-(3-trifluoromethyldiazirino)benzoic acid-lysine residues were incorporated into nascent secretory preprolactin and a nonsecretory preprolactin mutant. Trigger factor crosslinked to both types of nascent chains, provided they were ribosome bound. Trigger factor contains key residues of the substrate-binding pocket of FK506-binding protein-type peptidyl-prolyl-cis/trans-isomerases and has prolyl isomerase activity in vitro. We propose that trigger factor is a folding catalyst acting cotranslationally.

Least activation path for protein folding: investigation of staphylococcal nuclease folding by stopped-flow circular dichroism.

Is the pathway of protein folding determined by the relative stability of folding intermediates, or by the relative height of the activation barriers leading to these intermediates? This is a fundamental question for resolving the Levinthal paradox, which stated that protein folding by a random search mechanism would require a time too long to be plausible. To answer this question, we have studied the guanidinium chloride (GdmCl)-induced folding/unfolding of staphylococcal nuclease [(SNase, formerly EC 3.1.4.7; now called microbial nuclease or endonuclease, EC 3.1.31.1] by stopped-flow circular dichroism (CD) and differential scanning microcalorimetry (DSC). The data show that while the equilibrium transition is a quasi-two-state process, kinetics in the 2-ms to 500-s time range are triphasic. Data support the sequential mechanism for SNase folding: U3 <--> U2 <--> U1 <--> N0, where U1, U2, and U3 are substates of the unfolded protein and N0 is the native state. Analysis of the relative population of the U1, U2, and U3 species in 2.0 M GdmCl gives delta-G values for the U3 --> U2 reaction of +0.1 kcal/mol and for the U2 --> U1 reaction of -0.49 kcal/mol. The delta-G value for the U1 --> N0 reaction is calculated to be -4.5 kcal/mol from DSC data. The activation energy, enthalpy, and entropy for each kinetic step are also determined. These results allow us to make the following four conclusions. (i) Although the U1, U2, and U3 states are nearly isoenergetic, no random walk occurs among them during the folding. The pathway of folding is unique and sequential. In other words, the relative stability of the folding intermediates does not dictate the folding pathway. Instead, the folding is a descent toward the global free-energy minimum of the native state via the least activation path in the vast energy landscape. Barrier avoidance leads the way, and barrier height limits the rate. Thus, the Levinthal paradox is not applicable to the protein-folding problem. (ii) The main folding reaction (U1 --> N0), in which the peptide chain acquires most of its free energy (via van der Waals' contacts, hydrogen bonding, and electrostatic interactions), is a highly concerted process. These energy-acquiring events take place in a single kinetic phase. (iii) U1 appears to be a compact unfolded species; the rate of conversion of U2 to U1 depends on the viscosity of solution. (iv) All four relaxation times reported here depend on GdmCl concentrations: it is likely that none involve the cis/trans isomerization of prolines. Finally, a mechanism is presented in which formation of sheet-like chain conformations and a hydrophobic condensation event precede the main-chain folding reaction.

The Drosophila melanogaster dodo (dod) gene, conserved in humans, is functionally interchangeable with the ESS1 cell division gene of Saccharomyces cerevisiae.

We have sequenced the region of DNA adjacent to and including the flightless (fli) gene of Drosophila melanogaster and molecularly characterized four transcription units within it, which we have named tweety (twe), flightless (fli), dodo (dod), and penguin (pen). We have performed deletion and transgenic analysis to determine the consequences of the quadruple gene removal. Only the flightless gene is vital to the organism; the simultaneous absence of the other three allows the overriding majority of individuals to develop to adulthood and to fly normally. These gene deletion results are evaluated in the context of the redundancy and degeneracy inherent in many genetic networks. Our cDNA analyses and data-base searches reveal that the predicted dodo protein has homologs in other eukaryotes and that it is made up of two different domains. The first, designated WW, is involved in protein-protein interactions and is found in functionally diverse proteins including human dystrophin. The second is involved in accelerating protein folding and unfolding and is found in Escherichia coli in a new family of peptidylprolyl cis-trans isomerases (PPIases; EC 5.2.1.8). In eukaryotes, PPIases occur in the nucleus and the cytoplasm and can form stable associations with transcription factors, receptors, and kinases. Given this particular combination of domains, the dodo protein may well participate in a multisubunit complex involved in the folding and activation of signaling molecules. When we expressed the dodo gene product in Saccharomyces cerevisiae, it rescued the lethal phenotype of the ESS1 cell division gene.

Cyclophilin catalyzes protein folding in yeast mitochondria.

Cyclophilins are a family of ubiquitous proteins that are the intracellular target of the immunosuppressant drug cyclosporin A. Although cyclophilins catalyze peptidylprolyl cis-trans isomerization in vitro, it has remained open whether they also perform this function in vivo. Here we show that Cpr3p, a cyclophilin in the matrix of yeast mitochondria, accelerates the refolding of a fusion protein that was synthesized in a reticulocyte lysate and imported into the matrix of isolated yeast mitochondria. The fusion protein consisted of the matrix-targeting sequence of subunit 9 of F1F0-ATPase fused to mouse dihydrofolate reductase. Refolding of the dihydrofolate reductase moiety in the matrix was monitored by acquisition of resistance to proteinase K. The rate of refolding was reduced by a factor of 2-6 by 2.5 microM cyclosporin A. This reduced rate of folding was also observed with mitochondria lacking Cpr3p. In these mitochondria, protein folding was insensitive to cyclosporin A. The rate of protein import was not affected by cyclosporin A or by deletion of Cpr3p.

A cyclophilin from the polycentric anaerobic rumen fungus Orpinomyces sp. strain PC-2 is highly homologous to vertebrate cyclophilin B.

A cyclophilin (CyP) purified to homogeneity from the polycentric anaerobic rumen fungus Orpinomyces sp. strain PC-2 had a molecular mass of 20.5 kDa and a pI of 8.1. The protein catalyzed the isomerization of the prolyl peptide bond of N-succinyl-Ala-Ala-(cis,trans)-Pro-Phe p-nitroanilide with a kcat/Km value of 9.3 x 10(6) M-1.s-1 at 10 degrees C and pH 7.8. Cyclosporin A strongly inhibited this peptidylprolyl cis-trans isomerase activity with an IC50 of 19.6 nM. The sequence of the first 30 N-terminal amino acids of this CyP had high homology with the N-terminal sequences of other eukaryotic CyPs. By use of a DNA hybridization probe amplified by PCR with degenerate oligonucleotide primers designed based on the amino acid sequences of the N terminus of this CyP and highly conserved internal regions of other CyPs, a full-length cDNA clone was isolated. It possessed an open reading frame encoding a polypeptide of 203 amino acids with a calculated molecular weight of 21,969, containing a putative hydrophobic signal peptide sequence of 22 amino acids preceding the N terminus of the mature enzyme and a C-terminal sequence, Lys-Ala-Glu-Leu, characteristic of an endoplasmic reticulum retention signal. The Orpinomyces PC-2 CyP is a typical type B CyP. The amino acid sequence of the Orpinomyces CyP exhibits striking degrees of identity with the corresponding human (70%), bovine (69%), mouse (68%), chicken (66%), maize (61%), and yeast (54%) proteins. Phylogenetic analysis based on the CyP sequences indicated that the evolutionary origin of the Orpinomyces CyP was closely related with CyPs of animals.