Crystal structure analysis,covalent docking,and molecular dynamics calculations reveal a conformational switch in PhaZ7 PHB depolymerase

An open and a closed conformation of a surface loop in PhaZ7 extracellular poly(3‐hydroxybutyrate) depolymerase were identified in two high‐resolution crystal structures of a PhaZ7 Y105E mutant. Molecular dynamics (MD) simulations revealed high root mean square fluctuations (RMSF) of the 281–295 loop, in particular at residue Asp289 (RMSF 7.62 Å). Covalent docking between a 3‐hydroxybutyric acid trimer and the catalytic residue Ser136 showed that the binding energy of the substrate is significantly more favorable in the open loop conformation compared to that in the closed loop conformation. MD simulations with the substrate covalently bound depicted 1 Å RMSF higher values for the residues 281–295 in comparison to the apo (substrate‐free) form. In addition, the presence of the substrate in the active site enhanced the ability of the loop to adopt a closed form. Taken together, the analysis suggests that the flexible loop 281–295 of PhaZ7 depolymerase can act as a lid domain to control substrate access to the active site of the enzyme. Proteins 2017; 85:1351–1361. © 2017 Wiley Periodicals, Inc.     

The primary structure of the β-subunit of protocatechuate 3,4-dioxygenase from Pseudomonas aeruginosa

The complete amino acid sequence of the β-subunit of protocatechuate 3,4-dioxygenase was determined. The β-subunit contained four methionine residues. Thus, five peptides were obtained after cleavage of the carboxymethylated β-subunit with cyanogen bromide, and were isolated on Sephadex G-75 column chromatography. The amino acid sequences of the cyanogen bromide peptides were established by characterization of the peptides obtained after digestion with trypsin, chymotrypsin, thermolysin, or Staphylococcus aureus protease. The major sequencing techniques used were automated and manual Edman degradations. The five cyanogen bromide peptides were aligned by means of the amino acid sequences of the peptides containing methionine purified from the tryptic hydrolysate of the carboxymethylated β-subunit. The amino acid sequence of all the 238 residues was as follows: ProAlaGlnAspAsnSerArgPheValIleArgAsp ArgAsnTrpHis ProLysAlaLeuThrPro-Asp — TyrLysThrSerIleAlaArg SerProArgGlnAla LeuValSerIleProGlnSer — IleSerGluThrThrGly ProAsnPheSerHisLeu GlyPheGlyAlaHisAsp-His — AspLeuLeuLeuAsnPheAsn AsnGlyGlyLeu ProIleGlyGluArgIle-Ile — ValAlaGlyArgValValAsp GlnTyrGlyLysPro ValProAsnThrLeuValGluMet — TrpGlnAlaAsnAla GlyGlyArgTyrArg HisLysAsnAspArgTyrLeuAlaPro — LeuAspProAsn PheGlyGlyValGly ArgCysLeuThrAspSerAspGlyTyrTyr — SerPheArg ThrIleLysProGlyPro TyrProTrpArgAsnGlyProAsnAsp — TrpArgProAla HisIleHisPheGlyIle SerGlyProSerIleAlaThr-Lys — LeuIleThrGlnLeuTyr PheGluGlyAspPro LeuIleProMetCysProIleVal — LysSerIleAlaAsn ProGluAlaValGlnGln LeuIleAlaLysLeuAspMetAsnAsn — AlaAsnProMet AsnCysLeuAlaTyr ArgPheAspIleValLeuArgGlyGlnArgLysThrHis PheGluAsnCys. The sequence published earlier in summary form (Iwaki et al., 1979, J. Biochem.86, 1159–1162) contained a few errors which are pointed out in this paper.     

The primary structure of the COOH-terminal half of cholera toxin subunit A1 containing the ADP-ribosylation site

