Structural and kinetic analysis of caspase-3 reveals role for s5 binding site in substrate recognition

The molecular basis for the substrate specificity of human caspase-3 has been investigated using peptide analog inhibitors and substrates that vary at the P2, P3, and P5 positions. Crystal structures were determined of caspase-3 complexes with the substrate analogs at resolutions of 1.7 A to 2.3 A. Differences in the interactions of caspase-3 with the analogs are consistent with the Ki values of 1.3 nM, 6.5 nM, and 12.4 nM for Ac-DEVD-Cho, Ac-VDVAD-Cho and Ac-DMQD-Cho, respectively, and relative kcat/Km values of 100%, 37% and 17% for the corresponding peptide substrates. The bound peptide analogs show very similar interactions for the main-chain atoms and the conserved P1 Asp and P4 Asp, while interactions vary for P2 and P3. P2 lies in a hydrophobic S2 groove, consistent with the weaker inhibition of Ac-DMQD-Cho with polar P2 Gln. S3 is a surface hydrophilic site with favorable polar interactions with P3 Glu in Ac-DEVD-Cho. Ac-DMQD-Cho and Ac-VDVAD-Cho have hydrophobic P3 residues that are not optimal in the polar S3 site, consistent with their weaker inhibition. A hydrophobic S5 site was identified for caspase-3, where the side-chains of Phe250 and Phe252 interact with P5 Val of Ac-VDVAD-Cho, and enclose the substrate-binding site by conformational change. The kinetic importance of hydrophobic P5 residues was confirmed by more efficient hydrolysis of caspase-3 substrates Ac-VDVAD-pNA and Ac-LDVAD-pNA compared with Ac-DVAD-pNA. In contrast, caspase-7 showed less efficient hydrolysis of the substrates with P5 Val or Leu compared with Ac-DVAD-pNA. Caspase-3 and caspase-2 share similar hydrophobic S5 sites, while caspases 1, 7, 8 and 9 do not have structurally equivalent hydrophobic residues; these caspases are likely to differ in their selectivity for the P5 position of substrates. The distinct selectivity for P5 will help define the particular substrates and signaling pathways associated with each caspase.     

Substrate binding mode and its implication on drug design for botulinum neurotoxin A

The seven antigenically distinct serotypes of Clostridium botulinum neurotoxins, the causative agents of botulism, block the neurotransmitter release by specifically cleaving one of the three SNARE proteins and induce flaccid paralysis. The Centers for Disease Control and Prevention (CDC) has declared them as Category A biowarfare agents. The most potent among them, botulinum neurotoxin type A (BoNT/A), cleaves its substrate synaptosome-associated protein of 25 kDa (SNAP-25). An efficient drug for botulism can be developed only with the knowledge of interactions between the substrate and enzyme at the active site. Here, we report the crystal structures of the catalytic domain of BoNT/A with its uncleavable SNAP-25 peptide (197)QRATKM(202) and its variant (197)RRATKM(202) to 1.5 A and 1.6 A, respectively. This is the first time the structure of an uncleavable substrate bound to an active botulinum neurotoxin is reported and it has helped in unequivocally defining S1 to S5' sites. These substrate peptides make interactions with the enzyme predominantly by the residues from 160, 200, 250 and 370 loops. Most notably, the amino nitrogen and carbonyl oxygen of P1 residue (Gln197) chelate the zinc ion and replace the nucleophilic water. The P1'-Arg198, occupies the S1' site formed by Arg363, Thr220, Asp370, Thr215, Ile161, Phe163 and Phe194. The S2' subsite is formed by Arg363, Asn368 and Asp370, while S3' subsite is formed by Tyr251, Leu256, Val258, Tyr366, Phe369 and Asn388. P4'-Lys201 makes hydrogen bond with Gln162. P5'-Met202 binds in the hydrophobic pocket formed by the residues from the 250 and 200 loop. Knowledge of interactions between the enzyme and substrate peptide from these complex structures should form the basis for design of potent inhibitors for this neurotoxin.     

Carbohydrate recognition of gramicidin S analogues in aqueous medium.

