Crystallographic and calorimetric analysis of peptide binding to OppA protein

Identification of Ni-(L-His)₂ as a substrate for NikABCDE-dependent nickel uptake in Escherichia coli

Nickel is an important cofactor for several microbial enzymes. The ATP-dependent NikABCDE transporter is one of several types of uptake pathways known to be important for nickel acquisition in microbes. The Escherichia coli NikA periplasmic binding protein is structurally homologous to the di- and oligopeptide binding proteins, DppA and OppA. This structural similarity raises interesting questions regarding the evolutionary relationships between the recognition of nickel ions and short peptides. We find that in defined minimal growth medium NikABCDE transports nickel ions in the presence of exogenously added L-histidine (L-His), but not D-histidine. Both nickel uptake in cells and nickel binding to purified NikA showed an L-His concentration dependence consistent with recognition of a Ni-(L-His)₂ complex. This discovery reveals parallels to the transport of other metal complexes, notably iron, and suggests the structural diversity of nickel transporters may arise from the need to recognize extracellular nickel complexed with different organic ligands, whether they be exogenously or endogenously produced. Further, these results suggest that experiments examining the physiology and ecology of nickel-requiring microbes should account for the possibility that the growth medium may not support nickel uptake.

Rapid and accurate prediction and scoring of water molecules in protein binding sites

Water plays a critical role in ligand-protein interactions. However, it is still challenging to predict accurately not only where water molecules prefer to bind, but also which of those water molecules might be displaceable. The latter is often seen as a route to optimizing affinity of potential drug candidates. Using a protocol we call WaterDock, we show that the freely available AutoDock Vina tool can be used to predict accurately the binding sites of water molecules. WaterDock was validated using data from X-ray crystallography, neutron diffraction and molecular dynamics simulations and correctly predicted 97% of the water molecules in the test set. In addition, we combined data-mining, heuristic and machine learning techniques to develop probabilistic water molecule classifiers. When applied to WaterDock predictions in the Astex Diverse Set of protein ligand complexes, we could identify whether a water molecule was conserved or displaced to an accuracy of 75%. A second model predicted whether water molecules were displaced by polar groups or by non-polar groups to an accuracy of 80%. These results should prove useful for anyone wishing to undertake rational design of new compounds where the displacement of water molecules is being considered as a route to improved affinity.

Impact of ligand and protein desolvation on ligand binding to the S1 pocket of thrombin

In the present study, we investigate the impact of a tightly bound water molecule on ligand binding in the S1 pocket of thrombin. The S1 pocket contains a deeply buried deprotonated aspartate residue (Asp189) that is, due to its charged state, well hydrated in the uncomplexed state. We systematically studied the importance of this water molecule by evaluating a series of ligands that contains pyridine-type P1 side chains that could potentially alter the binding properties of this water molecule. All of the pyridine derivatives retain the original hydration state albeit sometimes with a slight perturbance. In order to prevent a direct H-bond formation with Asp189, and to create a permanent positive charge on the P1 side chain that is positioned adjacent to the Asp189 carboxylate anion, we methylated the pyridine nitrogen. This methylation resulted in displacement of water but was accompanied by a loss in binding affinity. Quantum chemical calculations of the ligand solvation free energy showed that the positively charged methylpyridinium derivatives suffer a large penalty of desolvation upon binding. Consequently, they have a substantially less favorable enthalpy of binding. In addition to the ligand desolvation penalty, the hydration shell around Asp189 has to be overcome, which is achieved in nearly all pyridinium derivatives. Only for the ortho derivative is a partial population of a water next to Asp189 found. Possibly, the gain of electrostatic interactions between the charged P1 side chain and Asp189 helps to compensate for the desolvation penalty. In all uncharged pyridine derivatives, the solvation shell remains next to Asp189, partly mediating interactions between ligand and protein. In the case of the para-pyridine derivative, a strongly disordered cluster of water sites is observed between ligand and Asp189.

Interfacial water molecules in SH3 interactions: a revised paradigm for polyproline recognition

In spite of its biomedical relevance, polyproline recognition is still not fully understood. The disagreement between the current description of SH3 (Src homology 3) complexes and their thermodynamic behaviour calls for a revision of the SH3-binding paradigm. Recently, Abl-SH3 was demonstrated to recognize its ligands by a dual binding mechanism involving a robust network of water-mediated hydrogen bonds that complements the canonical hydrophobic interactions. The systematic analysis of the SH3 structural database in the present study reveals that this dual binding mode is universal to SH3 domains. Tightly bound buried-interfacial water molecules were found in all SH3 complexes studied mediating the interaction between the peptide ligand and the domain. Moreover, structural waters were also identified in a high percentage of the free SH3 domains. A detailed analysis of the pattern of water-mediated interactions enabled the identification of conserved hydration sites in the polyproline-recognition region and the establishment of relationships between hydration profiles and the sequence of both ligands and SH3 domains. Water-mediated interactions were also systematically observed in WW (protein–protein interaction domain containing two conserved tryptophan residues), UEV (ubiquitin-conjugating enzyme E2 variant) and EVH-1 [Ena/VASP (vasodilator-stimulated phosphoprotein) homology 1] structures. The results of the present study clearly indicate that the current description of proline-rich sequence recognition by protein–protein interaction modules is incomplete and insufficient for a correct understanding of these systems. A new binding paradigm is required that includes interfacial water molecules as relevant elements in polyproline recognition.

Systematic placement of structural water molecules for improved scoring of protein-ligand interactions

Structural water molecules are found in many protein-ligand complexes. They are known to be vital in mediating hydrogen-bonding interactions and, in some cases, key for facilitating tight binding. It is thus very important to consider water molecules when attempting to model protein-ligand interactions for cognate ligand identification, virtual screening and drug design. While the rigid treatment of water molecules present in structures is feasible, the more relevant task of treating all possible positions and orientations of water molecules with each possible ligand pose is computationally daunting. Current methods in molecular docking provide partial treatment for such water molecules, with modest success. Here we describe a new method employing dead-end elimination to place water molecules within a binding site, bridging interactions between protein and ligand. Dead-end elimination permits a thorough, though still incomplete, treatment of water placement. The results show that this method is able to place water molecules correctly within known complexes and to create physically reasonable hydrogen bonds. The approach has also been incorporated within an inverse molecular design approach, to model a variety of compounds in the process of de novo ligand design. The inclusion of structural water molecules, combined with ranking based on the electrostatic contribution to binding affinity, improves a number of otherwise poor energetic predictions.

