Substrate stereoselectivity of mammalian D-aspartyl endopeptidase

The formation and accumulation of D-aspartate residue (D-Asp) in proteins caused by oxidative stress leads to dysfunction and/or denaturation of proteins, and is consequently responsible for aging-related misfolding diseases such as cataracts, prion disease, and Alzheimer's disease. We sought to identify that an unknown protease selectively degrades the noxious D-Asp-containing protein, namely D-aspartyl endopeptidase (DAEP), and finally purified it from the inner mitochondrial membrane of mouse liver. In order to analyze the substrate stereoselectivity of DAEP, we synthesized a peptide corresponding to 55-65 (Thr-Val-Leu-Asp-Ser-Gly-Ile-Ser-Glu-Val-Arg) of human αA-crystallin and its corresponding diastereoisomers in which L-α-Asp was replaced with L-β-, D-α- or D-β-Asp residue at position 58. Following incubation of that peptide with purified DAEP, it was only degraded at D-α-Asp(58), independent of ATP or NAD. This result indicates that DAEP stereoselectively recognizes and degrades its substrate at the internal D-α-Asp residue. DAEP therefore seems to physiologically serve as the quality control system against the noxious D-Asp-containing protein in the long life span of mammals.     

D-Amino acids in aged proteins: analysis and biological relevance

Homochirality is essential for life. L-Amino acids are exclusively used as substrates for the polymerization and formation of peptides and proteins in living systems. However, d-amino acids, which are enantiomers of L-amino acids, were recently detected in various living organisms in the form of free D-amino acids and D-amino acid residues in peptides and proteins. In particular, D-aspartyl (Asp) residues have been detected in various proteins from diverse tissues of elderly individuals. Here, we describe three important aspects of our research: (i) a method for detecting D-β-Asp at specific sites in particular proteins, (ii) a likely spontaneous mechanism by which Asp residues in proteins invert and isomerize to the D-β-form with age under physiological conditions, (iii) a discussion of factors that favor such a reaction.     

Age-related isomerization of Asp in human immunoglobulin G kappa chain

Isomerization of aspartate (Asp) is a common non-enzymatic posttranslational modification. Isomerized residues accumulate in proteins associated with age-related human disorders such as cataract and are well known to affect protein structure and function. We previously detected d-Asp-containing peptides in human serum. In this study, we investigated whether isomerized Asp residues are present in human immunoglobulin G (IgG) kappa chain by a qualitative d-amino acid analysis based on diastereomer formation and liquid chromatography tandem mass spectrometry (LC-MS/MS). We also investigated the d/l ratio of Asp residues in the IgG kappa chain in serum from donors aged 25, 37, 41, 54 and 67 years. As a result, two isomerized Asp residues, Asp151 and Asp170, were detected in the IgG kappa chain, and the d/l ratio of these residues was found to increase with aging. To assess the effects of this isomerization, we synthesized four isomeric peptides of IgG kappa chain containing lα-, lβ-, dα-, or dβ-Asp at position 170, and compared their secondary structures by CD spectroscopy. Peptide containing normal lα-Asp170 showed type II β-turn structure, while the other isomeric peptides showed random structure, clearly indicating that substitution of a single Asp isomer alters the secondary structure of the peptide. Because IgG is a main component of humoral immunity, Asp isomerization in IgG may reflect changes of structure and decrease in immune function. Proteome research on serum from the standpoint of racemization might enable us to develop new kinds of biomarker and new directions to study the aging process.     

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.     

Simultaneous stereoinversion and isomerization at the Asp-4 residue in βB2-crystallin from the aged human eye lenses