We have designed and synthesized of carbohydrate-binding peptides, gramicidin S analogues. Asn/Asp/Gln and Trp residues in the peptides were employed as the binding sites for carbohydrates by hydrogen-bonding interaction and the creation units for hydrophobic pocket to promote the interaction, respectively. The data of fluorescence spectroscopy and affinity column chromatography indicated that the peptides possessed the binding ability for some carbohydrates in aqueous medium. As a result of 1H NMR study, nuclear Overhauser effects between aromatic side chains of a peptide, [Gln(1,1'),Trp(3,3')]-gramisidin S and mannose were observed, indicating that the interaction of the peptide with the sugar occurred in the hydrophobic environment formed by Trp and Phe residues.     

The effect of changing the hydrophobic S1' subsite of thermolysin-like proteases on substrate specificity.

The hydrophobic S1' subsite is one of the major determinants of the substrate specificity of thermolysin and related M4 family proteases. In the thermolysin-like protease (TLP) produced by Bacillus stearothermophilus (TLP-ste), the hydrophobic S1' subsite is mainly formed by Phe130, Phe133, Val139 and Leu202. In the present study, we have examined the effects of replacing Leu202 by smaller (Gly, Ala, Val) and larger (Phe, Tyr) hydrophobic residues. The mutational effects showed that the wild-type S1' pocket is optimal for binding leucine side chains. Reduction of the size of residue 202 resulted in a higher efficiency towards substrates with Phe in the P1' position. Rather unexpectedly, the Leu202-->Phe and Leu202-->Tyr mutations, which were expected to decrease the size of the S1' subsite, resulted in a large increase in activity towards dipeptide substrates with Phe in the P1' position. This is probably due to the fact that 202Phe and 202Tyr adopt a second possible rotamer that opens up the subsite compared to Leu202, and also favours interactions with the substrate. To validate these results, we constructed variants of thermolysin with changes in the S1' subsite. Thermolysin and TLP-ste variants with identical S1' subsites were highly similar in terms of their preference for Phe vs. Leu in the P1' position.     

Substrate Specificity of the Highly Alkalophilic Bacterial Proteinase Esperase: Relation to the X-Ray Structure

Esperase is a highly alkalophilic bacterial proteinase produced by Bacillus lentus. The enzyme hydrolyzes peptide bonds comprising the carboxylic groups of hydrophobic as well as hydrophilic residues in the oxidized insulin B chain. Some of these bonds are not attacked by other alkaline microbial proteinases. P1-P4 specificity was determined by a series of peptide nitroanilides. The S1 recognition loop exhibits a preference for Phe. The "cleft" of the smallest subsite S2 prefers Ala and exhibits low affinity for the larger chain of Leu. S3 is more open than the other subsites and can accept a variety of residues. Hydrophobic interactions predominate in the S4-P4 interactions because S4 can accommodate Phe very well. The results characterize Esperase as an endopeptidase with a broader specificity in comparison with other microbial serine proteinases. This is probably owing to a more flexible substrate binding site.     

Multiple pocket recognition of SNAP25 by botulinum neurotoxin serotype E

Botulinum neurotoxins (BoNTs) are zinc proteases that cleave SNARE proteins to elicit flaccid paralysis by inhibiting the fusion of neurotransmitter-carrying vesicles to the plasma membrane of peripheral neurons. There are seven serotypes of BoNT, termed A-G. The molecular basis for SNAP25 recognition and cleavage by BoNT serotype E is currently unclear. Here we define the multiple pocket recognition of SNAP25 by LC/E. The initial recognition of SNAP25 is mediated by the binding of the B region of SNAP25 to the substrate-binding (B) region of LC/E comprising Leu166, Arg167, Asp127, Ala128, Ser129, and Ala130. The mutations at these residues affected substrate binding and catalysis. Three additional residues participate in scissile bond cleavage of SNAP25 by LC/E. The P3 site residues, Ile178, of SNAP25 interacted with the S3 pocket in LC/E through hydrophobic interactions. The S3 pocket included Ile47, Ile164, and Ile182 and appeared to align the P1' and P2 residues of SNAP25 with the S1' and S2 pockets of LC/E. The S1' pocket of LC/E included three residues, Phe191, Thr159, and Thr208, which contribute hydrophobic and steric interactions with the SNAP25 P1' residue Ile181. The S2 pocket residue of LC/E, Lys224, binds the P2 residue of SNAP25, Asp179, through ionic interactions. Deletion mapping indicates that main chain interaction(s) of residues 182-186 of SNAP25 contribute to substrate recognition by LC/E. Understanding the mechanism for substrate specificity provides insight for the development of inhibitors against the botulinum neurotoxins.     