Compensating stereochemical changes allow murein tripeptide to be accommodated in a conventional peptide-binding protein

The oligopeptide permease (Opp) of Escherichia coli is an ATP-binding cassette transporter that uses the substrate-binding protein (SBP) OppA to bind peptides and deliver them to the membrane components (OppBCDF) for transport. OppA binds conventional peptides 2-5 residues in length regardless of their sequence, but does not facilitate transport of the cell wall component murein tripeptide (Mtp, L-Ala-γ-D-Glu-meso-Dap), which contains a D-amino acid and a γ-peptide linkage. Instead, MppA, a homologous substrate-binding protein, forms a functional transporter with OppBCDF for uptake of this unusual tripeptide. Here we have purified MppA and demonstrated biochemically that it binds Mtp with high affinity (K(D) ∼ 250 nM). The crystal structure of MppA in complex with Mtp has revealed that Mtp is bound in a relatively extended conformation with its three carboxylates projecting from one side of the molecule and its two amino groups projecting from the opposite face. Specificity for Mtp is conferred by charge-charge and dipole-charge interactions with ionic and polar residues of MppA. Comparison of the structure of MppA-Mtp with structures of conventional tripeptides bound to OppA, reveals that the peptide ligands superimpose remarkably closely given the profound differences in their structures. Strikingly, the effect of the D-stereochemistry, which projects the side chain of the D-Glu residue at position 2 in the direction of the main chain in a conventional tripeptide, is compensated by the formation of a γ-linkage to the amino group of diaminopimelic acid, mimicking the peptide bond between residues 2 and 3 of a conventional tripeptide.

Role of bridging water molecules in GSK3β-inhibitor complexes: insights from QM/MM, MD, and molecular docking studies

The role of water molecules is increasingly gaining interest in drug design, and several studies have highlighted their paramount contributions to the specificity and the affinity of ligand binding. In this study, we employ the two-layer ONIOM-based quantum mechanics/molecular mechanics (QM/MM) calculations, molecular dynamics (MD) simulations, and molecular docking studies to investigate the effect of bridging water molecules at the GSK3β-inhibitors interfaces. The results obtained from the ONIOM geometry optimization and AIM analysis corroborated the presence of bridging water molecules that form hydrogen bonds with protein side chain of Thr138 and/or backbone of Gln185, and mediate interactions with inhibitors in the 10 selected GSK3β-inhibitor complexes. Subsequently, MD simulations carried out on a representative system of 1R0E demonstrated that the bridging water molecule is stable at the GSK3β-inhibitor interface and appears to contribute to the stability of the protein-inhibitor interactions. Furthermore, molecular docking studies of GSK3β-inhibitor complexes indicated that the inhibitors can increase binding affinities and the better docked conformation of inhibitors can be obtained by inclusion of the bridging water molecules, especially for the flexible inhibitors, in docking experiments into individual protein conformations. Our results elucidate the importance of bridging water molecules at the GSK3β-inhibitor interfaces and suggest that they might prove useful in rational drug design.

Structural analysis of semi-specific oligosaccharide recognition by a cellulose-binding protein of thermotoga maritima reveals adaptations for functional diversification of the oligopeptide periplasmic binding protein fold

Periplasmic binding proteins (PBPs) constitute a protein superfamily that binds a wide variety of ligands. In prokaryotes, PBPs function as receptors for ATP-binding cassette or tripartite ATP-independent transporters and chemotaxis systems. In many instances, PBPs bind their cognate ligands with exquisite specificity, distinguishing, for example, between sugar epimers or structurally similar anions. By contrast, oligopeptide-binding proteins bind their ligands through interactions with the peptide backbone but do not distinguish between different side chains. The extremophile Thermotoga maritima possesses a remarkable array of carbohydrate-processing metabolic systems, including the hydrolysis of cellulosic polymers. Here, we present the crystal structure of a T. maritima cellobiose-binding protein (tm0031) that is homologous to oligopeptide-binding proteins. T. maritima cellobiose-binding protein binds a variety of lengths of beta(1-->4)-linked glucose oligomers, ranging from two rings (cellobiose) to five (cellopentaose). The structure reveals that binding is semi-specific. The disaccharide at the nonreducing end binds specifically; the other rings are located in a large solvent-filled groove, where the reducing end makes several contacts with the protein, thereby imposing an upper limit of the oligosaccharides that are recognized. Semi-specific recognition, in which a molecular class rather than individual species is selected, provides an efficient solution for the uptake of complex mixtures.

A balancing act between net uptake of water during dihydrofolate binding and net release of water upon NADPH binding in R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. This enzyme is a homotetramer possessing 222 symmetry, and a single active site pore traverses the length of the protein. A promiscuous binding surface can accommodate either DHF or NADPH, thus two nonproductive complexes can form (2NADPH or 2DHF) as well as a productive complex (NADPH.DHF). The role of water in binding was monitored using a number of different osmolytes. From isothermal titration calorimetry (ITC) studies, binding of NADPH is accompanied by the net release of 38 water molecules. In contrast, from both steady state kinetics and ITC studies, binding of DHF is accompanied by the net uptake of water. Although different osmolytes have similar effects on NADPH binding, variable results are observed when DHF binding is probed. Sensitivity to water activity can also be probed by an in vivo selection using the antibacterial drug, trimethoprim, where the water content of the media is decreased by increasing concentrations of sorbitol. The ability of wild type and mutant clones of R67 DHFR to allow host Escherichia coli to grow in the presence of trimethoprim plus added sorbitol parallels the catalytic efficiency of the DHFR clones, indicating water content strongly correlates with the in vivo function of R67 DHFR.

