Structure- and Ligand-Based Structure-Activity Relationships for a Series of Inhibitors of Aldolase

Aldolase has emerged as a promising molecular target for the treatment of human African trypanosomiasis. Over the last years, due to the increasing number of patients infected with Trypanosoma brucei, there is an urgent need for new drugs to treat this neglected disease. In the present study, two-dimensional fragment-based quantitative-structure activity relationship (QSAR) models were generated for a series of inhibitors of aldolase. Through the application of leave-one-out and leave-many-out cross-validation procedures, significant correlation coefficients were obtained (r(2) = 0.98 and q(2) = 0.77) as an indication of the statistical internal and external consistency of the models. The best model was employed to predict pK(i) values for a series of test set compounds, and the predicted values were in good agreement with the experimental results, showing the power of the model for untested compounds. Moreover, structure-based molecular modeling studies were performed to investigate the binding mode of the inhibitors in the active site of the parasitic target enzyme. The structural and QSAR results provided useful molecular information for the design of new aldolase inhibitors within this structural class.

A fragment-based approach for ligand binding affinity and selectivity for the liver X receptor beta

Selective modulation of liver X receptor beta (LXR beta) has been recognized as an important approach to prevent or reverse the atherosclerotic process. In the present work, we have developed robust conformation-independent fragment-based quantitative structure-activity and structure-selectivity relationship models for a series of quinolines and cinnolines as potent modulators of the two LXR sub-types. The generated models were then used to predict the potency of an external test set and the predicted values were in good agreement with the experimental results, indicating the potential of the models for untested compounds. The final 2D molecular recognition patterns obtained were integrated to 3D structure-based molecular modeling studies to provide useful insights into the chemical and structural determinants for increased LXR beta binding affinity and selectivity. (C) 2011 Elsevier Inc. All rights reserved.

Expression, purification and structural analysis of the Pyrococcus abyssi RNA binding protein PAB1135

Abstract Background The gene coding for the uncharacterized protein PAB1135 in the archaeon Pyrococcus abyssi is in the same operon as the ribonuclease P (RNase P) subunit Rpp30. Findings Here we report the expression, purification and structural analysis of PAB1135. We analyzed the interaction of PAB1135 with RNA and show that it binds efficiently double-stranded RNAs in a non-sequence specific manner. We also performed molecular modeling of the PAB1135 structure using the crystal structure of the protein Af2318 from Archaeoglobus fulgidus (2OGK) as the template. Conclusions Comparison of this model has lead to the identification of a region in PAB1135 that could be involved in recognizing double-stranded RNA.

Structure-based drug design studies on a series of aldolase inhibitors

Human African trypanosomiasis, also known as sleeping sickness, is a major cause of death in Africa, and for which there are no safe and effective treatments available. The enzyme aldolase from Trypanosoma brucei is an attractive, validated target for drug development. A series of alkyl‑glycolamido and alkyl-monoglycolate derivatives was studied employing a combination of drug design approaches. Three-dimensional quantitative structure-activity relationships (3D QSAR) models were generated using the comparative molecular field analysis (CoMFA). Significant results were obtained for the best QSAR model (r2 = 0.95, non-cross-validated correlation coefficient, and q2 = 0.80, cross-validated correlation coefficient), indicating its predictive ability for untested compounds. The model was then used to predict values of the dependent variables (pKi) of an external test set,the predicted values were in good agreement with the experimental results. The integration of 3D QSAR, molecular docking and molecular dynamics simulations provided further insight into the structural basis for selective inhibition of the target enzyme.

Theoretical and experimental studies on the electronic, optical, and structural properties of poly-pyrrole-2-carboxylic acid films

A theoretical approach is used here to explain experimental results obtained from the electrosynthesis of polypyrrole-2-carboxylic acid (PPY-2-COOH) films in nonaqueous medium. An analysis of the Fukui function (reactivity index) indicates that the monomer (pyrrole-2-carboxylic acid, PY-2-COOH), and dimers and trimers are oxidized in the C4 or C5 positions of the heterocyclic ring of the PY-2-COOH structure. After calculating the heat of formation using semiempirical Austin Model 1 post-Hartree-Fock parameterization for dimer species, both C4 and C5 positions adjacent to the aromatic rings of PPY-2-COOH were considered the most susceptible ones to oxidative coupling reactions. The ZINDO-S/CI semiempirical method was used to simulate the electronic transitions typically seen in the UV-VIS-NIR range in monomer and oligomers with different conjugation lengths. The use of an electrochemical quartz crystal microbalance provides sufficient information to propose a polymerization mechanism of PY-2-COOH based on molecular modeling and experimental results.