The lens proteins are composed of α-, β-, and γ-crystallins that interact with each other to maintain the transparency and refractive power of the lens. Because the lens crystallins are long-lived proteins, they undergo various post-translational modifications including racemization, isomerization, deamidation, oxidation, glycation, and truncation. In βB2-crystallin, which is the most abundant β-crystallin, the deamidation of asparagine and glutamine residues has been reported. Here, we found that the aspartyl (Asp) residue at position 4 of βB2-crystallin in the lenses of elderly human individuals undergoes a significant degree of inversion and isomerization to the biologically uncommon residue D-β-Asp. Surprisingly, the D/L ratio of β-Asp at position 4 in βB2-crystallin from elderly donors (67-77 year old) was 0.88-3.21. A D/L ratio of amino acids greater than 1.0 is defined as an inversion of configuration from the L- to D-form, rather than a racemization. These extremely high D/L ratios are equivalent to those of Asp-58 and Asp-151 (D/L ratio: 3.1 for Asp-58 and 5.7 for Asp-151) in αA-crystallin from elderly donors (~80 year old) as reported previously. Initially, we identified specific Asp residues in the β-crystallin family of proteins that undergo a high degree of inversion. These results show that the isomerization and inversion of Asp residues occurs both in the α- and β-crystallins of the lens. Inversion of these Asp residues directly affects the higher order structure of the protein. Hence, this modification may change crystallin-crystallin interactions and disrupt the function of crystallins in the lens.     

Side-chain conformers to allow conversion from normal to isoaspartate in age-related proteins and peptides

Dβ (or D-iso)- and Lβ- (or iso)- aspartyl (Asp) residues are accumulated in aged lens crystallins and amyloid beta (Aβ) proteins, respectively, as a result of spontaneous, nonenzymatic isomerization of normal Lα-Asp. To explore why such uncommon Asp isomers are accumulated, the stability of Lα-, Lβ-, and Dβ-Asp was compared in view of the staggered side-chain conformers. By using cylindrin (KVKVLGD⁷VIEV) from αB-crystallin and Aβ17-25 (L¹⁷VFF²⁰AED²³)VG²⁵) containing Asp isomers, the vicinal spin-spin coupling constants of Asp Hα-H_β1 and Hα-H_β2 were quantified by high-resolution solution ¹H NMR. It was found that the trans conformer was extremely preferred in Dβ-Asp7 side-chain of cylindrin. In Aβ17–25, the side chain of Lβ-Asp23 was likely to adopt trans conformer, while gauche conformers were rather rich in Lα-Asp23. In gauche conformers, the close distance between Asp carboxylate carbon (C_COO-) and backbone nitrogen (N) next to Asp is advantageous to the intramolecular cyclization to form succinimide intermediate, followed by the conversion from α- to β-Asp. The cyclization is limited in the trans conformer because of the long distance between C_COO- and N, to keep Dβ- or Lβ-Asp stable. This would be the reason for the site specificity of Asp isomerization in proteins. The higher population of trans conformer in Asp side chain, the less isomerization of Asp as shown as Asp76 in αA-crystallin. The stability and less reactivity of normal Asp and its isomers are the potential factors to determine whether or not the abnormal accumulation is permitted in aged crystallins and Aβ.     

Comparison of the activation energy barrier for succinimide formation from α- and β-aspartic acid residues obtained from density functional theory calculations

The l-α-Asp residues in peptides or proteins are prone to undergo nonenzymatic reactions to form l-β-Asp, d-α-Asp, and d-β-Asp residues via a succinimide five-membered ring intermediate. From these three types of isomerized aspartic acid residues, particularly d-β-Asp has been widely detected in aging tissue. In this study, we computationally investigated the cyclization of α- and β-Asp residues to form succinimide with dihydrogen phosphate ion as a catalyst (H₂PO₄⁻). We performed the study using B3LYP/6-31 + G(d,p) density functional theory calculations. The comparison of the activation barriers of both residues is discussed. All the calculations were performed using model compounds in which an α/β-Asp-Gly sequence is capped with acetyl and methylamino groups on the N- and C-termini, respectively. Moreover, H₂PO₄⁻ catalyzes all the steps of the succinimide formation (cyclization-dehydration) acting as a proton-relay mediator. The calculated activation energy barriers for succinimide formation of α- and β-Asp residues are 26.9 and 26.0 kcal mol^− 1, respectively. Although it was experimentally confirmed that β-Asp has higher stability than α-Asp, there was no clear difference between the activation barriers. Therefore, the higher stability of β-Asp residue than α-Asp residue may be caused by an entropic effect associated with the succinimide formation.     

An S188V Mutation Alters Substrate Specificity of Non-Stereospecific α-Haloalkanoic Acid Dehalogenase E (DehE)