Evidence that bacterial ABC-type transporter imports free EDTA for metabolism

EDTA, a common chelating agent, is becoming a major organic pollutant in the form of metal-EDTA complexes in surface waters, partly due to its recalcitrance to biodegradation. Even an EDTA-degrading bacterium, BNC1, does not degrade stable metal-EDTA complexes. In the present study, an ABC-type transporter was identified for possible uptake of EDTA because the transporter genes and the EDTA monooxygenase gene were expressed from a single operon in BNC1. The ABC-type transporter had a periplasmic-binding protein (EppA) that should confer the substrate specificity for the transporter; therefore, EppA was produced in Escherichia coli, purified, and characterized. EppA was shown to bind free EDTA with a dissociation constant as low as 25 nM by using isothermal titration calorimetry. When unstable metal-EDTA complexes, e.g., (Mg-EDTA)(2-), were added to the EppA solution, binding was also observed. However, experimental data and theoretical analysis supported EppA binding only of free EDTA. When stable metal-EDTA complexes, e.g., (Cu-EDTA)(2-), were titrated into the EppA solution, no binding was observed. Since EDTA monooxygenase in the cytoplasm uses some of the stable metal-EDTA complexes as substrates, we suggest that the lack of EppA binding and EDTA uptake are responsible for the failure of BNC1 cells to degrade the stable complexes.

MALISAM: a database of structurally analogous motifs in proteins

MALISAM (manual alignments for structurally analogous motifs) represents the first database containing pairs of structural analogs and their alignments. To find reliable analogs, we developed an approach based on three ideas. First, an insertion together with a part of the evolutionary core of one domain family (a hybrid motif) is analogous to a similar motif contained within the core of another domain family. Second, a motif at an interface, formed by secondary structural elements (SSEs) contributed by two or more domains or subunits contacting along that interface, is analogous to a similar motif present in the core of a single domain. Third, an artificial protein obtained through selection from random peptides or in sequence design experiments not biased by sequences of a particular homologous family, is analogous to a structurally similar natural protein. Each analogous pair is superimposed and aligned manually, as well as by several commonly used programs. Applications of this database may range from protein evolution studies, e.g. development of remote homology inference tools and discriminators between homologs and analogs, to protein-folding research, since in the absence of evolutionary reasons, similarity between proteins is caused by structural and folding constraints. The database is publicly available at http://prodata.swmed.edu/malisam.

Synthesis of a pyridoxine-peptide based delivery system for nucleotides

As a first step in the development of a nucleotide delivery system making use of oligopeptide permease, we have synthesized pyridoxine-peptide-nucleotide conjugates. The nucleotides are bonded on serine residues. The peptides terminate with a pyroglutamate residue. In the first example, the pyridoxine moiety is connected at the end of the peptides, while, in the second example, the pyridoxine moiety is bonded at an aspartic acid residue in a middle position of the peptide. Nucleotides are introduced as phosphoramidites. The synthetic strategy involves a series of protection, deprotection, and coupling reactions, and integrates peptide, nucleotide, and pyridoxine chemistry. The final deprotection step was carried out in basic conditions using 10% K2CO3 in MeOH.

Classification of water molecules in protein binding sites

Water molecules play a crucial role in mediating the interaction between a ligand and a macromolecular receptor. An understanding of the nature and role of each water molecule in the active site of a protein could greatly increase the efficiency of rational drug design approaches: if the propensity of a water molecule for displacement can be determined, then synthetic effort may be most profitably applied to the design of specific ligands with the displacement of this water molecule in mind. In this paper, a thermodynamic analysis of water molecules in the binding sites of six proteins, each complexed with a number of inhibitors, is presented. Two classes of water molecules were identified: those conserved and not displaced by any of the ligands, and those that are displaced by some ligands. The absolute binding free energies of 54 water molecules were calculated using the double decoupling method, with replica exchange thermodynamic integration in Monte Carlo simulations. It was found that conserved water molecules are on average more tightly bound than displaced water molecules. In addition, Bayesian statistics is used to calculate the probability that a particular water molecule may be displaced by an appropriately designed ligand, given the calculated binding free energy of the water molecule. This approach therefore allows the numerical assessment of whether or not a given water molecule should be targeted for displacement as part of a rational drug design strategy.

Analysis of ligand-bound water molecules in high-resolution crystal structures of protein-ligand complexes

We have performed a comprehensive analysis of water molecules at the protein-ligand interfaces observed in 392 high-resolution crystal structures. There are a total of 1829 ligand-bound water molecules in these 392 complexes; 18% are surface water molecules, and 72% are interfacial water molecules. The number of ligand-bound water molecules in each complex structure ranges from 0 to 21 and has an average of 4.6. Of these interfacial water molecules, 76% are considered to be bridging water molecules, characterized by having polar interactions with both ligand and protein atoms. Among a number of factors that may influence the number of ligand-bound water molecules, the polar van der Waals (vdw) surface area of ligands has the highest Pearson linear correlation coefficient of 0.63. Our regression analysis predicted that one more ligand-bound water molecule is expected for every additional 24 A2 in the polar vdw surface area of the ligand. In contrast to the observation that the resolution is the primary factor influencing the number of water molecules in crystallographic models of proteins, we found that there is only a weak relationship between the number of ligand-bound water molecules and the resolution of the crystal structures. An analysis of the isotropic B factors of buried ligand-bound water molecules suggested that, when water molecules have fewer than two polar interactions with the protein-ligand complex, they are more mobile than protein atoms in the crystal structures; when they have more than three polar interactions, they are significantly less mobile than protein atoms.

Do Structurally Similar Ligands Bind in a Similar Fashion?