The sequence of 96 amino acid residues from the COOH-terminus of the active subunit of cholera toxin, A₁, has been determined as PheAsnValAsnAspVal LeuGlyAlaTyrAlaProHisProAsxGluGlu GluValSerAlaLeuGlyGly IleProTyrSerGluIleTyrGlyTrpTyrArg ValHisPheGlyValLeuAsp GluGluLeuHisArgGlyTyrArgAspArgTyr TyrSerAsnLeuAspIleAla ProAlaAlaAspGlyTyrGlyLeuAlaGlyPhe ProProGluHisArgAlaTrp ArgGluGluProTrpIleHisHisAlaPro ProGlyCysGlyAsnAlaProArg(OH). This is the largest fragment obtained by BrCN cleavage of the subunit A₁ (M_r 23,000), and has previously been indicated to contain the active site for the adenylate cyclase-stimulating activity. Unequivocal identification of the COOH-terminal structure was achieved by separation and analysis of the terminal peptide after the specific chemical cleavage at the only cysteine residue in A₁ polypeptide. The site of self ADP-ribosylation in the A₁ subunit [C. Y. Lai, Q.-C. Xia, and P. T. Salotra (1983) Biochem. Biophys. Res. Commun.116, 341–348] has now been identified as Arg-50 of this peptide, 46 residues removed from the COOH-terminus. The cysteine that forms disulfide bridge to A₂ subunit in the holotoxin is at position 91.     

Sequences of tryptic peptides containing the five cysteinyl residues of ribulosebisphosphate carboxylase/oxygenase from Rhodospirillum rubrum

Tryptic peptides which account for all five cysteinyl residues in ribulosebisphosphate carboxylase/oxygenase from Rhodospirillum rubrum have been purified and sequenced. Collectively, these peptides contain 94 of the approximately 500 amino acid residues per molecule of subunit. Due to one incomplete cleavage at a site for trypsin and two incomplete chymotryptic-like cleavages, eight major radioactive peptides (rather than five as predicted) were recovered from tryptic digests of the enzyme that had been carboxymethylated with [³H]iodoacetate. The established sequences are: GlyTyrThrAlaPheValHisCys¹Lys TyrValAspLeuAlaLeuLysGluGluAspLeuIleAla GlyGlyGluHisValLeuCys¹AlaTyr AlaGlyTyrGlyTyrValAlaThrAlaAlaHisPheAla AlaGluSerSerThrGlyThrAspValGluValCys¹ ThrThrAsxAsxPheThrArg AlaCys¹ThrProIleIleSerGlyGlyMetAsnAla LeuArg ProPheAlaGluAlaCys¹HisAlaPheTrpLeuGly GlyAsnPheIleLys In these peptides, radioactive carboxymethylcysteinyl residues are denoted with asterisks and the sites of incomplete cleavage with vertical wavy lines. None of the peptides appear homologous with either of two cysteinyl-containing, active-site peptides previously isolated from spinach ribulosebisphosphate carboxylase/oxygenase.     

Interaction Between Umami Peptide and Taste Receptor T1R1/T1R3

The umami taste receptor is a heterodimer composed of two members of the T1R taste receptor family: T1R1 and T1R3. The homology models of the ligand binding domains of the human umami receptor have been constructed based on crystallographic structures of the taste receptor of the central nervous system. Furthermore, the molecular simulations of the ligand binding domain show that the likely conformation was that T1R1 protein exists in the closed conformation, and T1R3 in the open conformation in the heterodimer. The molecular docking study of T1R1 and T1R3 in complex with four peptides, including Lys–Gly–Asp–Glu–Ser–Leu–Leu–Ala, Ser–Glu–Glu, G1u–Ser, and Asp–Glu–Ser, displayed that the amino acid residue of SER146 and Glu277 in T1R3 may play great roles in the synergism of umami taste. This docking result further validated the robustness of the model. In the paper, binding of umami peptide and the T1R1/T1R3 receptor was first described and the interaction is the base of umami activity theory.     

The Influence of Salinity Acclimation on Free Amino Acids and Enzyme Activities in the Intestinal Mucosa of Rainbow Trout,Oncorhynchus mykiss (Walbaum)

Postprandial changes of Arg, Leu, Val, Ala, Asp, Glu, Gly, Pro and Tau as well as activities of three enzymes of the transdeamination system in the midgut mucosa and, for comparison, in the liver of freshwater and seawater acclimated Oncorhynchus mykiss were studied. In the mucosa a postprandial increase of Arg, Leu, Val, Ala, Asp, Glu, Gly and Pro occurred. In contrast, only the postprandial Arg level increased strongly in the liver. Levels of Leu, Val, Ala, Asp, Glu, Gly, Pro and Tau remained stable. Concentrations of Ala, Asp, Glu and Pro are higher in the liver than the mucosa. Tau is the most important osmotic effector in both organs, but its concentration is much lower in the liver. Its postprandial concentrations remained stable in both tissues but were significantly higher in seawater trout. The trend of a stronger postprandial rise of Arg, Leu, Val, Ala, Asp, Glu, Gly and Pro levels in seawater trout than in freshwater trout was shown. In mucosa tissue aspartate aminotransferase activities were higher in seawater trout. Ratios of aspartate aminotransferase, alanine aminotransferase and glutamate dehydrogenase are similar to those of the gills.     