Computational analyses on the structure-function relation in ion channels

Ion channels are protein molecules, embedded in the lipid bilayer of the cell membranes. They act as powerful sensing elements switching chemicalphysical stimuli into ion-fluxes. At a glance, ion channels are water-filled pores, which can open and close in response to different stimuli (gating), and one once open select the permeating ion species (selectivity). They play a crucial role in several physiological functions, like nerve transmission, muscular contraction, and secretion. Besides, ion channels can be used in technological applications for different purpose (sensing of organic molecules, DNA sequencing). As a result, there is remarkable interest in understanding the molecular determinants of the channel functioning. Nowadays, both the functional and the structural characteristics of ion channels can be experimentally solved. The purpose of this thesis was to investigate the structure-function relation in ion channels, by computational techniques. Most of the analyses focused on the mechanisms of ion conduction, and the numerical methodologies to compute the channel conductance. The standard techniques for atomistic simulation of complex molecular systems (Molecular Dynamics) cannot be routinely used to calculate ion fluxes in membrane channels, because of the high computational resources needed. The main step forward of the PhD research activity was the development of a computational algorithm for the calculation of ion fluxes in protein channels. The algorithm - based on the electrodiffusion theory - is computational inexpensive, and was used for an extensive analysis on the molecular determinants of the channel conductance. The first record of ion-fluxes through a single protein channel dates back to 1976, and since then measuring the single channel conductance has become a standard experimental procedure. Chapter 1 introduces ion channels, and the experimental techniques used to measure the channel currents. The abundance of functional data (channel currents) does not match with an equal abundance of structural data. The bacterial potassium channel KcsA was the first selective ion channels to be experimentally solved (1998), and after KcsA the structures of four different potassium channels were revealed. These experimental data inspired a new era in ion channel modeling. Once the atomic structures of channels are known, it is possible to define mathematical models based on physical descriptions of the molecular systems. These physically based models can provide an atomic description of ion channel functioning, and predict the effect of structural changes. Chapter 2 introduces the computation methods used throughout the thesis to model ion channels functioning at the atomic level. In Chapter 3 and Chapter 4 the ion conduction through potassium channels is analyzed, by an approach based on the Poisson-Nernst-Planck electrodiffusion theory. In the electrodiffusion theory ion conduction is modeled by the drift-diffusion equations, thus describing the ion distributions by continuum functions. The numerical solver of the Poisson- Nernst-Planck equations was tested in the KcsA potassium channel (Chapter 3), and then used to analyze how the atomic structure of the intracellular vestibule of potassium channels affects the conductance (Chapter 4). As a major result, a correlation between the channel conductance and the potassium concentration in the intracellular vestibule emerged. The atomic structure of the channel modulates the potassium concentration in the vestibule, thus its conductance. This mechanism explains the phenotype of the BK potassium channels, a sub-family of potassium channels with high single channel conductance. The functional role of the intracellular vestibule is also the subject of Chapter 5, where the affinity of the potassium channels hEag1 (involved in tumour-cell proliferation) and hErg (important in the cardiac cycle) for several pharmaceutical drugs was compared. Both experimental measurements and molecular modeling were used in order to identify differences in the blocking mechanism of the two channels, which could be exploited in the synthesis of selective blockers. The experimental data pointed out the different role of residue mutations in the blockage of hEag1 and hErg, and the molecular modeling provided a possible explanation based on different binding sites in the intracellular vestibule. Modeling ion channels at the molecular levels relates the functioning of a channel to its atomic structure (Chapters 3-5), and can also be useful to predict the structure of ion channels (Chapter 6-7). In Chapter 6 the structure of the KcsA potassium channel depleted from potassium ions is analyzed by molecular dynamics simulations. Recently, a surprisingly high osmotic permeability of the KcsA channel was experimentally measured. All the available crystallographic structure of KcsA refers to a channel occupied by potassium ions. To conduct water molecules potassium ions must be expelled from KcsA. The structure of the potassium-depleted KcsA channel and the mechanism of water permeation are still unknown, and have been investigated by numerical simulations. Molecular dynamics of KcsA identified a possible atomic structure of the potassium-depleted KcsA channel, and a mechanism for water permeation. The depletion from potassium ions is an extreme situation for potassium channels, unlikely in physiological conditions. However, the simulation of such an extreme condition could help to identify the structural conformations, so the functional states, accessible to potassium ion channels. The last chapter of the thesis deals with the atomic structure of the !- Hemolysin channel. !-Hemolysin is the major determinant of the Staphylococcus Aureus toxicity, and is also the prototype channel for a possible usage in technological applications. The atomic structure of !- Hemolysin was revealed by X-Ray crystallography, but several experimental evidences suggest the presence of an alternative atomic structure. This alternative structure was predicted, combining experimental measurements of single channel currents and numerical simulations. This thesis is organized in two parts, in the first part an overview on ion channels and on the numerical methods adopted throughout the thesis is provided, while the second part describes the research projects tackled in the course of the PhD programme. The aim of the research activity was to relate the functional characteristics of ion channels to their atomic structure. In presenting the different research projects, the role of numerical simulations to analyze the structure-function relation in ion channels is highlighted.