The scope of the current work is to investigate whether structurally similar ligands bind in a similar fashion by exhaustively analyzing experimental data from the protein database (PDB). The complete PDB was searched for pairs of structurally similar ligands binding to the same biological target. The binding sites of the pairs of proteins complexing structurally similar ligands were found to differ in 83% of the cases. The most recurrent structural change among the pairs involves different water molecule architecture. Side-chain movements are observed in half of the pairs, whereas backbone movements rarely occurred. However, two structurally similar ligands generally confirm a high degree of structural conservation. That is, a majority of the ligand pairs occupy the same region in the binding sites, providing support for the use of shape matching in the drug design process. We allow ourselves to draw general conclusions because our data set consists of ligands with drug-like physicochemical properties complexed to a broad spectrum of different protein classes.

WATER MEDIATION IN PROTEIN FOLDING AND MOLECULAR RECOGNITION

Water is essential for life in many ways, and without it biomolecules might no longer truly be biomolecules. In particular, water is important to the structure, stability, dynamics, and function of biological macromolecules. In protein folding, water mediates the collapse of the chain and the search for the native topology through a funneled energy landscape. Water actively participates in molecular recognition by mediating the interactions between binding partners and contributes to either enthalpic or entropic stabilization. Accordingly, water must be included in recognition and structure prediction codes to capture specificity. Thus water should not be treated as an inert environment, but rather as an integral and active component of biomolecular systems, where it has both dynamic and structural roles. Focusing on water sheds light on the physics and function of biological machinery and self-assembly and may advance our understanding of the natural design of proteins and nucleic acids.

Protein translocation by the Sec61/SecY channel

The conserved protein-conducting channel, referred to as the Sec61 channel in eukaryotes or the SecY channel in eubacteria and archaea, translocates proteins across cellular membranes and integrates proteins containing hydrophobic transmembrane segments into lipid bilayers. Structural studies illustrate how the protein-conducting channel accomplishes these tasks. Three different mechanisms, each requiring a different set of channel binding partners, are employed to move polypeptide substrates: The ribosome feeds the polypeptide chain directly into the channel, a ratcheting mechanism is used by the eukaryotic endoplasmic reticulum chaperone BiP, and a pushing mechanism is utilized by the bacterial ATPase SecA. We review these translocation mechanisms, relating biochemical and genetic observations to the structures of the protein-conducting channel and its binding partners.

LigScore: a novel scoring function for predicting binding affinities

We present two new empirical scoring functions, LigScore1 and LigScore2, that attempt to accurately predict the binding affinity between ligand molecules and their protein receptors. The LigScore functions consist of three distinct terms that describe the van der Waals interaction, the polar attraction between the ligand and protein, and the desolvation penalty attributed to the binding of the polar ligand atoms to the protein and vice versa. Utilizing a regression approach on a data set of 118 protein-ligand complexes we have obtained a linear equation, LigScore2, using these three descriptors. LigScore2 has good predictability with regard to experimental pKi values yielding a correlation coefficient, r2), of 0.75 and a standard deviation of 1.04 over the training data set, which consists of a diverse set of proteins that span more than seven protein families.

The Structure of the Oligopeptide-binding Protein, AppA, from Bacillus subtilis in Complex with a Nonapeptide

Besides their role as a source of amino acids for Bacillus subtilis, exogenous peptides play important roles in the signalling pathways leading to the development of competence and sporulation. B.subtilis has three peptide transport systems all belonging to the ATP-binding cassette family, a dipeptide permease (Dpp) and two oligopeptide permeases (Opp and App) with overlapping specificity. These comprise a membrane-spanning channel through which the peptide passes, a pair of ATPases which couple ATP hydrolysis to peptide translocation and a lipid-modified, membrane-anchored extracellular "binding-protein" that serves as the receptor for the system. Here, we present the crystal structure of a soluble form of the peptide-binding protein AppA, which has been solved to 1.6 A spacing by anomalous scattering and molecular replacement methods. The structure reveals a protein made of two distinct lobes with a topology similar to those of DppA from Escherichia coli and OppA from Salmonella typhimurium. Examination of the interlobe region reveals an enlarged pocket, containing electron density defining a nonapeptide ligand. The main-chain of the peptide is well defined and makes a series of polar contacts with the protein including salt-bridges at both its termini. The side-chain density is ambiguous in places, consistent with the interpretation that a population of peptides is bound, whose average electron density resembles the amino acid sequence N-VDSKNTSSW-C.

Thermodynamic consequences of disrupting a water-mediated hydrogen bond network in a protein:pheromone complex

The mouse pheromones (+/-)-2-sec-butyl-4,5-dihydrothiazole (SBT) and 6-hydroxy-6-methyl-3-heptanone (HMH) bind into an occluded hydrophobic cavity in the mouse major urinary protein (MUP-1). Although the ligands are structurally unrelated, in both cases binding is accompanied by formation of a similar buried, water-mediated hydrogen bond network between the ligand and several backbone and side chain groups on the protein. To investigate the energetic contribution of this hydrogen bond network to ligand binding, we have applied isothermal titration calorimetry to measure the binding thermodynamics using several MUP mutants and ligand analogs. Mutation of Tyr-120 to Phe, which disrupts a hydrogen bond from the phenolic hydroxyl group of Tyr-120 to one of the bound water molecules, results in a substantial loss of favorable binding enthalpy, which is partially compensated by a favorable change in binding entropy. A similar thermodynamic effect was observed when the hydrogen bonded nitrogen atom of the heterocyclic ligand was replaced by a methyne group. Several other modifications of the protein or ligand had smaller effects on the binding thermodynamics. The data provide supporting evidence for the role of the hydrogen bond network in stabilizing the complex.

Evaluation of the relative stability of liganded versus ligand-free protein conformations using Simplicial Neighborhood Analysis of Protein Packing (SNAPP) method