Mutational replacements in subtilisin 309. Val104 has a modulating effect on the P4 substrate preference.

The previous notion that the amino acid side chain at position 104 of subtilisins is involved in the binding of the side chain at position P4 of the substrate has been investigated. The amino acid residue Val104 in subtilisin 309 has been replaced by Ala, Arg, Asp, Phe, Ser, Trp and Tyr by site-directed mutagenesis. It is shown that the P4 specificity of this enzyme is not determined solely by the amino acid residue occupying position 104, as the enzyme exhibits a marked preference for aromatic groups in P4, regardless of the nature of the position-104 residue. With hydrophilic amino acid residues at this position, no involvement is seen in binding of either hydrophobic or hydrophilic amino acid residues at position P4 of the substrates. The substrate with Asp in P4 is an exception, as the preference for this substrate is increased dramatically by introduction of an arginine residue at position 104 in the enzyme, presumably due to a substrate-induced conformational change. However, when position 104 is occupied by hydrophobic residues, it is highly involved in binding of hydrophobic amino acid residues, either by increasing the hydrophobicity of S4 or by determining the size of the pocket. The results suggest that the amino acid residue at position 104 is mobile such that it is positioned in the S4 binding site only when it can interact favourably with the substrate's side chain at position P4.     

Homology modeling of nematode Caenorhabditis elegans CED3 protein-inhibitor complex

CED3 protein, the product of a gene necessary for programmed cell death in the nematode Caenorhabditis elegans, is related to a highly specific cysteine protease family i.e., caspases. A tertiary-structural model has been constructed of a complex of the CED3 protein with tetrapeptide-aldehyde inhibitor, Ac-DEVD-CHO. The conformation of CED3 protein active site and the general binding features of inhibitor residues are similar to those observed in other caspases. The loop segment (Phe380-Pro387) binds with the P4 Asp in a different fashion compared to caspase-3. The comparative modeling of active sites from caspase-3 and CED3 protein indicated that although these enzymes require Asp at the position P4, variation could occur in the binding of this residue at the S4 subsite. This model allowed the definition of substrate specificity of CED3 protein from the structural standpoint and provided insight in designing of mutants for structure-function studies of this classical caspase homologue.     

Subsite specificity of memapsin 2 (beta-secretase): implications for inhibitor design

Memapsin 2 is the protease known as beta-secretase whose action on beta-amyloid precursor protein leads to the production of the beta-amyloid (Abeta) peptide. Since the accumulation of Abeta in the brain is a key event in the pathogenesis of Alzheimer's disease, memapsin 2 is an important target for the design of inhibitory drugs. Here we describe the residue preference for the subsites of memapsin 2. The relative k(cat)/K(M) values of residues in each of the eight subsites were determined by the relative initial cleavage rates of substrate mixtures as quantified by MALDI-TOF mass spectrometry. We found that each subsite can accommodate multiple residues. The S(1) subsite is the most stringent, preferring residues in the order of Leu > Phe > Met > Tyr. The preferences of other subsites are the following: S(2), Asp > Asn > Met; S(3), Ile > Val > Leu; S(4), Glu > Gln > Asp; S(1)', Met > Glu > Gln > Ala; S(2)', Val > Ile > Ala; S(3)', Leu > Trp > Ala; S(4)', Asp > Glu > Trp. In general, S subsites are more specific than the S' subsites. A peptide comprising the eight most favored residues (Glu-Ile-Asp-Leu-Met-Val-Leu-Asp) was found to be hydrolyzed with the highest k(cat)/K(M) value so far observed for memapsin 2. Residue preferences at four subsites were also studied by binding of memapsin 2 to a combinatorial inhibitor library. From 10 tight binding inhibitors, the consensus preferences were as follows: S(2), Asp and Glu; S(3), Leu and Ile; S(2)', Val; and S(3)', Glu and Gln. An inhibitor, OM00-3, Glu-Leu-Asp-LeuAla-Val-Glu-Phe (where the asterisk represents the hydroxyethylene tansition-state isostere), designed from the consensus residues, was found to be the most potent inhibitor of memapsin 2 so far reported (K(i) of 3.1 x 10(-10) M). A molecular model of OM00-3 binding to memapsin 2 revealed critical improvement of the interactions between inhibitor side chains with enzyme over a previous inhibitor, OM99-2 [Ghosh, A. K., et al. (2000) J. Am. Chem. Soc. 14, 3522-3523].     

