Crystal Structures of an Oligopeptide-Binding Protein from the Biosynthetic Pathway of the β-Lactamase Inhibitor Clavulanic Acid

Clavulanic acid (CA) is a clinically important β-lactamase inhibitor that is produced by fermentation of Streptomyces clavuligerus. The CA biosynthesis pathway starts from arginine and glyceraldehyde-3-phosphate and proceeds via (3S,5S)-clavaminic acid, which is converted to (3R,5R)-clavaldehyde, the immediate precursor of (3R,5R)-CA. Open reading frames 7 (orf7) and 15 (orf15) of the CA biosynthesis cluster encode oligopeptide-binding proteins (OppA1 and OppA2), which are essential for CA biosynthesis. OppA1/2 are proposed to be involved in the binding and/or transport of peptides across the S. clavuligerus cell membrane. Peptide binding assays reveal that recombinant OppA1 and OppA2 bind di-/tripeptides containing arginine and certain nonapeptides including bradykinin. Crystal structures of OppA2 in its apo form and in complex with arginine or bradykinin were solved to 1.45, 1.7, and 1.7 Å resolution, respectively. The overall fold of OppA2 consists of two lobes with a deep cavity in the center, as observed for other oligopeptide-binding proteins. The large cavity creates a peptide/arginine binding cleft. The crystal structures of OppA2 in complex with arginine or bradykinin reveal that the C-terminal arginine of bradykinin binds similarly to arginine. The results are discussed in terms of the possible roles of OppA1/2 in CA biosynthesis.     

Single-Stranded DNA-Binding Protein Complex from Helicobacter pylori Suggests an ssDNA-Binding Surface

Single-stranded DNA (ssDNA)-binding protein (SSB) plays an important role in DNA replication, recombination, and repair. SSB consists of an N-terminal ssDNA-binding domain with an oligonucleotide/oligosaccharide binding fold and a flexible C-terminal tail involved in protein-protein interactions. SSB from Helicobacter pylori (HpSSB) was isolated, and the ssDNA-binding characteristics of HpSSB were analyzed by fluorescence titration and electrophoretic mobility shift assay. Tryptophan fluorescence quenching was measured as 61%, and the calculated cooperative affinity was 5.4 × 10⁷ M^− 1 with an ssDNA-binding length of 25-30 nt. The crystal structure of the C-terminally truncated protein (HpSSBc) in complex with 35-mer ssDNA [HpSSBc-(dT)₃₅] was determined at a resolution of 2.3 Å. The HpSSBc monomer folds as an oligonucleotide/oligosaccharide binding fold with a Y-shaped conformation. The ssDNA wrapped around the HpSSBc tetramer through a continuous binding path comprising five essential aromatic residues and a positively charged surface formed by numerous basic residues.     

The crystal structure of the substrate-binding protein OpuBC from Bacillus subtilis in complex with choline

Bacillus subtilis can synthesize the compatible solute glycine betaine as an osmoprotectant from an exogenous supply of the precursor choline. Import of choline is mediated by two osmotically inducible ABC transport systems: OpuB and OpuC. OpuC catalyzes the import of various osmoprotectants, whereas OpuB is highly specific for choline. OpuBC is the substrate-binding protein of the OpuB transporter, and we have analyzed the affinity of the OpuBC/choline complex by intrinsic tryptophan fluorescence and determined a K_d value of about 30 μM. The X-ray crystal structure of the OpuBC/choline complex was solved at a resolution of 1.6 Å and revealed a fold typical of class II substrate-binding proteins. The positively charged trimethylammonium head group of choline is wedged into an aromatic cage formed by four tyrosine residues and is bound via cation-pi interactions. The hydroxyl group of choline protrudes out of this aromatic cage and makes a single interaction with residue Gln19. The substitution of this residue by Ala decreases choline binding affinity by approximately 15-fold. A water network stabilizes choline within its substrate-binding site and promotes indirect interactions between the two lobes of the OpuBC protein. Disruption of this intricate water network by site-directed mutagenesis of selected residues in OpuBC either strongly reduces choline binding affinity (between 18-fold and 25-fold) or abrogates ligand binding. The crystal structure of the OpuBC/choline complex provides a rational for the observed choline specificity of the OpuB ABC importer in vivo and explains its inability to catalyze the import of glycine betaine into osmotically stressed B. subtilis cells.     

Thermodynamics of the interactions of tryptophan-rich cathelicidin antimicrobial peptides with model and natural membranes

