Side-chain conformational changes upon Protein-Protein Association

Conformational changes upon protein-protein association are the key element of the binding mechanism. The study presents a systematic large-scale analysis of such conformational changes in the side chains. The results indicate that short and long side chains have different propensities for the conformational changes. Long side chains with three or more dihedral angles are often subject to large conformational transition. Shorter residues with one or two dihedral angles typically undergo local conformational changes not leading to a conformational transition. A relationship between the local readjustments and the equilibrium fluctuations of a side chain around its unbound conformation is suggested. Most of the side chains undergo larger changes in the dihedral angle most distant from the backbone. The frequencies of the core-to-surface interface transitions of six nonpolar residues and Tyr are larger than the frequencies of the opposite surface-to-core transitions. The binding increases both polar and nonpolar interface areas. However, the increase of the nonpolar area is larger for all considered classes of protein complexes, suggesting that the protein association perturbs the unbound interfaces to increase the hydrophobic contribution to the binding free energy. To test modeling approaches to side-chain flexibility in protein docking, conformational changes in the X-ray set were compared with those in the docking decoy sets. The results lead to a better understanding of the conformational changes in proteins and suggest directions for efficient conformational sampling in docking protocols.     

Structure of the conserved core of the yeast Dot1p, a nucleosomal histone H3 lysine 79 methyltransferase

Methylation of Lys79 on histone H3 by Dot1p is important for gene silencing. The elongated structure of the conserved core of yeast Dot1p contains an N-terminal helical domain and a seven-stranded catalytic domain that harbors the binding site for the methyl-donor and an active site pocket sided with conserved hydrophobic residues. The S-adenosyl-L-homocysteine exhibits an extended conformation distinct from the folded conformation observed in structures of SET domain histone lysine methyltransferases. A catalytic asparagine (Asn479), located at the bottom of the active site pocket, suggests a mechanism similar to that employed for amino methylation in DNA and protein glutamine methylation. The acidic, concave cleft between the two domains contains two basic residue binding pockets that could accommodate the outwardly protruding basic side chains around Lys79 of histone H3 on the disk-like nucleosome surface. Biochemical studies suggest that recombinant Dot1 proteins are active on recombinant nucleosomes, free of any modifications.     

Molecular Dynamic Simulations Reveal the Structural Determinants of Fatty Acid Binding to Oxy-Myoglobin

The mechanism(s) by which fatty acids are sequestered and transported in muscle have not been fully elucidated. A potential key player in this process is the protein myoglobin (Mb). Indeed, there is a catalogue of empirical evidence supporting direct interaction of globins with fatty acid metabolites; however, the binding pocket and regulation of the interaction remains to be established. In this study, we employed a computational strategy to elucidate the structural determinants of fatty acids (palmitic & oleic acid) binding to Mb. Sequence analysis and docking simulations with a horse (Equus caballus) structural Mb reference reveals a fatty acid-binding site in the hydrophobic cleft near the heme region in Mb. Both palmitic acid and oleic acid attain a “U” shaped structure similar to their conformation in pockets of other fatty acid-binding proteins. Specifically, we found that the carboxyl head group of palmitic acid coordinates with the amino group of Lys45, whereas the carboxyl group of oleic acid coordinates with both the amino groups of Lys45 and Lys63. The alkyl tails of both fatty acids are supported by surrounding hydrophobic residues Leu29, Leu32, Phe33, Phe43, Phe46, Val67, Val68 and Ile107. In the saturated palmitic acid, the hydrophobic tail moves freely and occasionally penetrates deeper inside the hydrophobic cleft, making additional contacts with Val28, Leu69, Leu72 and Ile111. Our simulations reveal a dynamic and stable binding pocket in which the oxygen molecule and heme group in Mb are required for additional hydrophobic interactions. Taken together, these findings support a mechanism in which Mb acts as a muscle transporter for fatty acid when it is in the oxygenated state and releases fatty acid when Mb converts to deoxygenated state.     

Dynamic structure based pharmacophore modeling of the Acetylcholinesterase reveals several potential inhibitors