The amino acid sequence of human chorionic gonadotropin. The alpha subunit and beta subunit.

The amino acid sequences of both the alpha and beta subunits of human chorionic gonadotropin have been determined. The amino acid sequence of the alpha subunit is: Ala - Asp - Val - Gln - Asp - Cys - Pro - Glu - Cys-10 - Thr - Leu - Gln - Asp - Pro - Phe - Ser - Gln-20 - Pro - Gly - Ala - Pro - Ile - Leu - Gln - Cys - Met - Gly-30 - Cys - Cys - Phe - Ser - Arg - Ala - Tyr - Pro - Thr - Pro-40 - Leu - Arg - Ser - Lys - Lys - Thr - Met - Leu - Val - Gln-50 - Lys - Asn - Val - Thr - Ser - Glu - Ser - Thr - Cys - Cys-60 - Val - Ala - Lys - Ser - Thr - Asn - Arg - Val - Thr - Val-70 - Met - Gly - Gly - Phe - Lys - Val - Glu - Asn - His - Thr-80 - Ala - Cys - His - Cys - Ser - Thr - Cys - Tyr - Tyr - His-90 - Lys - Ser. Oligosaccharide side chains are attached at residues 52 and 78. In the preparations studied approximately 10 and 30% of the chains lack the initial 2 and 3 NH2-terminal residues, respectively. This sequence is almost identical with that of human luteinizing hormone (Sairam, M. R., Papkoff, H., and Li, C. H. (1972) Biochem. Biophys. Res. Commun. 48, 530-537). The amino acid sequence of the beta subunit is: Ser - Lys - Glu - Pro - Leu - Arg - Pro - Arg - Cys - Arg-10 - Pro - Ile - Asn - Ala - Thr - Leu - Ala - Val - Glu - Lys-20 - Glu - Gly - Cys - Pro - Val - Cys - Ile - Thr - Val - Asn-30 - Thr - Thr - Ile - Cys - Ala - Gly - Tyr - Cys - Pro - Thr-40 - Met - Thr - Arg - Val - Leu - Gln - Gly - Val - Leu - Pro-50 - Ala - Leu - Pro - Gin - Val - Val - Cys - Asn - Tyr - Arg-60 - Asp - Val - Arg - Phe - Glu - Ser - Ile - Arg - Leu - Pro-70 - Gly - Cys - Pro - Arg - Gly - Val - Asn - Pro - Val - Val-80 - Ser - Tyr - Ala - Val - Ala - Leu - Ser - Cys - Gln - Cys-90 - Ala - Leu - Cys - Arg - Arg - Ser - Thr - Thr - Asp - Cys-100 - Gly - Gly - Pro - Lys - Asp - His - Pro - Leu - Thr - Cys-110 - Asp - Asp - Pro - Arg - Phe - Gln - Asp - Ser - Ser - Ser - Ser - Lys - Ala - Pro - Pro - Pro - Ser - Leu - Pro - Ser-130 - Pro - Ser - Arg - Leu - Pro - Gly - Pro - Ser - Asp - Thr-140 - Pro - Ile - Leu - Pro - Gln. Oligosaccharide side chains are found at residues 13, 30, 121, 127, 132, and 138. The proteolytic enzyme, thrombin, which appears to cleave a limited number of arginyl bonds, proved helpful in the determination of the beta sequence.     

Amino acid sequence of a protease inhibitor isolated from Sarcophaga bullata determined by mass spectrometry.

