Interaction of the factor XIII activation peptide with alpha -thrombin. Crystal structure of its enzyme-substrate analog complex

The serine protease thrombin proteolytically activates blood coagulation factor XIII by cleavage at residue Arg(37); factor XIII in turn cross-links fibrin molecules and gives mechanical stability to the blood clot. The 2.0-A resolution x-ray crystal structure of human alpha-thrombin bound to the factor XIII-(28-37) decapeptide has been determined. This structure reveals the detailed atomic level interactions between the factor XIII activation peptide and thrombin and provides the first high resolution view of this functionally important part of the factor XIII molecule. A comparison of this structure with the crystal structure of fibrinopeptide A complexed with thrombin highlights several important determinants of thrombin substrate interaction. First, the P1 and P2 residues must be compatible with the geometry and chemistry of the S1 and S2 specificity sites in thrombin. Second, a glycine in the P5 position is necessary for the conserved substrate conformation seen in both factor XIII-(28-37) and fibrinopeptide A. Finally, the hydrophobic residues, which occupy the aryl binding site of thrombin determine the substrate conformation further away from the catalytic residues. In the case of factor XIII-(28-37), the aryl binding site is shared by hydrophobic residues P4 (Val(34)) and P9 (Val(29)). A bulkier residue in either of these sites may alter the substrate peptide conformation.     

Docking and binding mode analysis of aryl diketoacids (ADK) at the active site of HCV RNA-dependent RNA polymerase

The pharmacophore-guided docking study of aryl diketoacid (ADK) analogues revealed two distinctive hydrophobic binding sites (a pocket and a groove) around the UTP-binding site of hepatitis C virus (HCV) RNA-dependent RNA polymerase (RdRp). Interestingly, the hydrophobic binding sites have appropriate shape and size to specifically substituted aromatic rings, which suggests the specific role of substituents on the aromatic ring in determining the binding affinity of the ADK analogue to the active site of the target enzyme. Binding mode analysis of ADK analogues with potent antiviral activity shows highly substituted aromatic rings map well onto the hydrophobic binding sites. For less active compounds, their lack of aromatic substitution and thereby insufficient size can be primarily ascribed to their inability to bind to the hydrophobic binding site. The characteristic binding mode of ADK analogues proposed in this study provides a useful tool in designing a structure–activity relationship study of novel ADK analogues based on various aromatic substituents.     

The refined 1.9 A crystal structure of human alpha-thrombin: interaction with D-Phe-Pro-Arg chloromethylketone and significance of the Tyr-Pro-Pro-Trp insertion segment.

A stoichiometric complex formed between human alpha-thrombin and D-Phe-Pro-Arg chloromethylketone was crystallized in an orthorhombic crystal form. Orientation and position of a starting model derived from homologous modelling were determined by Patterson search methods. The thrombin model was completed in a cyclic modelling-crystallographic refinement procedure to a final R-value of 0.171 for X-ray data to 1.92 A. The structure is in full agreement with published cDNA sequence data. The A-chain, ordered only in its central part, is positioned along the molecular surface opposite to the active site. The B-chain exhibits the characteristic polypeptide fold of trypsin-like proteinases. Several extended insertions form, however, large protuberances; most important for interaction with macromolecular substrates is the characteristic thrombin loop around Tyr60A-Pro60B-Pro60C-Trp60D (chymotrypsinogen numbering) and the enlarged loop around the unique Trp148. The former considerably restricts the active site cleft and seems likely to be responsible for poor binding of most natural proteinase inhibitors to thrombin. The exceptional specificity of D-Phe-Pro-Arg chloromethylketone can be explained by a hydrophobic cage formed by Ile174, Trp215, Leu99, His57, Tyr60A and Trp60D. The narrow active site cleft, with a more polar base and hydrophobic rims, extends towards the arginine-rich surface of loop Lys70-Glu80 that probably represents part of the anionic binding region for hirudin and fibrinogen.     

Beta-ketophosphonates as beta-lactamase inhibitors: Intramolecular cooperativity between the hydrophobic subsites of a class D beta-lactamase