Azobenzol- und Perylendiimid-funktionalisierte Polyphenylen-Dendrimere

Die Dissertation 'Azobenzol- und Perylendiimid-funktionalisierte Polyphenylen-Dendrimere - Synthese, Charakterisierung und Eigenschaften' gliedert sich in vier Themengebiete. Der erste Abschnitt beschäftigt sich mit der Synthese unterschiedlich dichter Dendrimere um einen Azobenzol-Kern. Einkristallstrukturen und Molekülvisualisierungen verdeutlichen die dreidimensionale Gestalt der Dendrimere. Die Dendrimere zeigen erstmalig eine Abhängigkeit des Isomerisationsverhaltens von der das Chromophor umgebenden Struktur. Der zweite Abschnitt hat Interaktionen von Chromophoren, deren Distanz und Orientierung zueinander gezielt durch einen äußeren Impuls geändert werden können, zum Thema. Die Verbindung von Azobenzol und PMI führt durch deren gegenseitige Beeinflussung zu einem Verlust der charakteristischen Eigenschaften der Chromophore. Eine Oligo-L-Lysinkette, deren Enden mit NMI und PMI funktionalisiert sind, stellt ein FRET-System dar. Distanz und Orientierung der Chromophore zueinander werden durch den mittels TFE induzierten Übergang des Peptids vom Knäuel zur Helix verändert. Der dritte Abschnitt führt die Synthese von PDI-gekernten Dendrimeren durch Substitution in der bay-Region des Chromophors ein. Die Eignenschaften der Verbindungen wurden mittels optischer Methoden und cyclovoltammetrischen Studien untersucht. Weiter wurde die Oberflächenfunktionalisierung mit Aminosäuren und Oligopeptiden zu wasserlöslichen Dendrimeren mit hoher Oberflächenladung verfolgt. Das letzte Kapitel stellt Untersuchungen zur Organisation von Polyphenylen-Dendrimeren auf HOPG vor. Es lassen sich einerseits Nanofasern formieren, andererseits können auch geordnete Mono- und Multilagen erzeugt werden.

First principle study of ions in water and ionic liquids

Die beiden in dieser Arbeit betrachteten Systeme, wässrige Lösungen von Ionen und ionische Flüssigkeiten, zeigen vielfältige Eigenschaften und Anwendungsmöglichkeiten, im Gegensatz zu anderen Systemen. Man findet sie beinahe überall im normalen Leben (Wasser), oder ihre Bedeutung wächst (ioinische Flüssigkeiten). Der elektronische Anteil und der atomare Anteil wurden getrennt voneinander untersucht und im Zusammenhang analysiert. Mittels dieser Methode konnten die in dem jeweiligen System auftretenden Mechanismen genauer untersucht werden. Diese Methode wird "Multiscale Modeling" genannt, dabei werden die Untereinheiten eines Systems genauer betrachtet, wie in diesem Fall die elektronischen and atomaren Teilsystem. Die Ergebnisse, die aus den jeweiligen Betrachtungen hervorgehen, zeigen, dass, im Falle von hydratisierten Ionen die Wasser-Wasser Wechselwirkungen wesentlich stärker sind als die elektrostatischen Wechselwirkung zwischen Wasser und dem Ion. Anhand der Ergebnisse ergibt sich, dass normale nicht-polarisierbare Modelle ausreichen, um Ionen-Wasser Lösungen zu beschreiben. Im Falle der ionischen Flüssigkeiten betrachten wir die elektronische Ebene mittels sehr genauer post-Hartree-Fock Methoden und DFT, deren Ergebnisse dann mit denen auf molekularer Ebene (mithilfe von CPMD/klassischer MD) in Beziehung gesetzt werden. Die bisherigen Ergebnisse zeigen, dass die Wasserstoff-Brückenbindungen im Fall der ionischen Flüssigkeiten nicht vernachässigt werden können. Weiterhin hat diese Studie herausgefunden, dass die klassischen Kraftfelder die Elektrostatik (Dipol- und Quadrupolmomente) nicht genau genug beschreibt. Die Kombination des mikroskopischen Mechanismus und der molekularen Eigenschaften ist besonders sinnvoll um verschiedene Anhaltspunkte von Simualtionen (z.B. mit klassische Molekular-Dynamik) oder Experimenten zu liefern oder solche zu erklären.