Acetylcholinesterase is a critical enzyme that regulates neurotransmission by catalyzing the breakdown of neurotransmitter acetylcholine in synapses of the nervous system. It is an important target for therapeutic drugs that treat Alzheimer’s disease. Since, the degree of flexibility of the side chains of the residues in the active-site gorge of Acetylcholinesterase is diverse it results in different bound ligand conformations. The side-chain conformations of Ser293, Tyr341, Leu76, and Val73 are flexible, while the side-chain conformations of Tyr72, Tyr 124, Ser125, Phe295, and Arg296 appear to be fixed. In this study, multi-conformation dynamic pharmacophore models from the donepezyl-binding pocket based on highly populated structures chosen from molecular dynamics simulations were used for screening compounds that can properly bind acetylcholinesterase. Based on these structures, three pharmacophore models were generated. Consequently, 14 hits were retrieved as final candidates by utilizing virtual screening of ZINC database and molecular docking.     

NMR studies of 4-19F-phenylalanine-labeled intestinal fatty acid binding protein: evidence for conformational heterogeneity in the native state

(19)F-Nuclear magnetic resonance (NMR) studies have been carried out after incorporation of 4-(19)F-phenylalanine into the intestinal fatty acid binding protein (IFABP), a protein composed of two beta-sheets containing a large hydrophobic cavity into which ligands bind. NMR spectra have been obtained with both the ligand-free and ligand-bound (oleate) forms. There are 29 residues involved in van der Waals or hydrophobic interactions or both to form a U-shaped ligand binding pocket (Sacchettni J. C., Scapin G., Gopaul D., and Gordon J. I. (1992) J. Biol. Chem. 267, 23534-23545). The protein contains eight phenylalanines, and all are included in those residues that line the pocket. Peak assignments were made using site-specific incorporation of 4-(19)F-phenylalanine. Fluorine is a highly sensitive probe to monitor the conformation and dynamics of the side chains in native state. We find that chemical exchange in the binding pocket exists in the native apo- and holo-state. Of the eight phenylalanine residues, Phe2, Phe47, Phe62, Phe68, and Phe93 are arranged on one side of the binding pocket, and all exist in two conformations with Phe2, Phe47, and Phe62 showing exchange cross-peaks with minor conformation in (19)F-(19)F nuclear Overhauser effect (NOESY) spectra. The line widths of Phe68 and Phe93 are broader than those of other phenylalanine residues and can be deconvoluted into two peaks. Phe47, Phe62, Phe68, Phe93, and Trp82 have been proposed to be involved in the early stage of collapse (Ropson, I. J., and Frieden, C. (1992) Proc. Natl. Acad. Sci U.S.A. 89, 7222-7226), but a temperature study suggests that Phe47 behaves differently than other residues and may be more involved in a later stage of folding, for example, side chain stabilization. In the holo-form, Phe17 shows an extra exchange cross-peak in addition to those exchange cross-peaks observed in apo-form. Holo-IFABP exhibits broader line width than the apo-form, suggesting more flexibility of the binding cavity upon ligand binding.     

An activation switch in the ligand binding pocket of the C5a receptor

Although agonists are thought to occupy binding pockets within the seven-helix core of serpentine receptors, the topography of these binding pockets and the conformational changes responsible for receptor activation are poorly understood. To identify the ligand binding pocket in the receptor for complement factor 5a (C5aR), we assessed binding affinities of hexapeptide ligands, each mutated at a single position, for seven mutant C5aRs, each mutated at a single position in the putative ligand binding site. In ChaW (an antagonist) and W5Cha (an agonist), the side chains at position 5 are tryptophan and cyclohexylalanine, respectively. Comparisons of binding affinities indicated that the hexapeptide residue at this position interacts with two C5aR residues, Ile-116 (helix III) and Val-286 (helix VII); in a C5aR model these two side chains point toward one another. Both the I116A and the V286A mutations markedly increased binding affinity of W5Cha but not that of ChaW. Moreover, ChaW, the antagonist hexapeptide, acted as a full agonist on the I116A mutant. These results argue that C5aR residues Ile-116 and Val-286 interact with the side chain at position 5 of the hexapeptide ligand to form an activation switch. Based on this and previous work, we present a docking model for the hexapeptide within the C5aR binding pocket. We propose that agonists induce a small change in the relative orientations of helices III and VII and that these helices work together to allow movement of helix VI away from the receptor core, thereby triggering G protein activation.     