A series of aryl and arylmethyl beta-aryl-beta-ketophosphonates have been prepared as potential beta-lactamase inhibitors. These compounds, as fast, reversible, competitive inhibitors, were most effective (micromolar K(i) values) against the class D OXA-1 beta-lactamase but had less activity against the OXA-10 enzyme. They were also quite effective against the class C beta-lactamase of Enterobacter cloacae P99 but less so against the class A TEM-2 enzyme. Reduction of the keto group to form the corresponding beta-hydroxyphosphonates led to reduced inhibitory activity. Molecular modeling, based on the OXA-1 crystal structure, suggested interaction of the aryl groups with the hydrophobic elements of the enzyme's active site and polar interaction of the keto and phosphonate groups with the active site residues Ser 115, Lys 212 and Thr 213 and with the non-conserved Ser 258. Analysis of binding free energies showed that the beta-aryl and phosphonate ester aryl groups interacted cooperatively within the OXA-1 active site. Overall, the results suggest that quite effective inhibitors of class C and some class D beta-lactamases could be designed, based on the beta-ketophosphonate platform.     

Asymmetric sulfoxidations mediated by alpha-chymotrypsin

The oxidation of aryl alkyl sulfides with H(2)O(2) in aqueous solution is a reasonably facile reaction producing racemic sulfoxides. We show that in the presence of the hydrolytic enzyme alpha-chymotrypsin such a sulfoxidation is accelerated and, more importantly, becomes stereoselective. With phenyl isobutyl sulfide as a model, the chymotrypsin-mediated, highly asymmetric oxidation is shown to occur in the hydrophobic binding pocket of the enzyme active site. The stereoselectivity of the chymotrypsin-mediated sulfoxidations is correctly explained by means of structure-based molecular modeling of the enzyme-sulfide complexes.     

Isoprenoid flavonoids are new leads in the modulation of chemoresistance

Flavonoid compounds are able to bind to P-glycoprotein (P-gp), a transporter involved in chemoresistance of cancer cells. This interaction involves, at least in part, two overlapping sites in the cytosolic domains of P-gp, the ATP site and a hydrophobic steroid-binding site. We have studied the structure-activity relationships toward binding to P-gp. Modification of the substitution pattern of the flavonoid ring by hydroxylation, methoxylation or introduction of nitrogen-containing substituents had little effect. On the contrary, the presence of a 3-hydroxyl (flavonols), and especially of a C-isoprenoid chain increased the affinity of flavonoids towards P-gp. More detailed examination of the interaction with the ATP site was conducted through inhibition by flavonoids of the photolabeling by radioactive 8-azido-TNP-ATP. Only simple flavonols were demonstrated to bind to the ATP site. When position 3 was free (flavones) or when a hydrophobic C-prenyl substituent was present, interaction was rather directed to the hydrophobic site. A number of flavonoid compounds were tested for their ability to modulate multidrug resistance in resistant leukemic K562/R7 cells. Again, prenyl flavonoids were potent modulators. Simple flavonoids were ineffective in this model. The beneficial effect of prenylation was lower in polyhydroxylated compounds, suggesting a crucial role of hydrophobicity in P-gp modulation. This revised version was published online in June 2006 with corrections to the Cover Date.     

Three-dimensional structure prediction of the NAD binding site of proton-pumping transhydrogenase from Escherichia coli

A three-dimensional structure of the NAD site of Escerichia coli transhydrogenase has been predicted. The model is based on analysis of conserved residues among the transhydrogenases from five different sources, homologies with enzymes using NAD as cofactors or substrates, hydrophilicity profiles, and secondary structure predictions. The present model supports the hypothesis that there is one binding site, located relatively close to the N-terminus of the α-subunit. The proposed structure spans residues α145 to α287, and it includes five β-strands and five α-helices oriented in a typical open twisted α/β conformation. The amino acid sequence following the GXGXXG dinucleotide binding consensus sequence (residues α172 to α177) correlates exactly to a typical fingerprint region for ADP binding βαβ folds in dinucleotide binding enzymes. In the model, aspartic acid α195 forms hydrogen bonds to one or both hydroxyl groups on the adenosine ribose sugar moiety. Threonine α196 and alanine α256, located at the end of βB and βD, respectively, create a hydrophobic sandwich with the adenine part of NAD buried inside. The nicotinamide part is located in a hydrophobic cleft between αA and βE. Mutagenesis work has been carried out in order to test the predicted model and to determine whether residues within this domain are important for proton pumping directly. All data support the predicted structure, and no residue crucial for proton pumping Was detected. Since no three-dimensional structure of transhydrogenase has been solved, a well based tertiary structure prediction is of great value for further experimental design in trying to elucidate the mechanism of the energy-linked proton pump. © 1995 Wiley-Liss, Inc.     