Zur Struktur und Funktion der Typ-3-Kupferproteine sowie der Lichtsammlerkomplexe LHCI und LHCII

In der vorliegenden Arbeit wurden die bioinformatischen Methoden der Homologie-Modellierung und Molekularen Modellierung dazu benutzt, die dreidimensionalen Strukturen verschiedenster Proteine vorherzusagen und zu analysieren. Experimentelle Befunde aus Laborversuchen wurden dazu benutzt, die Genauigkeit der Homologie-Modelle zu erhöhen. Die Ergebnisse aus den Modellierungen wurden wiederum dazu benutzt, um neue experimentelle Versuche vorzuschlagen. Anhand der erstellten Modelle und bekannten Kristallstrukturen aus der Protein-Datenbank PDB wurde die Struktur-Funktionsbeziehung verschiedener Tyrosinasen untersucht. Dazu gehörten sowohl die Tyrosinase des Bakteriums Streptomyces als auch die Tyrosinase der Hausmaus. Aus den vergleichenden Strukturanalysen der Tyrosinasen resultierten Mechanismen für die Monophenolhydroxylase-Aktivität der Tyrosinasen sowie für den Import der Kupferionen ins aktive Zentrum. Es konnte der Beweis geführt werden, daß die Blockade des CuA-Zentrums tatsächlich der Grund für die unterschiedliche Aktivität von Tyrosinasen und Catecholoxidasen ist. Zum ersten Mal konnte mit der Maus-Tyrosinase ein vollständiges Strukturmodell einer Säugetier-Tyrosinase erstellt werden, das dazu in der Lage ist, die Mechanismen bekannter Albino-Mutationen auf molekularer Ebene zu erklären. Die auf der Basis des ermittelten 3D-Modells gewonnenen Erkenntnisse über die Wichtigkeit bestimmter Aminosäuren für die Funktion wurde durch gerichtete Mutagenese an der rekombinant hergestellten Maus-Tyrosinase getestet und bestätigt. Weiterhin wurde die Struktur der Tyrosinase des Krebses Palinurus elephas durch eine niedrigaufgelöste 3D-Rekonstruktion aus elektronenmikroskopischen Bildern aufgeklärt. Der zweite große Themenkomplex umfasst die Strukturanalyse der Lichtsammlerkomplexe LHCI-730 und LHCII. Im Falle des LHCII konnte der Oligomerisierungszustand der LHCMoleküle mit diskreten Konformationen des N-Terminus korreliert werden. Auch hier kam eine Kombination von Homologie-Modellierung und einer experimentellen Methode, der Elektronen-Spin-Resonanz-Messung, zum Einsatz. Die Änderung des Oligomerisierungszustands des LHCII kontrolliert den Energiezufluß zu den Photosystemen PS I und PS II. Des Weiteren wurde ein vollständiges Modell des LHCI-730 erstellt, um die Auswirkungen gerichteter Mutagenese auf das Dimerisierungsverhalten zu untersuchen. Auf Basis dieses Modells wurden die Wechselwirkungen zwischen den Monomeren Lhca1 und Lhca4 evaluiert und potentielle Bindungspartner identifiziert.

Sensitive and selective detection of free FXIII activation peptide: a potential marker of acute thrombotic events