Tritrpticin and indolicidin are short 13-residue tryptophan-rich antimicrobial peptides that hold potential as future alternatives for antibiotics. Isothermal titration calorimetry (ITC) has been applied as the main tool in this study to investigate the thermodynamics of the interaction of these two cathelicidin peptides as well as five tritrpticin analogs with large unilamellar vesicles (LUVs), representing model and natural anionic membranes. The anionic LUVs were composed of (a) 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPE/POPG) (7:3) and (b) natural E. coli polar lipid extract. 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) was used to make model zwitterionic membranes. Binding isotherms were obtained to characterize the antimicrobial peptide binding to the LUVs, which then allowed for calculation of the thermodynamic parameters of the interaction. All peptides exhibited substantially stronger binding to anionic POPE/POPG and E. coli membrane systems than to the zwitterionic POPC system due to strong electrostatic attractions between the highly positively charged peptides and the negatively charged membrane surface, and results with tritrpticin derivatives further revealed the effects of various amino acid substitutions on membrane binding. No significant improvement was observed upon increasing the Tritrp peptide charge from + 4 to + 5. Replacement of Arg residues with Lys did not substantially change peptide binding to anionic vesicles but moderately decreased the binding to zwitterionic LUVs. Pro to Ala substitutions in tritrpticin, allowing the peptide to adopt an α-helical structure, resulted in a significant increase of the binding to both anionic and zwitterionic vesicles and therefore reduced the selectivity for bacterial and mammalian membranes. In contrast, substitution of Trp with other aromatic amino acids significantly decreased the peptide's ability to bind to anionic LUVs and essentially eliminated binding to zwitterionic LUVs. The ITC results were consistent with the outcome of fluorescence spectroscopy membrane binding and perturbation studies. Overall, our work showed that a natural E. coli polar lipid extract as a bacterial membrane model was advantageous compared to the simpler and more widely used POPE/POPG lipid system.     

Functional testing of putative oligopeptide permease (Opp) proteins of Borrelia burgdorferi: a complementation model in opp− Escherichia coli

Studies of the protein function of Borrelia burgdorferi have been limited by a lack of tools for manipulating borrelial DNA. We devised a system to study the function of a B. burgdorferi oligopeptide permease (Opp) orthologue by complementation with Escherichia coli Opp proteins. The Opp system of E. coli has been extensively studied and has well defined substrate specificities. The system is of interest in B. burgdorferi because analysis of its genome has revealed little identifiable machinery for synthesis or transport of amino acids and only a single intact peptide transporter operon. As such, peptide uptake may play a major role in nutrition for the organism. Substrate specificity for ABC peptide transporters in other organisms is determined by their substrate binding protein. The B. burgdorferi Opp operon differs from the E. coli Opp operon in that it has three separate substrate binding proteins, OppA-1, -2 and -3. In addition, B. burgdorferi has two OppA orthologues, OppA-4 and -5, encoded on separate plasmids. The substrate binding proteins interact with integral membrane proteins, OppB and OppC, to transport peptides into the cell. The process is driven by two ATP binding proteins, OppD and OppF. Using opp-deleted E. coli mutants, we transformed cells with B. burgdorferi oppA-1, -2, -4 or -5 and E. coli oppBCDF. All of the B. burgdorferi OppA proteins are able to complement E. coli OppBCDF to form a functional Opp transport system capable of transporting peptides for nutritional use. Although there is overlap in substrate specificities, the substrate specificities for B. burgdorferi OppAs are not identical to that of E. coli OppA. Transport of toxic peptides by B. burgdorferi grown in nutrient-rich medium parallels borrelial OppA substrate specificity in the complementation system. Use of this complementation system will pave the way for more detailed studies of B. burgdorferi peptide transport than currently available tools for manipulating borrelial DNA will allow.     

Structural and thermodynamic signatures of DNA recognition by Mycobacterium tuberculosis DnaA