The impact of protein flexibility on protein-protein docking

Accounting for protein flexibility in protein-protein docking algorithms is challenging, and most algorithms therefore treat proteins as rigid bodies or permit side-chain motion only. While the consequences are obvious when there are large conformational changes upon binding, the situation is less clear for the modest conformational changes that occur upon formation of most protein-protein complexes. We have therefore studied the impact of local protein flexibility on protein-protein association by means of rigid body and torsion angle dynamics simulation. The binding of barnase and barstar was chosen as a model system for this study, because the complexation of these 2 proteins is well-characterized experimentally, and the conformational changes accompanying binding are modest. On the side-chain level, we show that the orientation of particular residues at the interface (so-called hotspot residues) have a crucial influence on the way contacts are established during docking from short protein separations of approximately 5 A. However, side-chain torsion angle dynamics simulations did not result in satisfactory docking of the proteins when using the unbound protein structures. This can be explained by our observations that, on the backbone level, even small (2 A) local loop deformations affect the dynamics of contact formation upon docking. Complementary shape-based docking calculations confirm this result, which indicates that both side-chain and backbone levels of flexibility influence short-range protein-protein association and should be treated simultaneously for atomic-detail computational docking of proteins.     

Side-chain flexibility in protein-ligand binding: the minimal rotation hypothesis

The goal of this work is to learn from nature about the magnitudes of side-chain motions that occur when proteins bind small organic molecules, and model these motions to improve the prediction of protein-ligand complexes. Following analysis of protein side-chain motions upon ligand binding in 63 complexes, we tested the ability of the docking tool SLIDE to model these motions without being restricted to rotameric transitions or deciding which side chains should be considered as flexible. The model tested is that side-chain conformational changes involving more atoms or larger rotations are likely to be more costly and less prevalent than small motions due to energy barriers between rotamers and the potential of large motions to cause new steric clashes. Accordingly, SLIDE adjusts the protein and ligand side groups as little as necessary to achieve steric complementarity. We tested the hypothesis that small motions are sufficient to achieve good dockings using 63 ligands and the apo structures of 20 different proteins and compared SLIDE side-chain rotations to those experimentally observed. None of these proteins undergoes major main-chain conformational change upon ligand binding, ensuring that side-chain flexibility modeling is not required to compensate for main-chain motions. Although more frugal in the number of side-chain rotations performed, this model substantially mimics the experimentally observed motions. Most side chains do not shift to a new rotamer, and small motions are both necessary and sufficient to predict the correct binding orientation and most protein-ligand interactions for the 20 proteins analyzed.     

Flexibility of metal binding sites in proteins on a database scale

Protein metal binding sites in the pre-bound (apo) state, and their rearrangements upon metal binding were not analyzed previously at a database scale. Such a study may provide valuable information for metal binding site prediction and design. A high resolution, nonredundant dataset of 210 metal binding sites was created, containing all available representatives of apo-holo pairs for the most populated metals in the PDB. More than 40% of the sites underwent rearrangements upon metal binding. In 30 cases rearrangements involved the backbone. The tendency for side-chain rearrangement inversely correlates with the number of first-shell residues. Analysis of side-chain reorientations as a result of metal binding showed that in 95% of the rigid-backbone binding sites at most one side chain moved. Thus, in general, part of the first coordination shell is already in place in the pre-bound form. The frequencies of side-chain reorientation directly correlated with metal ligand flexibility and solvent accessibility in the apo state.     

Monitoring aromatic picosecond to nanosecond dynamics in proteins via 13C relaxation: expanding perturbation mapping of the rigidifying core mutation, V54A, in eglin c