The non-stereospecific α-haloalkanoic acid dehalogenase E (DehE) degrades many halogenated compounds but is ineffective against β-halogenated compounds such as 3-chloropropionic acid (3CP). Using molecular dynamics (MD) simulations and site-directed mutagenesis we show here that introducing the mutation S188V into DehE improves substrate specificity towards 3CP. MD simulations showed that residues W34, F37, and S188 of DehE were crucial for substrate binding. DehE showed strong binding ability for D-2-chloropropionic acid (D-2CP) and L-2-chloropropionic acid (L-2CP) but less affinity for 3CP. This reduced affinity was attributed to weak hydrogen bonding between 3CP and residue S188, as the carboxylate of 3CP forms rapidly interconverting hydrogen bonds with the backbone amide and side chain hydroxyl group of S188. By replacing S188 with a valine residue, we reduced the inter-molecular distance and stabilised bonding of the carboxylate of 3CP to hydrogens of the substrate-binding residues. Therefore, the S188V can act on 3CP, although its affinity is less strong than for D-2CP and L-2CP as assessed by K_m. This successful alteration of DehE substrate specificity may promote the application of protein engineering strategies to other dehalogenases, thereby generating valuable tools for future bioremediation technologies.     

Brownian dynamics simulations of the recognition of the scorpion toxin maurotoxin with the voltage-gated potassium ion channels