The amino acid sequence of a protease inhibitor isolated from the hemolymph of Sarcophaga bullata larvae was determined by tandem mass spectrometry. Homology considerations with respect to other protease inhibitors with known primary structures assisted in the choice of the procedure followed in the sequence determination and in the alignment of the various peptides obtained from specific chemical cleavage at cysteines and enzyme digests of the S. bullata protease inhibitor. The resulting sequence of 57 residues is as follows: Val Asp Lys Ser Ala Cys Leu Gln Pro Lys Glu Val Gly Pro Cys Arg Lys Ser Asp Phe Val Phe Phe Tyr Asn Ala Asp Thr Lys Ala Cys Glu Glu Phe Leu Tyr Gly Gly Cys Arg Gly Asn Asp Asn Arg Phe Asn Thr Lys Glu Glu Cys Glu Lys Leu Cys Leu.     

DNA polymerase beta catalytic efficiency mirrors the Asn279-dCTP H-bonding strength

Ternary complexes of wild type or mutant form of human DNA polymerase beta (pol beta) bound to DNA and dCTP substrates were studied by molecular dynamics (MD) simulations. The occurrences of contact configurations (CC) of structurally important atom pairs were sampled along the MD trajectories, and converted into free-energy differences, DeltaG(CC). DeltaG(CC) values were correlated with the experimental binding and catalytic free energies for the wild type pol beta and its Arg183Ala, Tyr271Ala, Asp276Val, Lys280Gly, Arg283Ala, and Glu295Ala mutants. The correlation coefficients show that the strength of the H-bond between dCTP and Asn279 is a strong predictor of the mutation-induced changes in the catalytic efficiency of pol beta. This finding is consistent with the view that enzyme preorganization plays a major role in controlling DNA polymerase specific activity.     

The covalent structure of pig kidney fructose 1,6-bisphosphatase: sequence of the 60-residue NH2-terminal peptide produced by digestion with subtilisin

Digestion of the native pig kidney fructose 1,6-bisphosphatase tetramer with subtilisin cleaves each of the 35,000-molecular-weight subunits to yield two major fragments: the S-subunit (M_{r ca.} 29,000), and the S-peptide (M_r 6,500). The following amino acid sequence has been determined for the S peptide: AcThrAspGlnAlaAlaPheAspThrAsnIle Val ThrLeuThrArgPheValMetGluGlnGlyArgLysAla ArgGlyThrGlyGlu MetThrGlnLeuLeuAsnSerLeuCysThrAlaValLys AlaIleSerThrAla z.sbnd;ValArgLysAlaGlyIleAlaHisLeuTyrGlyIleAla. Comparison of this sequence with that of the NH₂-terminal 60 residues of the enzyme from rabbit liver (El-Dorry et al., 1977, Arch. Biochem. Biophys.182, 763) reveals strong homology with 52 identical positions and absolute identity in sequence from residues 26 to 60.Although subtilisin cleavage of fructose 1,6-bisphosphatase results in diminished sensitivity of the enzyme to AMP inhibition, we have found no AMP inhibition-related amino acid residues in the sequenced S-peptide. The loss of AMP sensitivity that occurs upon pyridoxal-P modification of the enzyme does not result in the modification of lysyl residues in the S-peptide. Neither photoaffinity labeling of fructose 1,6-bisphosphatase with 8-azido-AMP nor modification of the cysteinyl residue proximal to the AMP allosteric site resulted in the modification of residues located in the NH₂-terminal 60-amino acid peptide.     

Probing the affinity and specificity of yeast alcohol dehydrogenase I for coenzymes.

Yeast (Saccharomyces cerevisiae) alcohol dehydrogenase I (SceADH) binds NAD+ and NADH less tightly and turns over substrates more rapidly than does horse (Equus caballus) liver alcohol dehydrogenase E isoenzyme (EcaADH), and neither enzyme uses NADP efficiently. Amino acid residues in the proposed adenylate binding pocket of SceADH were substituted in attempts to improve affinity for coenzymes or reactivity with NADP. Substitutions in SceADH (Gly202Ile or Ser246Ile) with the corresponding residues in the adenine binding site of the homologous EcaADH have modest effects on coenzyme binding and other kinetic constants, but the Ser246Ile substitution decreases turnover numbers by 350-fold. The Ser176Phe substitution (also near adenine site) significantly decreases affinity for coenzymes and turnover numbers. In the consensus nucleotide-binding betaalphabeta fold sequence, SceADH has two alanine residues (177-GAAGGLG-183) instead of the Leu200 in EcaADH (199-GLGGVG-204); the Ala178-Ala179 to Leu substitution significantly decreases affinity for coenzymes and turnover numbers. Some NADP-dependent enzymes have an Ala corresponding to Gly183 in SceADH; the Gly183Ala substitution significantly decreases affinity for coenzymes and turnover numbers. NADP-dependent enzymes usually have a neutral residue instead of the Asp (Asp201 in SceADH) that interacts with the hydroxyl groups of the adenosine ribose, along with a basic residue (at position 202 or 203) to stabilize the 2'-phosphate of NADP. The Gly203Arg change in SceADH does not significantly affect the kinetics. The Gly183Ala or Gly203Arg substitutions do not enable SceADH to use NADP+ as coenzyme. SceADH with the single Asp201Gly or double Asp201Gly:Gly203Arg substitutions have similar, low activity with NADP+. The results suggest that several of the amino acid residues participate in coenzyme binding and that conversion of specificity for coenzyme requires multiple substitutions.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.     