Investigation into new anticonvulsant derivatives of alpha-substituted N-benzylamides of gamma-hydroxy- and gamma-acetoxybutyric acid. Part 5: search for new anticonvulsant compounds

A series of four N-benzylamides of gamma-hydroxybutyric acid (GHB), that contain N-(4-phenylpiperazine)-, N-(4-benzylpiperazine)rings, N-benzylamino-, or N-(2-phenylethylamine)-groups in the alpha-position of GHB were selected as model compounds, for determining the structural elements responsible for their potential anticonvulsant action. Based on the results of pharmacological, physicochemical, and molecular modelling investigations, the pharmacophore model for anticonvulsant N-substituted amides of GHB was defined. In this model, the presence of the N-benzylamide fragment is essential for activity. In addition, all of the amides contained another hydrophobic unit (aryl ring) as a distal binding site and H-bond donor. In consideration of these model parameters, a number of N-substituted amides of GHB, containing a hydrophobic moiety such as: N-benzylamino or N-(4-chlorobenzylamino) group in the alpha-position of GHB, and a lipophilic substituent in the amide portion, were prepared. It has been shown that the anticonvulsant activities of the newly synthesized compounds might partially be explained on the basis of their lipophilicity (calculated log P values) and the presence of a hydroxyl group in the molecule.     

Specificity of Prodan for the self-associating domain of spectrin: a molecular docking study

The hydrophobic fluorescent probe Prodan binds to the self-associating domain of spectrin with 1:1 stoichiometry. A model of the self-associating domain was generated based on its homology with other domains of spectrin. Prodan was then docked onto the model, and several sites with low interaction energy were identified. To verify whether the binding of Prodan is specific towards the self-associating domain of spectrin, it was docked on to several other domains of spectrin, having a known three-dimensional structure. Analysis of the docking results suggests that the binding of Prodan to the self-associating domain of spectrin will involve hydrophobic and hydrophilic groups of Prodan. The results clearly indicate the preference of Prodan for a particular binding site of the self-associating domain.     

4-Aminoethylpiperazinyl aryl ketones with 5-HT₁A/5-HT₇ selectivity

The well-known 5-HT(1A)/5-HT(7) selectivity issue was tackled by a new series of 4-aminoethylpiperazinyl aryl ketones (1a-1l) specifically designed to distinguish the two hydrophobic sites centered at the anchoring salt bridge. The 4-aminoethylpiperazinyl aryl ketones showed a wide spectrum of activity and selectivity for the 5-HT receptors depending on the type of the hydrophobic groups attached at the aryl piperazinyl ketone scaffold. Docking study of the most active compounds against 5-HT(7)R and 5-HT(1A)R revealed that both receptors have two hydrophobic pockets around the anchoring salt bridge. These two binding sites are perpendicular to each other in 5-HT(7)R but parallel in 5-HT(1A)R, and this observation is well matched with the previous report which claimed that 5-HT(7)R affinity arises from bent conformation of the bound ligand whereas an extended one is best suited for 5-HT(1A)R selectivity. Also, as these pockets have different size and shape, inhibitory activity as well as selectivity of the 4-aminoethylpiperazinyl aryl ketones against 5-HT(7)R and 5-HT(1A)R seemed to be determined by combination of two hydrophobic substituents attached at both ends of the title compounds.     

Specificity of Prodan for the Self-associating Domain of Spectrin: A Molecular Docking Study

Abstract

The hydrophobic fluorescent probe Prodan binds to the self-associating domain of spectrin with 1:1 stoichiometry. A model of the self-associating domain was generated based on its homology with other domains of spectrin. Prodan was then docked onto the model, and several sites with low interaction energy were identified. To verify whether the binding of Prodan is specific towards the self-associating domain of spectrin, it was docked on to several other domains of spectrin, having a known three-dimensional structure. Analysis of the docking results suggests that the binding of Prodan to the self-associating domain of spectrin will involve hydrophobic and hydrophilic groups of Prodan. The results clearly indicate the preference of Prodan for a particular binding site of the self-associating domain.     