Amphipathic property of free thiol group contributes to an increase in the catalytic efficiency of carboxypeptidase Y.

Cys341 of carboxypeptidase Y, which constitutes one side of the solvent-accessible surface of the S1 binding pocket, was replaced with Gly, Ser, Asp, Val, Phe or His by site-directed mutagenesis. Kinetic analysis, using Cbz-dipeptide substrates, revealed that polar amino acids at the 341 position increased K(m) whereas hydrophobic amino acids in this position tended to decrease K(m). This suggests the involvement of Cys341 in the formation of the Michaelis complex in which Cys341 favors the formation of hydrophobic interactions with the P1 side chain of the substrate as well as with residues comprising the surface of the S1 binding pocket. Furthermore, C341G and C341S mutants had significantly higher k(cat) values with substrates containing the hydrophobic P1 side chain than C341V or C341F. This indicates that the nonhydrophobic property conferred by Gly or Ser gives flexibility or instability to the S1 pocket, which contributes to the increased k(cat) values of C341G or C341S. The results suggest that Cys341 may interact with His397 during catalysis. Therefore, we propose a dual role for Cys341: (a) its hydrophobicity allows it to participate in the formation of the Michaelis complex with hydrophobic substrates, where it maintains an unfavorable steric constraint in the S1 subsite; (b) its interaction with the imidazole ring of His397 contributes to the rate enhancement by stabilizing the tetrahedral intermediate in the transition state.     

Caspase-8 specificity probed at subsite S(4): crystal structure of the caspase-8-Z-DEVD-cho complex

Caspase-8 is an initiator enzyme in the Fas-mediated pathway of which the downstream executioner caspase-3 is a physiological target. Caspases are cysteine proteases that are specific for substrates with an aspartic acid residue at the P(1) position and have an optimal recognition motif that incorporates four amino acid residues N-terminal to the cleavage site. Caspase-8 has been classified as a group III caspase member because it shows a preference for a small hydrophobic residue at the P(4) substrate position. We report the X-ray crystallographic structure of caspase-8 in complex with benzyloxycarbonyl-Asp-Glu-Val-Asp-aldehyde (Z-DEVD), a specific group II caspase inhibitor. The structure shows that the inhibitor interacts favourably with the enzyme in subsite S(4). Kinetic data reveal that Z-DEVD (K(i) 2 nM) is an almost equally potent inhibitor of caspase-8 as the specific group III inhibitor Boc-IETD-aldehyde (K(i) 1 nM). In view of this finding, the original classification of caspases into three specificity groups needs to be modified, at least for caspase-8, which tolerates small hydrophobic residues as well as the acidic residue Asp in subsite S(4). We propose that the subsite S(3) must be considered as an important specificity-determining factor.     

Human leukocyte cathepsin G. Subsite mapping with 4-nitroanilides, chemical modification, and effect of possible cofactors

The extended substrate binding site of cathepsin G from human leukocytes has been mapped by using a series of peptide 4-nitroanilide substrates. The enzyme has a significant preference for substrates with a P1 Phe over those with the other aromatic amino acids Tyr and Trp. The S2 subsite was mapped with the substrates Suc-Phe-AA-Phe-NA where AA was 13 of the 20 amino acid residues commonly found in proteins. The best residues were Pro and Met. The S3 subsite was mapped with the sequence Suc-AA-Pro-Phe-NA by using 14 different amino acid residues for AA. The two best residues were the isosteric Val and Thr. No significant improvement in reactivity was obtained by extending the substrate to include seven different P4 residues. The kinetic parameters for cathepsin G are significantly slower than those for many other serine proteases. Changes in the reaction conditions and addition of possible cofactors or ligands were in general found to have little effect on the enzymatic activity, while chemical modifications and proteolysis destroyed the activity of cathepsin G. Cathepsin G hydrolyzed peptides containing model desmosine residues and prefers the hydrophobic picolinoyllysine derivative over lysine by substantial margins at both the S4 and S2 subsites but will not tolerate it at S3. Substrates with sequences related to the cathepsin G cleavage site in angiotensin I and angiotensinogen, and the reactive site of alpha 1-antichymotrypsin, were hydrolyzed effectively by enzyme, but with unexceptional rates. Our results indicate that the natural substrate(s) and function(s) of cathepsin G still remain to be discovered.     