Structural insight into epothilones antitumor activity based on the conformational preferences and tubulin binding modes of epothilones A and B obtained from molecular dynamics simulations

Molecular dynamics simulations were employed to analyze the conformational preferences and binding modes of epothilones A and B as a source of structural information regarding the antitumor properties of these species. Our results suggest that the conformation of free and tubulin-bound epothilones is strongly influenced by the presence of a methyl group at C12 and that epothilones A and B exploit the binding cavity in a unique and different way. The binding sites of epothilones A and B share a common region of association (Leu215, Leu217, His227, Leu228, Ala231, Phe270, Gly360, and Leu361), but lead to different ligand–residue interactions. Average interaction energies predict a larger stabilization for the epothilone B–tubulin complex, which is mainly driven by the enhancement of the electrostatic component of ligand–residue interactions compared to the epothilone A–tubulin complex. These structural and energetic results can be useful to account for the activity difference between epothilones A and B, and to design more active and potent analogs that resemble the mechanism of action of epothilones against cancer cells.     

Loop propensity of the sequence YKGQP from staphylococcal nuclease: implications for the folding of nuclease.

Recently we performed molecular dynamics (MD) simulations on the folding of the hairpin peptide DTVKLMYKGQPMTFR from staphylococcal nuclease in explicit water. We found that the peptide folds into a hairpin conformation with native and nonnative hydrogen-bonding patterns. In all the folding events observed in the folding of the hairpin peptide, loop formation involving the region YKGQP was an important event. In order to trace the origins of the loop propensity of the sequence YKGQP, we performed MD simulations on the sequence starting from extended, polyproline II and native type I' turn conformations for a total simulation length of 300 ns, using the GROMOS96 force field under constant volume and temperature (NVT) conditions. The free-energy landscape of the peptide YKGQP shows minima corresponding to loop conformation with Tyr and Pro side-chain association, turn and extended conformational forms, with modest free-energy barriers separating the minima. To elucidate the role of Gly in facilitating loop formation, we also performed MD simulations of the mutated peptide YKAQP (Gly --> Ala mutation) under similar conditions starting from polyproline II conformation for 100 ns. Two minima corresponding to bend/turn and extended conformations were observed in the free-energy landscape for the peptide YKAQP. The free-energy barrier between the minima in the free-energy landscape of the peptide YKAQP was also modest. Loop conformation is largely sampled by the YKGQP peptide, while extended conformation is largely sampled by the YKAQP peptide. We also explain why the YKGQP sequence samples type II turn conformation in these simulations, whereas the sequence as part of the hairpin peptide DTVKLMYKGQPMTFR samples type I' turn conformation both in the X-ray crystal structure and in our earlier simulations on the folding of the hairpin peptide. We discuss the implications of our results to the folding of the staphylococcal nuclease.     

Importance of stalk segment S5 for intramolecular communication in the sarcoplasmic reticulum Ca2+-ATPase