The recognition of the scorpion toxin maurotoxin (MTX) by the voltage-gated potassium (Kv1) channels, Kv1.1, Kv1.2, and Kv1.3, has been studied by means of Brownian dynamics (BD) simulations. All of the 35 available structures of MTX in the Protein Data Bank (http://www.rcsb.org/pdb) determined by nuclear magnetic resonance were considered during the simulations, which indicated that the conformation of MTX significantly affected both the recognition and the binding between MTX and the Kv1 channels. Comparing the top five highest-frequency structures of MTX binding to the Kv1 channels, we found that the Kv1.2 channel, with the highest docking frequencies and the lowest electrostatic interaction energies, was the most favorable for MTX binding, whereas Kv1.1 was intermediate, and Kv1.3 was the least favorable one. Among the 35 structures of MTX, the 10th structure docked into the binding site of the Kv1.2 channel with the highest probability and the most favorable electrostatic interactions. From the MTX-Kv1.2 binding model, we identified the critical residues for the recognition of these two proteins through triplet contact analyses. MTX locates around the extracellular mouth of the Kv1 channels, making contacts with its beta-sheets. Lys23, a conserved amino acid in the scorpion toxins, protrudes into the pore of the Kv1.2 channel and forms two hydrogen bonds with the conserved residues Gly401(D) and Tyr400(C) and one hydrophobic contact with Gly401(C) of the Kv1.2 channel. The critical triplet contacts for recognition between MTX and the Kv1.2 channel are Lys23(MTX)-Asp402(C)(Kv1), Lys27(MTX)-Asp378(D)(Kv1), and Lys30(MTX)-Asp402(A)(Kv1). In addition, six hydrogen-bonding interactions are formed between residues Lys23, Lys27, Lys30, and Tyr32 of MTX and residues Gly401, Tyr400, Asp402, Asp378, and Thr406 of Kv1.2. Many of them are formed by side chains of residues of MTX and backbone atoms of the Kv1.2 channel. Five hydrophobic contacts exist between residues Pro20, Lys23, Lys30 and Tyr32 of MTX and residues Asp402, Val404, Gly401, and Arg377 of the Kv1.2 channel. The simulation results are in agreement with the previous molecular biology experiments and explain the binding phenomena between MTX and Kv1 channels at the molecular level. The consistency between the results of the BD simulations and the experimental data indicated that our three-dimensional model of the MTX-Kv1.2 channel complex is reasonable and can be used in additional biological studies, such as rational design of novel therapeutic agents blocking the voltage-gated channels and in mutagenesis studies in both the toxins and the Kv1 channels. In particular, both the BD simulations and the molecular mechanics refinements indicate that residue Asp378 of the Kv1.2 channel is critical for its recognition and binding functionality toward MTX. This phenomenon has not been appreciated in the previous mutagenesis experiments, indicating this might be a new clue for additional functional study of Kv1 channels.     

D-Amino acids in protein: The mirror of life as a molecular index of aging

Proteins are composed exclusively of l-amino acids. Among elderly individuals, however, d-aspartic acid (d-Asp) residues have been found in eye lens and brain, as well as in other tissues. The presence of d-Asp may change the higher-order structure of a protein, which in turn may have a role in age-related disorders such as cataract and Alzheimer's disease. d-Asp results from the spontaneous racemization of Asp residues in susceptible proteins. During aging, natural lα-Asp residues in proteins are non-enzymatically isomerized via a succinimidyl intermediate to l-β-, d-α- and d-β-isomers. This isomerization does not happen uniformly, but instead occurs at specific residues that are susceptible to isomerization due to their sequence or structural context. Thus, it is necessary to establish the nature of each individual Asp residue in susceptible proteins. Recently, a new method based on LC-MS/MS for the analysis of Asp isomerization at specific protein sites has been described. In this review, we first show that the homochirality of amino acids in proteins is not guaranteed throughout life. We then describe the development of a new method for protein-bound d-amino acid analysis, and discuss the negative influence that d-Asp has on protein structure and function.     

Brownian dynamics simulations of the recognition of the scorpion toxin P05 with the small-conductance calcium-activated potassium channels

The recognition of the scorpion toxin P05 and the small-conductance, calcium-activated potassium (SK) channels, rsk1, rsk2, and rsk3, has been studied by means of the Brownian dynamics (BD) method. All of the 25 available structures of P05 in the RCSB Protein Data Bank determined by NMR were considered during the simulation, which indicated that the conformation of P05 affects both the recognition and binding between the two proteins significantly. Comparing the top four high-frequency structures of P05 binding to the SK channels, we found that the rsk2 channel, with high frequencies and lowest electrostatic interaction energies (E (int)(ES)), is the most favorable for P05 binding, while rsk3 is intermediate, and rsk1 is the least favorable. Among the 25 structures of P05, the 13th structure docks into the binding site of the rsk2 channel with the highest probability and most favorable electrostatic interactions. From the P05-rsk2 channel binding model, we identified the residues critical for the recognition of these two proteins through triplet contact analyses. P05 locates around the extracellular mouth of the SK channels and contacts the SK channels using its alpha-helix rather than beta-sheets. The critical triplet contacts for recognition between P05 and the rsk2 channel are Arg6 (P05)-Asp364 (SK), Arg7 (P05)-Asn368 (SK), and Arg13 (P05)-Asp341 (SK). The structure of the P05-rsk2 complex with the most favorable electrostatic interaction energy was further refined by molecular mechanics, showing that six hydrogen bonding interactions exist between P05 and the rsk2 channel: one hydrogen bond is formed between Arg6 (P05) and Asp364(D) (rsk2); Arg7 (P05) forms three hydrogen bonds with Asp341(B) (rsk2)) and Asp364(C) (rsk2); two hydrogen bonds are formed by Arg13 (P05) with Asp341(A) (rsk2) and Asp364(B) (rsk2). The simulation results are in good agreement with the previous molecular biological experiments and can explain the binding phenomena between P05 and SK channels at the level of molecular structure. The consistency between the results of the BD simulations and the experimental data indicated that our 3D model of the P05-rsk2 channel complex is reasonable and can be employed in further biological studies, such as rational design of the novel therapeutic agents blocking the small-conductance, calcium-activated and apamin-sensitive potassium channels, and for mutagenesis studies in both toxins and SK channels. In particular, both the BD simulations and the molecular mechanics refinements indicate that residue Asp364 of the rsk2 channel is critical for its recognition and binding functionality towards P05. This phenomenon has not been appreciated in the previous mutagenesis experiments, indicating that this might be a new clue for further functional study of SK channels.     

Caspase-3 binds diverse P4 residues in peptides as revealed by crystallography and structural modeling