Long-range effects, such as allostery, have evolved in proteins as a means of regulating function via communication between distal sites. An NMR-based perturbation mapping approach was used to more completely probe the dynamic response of the core mutation V54A in the protein eglin c by monitoring changes in picosecond to nanosecond aromatic side-chain dynamics and H/D exchange stabilities. Previous side-chain dynamics studies on this mutant were limited to methyl-bearing residues, most of which were found to rigidify on the picosecond to nanosecond time scale in the form of a contiguous "network". Here, high precision (13)C relaxation data from 13 aromatic side chains were acquired by applying canonical relaxation experiments to a newly developed carbon labeling scheme [Teilum et al. (2006) J. Am. Chem. Soc. 128, 2506-2507]. The fitting of model-free parameters yielded S (2) variability which is intermediate with respect to backbone and methyl-bearing side-chain variability and tau e values that are approximately 1 ns. Inclusion of the aromatic dynamic response results in an expanded network of dynamically coupled residues, with some aromatics showing increases in flexibility, which partially offsets the rigidification in methyl side chains. Using amide hydrogen exchange, dynamic propagation on a slower time scale was probed in response to the V54A perturbation. Surprisingly, regional stabilization (slowed exchange) 10-12 A from the site of mutation was observed despite a global destabilization of 1.5 kcal x mol (-1). Furthermore, this unlikely pocket of stabilized residues colocalizes with increases in aromatic flexibility on the faster time scale. Because the converse is also true (destabilized residues colocalize with rigidification on the fast time scale), a plausible entropy-driven mechanism is discussed for relating colocalization of opposing dynamic trends on vastly different time scales.     

Cross Talk between KGF and KITLG Proteins Implicated with Ovarian Folliculogenesis in Buffalo Bubalus bubalis

Molecular interactions between mesenchymal-derived Keratinocyte growth factor (KGF) and Kit ligand (KITLG) are essential for follicular development. These factors are expressed by theca and granulosa cells. We determined full length coding sequence of buffalo KGF and KITLG proteins having 194 and 274 amino acids, respectively. The recombinant KGF and KITLG proteins were solubilized in 10 mM Tris, pH 7.5 and 50 mM Tris, pH 7.4 and purified using Ni-NTA column and GST affinity chromatography, respectively. The purity and molecular weight of His-KGF (~23 kDa) and GST-KITLG (~57 kDa) proteins were confirmed by SDS-PAGE and western blotting. The co-immunoprecipitation assay accompanied with computational analysis demonstrated the interaction between KGF and KITLG proteins. We deduced 3D structures of the candidate proteins and assessed their binding based on protein docking. In the process, KGF specific residues, Lys123, Glu135, Lys140, Lys155 and Trp156 and KITLG specific ones, Ser226, Phe233, Gly234, Ala235, Phe236, Trp238 and Lys239 involved in the formation of KGF-KITLG complex were detected. The hydrophobic interactions surrounding KGF-KITLG complex affirmed their binding affinity and stability to the interacting interface. Additionally, in-silico site directed mutagenesis enabled the assessment of changes that occurred in the binding energies of mutated KGF-KITLG protein complex. Our results demonstrate that in the presence of KITLG, KGF mimics its native binding mode suggesting all the KGF residues are specific to their binding complex. This study provides an insight on the critical amino acid residues participating in buffalo ovarian folliculogenesis.     

Structural analysis of protein‐ligand interactions: the binding of endogenous compounds and of synthetic drugs

The large number of macromolecular structures deposited with the Protein Data Bank (PDB) describing complexes between proteins and either physiological compounds or synthetic drugs made it possible a systematic analysis of the interactions occurring between proteins and their ligands. In this work, the binding pockets of about 4000 PDB protein‐ligand complexes were investigated and amino acid and interaction types were analyzed. The residues observed with lowest frequency in protein sequences, Trp, His, Met, Tyr, and Phe, turned out to be the most abundant in binding pockets. Significant differences between drug‐like and physiological compounds were found. On average, physiological compounds establish with respect to drugs about twice as many hydrogen bonds with protein atoms, whereas drugs rely more on hydrophobic interactions to establish target selectivity. The large number of PDB structures describing homologous proteins in complex with the same ligand made it possible to analyze the conservation of binding pocket residues among homologous protein structures bound to the same ligand, showing that Gly, Glu, Arg, Asp, His, and Thr are more conserved than other amino acids. Also in the cases in which the same ligand is bound to unrelated proteins, the binding pockets showed significant conservation in the residue types. In this case, the probability of co‐occurrence of the same amino acid type in the binding pockets could be up to thirteen times higher than that expected on a random basis. The trends identified in this study may provide an useful guideline in the process of drug design and lead optimization. Copyright © 2014 John Wiley & Sons, Ltd.     