Many proteins change their conformation upon ligand binding. For instance, bacterial periplasmic binding proteins (bPBPs), which transport nutrients into the cytoplasm, generally consist of two globular domains connected by strands, forming a hinge. During ligand binding, hinge motion changes the conformation from the open to the closed form. Both forms can be crystallized without a ligand, suggesting that the energy difference between them is small. We applied Simplicial Neighborhood Analysis of Protein Packing (SNAPP) as a method to evaluate the relative stability of open and closed forms in bPBPs. Using united residue representation of amino acids, SNAPP performs Delaunay tessellation of the protein, producing an aggregate of space-filling, irregular tetrahedra with nearest neighbor residues at the vertices. The SNAPP statistical scoring function is derived from log-likelihood scores for all possible quadruplet compositions of amino acids found in a representative subset of the Protein Data Bank, and the sum of the scores for a given protein provides the total SNAPP score. Results of scoring for bPBPs suggest that in most cases, the unliganded form is more stable than the liganded form, and this conclusion is corroborated by similar observations of other proteins undergoing conformation changes upon binding their ligands. The results of these studies suggest that the SNAPP method can be used to predict the relative stability of accessible protein conformations. Furthermore, the SNAPP method allows delineation of the role of individual residues in protein stabilization, thereby providing new testable hypotheses for rational site-directed mutagenesis in the context of protein engineering.

Enterococcal peptide sex pheromones: synthesis and control of biological activity

The enterococcal pheromone-inducible plasmids such as pCF10 represent a unique class of mobile genetic elements whose transfer functions are induced by peptide sex pheromones. These pheromones are excreted by potential recipient cells and detected by plasmid-containing donor cells at the cell surface, where the pheromone is imported and signals induction of the plasmid transfer system. Pheromone is processed from a chromosomally encoded lipoprotein and excreted by both the donor and recipient cells, but a carefully controlled detection system prevents a response to self-pheromone while still allowing an extremely sensitive response to exogenous pheromone.

A large conformational change of the translocation ATPase SecA

The ATPase SecA mediates the posttranslational translocation of a wide range of polypeptide substrates through the SecY channel in the cytoplasmic membrane of bacteria. We have determined the crystal structure of a monomeric form of Bacillus subtilis SecA at a 2.2-A resolution. A comparison with the previously determined structures of SecA reveals a nucleotide-independent, large conformational change that opens a deep groove similar to that in other proteins that interact with diverse polypeptides. We propose that the open form of SecA represents an activated state.

Analysis of differences in the functional properties of the substrate binding proteins of the Borrelia burgdorferi oligopeptide permease (Opp) operon

The Borrelia burgdorferi genome encodes five orthologues of the substrate binding protein oligopeptide permease A (OppA). It was previously shown that these genes are under the control of separate promoters and are differentially expressed under various environmental conditions. We were interested in determining whether there are also differences in substrate specificities among the proteins. The substrate specificities of recombinant proteins were determined by screening for high-affinity peptides by use of a combinatorial phage display heptapeptide library. Different heptapeptides with high affinities for OppA-1, OppA-2, and OppA-3 were identified. No heptapeptide binding OppA-4 or OppA-5 could be identified. Competitive binding assays were performed under various conditions to determine the substrate preferences of the OppA proteins. OppA-1 retained maximal activity over a broad range of pHs (5.5 to 7.5), whereas OppA-2 and OppA-3 showed peak activities at pHs below 5.5. OppA-1 and OppA-2 showed preferences for tripeptides over dipeptides and longer-chain peptides. Although a wide variety of amino acyl side chains were tolerated by all three OppA proteins, OppA-1 showed the broadest substrate specificity and was able to accommodate peptides composed of bulky hydrophobic residues; OppA-2 and OppA-3 showed preferences for peptides composed of small nonpolar amino acids. All three OppA proteins showed preferences for peptides composed of L- rather than D-amino acids. OppA-3 showed the greatest tolerance for changes in stereochemistry. Substantial differences in the substrate specificities of the OppA proteins of B. burgdorferi suggest that they may have distinct functions in the organism.

Peptide-protein interactions: photoinduced electron-transfer within the preformed and encounter complexes of a designed metallopeptide and cytochrome c