Caspase-3 recognition of various P4 residues in its numerous protein substrates was investigated by crystallography, kinetics, and calculations on model complexes. Asp is the most frequent P4 residue in peptide substrates, although a wide variety of P4 residues are found in the cellular proteins cleaved by caspase-3. The binding of peptidic inhibitors with hydrophobic P4 residues, or no P4 residue, is illustrated by crystal structures of caspase-3 complexes with Ac-IEPD-Cho, Ac-WEHD-Cho, Ac-YVAD-Cho, and Boc-D(OMe)-Fmk at resolutions of 1.9–2.6 Å. The P4 residues formed favorable hydrophobic interactions in two separate hydrophobic regions of the binding site. The side chains of P4 Ile and Tyr form hydrophobic interactions with caspase-3 residues Trp206 and Trp214 within a non-polar pocket of the S4 subsite, while P4 Trp interacts with Phe250 and Phe252 that can also form the S5 subsite. These interactions of hydrophobic P4 residues are distinct from those for polar P4 Asp, which indicates the adaptability of caspase-3 for binding diverse P4 residues. The predicted trends in peptide binding from molecular models had high correlation with experimental values for peptide inhibitors. Analysis of structural models for the binding of 20 different amino acids at P4 in the aldehyde peptide Ac-XEVD-Cho suggested that the majority of hydrophilic P4 residues interact with Phe250, while hydrophobic residues interact with Trp206, Phe250, and Trp214. Overall, the S4 pocket of caspase-3 exhibits flexible adaptation for different residues and the new structures and models, especially for hydrophobic P4 residues, will be helpful for the design of caspase-3 based drugs.     

Insight into the stereochemistry in the inhibition of carboxypeptidase A with N-(hydroxyaminocarbonyl)phenylalanine: binding modes of an enantiomeric pair of the inhibitor to carboxypeptidase A

Both D- and L-isomers of N-(hydroxyaminocarbonyl)phenylalanine () were shown to have strong binding affinity towards carboxypeptidase A (CPA) with D- being more potent than its enantiomer by 3-fold (Chung, S. J.; Kim, D. H. Bioorg. Med. Chem. 2001, 9, 185.). In order to understand the reversed stereochemical preference shown in the CPA inhibition, we have solved the crystal structures of CPA complexed with each enantiometer of up to 1.75 A resolution. Inhibitor L- whose stereochemistry belongs to the stereochemical series of substrate binds CPA like substrate does with its carbonyl oxygen coordinating to the active site zinc ion. Its hydroxyl is engaged in hydrogen bonding with the carboxylate of Glu-270. On the other hand, in binding of D- to CPA, its terminal hydroxyl group is involved in interactions with the active site zinc ion and the carboxylate of Glu-270. In both CPA small middle dot complexes, the phenyl ring in is fitted in the substrate recognition pocket at the S(1)' subsite, and the carboxylate of the inhibitors forms bifurcated hydrogen bonds with the guanidinium moiety of Arg-145 and a hydrogen bond with the guanidinium of Arg-127. In the complex of CPA small middle dotD-, the carboxylate of the inhibitor is engaged in hydrogen bonding with the phenolic hydroxyl of the down-positioned Tyr-248. While the L- binding induces a concerted movement of the backbone amino acid residues at the active site, only the downward movement of Tyr-248 was noted when D- binds to CPA.     

In silico characterization of a novel β-1,3-glucanase gene from Bacillus amyloliquefaciens—a bacterial endophyte of Hevea brasiliensis antagonistic to Phytophthora meadii

We report the molecular characterization of β-1,3-glucanase-producing Bacillus amyloliquefaciens—an endophyte of Hevea brasiliensis antagonistic to Phytophthora meadii. After cloning and sequencing, the β-1,3-glucanase gene was found to be 747 bp in length. A homology model of the β-1,3-glucanase protein was built from the amino acid sequence obtained upon translation of the gene. The target β-1,3-glucanase protein and the template protein, endo β-1,3-1,4-glucanase protein (PDB ID: 3o5s), were found to share 94 % sequence identity and to have similar secondary and tertiary structures. In the modeled structure, three residues in the active site region of the template—Asn52, Ile157 and Val158—were substituted with Asp, Leu and Ala, respectively. Computer-aided docking studies of the substrate disaccharide (β-1, 3-glucan) with the target as well as with the template proteins showed that the two protein-substrate complexes were stabilized by three hydrogen bonds and by many van der Waals interactions. Although the binding energies and the number of hydrogen bonds were the same in both complexes, the orientations of the substrate in the active sites of the two proteins were different. These variations might be due to the change in the three amino acids in the active site region of the two proteins. The difference in substrate orientation in the active site could also affect the catalytic potential of the β-1,3 glucanase enzyme.     