Crystal structure of the intrinsically flexible addiction antidote MazE

A specific camel VHH (variable domain of dromedary heavy chain antibody) fragment was used to crystallize the intrinsically flexible addiction antidote MazE. Only 45% of the polypeptide chain is found ordered in the crystal. The MazE monomer consisting of two beta-hairpins connected by a short alpha-helix has no hydrophobic core on its own and represents only one half of a typical protein domain. A complete domain structure is formed by the association of two chains, creating a hydrophobic core between two four-stranded beta-sheets. This hydrophobic core consists exclusively of short aliphatic residues. The folded part of MazE contains a novel DNA binding motif. A model for DNA binding that is consistent with the available biochemical data is presented.     

Non-thiol farnesyltransferase inhibitors: utilization of an aryl binding site by 5-arylacryloylaminobenzophenones

We recently described a novel aryl binding site of farnesyltransferase. The 2-naphthylacryloyl residue was developed as an appropriate substituent for our benzophenone-based AAX-peptidomimetic capable of occupying this binding site, resulting in a non-thiol farnesyltransferase inhibitor with nanomolar activity. The activity of this inhibitor is readily explained on the basis of docking studies which show the 2-naphthyl residue fitting into the aryl binding site.     

Hydrophilic residues at the apical domain of GroEL contribute to GroES binding but attenuate polypeptide binding

The GroES binding site at the apical domain of GroEL, mostly consisting of hydrophobic residues, overlaps largely with the substrate polypeptide binding site. Essential contribution of hydrophobic interaction to the binding of both GroES and polypeptide was exemplified by the mutant GroEL(L237Q) which lost the ability to bind either of them. The binding site, however, contains three hydrophilic residues, E238, T261, and N265. For GroES binding, N265 is essential since GroEL(N265A) is unable to bind GroES. E238 contributes to rapid GroES binding to GroEL because GroEL(E238A) is extremely sluggish in GroES binding. Polypeptide binding was not impaired by any mutations of E238A, T261A, and N265A. Rather, these mutants, especially GroEL(N265A), showed stronger polypeptide binding affinity than wild-type GroEL. Thus, these hydrophilic residues have a dual role; they help GroES binding on one hand but attenuate polypeptide binding on the other hand.     

Glycine-Rich Loop of Mitochondrial Processing Peptidase α-Subunit Is Responsible for Substrate Recognition by a Mechanism Analogous to Mitochondrial Receptor Tom20

Tryptophan fluorescence measurements were used to characterize the local dynamics of the highly conserved glycine-rich loop (GRL) of the mitochondrial processing peptidase (MPP) α-subunit in the presence of the substrate precursor. Reporter tryptophan residue was introduced into the GRL of the yeast α-MPP (Y299W) or at a proximal site (Y303W). Time-resolved and steady-state fluorescence spectroscopy demonstrated that for Trp299, the primary contact with the yeast malate dehydrogenase precursor evokes a change of the local GRL mobility. Moreover, time-resolved measurements showed that a functionless α-MPP with a single-residue deletion in the loop (Y303W/ΔG292) is defective particularly in the primary contact with substrate. Thus, the GRL was proved to be part of a contact site of the enzyme specifically recognizing the substrate. Regarding the surface exposure and presence of the hydrophobic patches within the GRL, we proposed a functional analogy between the presequence recognition by the hydrophobic binding groove of the Tom20 mitochondrial import receptor and the GRL of the α-MPP. A molecular dynamics (MD) simulation of the MPP-substrate peptide complex model was employed to test this hypothesis. The initial positioning and conformation of the substrate peptide in the model fitting were chosen based on the analogy of its interaction with the Tom20 binding groove. MD simulation confirmed the stability of the proposed interaction and showed also a decrease in GRL flexibility in the presence of substrate, in agreement with fluorescence measurements. Moreover, conserved substrate hydrophobic residues in positions + 1 and − 4 to the cleavage site remain in close contact with the side chains of the GRL during the entire production part of MD simulation as stabilizing points of the hydrophobic interaction. We conclude that the GRL of the MPP α-subunit is the crucial evolutional outcome of the presequence recognition by MPP and represents a functional parallel with Tom20 import receptor.     