Site-Directed Mutagenesis of Cytochrome P450scc (CYP11A1). Effect of Lysine Residue Substitution on Its Structural and Functional Properties

Our previous chemical modification and cross-linking studies identified some positively charged amino acid residues of cytochrome P450scc that may be important for its interaction with adrenodoxin and for its functional activity. The present study was undertaken to further evaluate the role of these residues in the interaction of cytochrome P450scc with adrenodoxin using site-directed mutagenesis. Six cytochrome P450scc mutants containing replacements of the surface-exposed positively charged residues (Lys103Gln, Lys110Gln, Lys145Gln, Lys394Gln, Lys403Gln, and Lys405Gln) were expressed in E. coli cells, purified as a substrate-bound high-spin form, and characterized as compared to the wild-type protein. The replacement of the surface Lys residues does not dramatically change the protein folding or the heme pocket environment as judged from limited proteolysis and spectral studies of the cytochrome P450 mutants. The replacement of Lys in the N-terminal sequence of P450scc does not dramatically affect the activity of the heme protein. However, mutant Lys405Gln revealed rather dramatic loss of cholesterol side-chain cleavage activity, efficiency of enzymatic reduction in a reconstituted system, and apparent dissociation constant for adrenodoxin binding. The present results, together with previous findings, suggest that the changes in functional activity of mutant Lys405Gln may reflect the direct participation of this amino acid residue in the electrostatic interaction of cytochrome P450scc with its physiological partner, adrenodoxin.     

Dynamics of sialyl Lewis(a) in aqueous solution and prediction of the structure of the sialyl Lewis(a)-SelectinE complex

In this article we investigate all possible three-dimensional structures for sialyl Lewis^a (SLe^a) in aqueous solution and we predict without a priori experimental information its conformation when bound to SelectinE by using a combination of long molecular dynamics (MD) simulations. Based on 10 ns MD studies, three structures differing in glycosidic conformations are proposed for SLe^a in aqueous solution. Based on a 4 ns MD study of the SLe^a-SelectinE complex with initial structures derived from our prediction tools, we find that, fucose and N-acetyl neuraminic acid are in close contact with SelectinE and therefore expect interactions of the protein with these two sugar rings to be significantly more important than in the case of galactose and N-acetyl glucosamine. Our predictions indicate that the N-acetyl glucosamine of SLe^a is positioned primarily in the aqueous phase. In order to be able to interact with SLe^a the side chains of amino acid residues Lys99 and Lys111 in SelectinE appear to undergo large conformational changes when contrasted with the positions of these residues in the X-ray crystal structure. Furthermore, amino acid residues Arg97, Glu98 and Lys99 are acting as a holding arm to position the NeuNAc of SLe^a in the binding pocket.     

The membrane-associated lipoprotein-9 GmpC from Staphylococcus aureus binds the dipeptide GlyMet via side chain interactions