An essential protein, DnaA, binds to 9-bp DNA sites within the origin of replication oriC. These binding events are prerequisite to forming an enigmatic nucleoprotein scaffold that initiates replication. The number, sequences, positions, and orientations of these short DNA sites, or DnaA boxes, within the oriCs of different bacteria vary considerably. To investigate features of DnaA boxes that are important for binding Mycobacterium tuberculosis DnaA (MtDnaA), we have determined the crystal structures of the DNA binding domain (DBD) of MtDnaA bound to a cognate MtDnaA-box (at 2.0 Å resolution) and to a consensus Escherichia coli DnaA-box (at 2.3 Å). These structures, complemented by calorimetric equilibrium binding studies of MtDnaA DBD in a series of DnaA-box variants, reveal the main determinants of DNA recognition and establish the [T/C][T/A][G/A]TCCACA sequence as a high-affinity MtDnaA-box. Bioinformatic and calorimetric analyses indicate that DnaA-box sequences in mycobacterial oriCs generally differ from the optimal binding sequence. This sequence variation occurs commonly at the first 2 bp, making an in vivo mycobacterial DnaA-box effectively a 7-mer and not a 9-mer. We demonstrate that the decrease in the affinity of these MtDnaA-box variants for MtDnaA DBD relative to that of the highest-affinity box TTGTCCACA is less than 10-fold. The understanding of DnaA-box recognition by MtDnaA and E. coli DnaA enables one to map DnaA-box sequences in the genomes of M. tuberculosis and other eubacteria.     

Relating structure to thermodynamics: the crystal structures and binding affinity of eight OppA-peptide complexes.

The oligopeptide-binding protein OppA provides a useful model system for studying the physical chemistry underlying noncovalent interactions since it binds a variety of readily synthesized ligands. We have studied the binding of eight closely related tripeptides of the type Lysine-X-Lysine, where X is an abnormal amino acid, by isothermal titration calorimetry (ITC) and X-ray crystallography. The tripeptides fall into three series of ligands, which have been designed to examine the effects of small changes to the central side chain. Three ligands have a primary amine as the second side chain, two have a straight alkane chain, and three have ring systems. The results have revealed a definite preference for the binding of hydrophobic residues over the positively charged side chains, the latter binding only weakly due to unfavorable enthalpic effects. Within the series of positively charged groups, a point of lowest affinity has been identified and this is proposed to arise from unfavorable electrostatic interactions in the pocket, including the disruption of a key salt bridge. Marked entropy-enthalpy compensation is found across the series, and some of the difficulties in designing tightly binding ligands have been highlighted.     

Crystal Structure of Methanocaldococcus jannaschii Trm4 Complexed with Sinefungin

tRNA:m⁵C methyltransferase Trm4 generates the modified nucleotide 5-methylcytidine in archaeal and eukaryotic tRNA molecules, using S-adenosyl-l-methionine (AdoMet) as methyl donor. Most archaea and eukaryotes possess several Trm4 homologs, including those related to diseases, while the archaeon Methanocaldococcus jannaschii has only one gene encoding a Trm4 homolog, MJ0026. The recombinant MJ0026 protein catalyzed AdoMet-dependent methyltransferase activity on tRNA in vitro and was shown to be the M. jannaschii Trm4. We determined the crystal structures of the substrate-free M. jannaschii Trm4 and its complex with sinefungin at 1.27 Å and 2.3 Å resolutions, respectively. This AdoMet analog is bound in a negatively charged pocket near helix α8. This helix can adopt two different conformations, thereby controlling the entry of AdoMet into the active site. Adjacent to the sinefungin-bound pocket, highly conserved residues form a large, positively charged surface, which seems to be suitable for tRNA binding. The structure explains the roles of several conserved residues that were reportedly involved in the enzymatic activity or stability of Trm4p from the yeast Saccharomyces cerevisiae. We also discuss previous genetic and biochemical data on human NSUN2/hTrm4/Misu and archaeal PAB1947 methyltransferase, based on the structure of M. jannaschii Trm4.     

Crystal Structure and Functional Analysis of the Protein Disulfide Isomerase-Related Protein ERp29

The protein disulfide isomerase-related protein ERp29 is a putative chaperone involved in processing and secretion of secretory proteins. Until now, however, both the structure and the exact nature of interacting substrates remained unclear. We provide for the first time a crystal structure of human ERp29, refined to 2.9 Å, and show that the protein has considerable structural homology to its Drosophila homolog Wind. We show that ERp29 binds directly not only to thyroglobulin and thyroglobulin-derived peptides in vitro but also to the Wind client protein Pipe and Pipe-derived peptides, although it fails to process Pipe in vivo. A monomeric mutant of ERp29 and a D domain mutant in which the second peptide binding site is inactivated also bind protein substrates, indicating that the monomeric thioredoxin domain is sufficient for client protein binding. Indeed, the b domains of ERp29 or Wind, expressed alone, are sufficient for binding proteins and peptides. Interacting peptides have in common two or more aromatic residues, with stronger binding for sequences with overall basic character. Thus, the data allow a view of the two putative peptide binding sites of ERp29 and indicate that the apparent, different processing activity of the human and Drosophila proteins in vivo does not stem from differences in peptide binding properties.     