Sixteen residues in stalk segment S5 of the Ca(2+)-ATPase of sarcoplasmic reticulum were studied by site-directed mutagenesis. The rate of the Ca(2+) binding transition, determined at 0 degrees C, was enhanced relative to wild type in mutants Ile(743) --> Ala, Val(747) --> Ala, Glu(748) --> Ala, Glu(749) --> Ala, Met(757) --> Gly, and Gln(759) --> Ala and reduced in mutants Asp(737) --> Ala, Asp(738) --> Ala, Ala(752) --> Leu, and Tyr(754) --> Ala. In mutant Arg(762) --> Ile, the rate of the Ca(2+) binding transition was wild type like at 0 degrees C, whereas it was 3.5-fold reduced relative to wild type at 25 degrees C. The rate of dephosphorylation of the ADP-insensitive phosphoenzyme was increased conspicuously in mutants Ile(743) --> Ala and Tyr(754) --> Ala (close to 20-fold in the absence of K(+)) and increased to a lesser extent in Asn(739) --> Ala, Glu(749) --> Ala, Gly(750) --> Ala, Ala(752) --> Gly, Met(757) --> Gly, and Arg(762) --> Ile, whereas it was reduced in mutants Asp(737) --> Ala, Val(744) --> Gly, Val(744) --> Ala, Val(747) --> Ala, and Ala(752) --> Leu. In mutants Ile(743) --> Ala, Tyr(754) --> Ala, and Arg(762) --> Ile, the apparent affinities for vanadate were enhanced 23-, 30-, and 18-fold, respectively, relative to wild type. The rate of Ca(2+) dissociation was 11-fold increased in Gly(750) --> Ala and 2-fold reduced in Val(747) --> Ala. Mutants with alterations to Arg(751) either were not expressed at a significant level or were completely nonfunctional. The findings show that S5 plays a crucial role in mediating communication between the Ca(2+) binding pocket and the catalytic domain and that Arg(751) is important for both structural and functional integrity of the enzyme.     

Effects of protonation state of Asp181 and position of active site water molecules on the conformation of PTP1B

In protein tyrosine phosphatase 1B (PTP1B), the flexible WPD loop adopts a closed conformation (WPD_closed) in the active state of PTP1B, bringing the catalytic Asp181 close to the active site pocket, while WPD loop is in an open conformation (WPD_open) in the inactive state. Previous studies showed that Asp181 may be protonated at physiological pH, and ordered water molecules exist in the active site. In the current study, molecular dynamics simulations are employed at different Asp181 protonation states and initial positions of active site water molecules, and compared with the existing crystallographic data of PTP1B. In WPD_closed conformation, the active site is found to maintain its conformation only in the protonated state of Asp181 in both free and liganded states, while Asp181 is likely to be deprotonated in WPD_open conformation. When the active site water molecule network that is a part of the free WPD_closed crystal structure is disrupted, intermediate WPD loop conformations, similar to that in the PTPRR crystal structure, are sampled in the MD simulations. In liganded PTP1B, one active site water molecule is found to be important for facilitating the orientation of Cys215 and the phosphate ion, thus may play a role in the reaction. In conclusion, conformational stability of WPD loop, and possibly catalytic activity of PTP1B, is significantly affected by the protonation state of Asp181 and position of active site water molecules, showing that these aspects should be taken into consideration both in MD simulations and inhibitor design. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

The NMR-derived conformation of orexin-A: an orphan G-protein coupled receptor agonist involved in appetite regulation and sleep

The conformation of orexin-A, an orphan G-protein coupled receptor agonist has been determined when bound to sodium dodecylsulphate-d(25) (SDS) micelles by (1)H and (13)C NMR and molecular modeling. Orexin-A has been implicated in sleep-wakefulness and feeding regulation. The conformational preference of orexin-A consists of a short helical section, involving Asp(5) to Gln(9) that makes up helix I, followed by a bend from Lys(10) to Ser(13). Residues Leu(16) to Gly(22) make up helix II. The conformation of orexin-A can now be used to explain the results of earlier Ala substitution mutagenesis experiments (J. G. Darker et al., Bioorg. Med. Chem. Lett. 11, 737-740 (2001); S. Ammoun, et al., J. Pharmacol. Expt. Ther. 305, 507-514 (2003)). Darker et al., working with orexin-A (15-33) amide, observed a significant drop in functional potency at the OX(1)R receptor when Leu(16), Leu(19), Leu(20), His(26), Gly(29), Ile(30), Leu(31), Thr(32), and Leu(33) were replaced by Ala. Ammoun et al. identified three areas of interest, which were the same for OX(1)R and OX(2)R receptors, as amino acids 15-17, 20 and 25-26 with the most marked reduction in activity being produced by the replacement of Leu(20) by Ala. We suggest that Leu(16), Leu(19), and Leu(20), which are in helix II, are likely responsible for binding orexin-A to the surface of the micelle.     