Bacterial dipeptide ABC transporters function to import a wide range of dipeptide substrates. This ability to transport a wide variety of dipeptides is conferred by the cognate substrate binding protein (SBP) of these transporters. SBPs bind dipeptides with little regard for their amino acid content. Here, we report the 1.7 A resolution structure of lipoprotein-9 (SA0422) of Staphylococcus aureus in complex with the dipeptide glycylmethionine. Experimental characterization of the subcellular location of the protein confirmed that SA0422 is an acylated, peripheral membrane protein. This is the first structure determined for an SBP of a Gram-positive dipeptide ABC transporter. Usually, binding of dipeptides occurs in a binding pocket that is largely hydrated and able to accommodate the side chains of several different amino acid residues. Unlike any other known SBP, lipoprotein-9 binds the side chains of the glycylmethionine dipeptide through very specific interactions. Lipoprotein-9 shares significant structural and sequence homology with the MetQ family of methionine SBP. Sequence comparisons between MetQ-like proteins and lipoprotein-9 suggest that the residues forming the tight interactions with the methionine side chains of the ligand are highly conserved between lipoprotein-9 and MetQ homologues, while the residues involved in coordinating the glycine residue are not. Modeling of the Vibrio cholerae MetQ and lipoprotein-9 binding pockets can account for lipoprotein-9 substrate specificity toward glycylmethionine. For this reason, we have designated lipoprotein-9 GmpC, for glycylmethionine binding protein.     

Investigation of the stereoselectivity of an anti-amino acid antibody using molecular modeling and ligand docking

The structure of the binding site of the stereoselective anti-D-amino acid antibody 67.36 was modeled utilizing web antibody modeling (WAM) and SWISS-MODEL. Although docking experiments performed with an aromatic amino acid as model ligand were unsuccessful with the WAM structure, ligand binding was achieved with the SWISS-MODEL structure. Incorporation of side-chain flexibility within the binding site resulted in a protein structure that stereoselectively binds to the D-enantiomer of the model ligand. In addition to four hydrogen bonds that are formed between amino acid residues in the binding site and the ligand, a number of hydrophobic interactions are involved in the formation of the antibody-ligand complex. The aromatic side chain of the ligand interacts with a tryptophan and a tyrosine residue in the binding site through pi-pi stacking. Fluorescence spectroscopic investigations also suggest the presence of tryptophan residues in the binding site, as ligand binding causes an enhancement of the antibody's intrinsic fluorescence at an emission wavelength of 350 nm. Based on the modeled antibody structure, the L-enantiomer of the model ligand cannot access the binding site due to steric hindrance. Additional docking experiments performed with D-phenylalanine and D-norvaline showed that these ligands are bound to the antibody in a way analogous to the D-enantiomer of the model ligand.     

Rotamer libraries and probabilities of transition between rotamers for the side chains in protein-protein binding

Conformational changes in the side chains are essential for protein-protein binding. Rotameric states and unbound- to-bound conformational changes in the surface residues were systematically studied on a representative set of protein complexes. The side-chain conformations were mapped onto dihedral angles space. The variable threshold algorithm was developed to cluster the dihedral angle distributions and to derive rotamers, defined as the most probable conformation in a cluster. Six rotamer libraries were generated: full surface, surface noninterface, and surface interface-each for bound and unbound states. The libraries were used to calculate the probabilities of the rotamer transitions upon binding. The stability of amino acids was quantified based on the transition maps. The noninterface residues' stability was higher than that of the interface. Long side chains with three or four dihedral angles were less stable than the shorter ones. The transitions between the rotamers at the interface occurred more frequently than on the noninterface surface. Most side chains changed conformation within the same rotamer or moved to an adjacent rotamer. The highest percentage of the transitions was observed primarily between the two most occupied rotamers. The probability of the transition between rotamers increased with the decrease of the rotamer stability. The analysis revealed characteristics of the surface side-chain conformational transitions that can be utilized in flexible docking protocols.     

A mining minima approach to exploring the docking pathways of p-nitrocatechol sulfate to YopH

Using the docking of p-nitrocatechol sulfate to Yersinia protein tyrosine phosphatase YopH as an example, we showed that an approach based on mining minima followed by cluster and similarity analysis could generate useful insights into docking pathways. Our simulation treated both the ligand and the protein as flexible molecules so that the coupling between their motion could be properly accounted for. Our simulation identified three docking poses; the one with the lowest energy agreed well with experimental structure. The model also predicted the side-chain conformations of the amino acids lying in the binding pocket correctly with the exception of three residues that appeared to be stabilized by two structural water molecules in the crystal structure. The implicit solvent model employed in the simulation could not capture such effects well. We also found four major pathways leading to these docking poses after the ligand entered the mouth of the binding pocket. In addition, the sulfate group of p-nitrocatechol sulfate was found to be important both in binding the ligand to the pocket and in guiding the ligand to dock into the pocket. The coupling of the motion between the protein and the ligand also played an important role in facilitating ligand loading and unloading.     