Interactions of the Escherichia coli Primosomal PriB Protein with the Single-stranded DNA. Stoichiometries, Intrinsic Affinities, Cooperativities, and Base Specificities

Quantitative analysis of the interactions of the Escherichia coli primosomal PriB protein with a single-stranded DNA was done using quantitative fluorescence titration, photocrosslinking, and analytical ultracentrifugation techniques. Stoichiometry studies were done with a series of etheno-derivatives of single-stranded (ss) DNA oligomers. Interactions with the unmodified nucleic acids were studied, using the macromolecular competition titration (MCT) method. The total site-size of the PriB dimer-ssDNA complex, i.e. the maximum number of nucleotides occluded by the PriB dimer in the complex, is 12 ± 1 nt. The protein has a single DNA-binding site, which is located centrally within the dimer and has a functionally homogeneous structure. The stoichiometry and photocrosslinking data show that only a single monomer of the PriB dimer engages in interactions with the nucleic acid. The analysis of the PriB binding to long oligomers was done using a statistical thermodynamic model that takes into account the overlap of potential binding sites and cooperative interactions. The PriB dimer binds the ssDNA with strong positive cooperativity. Both the intrinsic affinity and cooperative interactions are accompanied by a net ion release, with anions participating in the ion exchange process. The intrinsic binding process is an entropy-driven reaction, suggesting strongly that the DNA association induces a large conformational change in the protein. The PriB protein shows a dramatically strong preference for the homo-pyrimidine oligomers with an intrinsic affinity higher by about three orders of magnitude, as compared to the homo-purine oligomers. The significance of these results for PriB protein activity is discussed.     

MacA,a Periplasmic Membrane Fusion Protein of the Macrolide Transporter MacAB-TolC,Binds Lipopolysaccharide Core Specifically and with High Affinity

The Escherichia coli MacAB-TolC transporter has been implicated in efflux of macrolide antibiotics and secretion of enterotoxin STII. In this study, we found that purified MacA, a periplasmic membrane fusion protein, contains one tightly bound rough core lipopolysaccharide (R-LPS) molecule per MacA molecule. R-LPS was bound specifically to MacA protein with affinity exceeding that of polymyxin B. Sequence analyses showed that MacA contains two high-density clusters of positively charged amino acid residues located in the cytoplasmic N-terminal domain and the periplasmic C-terminal domain. Substitutions in the C-terminal cluster reducing the positive-charge density completely abolished binding of R-LPS. At the same time, these substitutions significantly reduced the functionality of MacA in the protection of E. coli against macrolides in vivo and in the in vitro MacB ATPase stimulation assays. Taken together, our results suggest that R-LPS or a similar glycolipid is a physiological substrate of MacAB-TolC.     

In Human Pseudouridine Synthase 1 (hPus1), a C-Terminal Helical Insert Blocks tRNA from Binding in the Same Orientation as in the Pus1 Bacterial Homologue TruA,Consistent with Their Different Target Selectivities

Human pseudouridine (Ψ) synthase Pus1 (hPus1) modifies specific uridine residues in several non-coding RNAs: tRNA, U2 spliceosomal RNA, and steroid receptor activator RNA. We report three structures of the catalytic core domain of hPus1 from two crystal forms, at 1.8 Å resolution. The structures are the first of a mammalian Ψ synthase from the set of five Ψ synthase families common to all kingdoms of life. hPus1 adopts a fold similar to bacterial Ψ synthases, with a central antiparallel β-sheet flanked by helices and loops. A flexible hinge at the base of the sheet allows the enzyme to open and close around an electropositive active-site cleft. In one crystal form, a molecule of Mes [2-(N-morpholino)ethane sulfonic acid] mimics the target uridine of an RNA substrate. A positively charged electrostatic surface extends from the active site towards the N-terminus of the catalytic domain, suggesting an extensive binding site specific for target RNAs. Two α-helices C-terminal to the core domain, but unique to hPus1, extend along the back and top of the central β-sheet and form the walls of the RNA binding surface. Docking of tRNA to hPus1 in a productive orientation requires only minor conformational changes to enzyme and tRNA. The docked tRNA is bound by the electropositive surface of the protein employing a completely different binding mode than that seen for the tRNA complex of the Escherichia coli homologue TruA.     