Plasticity of S2-S4 specificity pockets of executioner caspase-7 revealed by structural and kinetic analysis

Many protein substrates of caspases are cleaved at noncanonical sites in comparison to the recognition motifs reported for the three caspase subgroups. To provide insight into the specificity and aid in the design of drugs to control cell death, crystal structures of caspase-7 were determined in complexes with six peptide analogs (Ac-DMQD-Cho, Ac-DQMD-Cho, Ac-DNLD-Cho, Ac-IEPD-Cho, Ac-ESMD-Cho, Ac-WEHD-Cho) that span the major recognition motifs of the three subgroups. The crystal structures show that the S2 pocket of caspase-7 can accommodate diverse residues. Glu is not required at the P3 position because Ac-DMQD-Cho, Ac-DQMD-Cho and Ac-DNLD-Cho with varied P3 residues are almost as potent as the canonical Ac-DEVD-Cho. P4 Asp was present in the better inhibitors of caspase-7. However, the S4 pocket of executioner caspase-7 has alternate regions for binding of small branched aliphatic or polar residues similar to those of initiator caspase-8. The observed plasticity of the caspase subsites agrees very well with the reported cleavage of many proteins at noncanonical sites. The results imply that factors other than the P4-P1 sequence, such as exosites, contribute to the in vivo substrate specificity of caspases. The novel peptide binding site identified on the molecular surface of the current structures is suggested to be an exosite of caspase-7. These results should be considered in the design of selective small molecule inhibitors of this pharmacologically important protease.     

Protein engineering of the high-alkaline serine protease PB92 from Bacillus alcalophilus: functional and structural consequences of mutation at the S4 substrate binding pocket.

Serine endoproteases such as trypsins and subtilisins are known to have an extended substrate binding region that interacts with residues P6 to P3' of a substrate. In order to investigate the structural and functional effects of replacing residues at the S4 substrate binding pocket, the serine protease from the alkalophilic Bacillus strain PB92, which shows homology with the subtilisins, was mutated at positions 102 and 126-128. Substitution of Val102 by Trp results in a 12-fold increase in activity towards succinyl-L-Ala-L-Ala-L-Pro-L-Phe-p-nitroanilide (sAAPFpNA). An X-ray structure analysis of the V102W mutant shows that the Trp side chain occupies a hydrophobic pocket at the surface of the molecule leaving a narrow crevice for the P4 residue of a substrate. Better binding of sAAPFpNA by the mutant compared with the wild type protein as indicated by the kinetic data might be due to the hydrophobic interaction of Ala P4 of the substrate with the introduced Trp102 side chain. The observed difference in binding of sAAPFpNA by protease PB92 and thermitase, both of which possess a Trp at position 102, is probably related to the amino acid substitutions at positions 105 and 126 (in the protease PB92 numbering). Kinetic data for the variants obtained by random mutation of residues Ser126, Pro127 and Ser128 reveal that the activity towards sAAPFpNA increases when a hydrophobic residue is introduced at position 126.(ABSTRACT TRUNCATED AT 250 WORDS)     

Energetic and mechanistic studies of glucoamylase using molecular recognition of maltose OH groups coupled with site-directed mutagenesis