Reduced-bond tight-binding inhibitors of HIV-1 protease. Fine tuning of the enzyme subsite specificity.

Truncation of a peptide substrate in the N-terminus and replacement of its scissile amide bond with a non-cleavable reduced bond results in a potent inhibitor of HIV-1 protease. A series of such inhibitors has been synthesized, and S2-S3' subsites of the protease binding cleft mapped. The S2 pocket requires bulky Boc or PIV groups, large aromatic Phe residues are preferred in P1 and P1' and Glu in P2'. The S3' pocket prefers Phe over small Ala or Val. Introduction of a Glu residue into the P2' position yields a tight-binding inhibitor of HIV-1 protease, Boc-Phe-[CH2-NH]-Phe-Glu-Phe-OMe, with a subnanomolar inhibition constant. The relevant peptide derived from the same amino acid sequence binds to the protease with a Ki of 110 nM, thus still demonstrating a good fit of the amino acid residues into the protease binding pockets and also the importance of the flexibility of P1-P1' linkage for proper binding. A new type of peptide bond mimetic, N-hydroxylamine -CH2-N(OH)-, has been synthesized. Binding of hydroxylamino inhibitor of HIV-1 protease is further improved with respect to reduced-bond inhibitor.     

Binding of amino acid side chains to preformed cavities: interaction of serine proteinases with turkey ovomucoid third domains with coded and noncoded P1 residues.

In the association of serine proteinases with their cognate substrates and inhibitors an important interaction is the fitting of the P1 side chain of the substrate or inhibitor into a preformed cavity of the enzyme called the S1 pocket. In turkey ovomucoid third domain, which is a canonical protein proteinase inhibitor, the P1 residue is Leu18. Here we report the values of equilibrium constants, Ka, for turkey ovomucoid third domain and 13 additional Leu18X variants with six serine proteinases: bovine alpha chymotrypsin A, porcine pancreatic elastase, subtilisin Carlsberg, Streptomyces griseus proteinases A and B, and human leukocyte elastase. Eight of the Xs are coded amino acids: Ala, Ser, Val, Met, Gln, Glu, Lys, and Phe, and five are noncoded: Abu, Ape, Ahx, Ahp, and Hse. They were chosen to simplify the interamino acid comparisons. In the homologous series of straight-chain side chains Ala, Abu, Ape, Ahx, Ahp, free energy of binding decreases monotonically with the side-chain length for chymotrypsin with large binding pocket, but even for this enzyme shows curvature. For the two S. griseus enzymes a minimum appears to be reached at Ahp. A minimum is clearly evident for the two elastases, where increasing the side-chain length from Ahx to Ahp greatly weakens binding, but much more so for the apparently more rigid pancreatic enzyme than for the more flexible leukocyte enzyme. beta-Branching (Ape/Val) is very deleterious for five of the six enzymes; it is only slightly deleterious for the more flexible human leukocyte elastase. The effect of gamma-branching (Ahx/Leu), of introduction of heteroatoms (Abu/Ser), (Ape/Hse), and (Ahx/Met), and of introduction of charge (Gln/Glu) and (Ahp/Lys) are tabulated and discussed. An important component of the free energy of interaction is the distortion of the binding pocket by bulky or branched side chains. Most of the variants studied were obtained by enzymatic semisynthesis. X18 variants of the 6-18 peptide GlyNH2 were synthesized and combined with natural reduced peptide 19-56. Disulfide bridges were formed. The GlyNH2 was removed and the reactive-site peptide bond X18-Glu19 was synthesized by complex formation with proteinase K. The resultant complexes were dissociated by sudden pH drop. This kinetically controlled dissociation afforded virgin, reactive-site-intact inhibitor variants.