Coagulation factor XIII (FXIII) stabilizes fibrin fibers and is therefore a major player in the maintenance of hemostasis. FXIII is activated by thrombin resulting in cleavage and release of the FXIII activation peptide (AP-FXIII). The objective of this study was to characterize the released AP-FXIII and determine specific features that may be used for its specific detection. We analyzed the structure of bound AP-FXIII within the FXIII A-subunit and interactions of AP-FXIII by hydrogen bonds with both FXIII A-subunit monomers. We optimized our previously developed AP-FXIII ELISA by using 2 monoclonal antibodies. We determined high binding affinities between the antibodies and free AP-FXIII and demonstrated specific binding by epitope mapping analyses with surface plasmon resonance and enzyme-linked immunosorbent assay. Because the structure of free AP-FXIII had been characterized so far by molecular modeling only, we performed structural analysis by nuclear magnetic resonance. Recombinant AP-FXIII was largely flexible both in plasma and water, differing significantly from the rigid structure in the bound state. We suggest that the recognized epitope is either occluded in the noncleaved form or possesses a structure that does not allow binding to the antibodies. On the basis of our findings, we propose AP-FXIII as a possible new marker for acute thrombotic events.

Clinical and biochemical description of a novel CYP21A2 gene mutation 962_963insA using a new 3D model for the P450c21 protein

OBJECTIVE: A severely virilized 46, XX newborn girl was referred to our center for evaluation and treatment of congenital adrenal hyperplasia (CAH) because of highly elevated 17alpha-hydroxyprogesterone levels at newborn screening; biochemical tests confirmed the diagnosis of salt-wasting CAH. Genetic analysis revealed that the girl was compound heterozygote for a previously reported Q318X mutation in exon 8 and a novel insertion of an adenine between nucleotides 962 and 963 in exon 4 of the CYP21A2 gene. This 962_963insA mutation created a frameshift leading to a stop codon at amino acid 161 of the P450c21 protein. AIM AND METHODS: To better understand structure-function relationships of mutant P450c21 proteins, we performed multiple sequence alignments of P450c21 with three mammalian P450s (P450 2C8, 2C9 and 2B4) with known structures as well as with human P450c17. Comparative molecular modeling of human P450c21 was then performed by MODELLER using the X-ray crystal structure of rabbit P450 2B4 as a template. RESULTS: The new three dimensional model of human P450c21 and the sequence alignment were found to be helpful in predicting the role of various amino acids in P450c21, especially those involved in heme binding and interaction with P450 oxidoreductase, the obligate electron donor. CONCLUSION: Our model will help in analyzing the genotype-phenotype relationship of P450c21 mutations which have not been tested for their functional activity in an in vitro assay.

Modulation of human CYP19A1 activity by mutant NADPH P450 oxidoreductase

Mutations in NADPH P450 oxidoreductase (POR) cause a broad spectrum of human disease with abnormalities in steroidogenesis. We have studied the impact of P450 reductase mutations on the activity of CYP19A1. POR supported CYP19A1 activity with a calculated Km of 126 nm for androstenedione and a Vmax of 1.7 pmol/min. Mutations R457H and V492E located in the FAD domain of POR that disrupt electron transfer caused a complete loss of CYP19A1 activity. The A287P mutation of POR decreased the activities of CYP17A1 by 60-80% but had normal CYP19A1 activity. Molecular modeling and protein docking studies suggested that A287P is involved in the interaction of POR:CYP17A1 but not in the POR:CYP19A1 interaction. Mutations C569Y and V608F in the NADPH binding domain of POR had 49 and 28% of activity of CYP19A1 compared with normal reductase and were more sensitive to the amount of NADPH available for supporting CYP19A1 activity. Substitution of NADH for NADPH had a higher impact on C569Y and V608F mutants of POR. Similar effects were obtained at low/high (5.5/8.5) pH, but using octanol to limit the flux of electrons from POR to CYP19A1 inhibited activity supported by all variants. High molar ratios of KCl also reduced the CYP19A1 supporting activities of C569Y and V608F mutants of POR to a greater extent compared to normal POR and A287P mutant. Because POR supports many P450s involved in steroidogenesis, bone formation, and drug metabolism, variations in the effects of POR mutations on specific enzyme activities may explain the broad clinical spectrum of POR deficiency.

The crystal structure of the platelet activator aggretin reveals a novel (alphabeta)2 dimeric structure