Molecular recognition using a series of deoxygenated maltose analogues was used to determine the substrate transition-state binding energy profiles of 10 single-residue mutants at the active site of glucoamylase from Aspergillus niger. The individual contribution of each substrate hydroxyl group to transition-state stabilization with the wild type and each mutant GA was determined from the relation Delta(DeltaG()) = -RT ln[(k(cat)/K(M))(x)/(k(cat)/K(M))(y)], where x represents either a mutant enzyme or substrate analogue and y the wild-type enzyme or parent substrate. The resulting binding energy profiles indicate that disrupting an active site hydrogen bond between enzyme and substrate, as identified in crystal structures, not only sharply reduces or eliminates the energy contributed from that particular hydrogen bond but also perturbs binding contributions from other substrate hydroxyl groups. Replacing the active site acidic groups, Asp55, Glu180, or Asp309, with the corresponding amides, and the neutral Trp178 with the basic Arg, all substantially reduced the binding energy contribution of the 4'- and 6'-OH groups of maltose at subsite -1, even though both Glu180 and Asp309 are localized at subsite 1. In contrast, the substitution, Asp176 --> Asn, located near subsites -1 and 1, did not substantially perturb any of the individual hydroxyl group binding energies. Similarly, the substitutions Tyr116 --> Ala, Ser119 --> Tyr, or Trp120 --> Phe also did not substantially alter the energy profiles even though Trp120 has a critical role in directing conformational changes necessary for activity. Since the mutations at Trp120 and Asp176 reduced k(cat) values by 50- and 12-fold, respectively, a large effect on k(cat) is not necessarily accompanied by changes in hydroxyl group binding energy contributions. Two substitutions, Asn182 --> Ala and Tyr306 --> Phe, had significant though small effects on interactions with 3- and 4'-OH, respectively. Binding interactions between the enzyme and the glucosyl group in subsite -1, particularly with the 4'- and 6'-OH groups, play an important role in substrate binding, while subsite 1 interactions may play a more important role in product release.     

Crystal structure of an antibody bound to an immunodominant peptide epitope: novel features in peptide-antibody recognition

The crystal structure of Fab of an Ab PC283 complexed with its corresponding peptide Ag, PS1 (HQLDPAFGANSTNPD), derived from the hepatitis B virus surface Ag was determined. The PS1 stretch Gln2P to Phe7P is present in the Ag binding site of the Ab, while the next three residues of the peptide are raised above the binding groove. The residues Ser11P, Thr12P, and Asn13P then loop back onto the Ag-binding site of the Ab. The last two residues, Pro14P and Asp15P, extend outside the binding site without forming any contacts with the Ab. The PC283-PS1 complex is among the few examples where the light chain complementarity-determining regions show more interactions than the heavy chain complementarity-determining regions, and a distal framework residue is involved in Ag binding. As seen from the crystal structure, most of the contacts between peptide and Ab are through the five residues, Leu3-Asp4-Pro5-Ala6-Phe7, of PS1. The paratope is predominantly hydrophobic with aromatic residues lining the binding pocket, although a salt bridge also contributes to stabilizing the Ag-Ab interaction. The molecular surface area buried upon PS1 binding is 756 A(2) for the peptide and 625 A(2) for the Fab, which is higher than what has been seen to date for Ab-peptide complexes. A comparison between PC283 structure and a homology model of its germline ancestor suggests that paratope optimization for PS1 occurs by improving both charge and shape complementarity.     

Structural basis for executioner caspase recognition of P5 position in substrates

Caspase-3, -6 and -7 cleave many proteins at specific sites to induce apoptosis. Their recognition of the P5 position in substrates has been investigated by kinetics, modeling and crystallography. Caspase-3 and -6 recognize P5 in pentapeptides as shown by enzyme activity data and interactions observed in the crystal structure of caspase-3/LDESD and in a model for caspase-6. In caspase-3 the P5 main-chain was anchored by interactions with Ser209 in loop-3 and the P5 Leu side-chain interacted with Phe250 and Phe252 in loop-4 consistent with 50% increased hydrolysis of LDEVD relative to DEVD. Caspase-6 formed similar interactions and showed a preference for polar P5 in QDEVD likely due to interactions with polar Lys265 and hydrophobic Phe263 in loop-4. Caspase-7 exhibited no preference for P5 residue in agreement with the absence of P5 interactions in the caspase-7/LDESD crystal structure. Initiator caspase-8, with Pro in the P5-anchoring position and no loop-4, had only 20% activity on tested pentapeptides relative to DEVD. Therefore, caspases-3 and -6 bind P5 using critical loop-3 anchoring Ser/Thr and loop-4 side-chain interactions, while caspase-7 and -8 lack P5-binding residues.     

Interactions of human melanocortin 4 receptor with nonpeptide and peptide agonists