Identification of secretin,vasoactive intestinal peptide and glucagon binding sites: from chimaeric receptors to point mutations

We have identified two basic residues that are important for the recognition of secretin and vasoactive intestinal peptide (VIP) by their respective receptors. These two peptides containing an Asp residue at position 3 interacted with an arginine residue in transmembrane helix 2 (TM2) of the receptor, and the lysine residue in extracellular loop 1 (ECL1) stabilized the active receptor conformation induced by the ligand. The glucagon receptor possesses a Lys instead of an Arg in TM2, and an Ile instead of Lys in ECL1; it markedly prefers a Gln side chain in position 3 of the ligand. Our results suggested that, in the wild-type receptor, the Ile side chain prevented access to the TM2 Lys side chain, but oriented the glucagon Gln(3) side chain to its proper binding site. In the double mutant, the ECL1 Lys allowed an interaction between negatively charged residues in position 3 of glucagon and the TM2 Arg, resulting in efficient receptor activation by [Asp(3)]glucagon as well as by glucagon.     

UV B-irradiation enhances the racemization and isomerizaiton of aspartyl residues and production of Nε-carboxymethyl lysine (CML) in keratin of skin

UV-B irradiation is one of the risk factors in age-related diseases. We have reported that biologically uncommon D-β-Asp residues accumulate in proteins from sun-exposed elderly human skin. A previous study also reported that carboxymethyl lysine (CML; one of the advanced glycation end products (AGEs)) which is produced by the oxidation of glucose and peroxidation of lipid, also increases upon UV B irradiation. The formation of D-β-Asp and CML were reported as the alteration of proteins in UV B irradiated skin, independently. In this study, in order to clarify the relationship between the formation of D-β-Asp and CML, immunohistochemical analysis using anti-D-β-Asp containing peptide antibodies and anti-CML antibodies was performed in UV B irradiated mice. Immunohistochemical analyses clearly indicated that an anti-D-β-Asp containing peptide antibody and anti-CML antibody reacted at a common area in UV B irradiated skin. Western blot analyses of the proteins isolated from UV B irradiated skin demonstrated that proteins of 50-70 kDa were immunoreactive towards antibodies for both D-β-Asp containing peptide and CML. These proteins were identified by proteomic analysis as members of the keratin families including keratin-1, keratin-6B, keratin-10, and keratin-14.     

Mutagenesis and molecular modeling reveal the importance of the 5-HT3 receptor F-loop

The 5-HT(3) receptor is a member of the Cys-loop family of ligand-gated ion channels. The extracellular domains of these proteins contain six amino acid loops (A-F) that converge to form the ligand binding site. In this study we have mutated 21 residues in or close to the 5-HT(3) receptor F-loop (Ile(192) to Gly(212)) to Ala or to a residue with similar chemical properties. Mutant receptors were expressed in HEK293 cells, and binding affinity was measured using [(3)H]granisetron. Two regions displayed decreases in binding affinity when mutated to Ala (Ile(192)-Arg(196) and Asp(204)-Ser(206)), but only one region was sensitive when mutated to chemically similar residues (Ile(192)-Val(201)). Homology modeling using acetylcholine-binding protein crystal structures with a variety of different bound ligands suggests there may be distinct movements of Trp(195) and Asp(204) upon ligand binding, indicating that these residues and their immediate neighbors have the ability to interact differently with different ligands. The models suggest predominantly lateral movement around Asp(204) and rotational movement around Trp(195), indicating the former is in a more flexible region. Overall our results are consistent with a flexible 5-HT(3) receptor F-loop with two regions that have specific but distinct roles in ligand binding.     

C‐H…O hydrogen bonds in FK506‐binding protein–ligand interactions