Aggretin is a C-type lectin purified from Calloselasma rhodostoma snake venom. It is a potent activator of platelets, resulting in a collagen-like response by binding and clustering platelet receptor CLEC-2. We present here the crystal structure of aggretin at 1.7 A which reveals a unique tetrameric quaternary structure. The two alphabeta heterodimers are arranged through 2-fold rotational symmetry, resulting in an antiparallel side-by-side arrangement. Aggretin thus presents two ligand binding sites on one surface and can therefore cluster ligands in a manner reminiscent of convulxin and flavocetin. To examine the molecular basis of the interaction with CLEC-2, we used a molecular modeling approach of docking the aggretin alphabeta structure with the CLEC-2 N-terminal domain (CLEC-2N). This model positions the CLEC-2N structure face down in the "saddle"-shaped binding site which lies between the aggretin alpha and beta lectin-like domains. A 2-fold rotation of this complex to generate the aggretin tetramer reveals dimer contacts for CLEC-2N which bring the N- and C-termini into the proximity of each other, and a series of contacts involving two interlocking beta-strands close to the N-terminus are described. A comparison with homologous lectin-like domains from the immunoreceptor family reveals a similar but not identical dimerization mode, suggesting this structure may represent the clustered form of CLEC-2 capable of signaling across the platelet membrane.

Alkylamides from Echinacea are a new class of cannabinomimetics. Cannabinoid type 2 receptor-dependent and -independent immunomodulatory effects

Alkylamides (alkamides) from Echinacea modulate tumor necrosis factor alpha mRNA expression in human monocytes/macrophages via the cannabinoid type 2 (CB2) receptor (Gertsch, J., Schoop, R., Kuenzle, U., and Suter, A. (2004) FEBS Lett. 577, 563-569). Here we show that the alkylamides dodeca-2E,4E,8Z,10Z-tetraenoic acid isobutylamide (A1) and dodeca-2E,4E-dienoic acid isobutylamide (A2) bind to the CB2 receptor more strongly than the endogenous cannabinoids. The Ki values of A1 and A2 (CB2 approximately 60 nM; CB1 >1500 nM) were determined by displacement of the synthetic high affinity cannabinoid ligand [3H]CP-55,940. Molecular modeling suggests that alkylamides bind in the solvent-accessible cavity in CB2, directed by H-bonding and pi-pi interactions. In a screen with 49 other pharmacologically relevant receptors, it could be shown that A1 and A2 specifically bind to CB2 and CB1. A1 and A2 elevated total intracellular Ca2+ in CB2-positive but not in CB2-negative promyelocytic HL60 cells, an effect that was inhibited by the CB2 antagonist SR144528. At 50 nM, A1, A2, and the endogenous cannabinoid anandamide (CB2 Ki >200 nM) up-regulated constitutive interleukin (IL)-6 expression in human whole blood in a seemingly CB2-dependent manner. A1, A2, anandamide, the CB2 antagonist SR144528 (Ki <10 nM), and also the non-CB2-binding alkylamide undeca-2E-ene,8,10-diynoic acid isobutylamide all significantly inhibited lipopolysaccharide-induced tumor necrosis factor alpha, IL-1beta, and IL-12p70 expression (5-500 nM) in a CB2-independent manner. Alkylamides and anandamide also showed weak differential effects on anti-CD3-versus anti-CD28-stimulated cytokine expression in human whole blood. Overall, alkylamides, anandamide, and SR144528 potently inhibited lipopolysaccharide-induced inflammation in human whole blood and exerted modulatory effects on cytokine expression, but these effects are not exclusively related to CB2 binding.

Structural model of the CopA copper ATPase of Enterococcus hirae based on chemical cross-linking

The CopA copper ATPase of Enterococcus hirae belongs to the family of heavy metal pumping CPx-type ATPases and shares 43% sequence similarity with the human Menkes and Wilson copper ATPases. Due to a lack of suitable protein crystals, only partial three-dimensional structures have so far been obtained for this family of ion pumps. We present a structural model of CopA derived by combining topological information obtained by intramolecular cross-linking with molecular modeling. Purified CopA was cross-linked with different bivalent reagents, followed by tryptic digestion and identification of cross-linked peptides by mass spectrometry. The structural proximity of tryptic fragments provided information about the structural arrangement of the hydrophilic protein domains, which was integrated into a three-dimensional model of CopA. Comparative modeling of CopA was guided by the sequence similarity to the calcium ATPase of the sarcoplasmic reticulum, Serca1, for which detailed structures are available. In addition, known partial structures of CPx-ATPase homologous to CopA were used as modeling templates. A docking approach was used to predict the orientation of the heavy metal binding domain of CopA relative to the core structure, which was verified by distance constraints derived from cross-links. The overall structural model of CopA resembles the Serca1 structure, but reveals distinctive features of CPx-type ATPases. A prominent feature is the positioning of the heavy metal binding domain. It features an orientation of the Cu binding ligands which is appropriate for the interaction with Cu-loaded metallochaperones in solution. Moreover, a novel model of the architecture of the intramembranous Cu binding sites could be derived.