Protein-Precursor tRNA Contact Leads to Sequence-Specific Recognition of 5′ Leaders by Bacterial Ribonuclease P

Bacterial ribonuclease P (RNase P) catalyzes the cleavage of 5′ leader sequences from precursor tRNAs (pre-tRNAs). Previously, all known substrate nucleotide specificities in this system are derived from RNA-RNA interactions with the RNase P RNA subunit. Here, we demonstrate that pre-tRNA binding affinities for Bacillus subtilis and Escherichia coli RNase P are enhanced by sequence-specific contacts between the fourth pre-tRNA nucleotide on the 5′ side of the cleavage site (N(− 4)) and the RNase P protein (P protein) subunit. B. subtilis RNase P has a higher affinity for pre-tRNA with adenosine at N(− 4), and this binding preference is amplified at physiological divalent ion concentrations. Measurements of pre-tRNA-containing adenosine analogs at N(− 4) indicate that specificity arises from a combination of hydrogen bonding to the N6 exocyclic amine of adenosine and steric exclusion of the N2 amine of guanosine. Mutagenesis of B. subtilis P protein indicates that F20 and Y34 contribute to selectivity at N(− 4). The hydroxyl group of Y34 enhances selectivity, likely by forming a hydrogen bond with the N(− 4) nucleotide. The sequence preference of E. coli RNase P is diminished, showing a weak preference for adenosine and cytosine at N(− 4), consistent with the substitution of Leu for Y34 in the E. coli P protein. This is the first identification of a sequence-specific contact between P protein and pre-tRNA that contributes to molecular recognition of RNase P. Additionally, sequence analyses reveal that a greater-than-expected fraction of pre-tRNAs from both E. coli and B. subtilis contains a nucleotide at N(− 4) that enhances RNase P affinity. This observation suggests that specificity at N(− 4) contributes to substrate recognition in vivo. Furthermore, bioinformatic analyses suggest that sequence-specific contacts between the protein subunit and the leader sequences of pre-tRNAs may be common in bacterial RNase P and may lead to species-specific substrate recognition.     

Crystal structure of the GluR0 ligand-binding core from Nostoc punctiforme in complex with L-glutamate: structural dissection of the ligand interaction and subunit interface

GluR0 from Nostoc punctiforme (NpGluR0) is a bacterial homologue of the ionotropic glutamate receptor (iGluR). We have solved the crystal structure of the ligand-binding core of NpGluR0 in complex with l-glutamate at a resolution of 2.1 Å. The structure exhibits a noncanonical ligand interaction and two distinct subunit interfaces. The side-chain guanidium group of Arg80 forms a salt bridge with the γ-carboxyl group of bound l-glutamate: in GluR0 from Synechocystis (SGluR0) and other iGluRs, the equivalent residues are Asn or Thr, which cannot form a similar interaction. We suggest that the local positively charged environment and the steric constraint created by Arg80 mediate the selectivity of l-glutamate binding by preventing the binding of positively charged and hydrophobic amino acids. In addition, the NpGluR0 ligand-binding core forms a new subunit interface in which the two protomers are arranged differently than the known iGluR and SGluR0 dimer interfaces. The significance of there being two different dimer interfaces was investigated using analytical ultracentrifugation analysis.     

Coarse-grained models for simulations of multiprotein complexes: application to ubiquitin binding

We develop coarse-grained models and effective energy functions for simulating thermodynamic and structural properties of multiprotein complexes with relatively low binding affinity (K_d > 1 μM) and apply them to binding of Vps27 to membrane-tethered ubiquitin. Folded protein domains are represented as rigid bodies. The interactions between the domains are treated at the residue level with amino-acid-dependent pair potentials and Debye-Hückel-type electrostatic interactions. Flexible linker peptides connecting rigid protein domains are represented as amino acid beads on a polymer with appropriate stretching, bending, and torsion-angle potentials. In simulations of membrane-attached protein complexes, interactions between amino acids and the membrane are described by residue-dependent short-range potentials and long-range electrostatics. We parameterize the energy functions by fitting the osmotic second virial coefficient of lysozyme and the binding affinity of the ubiquitin-CUE complex. For validation, extensive replica-exchange Monte Carlo simulations are performed of various protein complexes. Binding affinities for these complexes are in good agreement with the experimental data. The simulated structures are clustered on the basis of distance matrices between two proteins and ranked according to cluster population. In ∼ 70% of the complexes, the distance root-mean-square is less than 5 Å from the experimental structures. In ∼ 90% of the complexes, the binding interfaces on both proteins are predicted correctly, and in all other cases at least one interface is correct. Transient and nonspecifically bound structures are also observed. With the validated model, we simulate the interaction between the Vps27 multiprotein complex and a membrane-tethered ubiquitin. Ubiquitin is found to bind preferentially to the two UIM domains of Vps27, but transient interactions between ubiquitin and the VHS and FYVE domains are observed as well. These specific and nonspecific interactions are found to be positively cooperative, resulting in a substantial enhancement of the overall binding affinity beyond the ∼ 300 μM of the specific domains. We also find that the interactions between ubiquitin and Vps27 are highly dynamic, with conformational rearrangements enabling binding of Vps27 to diverse targets as part of the multivesicular-body protein-sorting pathway.     