Photoinduced electron-transfer (ET) occurs between a negatively charged metallopeptide, [Ru(bpy)(2)(phen-am)-Cys-(Glu)(5)-Gly](3-) = RuCE(5)G, and ferricytochrome c = Cyt c. In the presence of Cyt c, the triplet state lifetime of the ruthenium metallopeptide is shortened, and the emission decays via biexponential kinetics, which indicates the existence of two excited-state populations of ruthenium peptides. The faster decay component displays concentration-independent kinetics demonstrating the presence of a preformed peptide-protein complex that undergoes intra-complex electron-transfer. Values of K(b) = (3.5 +/- 0.2) x 10(4) M(-1) and k(obs)(ET)= (2.7 +/- 0.4) x 10(6) s(-1) were observed at ambient temperatures. The magnitude of k(obs)(ET) decreases with increasing solvent viscosity, and the behavior can be fit to the expression k(obs)(ET) proportional to eta(-alpha) to give alpha = 0.59 +/- 0.05. The electron-transfer process occurring in the preformed complex is therefore gated by a rate-limiting configurational change of the complex. The slower decay component displays concentration-dependent kinetics that saturate at high concentrations of Cyt c. Analysis according to rapid equilibrium formation of an encounter complex that undergoes unimolecular electron-transfer yields K(b)' = (2.5 +/- 0.7) x 10(4) M(-1) and k(obs')(ET)= (7 +/- 3) x 10(5) s(-1). The different values of k(obs)(ET) and k(obs')(ET) suggest that the peptide lies farther from the heme when in the encounter complex. The value of k(obs')(ET) is viscosity dependent indicating that the reaction occurring within the encounter complex is also configurationally gated. A value of alpha = 0.98 +/- 0.14 is observed for k(obs')(ET), which suggests that the rate-limiting gating processes in the encounter complex is different from that in the preformed complex.

Utilising structural knowledge in drug design strategies: applications using Relibase

The concept of structure-based drug design is based upon an in-depth understanding of the principles of molecular recognition. Despite our lack of a thorough comprehension of these principles, the wealth of protein structures available opens up unprecedented possibilities for new insights from the analysis of these data. Unravelling universal rules of molecular recognition is certainly one of the most appealing goals. But our knowledge is enhanced also when studying the specific determinants that characterise single targets or target families only, and the factors governing and discriminating their recognition properties.Here, we illustrate how the structure-based design process can benefit from the consequent incorporation of database query tools. We discuss representative examples to address issues such as protein flexibility, water molecules in binding pockets, and ligand specificity as some of the most critical aspects of drug design. All studies are carried out using our database system Relibase. We also show the application of Relibase in searching for preferred geometrical patterns between interacting molecular fragments.

The role of hydrogen bonding via interfacial water molecules in antigen-antibody complexation. The HyHEL-10-HEL interaction

To study the role of hydrogen bonding via interfacial water molecules in protein-protein interactions, we examined the interaction between hen egg white lysozyme (HEL) and its HyHEL-10 variable domain fragment (Fv) antibody. We constructed three antibody mutants (l-Y50F, l-S91A, and l-S93A) and investigated the interactions between the mutant Fvs and HEL. Isothermal titration calorimetry indicated that the mutations significantly decreased the negative enthalpy change (8-25 kJ mol(-1)), despite some offset by a favorable entropy change. X-ray crystallography demonstrated that the complexes had nearly identical structures, including the positions of the interfacial water molecules. Taken together, the isothermal titration calorimetric and x-ray crystallographic results indicate that hydrogen bonding via interfacial water enthalpically contributes to the Fv-HEL interaction despite the partial offset because of entropy loss, suggesting that hydrogen bonding stiffens the antigen-antibody complex.

Trapping of peptide-based surrogates in an artificially created channel of cytochrome c peroxidase

As recently described, the deliberate removal of the proposed electron transfer pathway from cytochrome c peroxidase resulted in the formation of an extended ligand-binding channel. The engineered channel formed a template for the removed peptide segment, suggesting that synthetic surrogates might be introduced to replace the native electron transfer pathway. This approach could be united with the recent development of sensitizer-linked substrates to initiate and study electron transfer, allowing access to unresolved issues about redox mechanism of the enzyme. Here, we present the design, synthesis, and screening of a peptide library containing natural and unnatural amino acids to identify the structural determinants for binding this channel mutant. Only one peptide, (benzimidazole-propionic acid)-Gly-Ala-Ala, appeared to interact, and gave evidence for both reversible and kinetically trapped binding, suggesting multiple conformations for the channel protein. Notably, this peptide was the most analogous to the removed electron transfer sequence, supporting the use of a cavity-template strategy for design of specific sensitizer-linked substrates as replacements for the native electron transfer pathway.

Comparative binding energy (COMBINE) analysis of OppA-peptide complexes to relate structure to binding thermodynamics

The periplasmic oligopeptide binding component (OppA) of the oligopeptide permease found in Gram-negative bacteria acts as a receptor for peptide transport across the cell membrane and is a potential target for antibacterial drug design. OppA exhibits broad specificity, binding to diverse peptides of 2-5 amino acid residues length. Crystallographic and calorimetric measurements have been carried out by Tame et al. of the binding of 28 peptides of sequence K-X-K to OppA, where X is a natural or nonnatural amino acid. Despite this extensive experimental characterization, a clear relationship between structural and thermodynamic parameters could not be readily identified, with a complicating factor being the observation of varying numbers of water molecules at the binding interface in the different complexes. Consequently, we have applied COMparative BINding Energy (COMBINE) analysis to derive quantitative structure-activity relationships (QSARs) for these 28 OppA-tripeptide complexes. This is the first application of COMBINE analysis to predict binding enthalpies and entropies, and predictive QSAR models were obtained for these quantities as well as for binding free energies. These QSAR models highlight several protein residues and bound water molecules in the binding site, as well as the electrostatic desolvation energies of the protein and the peptides, as responsible for most of the differences in binding thermodynamics between the peptides studied. The QSAR models aid rationalization of the determinants of binding affinity of the OppA:peptide complexes and provide guides for further ligand design. This study also points to the general applicability of COMBINE analysis to estimating thermodynamic parameters for protein-peptide complexes.

Mcg light chain dimer as a model system for ligand design: a docking study

Mcg light chain dimer has been extensively studied by crystallography and peptide binding studies to investigate its peptide cross-reactivity as well as to use it as a model system for designing space filling peptide ligands. Here we extend these investigations by utilizing automated docking. Mcg light chain dimer is an ideal model system for such study due to the availability of experimental data for both the native structure and the 14 complexes with various peptide ligands. We show the ability of the docking approach to reproduce the experimental structures and discuss the limitations associated with such outcomes. We demonstrate the usefulness of the docking approach in generating structural information otherwise not available from the experiment.

Isothermal titration calorimetry in drug discovery

Isothermal titration calorimetry (ITC) follows the heat change when a test compound binds to a target protein. It allows precise measurement of affinity. The method is direct, making interpretation facile, because there is no requirement for competing molecules. Titration in the presence of other ligands rapidly provides information on the mechanism of action of the test compound, identifying the intermolecular complexes that are relevant for structure-based design. Calorimetry allows measurement of stoichiometry and so evaluation of the proportion of the sample that is functional. ITC can characterize protein fragments and catalytically inactive mutant enzymes. It is the only technique which directly measures the enthalpy of binding (delta H degree). Interpretation of delta H degree and its temperature dependence (delta Cp) is usually qualitative, not quantitative. This is because of complicated contributions from linked equilibria and a single change in structure giving modification of several physicochemical properties. Measured delta H degree values allow characterization of proton movement linked to the association of protein and ligand, giving information on the ionization of groups involved in binding. Biochemical systems characteristically exhibit enthalpy-entropy compensation where increased bonding is offset by an entropic penalty, reducing the magnitude of change in affinity. This also causes a lack of correlation between the free energy of binding (delta G degree) and delta H degree. When characterizing structure-activity relationships (SAR), most groups involved in binding can be detected as contributing to delta H degree, but not to affinity. Large enthalpy changes may reflect a modified binding mode, or protein conformation changes. Thus, delta H degree values may highlight a potential discontinuity in SAR, so that experimental structural data are likely to be particularly valuable in molecular design.

Evaluation of the conformational propensities of peptide isosteres as a basis for selecting bioactive pseudopeptides

Our aim was to compare the repertoires of conformers formed by the model zwitterionic peptides AA and AAA in aqueous solution with the conformational profiles of a range of their peptide isosteres, so as to facilitate selection of isosteres for synthesis and testing as biologically stable surrogates of bioactive di- and tripeptides. Comparisons were based upon the results of conformational analysis using a random search approach implemented within the SYBYL molecular modelling package, using zwitterionic molecules, simulated aqueous solvation using a dielectric constant of 80 and allowing all torsions to vary. For each compound, individual conformers were grouped on the basis of specific combinations of psi, phi and omega torsions and, using their energies, the aggregated percentage for each group was calculated using a Boltzmann distribution and displayed using a 3D pseudo Ramachandran plot relating percentage conformer to psi and phi torsions. Retroamide, N-methylamide and thioamide isosteres showed the best match to natural peptides and to the molecular recognition parameters defined for substrates of peptide transporters. The results should aid rational design of therapeutic agents in various areas, e.g. oral delivery of drugs by peptide transporters and of peptidase inhibitors. This approach may usefully be applied to various biochemical and pharmaceutical topics.

One site fits both: a model for the ternary complex of folate + NADPH in R67 dihydrofolate reductase, a D2 symmetric enzyme

R67 dihydrofolate reductase (DHFR) is a novel enzyme that confers resistance to the antibiotic trimethoprim. The crystal structure of R67 DHFR displays a toroidal structure with a central active-site pore. This homotetrameric protein exhibits 222 symmetry, with only a few residues from each chain contributing to the active site, so related sites must be used to bind both substrate (dihydrofolate) and cofactor (NADPH) in the productive R67 DHFR.NADPH.dihydrofolate complex. Whereas the site of folate binding has been partially resolved crystallographically, an interesting question remains: how can the highly symmetrical active site also bind and orient NADPH for catalysis? To model this ternary complex, we employed DOCK and SLIDE, two methods for docking flexible ligands into proteins using quite different algorithms. The bound pteridine ring of folate (Fol I) from the crystal structure of R67 DHFR was used as the basis for docking the nicotinamide-ribose-Pi (NMN) moiety of NADPH. NMN was positioned by both DOCK and SLIDE on the opposite side of the pore from Fol I, where it interacts with Fol I at the pore's center. Numerous residues serve dual roles in binding. For example, Gln 67 from both the B and D subunits has several contacts with the pteridine ring, while the same residue from the A and C subunits has several contacts with the nicotinamide ring. The residues involved in dual roles are generally amphipathic, allowing them to make both hydrophobic and hydrophilic contacts with the ligands. The result is a 'hot spot' binding surface allowing the same residues to co-optimize the binding of two ligands, and orient them for catalysis.

Factorising ligand affinity: a combined thermodynamic and crystallographic study of trypsin and thrombin inhibition

The binding of a series of low molecular weight ligands towards trypsin and thrombin has been studied by isothermal titration calorimetry and protein crystallography. In a series of congeneric ligands, surprising changes of protonation states occur and are overlaid on the binding process. They result from induced pK(a) shifts depending on the local environment experienced by the ligand and protein functional groups in the complex (induced dielectric fit). They involve additional heat effects that must be corrected before any conclusion on the binding enthalpy (DeltaH) and entropy (DeltaS) can be drawn. After correction, trends in both contributions can be interpreted in structural terms with respect to the hydrogen bond inventory or residual ligand motions. For all inhibitors studied, a strong negative heat capacity change (DeltaC(p)) is detected, thus binding becomes more exothermic and entropically less favourable with increasing temperature. Due to a mutual compensation, Gibbs free energy remains virtually unchanged. The strong negative DeltaC(p) value cannot solely be explained by the removal of hydrophobic surface portions of the protein or ligand from water exposure. Additional contributions must be considered, presumably arising from modulations of the local water structure, changes in vibrational modes or other ordering parameters. For thrombin, smaller negative DeltaC(p) values are observed for ligand binding in the presence of sodium ions compared to the other alkali ions, probably due to stabilising effects on the protein or changes in the bound water structure.

Functional equality in the absence of structural similarity: an added dimension to molecular mimicry

The crystal structure of meso-tetrasulfonatophenylporphyrin complexed with concanavalin A (ConA) was determined at 1.9 A resolution. Comparison of this structure with that of ConA bound to methyl alpha-d-mannopyranoside provided direct structural evidence of molecular mimicry in the context of ligand receptor binding. The sulfonatophenyl group of meso-tetrasulfonatophenylporphyrin occupies the same binding site on ConA as that of methyl alpha-d-mannopyranoside, a natural ligand. A pair of stacked porphyrin molecules stabilizes the crystal structure by end-to-end cross-linking with ConA resulting in a network similar to that observed upon agglutination of cells by lectins. The porphyrin binds to ConA predominantly through hydrogen bonds and water-mediated interactions. The sandwiched water molecules in the complex play a cementing role, facilitating favorable binding of porphyrin. Seven of the eight hydrogen bonds observed between methyl alpha-d-mannopyranoside and ConA are mimicked by the sulfonatophenyl group of porphyrin after incorporating two water molecules. Thus, the similarity in chemical interactions was manifested in terms of functional mimicry despite the obvious structural dissimilarity between the sugar and the porphyrin.

Thermodynamic aspects of hydrophobicity and biological QSAR

A protein contains a large amount of water molecules, and the nature of the interactions of the water molecules with a protein play an important role in the thermodynamics of the ligand binding process. In this paper, thermodynamic aspects of drug-receptor interactions, enthalpy-entropy compensation or reinforcement, hydrophobicity, and biological 2D- and 3D-QSAR are discussed. Comparisons of the thermodynamic QSAR of phenyl esters of N-benzoyl L-alanine in phosphate buffer and pentanol provide useful insight for the ligand-enzyme interactions.

The intracellular function of extracellular signaling peptides

A novel class of extracellular signaling peptides has been identified in Gram-positive bacteria that are actively transported into the cell to interact with intracellular receptors. The defining members of this novel class of signaling peptides are the Phr peptides of Bacillus subtilis and the mating pheromones of Enterococcus faecalis. These peptides are small and unmodified, gene encoded, and secreted by the bacterium. Most of these peptides diffuse into the extracellular medium, and when their concentration is sufficiently high, they are then actively transported into the cell by an oligopeptide permease (Opp). Once inside the cell, these peptides interact with an array of intracellular receptors. In B. subtilis, the Phr peptides regulate development of environmentally resistant spores and genetically competent cells (i.e. the natural ability to take up exogenous DNA). In E. faecalis, the mating pheromones regulate cell-cell transfer of plasmids, many of which encode antibiotic resistance or virulence factors. At least one component of the signaling pathway for these peptides is conserved in many bacteria, Opp. Opp is a non-specific transporter that transports peptides for use as carbon and nitrogen sources. The possibility that other bacteria could possess similar intracellularly functioning signaling peptides is discussed.

Influence of Glu-376 --> Gln mutation on enthalpy and heat capacity changes for the binding of slightly altered ligands to medium chain acyl-CoA dehydrogenase

We showed that the alpha-CH(2) --> NH substitution in octanoyl-CoA alters the ground and transition state energies for the binding of the CoA ligands to medium-chain acyl-CoA dehydrogenase (MCAD), and such an effect is caused by a small electrostatic difference between the ligands. To ascertain the extent that the electrostatic contribution of the ligand structure and/or the enzyme site environment modulates the thermodynamics of the enzyme-ligand interaction, we undertook comparative microcalorimetric studies for the binding of 2-azaoctanoyl-CoA (alpha-CH(2) --> NH substituted octanoyl-CoA) and octenoyl-CoA to the wild-type and Glu-376 --> Gln mutant enzymes. The experimental data revealed that both enthalpy (DeltaH degrees ) and heat capacity changes (DeltaC(p) degrees ) for the binding of 2-azaoctanoyl-CoA (DeltaH degrees (298) = -21.7 +/- 0.8 kcal/mole, DeltaC(p) degrees = -0.627 +/- 0.04 kcal/mole/K) to the wild-type MCAD were more negative than those obtained for the binding of octenoyl-CoA (DeltaH degrees (298) = -17.2 +/- 1.6 kcal/mole, DeltaC(p) degrees = -0.526 +/- 0.03 kcal/mole/K). Of these, the decrease in the magnitude of DeltaC(p) degrees for the binding of 2-azaoctanoyl-CoA (vis-à-vis octenoyl-CoA) to the enzyme was unexpected, because the former ligand could be envisaged to be more polar than the latter. To our further surprise, the ligand-dependent discrimination in the above parameters was completely abolished on Glu-376 --> Gln mutation of the enzyme. Both DeltaH degrees and DeltaC(p) degrees values for the binding of 2-azaoctanoyl-CoA (DeltaH degrees (298) = -13.3 +/- 0.6 kcal/mole, DeltaC(p) degrees = -0.511 +/- 0.03 kcal/mole/K) to the E376Q mutant enzyme were found to be correspondingly identical to those obtained for the binding of octenoyl-CoA (DeltaH degrees (298) = -13.2 +/- 0.6 kcal/mole, DeltaC(p) degrees = -0.520 +/- 0.02 kcal/mole/K). However, in neither case could the experimentally determined DeltaC(p) degrees values be predicted on the basis of the changes in the water accessible surface areas of the enzyme and ligand species. Arguments are presented that the origin of the above thermodynamic differences lies in solvent reorganization and water-mediated electrostatic interaction between ligands and enzyme site groups, and such interactions are intrinsic to the molecular basis of the enzyme-ligand complementarity.

Plasticity in protein-peptide recognition: crystal structures of two different peptides bound to concanavalin A

The structures of concanavalin A (ConA) in complex with two carbohydrate-mimicking peptides, 10-mer (MYWYPYASGS) and 15-mer (RVWYPYGSYLTASGS) have been determined at 2.75 A resolution. In both crystal structures four independent peptide molecules bind to each of the crystallographically independent subunits of ConA tetramer. The peptides exhibit small but significant variability in conformations and interactions while binding to ConA. The crystal structure of another similar peptide, 12-mer (DVFYPYPYASGS), in complex with ConA has been determined (Jain, D., K. J. Kaur, B. Sundaravadivel, and D. M. Salunke. 2000. Structural and functional consequences of peptide-carbohydrate mimicry. J. Biol. Chem. 275:16098-16102). Comparison of the three complexes shows that the peptides bind to ConA at a common binding site, using different contacting residues and interactions depending on their sequence and the local environment at the binding site. The binding is also optimized by corresponding plasticity of the peptide binding site on ConA. The diversity in conformation and interactions observed here are in agreement with the structural leeway concerning plasticity of specific molecular recognition in biological processes. The adaptability of peptide-ConA interactions may also be correlated with the carbohydrate-mimicking property of these peptides.

Docking of combinatorial peptide libraries into a broadly cross-reactive human IgM

A monoclonal IgM cryoglobulin with diverse binding behavior was isolated from a patient (Mez) with Waldenström's macroglobulinemia. It gave very high titers in the binding of combinatorially synthesized libraries of peptides ranging in size from two to eight residues. The crystal structure of Mez Fv revealed that the binding site was divided into two cavities of unequal volumes with dimensions and chemical properties that were compatible with the binding of peptides. Access to this unique combination of structural information and peptide binding data led us to carry out Mez-peptide docking simulations to gain insight into the Mez binding propensities. In the present article, the results for docking of five peptide libraries are combined with discussions of the methods and approximations involved in the docking process. We analyze the origins of peptide binding affinity for Mez IgM in terms of its cross-reactivity and its structural preferences.