In silico mutations and molecular dynamics studies on a winged bean chymotrypsin inhibitor protein

Winged bean chymotrypsin inhibitor (WCI) has an intruding residue Asn14 that plays a crucial role in stabilizing the reactive site loop conformation. This residue is found to be conserved in the Kunitz (STI) family of serine protease inhibitors. To understand the contribution of this scaffolding residue on the stability of the reactive site loop, it was mutated in silico to Gly, Ala, Ser, Thr, Leu and Val and molecular dynamics (MD) simulations were carried out on the mutants. The results of MD simulations reveal the conformational variability and range of motions possible for the reactive site loop of different mutants. The N-terminus side of the scissile bond, which is close to a beta-barrel, is conformationally less variable, while the C-terminus side, which is relatively far from any such secondary structural element, is more variable and needs stability through hydrogen-bonding interactions. The simulated structures of WCI and the mutants were docked in the peptide-binding groove of the cognate enzyme chymotrypsin and the ability to form standard hydrogen-bonding interactions at P3, P1 and P2' residues were compared. The results of the MD simulations coupled with docking studies indicate that hydrophobic residues like Leu and Val at the 14th position are disruptive for the integrity of the reactive site loop, whereas a residue like Thr, which can stabilize the C-terminus side of the scissile bond, can be predicted at this position. However, the size and charge of the Asn residue made it most suitable for the best maintenance of the integrity of the reactive site loop, explaining its conserved nature in the family.     

Molecular modeling of SWR-0342SA, a beta3-selective agonist, with beta1- and beta3-adrenoceptor

The molecular dynamics (MD) simulations study in the formation of the complex between compound SWR-0342SA and beta-ARs suggested that upon binding SWR-0342SA stimulates receptor activation through residues network (Asp104, Leu335 in beta(1)-AR; Asp117, Ser209, Leu303, Ser191 in beta(3)-AR) in an active conformation state. The models suggest that the structural origin of the selectivity of SWR-0342SA to beta(3)-AR vs. beta(1)-AR comes from the following results: (a) the tight interaction between the agonist and the TMs 3, 5, 6 and 2 nd EC loop. Asp117 interacts with the cationic amino group of the agonist molecule. (b) Additional contacts are done with Ser209, Leu303 and Ser191. These results are in good agreement with the binding affinities (pKi values) of SWR-0342SA to beta-AR family expressed in recombinant mammalian cells.     

Vibrational circular dichroism and IR spectral analysis as a test of theoretical conformational modeling for a cyclic hexapeptide

A model cyclohexapeptide, cyclo-(Phe-(D)Pro-Gly-Arg-Gly-Asp) was synthesized and its IR and VCD spectra were used as a test of density functional theory (DFT) level predictions of spectral intensities for a peptide with a nonrepeating but partially constricted conformation. Peptide structure and flexibility was estimated by molecular dynamics (MD) simulations and the spectra were simulated using full quantum mechanical (QM) approaches for the complete peptide and for simplified models with truncated side chains. After simulated annealing, the backbone conformation of the ring structure is relatively stable, consisting of a normal beta-turn and a tight loop (no H-bond) which does not vary over short trajectories. Only in quite long MD runs at high temperatures do other conformations appear. MD simulations were carried out for the cyclic peptide in water and in TFE, which match experimental solvents, as well as with and without protonation of the Asp carboxyl group. DFT spectral simulations were made using the annealed structure and were extended to include basis set variation, to determine an optimal computational approach, and solvent simulation with a polarized continuum model (PCM). Stepwise full DFT simulation of spectra was done for various sequences with the same backbone geometry but based on (1) solely Gly residues, (2) Ala substitution except Gly and Pro, and (3) complete sequences with side chains. Additionally, a selection of structures was used to compute IR and VCD spectra with the optimal method to determine structural variation effects. The side chains, especially the Asp-COOH and Arg-NH(2) transitions, had an impact on the computed amide frequencies, IR intensities and VCD pattern. Since experimentally these groups would have little chirality, due to conformational variation, they do not impact the observed VCD spectra. Correcting for frequency shifts, the Ala model for the cyclopeptide gives the clearest representation of the amide VCD. The experimental sign pattern for the amide I' band in D(2)O and also the sharper, more intense amide I VCD band in TFE was seen to some degree in one conformer with Type II' turns, but the data favor a mix of structures.