A molecular model for the bilayer helices of the acetylcholine receptor including an acetylcholine binding site

A molecular model for the bilayer helices of the acetylcholine receptor is constructed from the 7 channel elements and the 17 hydrophobic helices of the 5 protein subunits. The acetylcholine binding site and the opening to the ion channel are included.     

Docking and molecular dynamics study on the inhibitory activity of <Emphasis Type="Italic">N,N</Emphasis>-disubstituted-trifluoro-3-amino-2-propanols-based inhibitors of cholesteryl ester transfer protein

Extensive studies suggest direct links between cholesteryl ester transfer protein (CETP), high-density lipoproteins-cholesterol level and cardiovascular diseases. Many therapeutic approaches are aimed at the CETP. A series of N, N-disubstituted-trifluoro-3-amino-2-propanol analogues are among the most highly potent and selective inhibitors of CETP described to date. For in-depth investigation into the structural and chemical features responsible for exploring the binding pocket of these compounds, as well as for the binding recognition mechanism concerned, we performed a series of automated molecular docking operations. Moreover, the docking results were quite robust as further validated by molecular dynamics. The docking results reveal that the binding site mainly consists of two hydrophobic regions (P1 and P2 site) which are able to accommodate the lipophilic arms of the compounds investigated. Val421 in P1 site and Met194 in P2 site could be considered to be two important residues in forming the two hydrophobic regions. The presence of residues Phe197 and Phe463 in P2 site may be responsible for the binding recognition through π-π stacking interactions. The hydrophobic 3-phenoxy substituent may be important in creating the preferable inhibitive capability for increasing the binding potency. The hydrophobic character of the tetrafluoroethoxybenzyl group at position 3 displays better hydrophobicity than a shorter hydrophobic substituent. An interaction model of CETP-inhibitors is derived that can be successfully used to explain the different biologic activities of these inhibitors. It is anticipated that the findings reported here may provide very useful information or clues for designing effective drugs for the therapeutic treatment of CETP-related cardiovascular diseases.     

Multisite Binding of a General Anesthetic to the Prokaryotic Pentameric Erwinia chrysanthemi Ligand-gated Ion Channel (ELIC)

Pentameric ligand-gated ion channels (pLGICs), such as nicotinic acetylcholine, glycine, γ-aminobutyric acid GABA_A/C receptors, and the Gloeobacter violaceus ligand-gated ion channel (GLIC), are receptors that contain multiple allosteric binding sites for a variety of therapeutics, including general anesthetics. Here, we report the x-ray crystal structure of the Erwinia chrysanthemi ligand-gated ion channel (ELIC) in complex with a derivative of chloroform, which reveals important features of anesthetic recognition, involving multiple binding at three different sites. One site is located in the channel pore and equates with a noncompetitive inhibitor site found in many pLGICs. A second transmembrane site is novel and is located in the lower part of the transmembrane domain, at an interface formed between adjacent subunits. A third site is also novel and is located in the extracellular domain in a hydrophobic pocket between the β7–β10 strands. Together, these results extend our understanding of pLGIC modulation and reveal several specific binding interactions that may contribute to modulator recognition, further substantiating a multisite model of allosteric modulation in this family of ion channels.     

Fluorescence labeling and computational analysis of the strut of myosin's 50kDa cleft

A new fluorescent labeling procedure specific for the strut sequence of myosin subfragment-1's 50kDa cleft was developed using CY3 N-hydroxy succinimidyl ester as a hydrophobic tag and hydrophobic interaction chromatography to purify the major labeled species which retained actin-activated ATPase activity. Stern-Volmer analysis suggests that the CY3 is in close proximity to basic residues, consistent with inspection of the mapped labeling site in the atomic model. Fluorescence polarization indicates that the CY3 becomes more mobile upon actin binding, supporting a location near the actomyosin interface. In contrast, nucleotide binding to myosin had little impact on the CY3. Molecular mechanics and stochastic dynamics simulations suggest that this labeling site is sensitive to forced cleft opening and closure, but the upper 50kDa cleft does not move easily. In addition, there appear to be some long-range effects of forced cleft opening and closing that could impact the lever arm position.