Discrimination between Closely Related Cellular Metabolites by the SAM-I Riboswitch

The SAM-I riboswitch is a cis-acting element of genetic control found in bacterial mRNAs that specifically binds S-adenosylmethionine (SAM). We previously determined the 2.9-Å X-ray crystal structure of the effector-binding domain of this RNA element, revealing details of RNA-ligand recognition. To improve this structure, variations were made to the RNA sequence to alter lattice contacts, resulting in a 0.5-Å improvement in crystallographic resolution and allowing for a more accurate refinement of the crystallographic model. The basis for SAM specificity was addressed by a structural analysis of the RNA complexed to S-adenosylhomocysteine (SAH) and sinefungin and by measuring the affinity of SAM and SAH for a series of mutants using isothermal titration calorimetry. These data illustrate the importance of two universally conserved base pairs in the RNA that form electrostatic interactions with the positively charged sulfonium group of SAM, thereby providing a basis for discrimination between SAM and SAH.     

Probing receptor-translocator interactions in the oligopeptide ABC transporter by fluorescence correlation spectroscopy

The oligopeptide transporter Opp is a five-component ABC uptake system. The extracytoplasmic lipid-anchored substrate-binding protein (or receptor) OppA delivers peptides to an integral membrane complex OppBCDF (or translocator), where, on ATP binding and hydrolysis, translocation across the membrane takes place. OppA and OppBCDF were labeled with fluorescent probes, reconstituted into giant unilamellar vesicles, and the receptor-translocator interactions were investigated by fluorescence correlation spectroscopy. Lateral mobility of OppA was reduced on incorporation of OppBCDF into giant unilamellar vesicles, and decreased even further on the addition of peptide. Fluorescence cross-correlation measurements revealed that OppBCDF distinguished liganded from unliganded OppA, binding only the former. Addition of ATP or its nonhydrolyzable analog AMP-PNP resulted in release of OppA from OppBCDF. In vanadate-trapped “transition state” conditions, OppA also was not bound by OppBCDF. A model is presented in which ATP-binding to OppDF results in donation of peptide to OppBC and simultaneous release of OppA. ATP-hydrolysis would complete the peptide translocation and reset the transporter for another catalytic cycle. Implications in terms of a general transport mechanism for ABC importers and exporters are discussed.     

More than a Simple Lipophilic Contact: A Detailed Thermodynamic Analysis of Nonbasic Residues in the S1 Pocket of Thrombin

The field of medicinal chemistry aims to design and optimize small molecule leads into drug candidates that may positively interfere with pathological disease situations in humans or combat the growth of infective pathogens. From the plethora of crystal structures of protein-inhibitor complexes we have learned how molecules recognize each other geometrically, but we still have rather superficial understanding of why they bind to each other. This contribution surveys a series of 26 thrombin inhibitors with small systematic structural differences to elucidate the rationale for their widely deviating binding affinity from 185 μM to 4 nM as recorded by enzyme kinetic measurements. Five well-resolved (resolution 2.30 - 1.47 Å) crystal structures of thrombin-inhibitor complexes and an apo-structure of the uncomplexed enzyme (1.50 Å) are correlated with thermodynamic data recorded by isothermal titration calorimetry with 12 selected inhibitors from the series. Taking solubility data into account, the variation in physicochemical properties allows conclusions to be reached about the relative importance of the enthalpic binding features as well as to estimate the importance of the parameters more difficult to capture, such as residual ligand entropy and desolvation properties. The collected data reveal a comprehensive picture of the thermodynamic signature that explains the so far poorly understood attractive force experienced by m-chloro-benzylamides to thrombin.