Hydrogen bonds are important interaction forces observed in protein structures. They can be classified as stronger or weaker depending on their energy, thereby reflecting on the type of donor. The contribution of weak hydrogen bonds is deemed as an important factor toward structure stability along with the stronger bonds. One such bond, the C‐H…O type hydrogen bond, is shown to make a contribution in maintaining three dimensional structures of proteins. Apart from their presence within protein structures, the role of these bonds in protein–ligand interactions is also noteworthy. In this study, we present a statistical analysis on the presence of C‐H…O hydrogen bonds observed between FKBPs and their cognate ligands. The FK506‐binding proteins (FKBPs) carry peptidyl cis–trans isomerase activity apart from the immunosuppressive property by binding to the immunosuppressive drugs FK506 or rapamycin. Because the active site of FKBPs is lined up by many hydrophobic residues, we speculated that the prevalence of C‐H…O hydrogen bonds will be considerable. In a total of 25 structures analyzed, a higher frequency of C‐H…O hydrogen bonds is observed in comparison with the stronger hydrogen bonds. These C‐H…O hydrogen bonds are dominated by a highly conserved donor, the C^α/β of Val55 and an acceptor, the backbone oxygen of Glu54. Both these residues are positioned in the β4‐α1 loop, whereas the other residues Tyr26, Phe36 and Phe99 with higher frequencies are lined up at the opposite face of the active site. These preferences could be implicated in FKBP pharmacophore models toward enhancing the ligand affinity. This study could be a prelude to studying other proteins with hydrophobic pockets to gain better insights into ligand recognition. Copyright © 2013 John Wiley & Sons, Ltd.     

Characterization of rat and mouse NAD+-dependent 3alpha/17beta/20alpha-hydroxysteroid dehydrogenases and identification of substrate specificity determinants by site-directed mutagenesis

In this study, we characterized rat and mouse aldo-keto reductases (AKR1C16 and AKR1C13, respectively) with 92% sequence identity. The recombinant enzymes oxidized non-steroidal alcohols using NAD⁺ as the preferred coenzyme, and showed low 3α/17β/20α-hydroxysteroid dehydrogenase (HSD) activities. The substrate specificity differs from that of rat NAD⁺-dependent 3α-HSD (AKR1C17) that shares 95% sequence identity with AKR1C16. To elucidate the residues determining the substrate specificity of the enzymes, we performed site-directed mutagenesis of Tyr24, Asp128 and Phe129 of AKR1C16 with the corresponding residues (Ser, Tyr and Leu, respectively) of AKR1C17. The double mutation (Asp128/Tyr-Phe129/Leu) had few effects on the substrate specificity, while the Tyr24/Ser mutant showed only 3α-HSD activity, and the triple mutation of the three residues produced an enzyme that had almost the same properties as AKR1C17. The importance of the residue 24 for substrate recognition was verified by the mutagenesis of Ser24/Tyr of AKR1C17 which resulted in a decrease in 3α-HSD activity and appearance of 17β- and 20α-HSD activities. AKR1C16 is also 92% identical with rat NAD⁺-dependent 17β-HSD (AKR1C24), which possesses Tyr24. The replacement of Asp128, Phe129 and Ser137 of AKR1C16 with the corresponding residues (Glu, Ser and Phe, respectively) of AKR1C24 increased the catalytic efficiency for 17β- and 20α-hydroxysteroids.     

Pharmacophore-directed Homology Modeling and Molecular Dynamics Simulation of G Protein-coupled Receptor： Study of Possible Binding Modes of 5-HT2c Receptor Agonists

A new pharmacophore-based modeling procedure, including homology modeling, pharmacophore study, flexible molecular docking, and long-time molecular dynamics (MD) simulations, was employed to construct the structure of the human 5-HT_(2C) receptor and determine the characteristics of binding modes of 5-HT_(2C) receptor agonists. An agonist-receptor complex has been constructed based on homology modeling and a pharmacophore hypothesis model based on some high active compounds. Then MD simulations of the ligand-receptor complex in an explicit membrane environment were carried out. The conformation of the 5- HT_(2C) receptor during MD simulation was explored, and the stable binding modes of the studied agonist were determined. Flexible molecular docking of several structurally diverse agonists of the human 5-HT_(2C) receptor was carried out, and the general binding modes of these agonists were investigated. According to the models presented in this work and the results of Flexi-Dock, the involvement of the amino acid residues Asp134, Ser138, Ash210, Asn331, Tyr358, Ile131, Ser132, Val135, Thr139, Ile189, Val202, Val208, Leu209, Phe214, Val215, Gly218, Ser219, Phe223, Trp324, Phe327, and Phe328 in agonist recognition was studied. The obtained binding modes of the human 5-HT_(2C) receptor agonists have good agreement with the site-directed mutagenesis data and other studies.