Specific interactions of human melanocortin-4 receptor (hMC4R) with its nonpeptide and peptide agonists were studied using alanine-scanning mutagenesis. The binding affinities and potencies of two synthetic, small-molecule agonists (THIQ, MB243) were strongly affected by substitutions in transmembrane alpha-helices (TM) 2, 3, 6, and 7 (residues Glu(100), Asp(122), Asp(126), Phe(261), His(264), Leu(265), and Leu(288)). In addition, a I129A mutation primarily affected the binding and potency of THIQ, while F262A, W258A, Y268A mutations impaired interactions with MB243. By contrast, binding affinity and potency of the linear peptide agonist NDP-MSH were substantially reduced only in D126A and H264A mutants. Three-dimensional models of receptor-ligand complexes with their agonists were generated by distance-geometry using the experimental, homology-based, and other structural constraints, including interhelical H-bonds and two disulfide bridges (Cys(40)-Cys(279), Cys(271)-Cys(277)) of hMC4R. In the models, all pharmacophore elements of small-molecule agonists are spatially overlapped with the corresponding key residues (His(6), d-Phe(7), Arg(8), and Trp(9)) of the linear peptide: their charged amine groups interact with acidic residues from TM2 and TM3, similar to His(6) and Arg(6) of NDP-MSH; their substituted piperidines mimic Trp(9) of the peptide and interact with TM5 and TM6, while the d-Phe aromatic rings of all three agonists contact with Leu(133), Trp(258), and Phe(261) residues.     

Investigate the Binding of Catechins to Trypsin Using Docking and Molecular Dynamics Simulation

To explore the inhibitory mechanism of catechins for digestive enzymes, we investigated the binding mode of catechins to a typical digestive enzyme-trypsin and analyzed the structure-activity relationship of catechins, using an integration of molecular docking, molecular dynamics simulation and binding free energy calculation. We found that catechins with different structures bound to a conservative pocket S1 of trypsin, which is comprised of residues 189–195, 214–220 and 225–228. In the trypsin-catechin complexes, Asp189 by forming strong hydrogen bonding, and Gln192, Trp215 and Gly216 through hydrophobic interactions, all significantly contribute to the binding of catechins. The number and the position of hydroxyl and aromatic groups, the structure of stereoisomers, and the orientation of catechins in the binding pocket S1 of trypsin all affect the binding affinity. The binding affinity is in the order of Epigallocatechin gallate (EGCG) > Epicatechin gallate (ECG) > Epicatechin (EC) > Epigallocatechin (EGC), and 2R-3R EGCG shows the strongest binding affinity out of other stereoisomers. Meanwhile, the synergic conformational changes of residues and catechins were also analyzed. These findings will be helpful in understanding the knowledge of interactions between catechins and trypsin and referable for the design of novel polyphenol based functional food and nutriceutical formulas.     

Molecular dynamics of NPY Y1 receptor activation

A three-dimensional model of the human neuropeptide Y(NPY)Y₁ receptor (hY₁) was constructed, energy refined and used to simulate molecular receptor interactions of the peptide ligands NPY, [L31, P34]NPY, peptide YY (PYY) and pancreatic polypeptide (PP), and of the nonpeptide antagonist R-N²-(diphenylacetyl)-N-(4-hydroxyphenyl)methyl-argininamide (BIBP3226) and its S-enantiomer BIBP3435. The best complementarity in charges between the receptor and the peptides, and the best structural accordance with experimental studies, was obtained with amino acid 1–4 of the peptides interacting with Asp194, Asp200, Gln201, Phe202 and Trp288 in the receptor. Arg33 and Arg35 of the peptides formed salt bridges with Asp104 and Asp287, respectively, while Tyr36 interacted in a binding pocket formed by Phe41, Thr42, Tyr100, Asn297, His298 and Phe302. Calculated electrostatic potentials around NPY and hY₁ molecules indicated that ligand binding is initiated by electrostatic interactions between a highly positive region in the N- and C-terminal parts of the peptides, and a negative region in the extracellular receptor domains. Molecular dynamics simulations of NPY and BIBP3226 interactions with the receptor indicated rigid body motions of TMH5 and TMH6 upon NPY binding as mechanisms of receptor activation, and that BIBP3226 may act as an antagonist by constraining these motions.