1
|
Park HS, Byun BJ, Kang YK. Exploring Conformational Preferences of Leu-enkephalin Using the Conformational Search and Double-Hybrid DFT Energy Calculations. ACS OMEGA 2022; 7:27755-27768. [PMID: 35967045 PMCID: PMC9366962 DOI: 10.1021/acsomega.2c03942] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/23/2022] [Accepted: 07/15/2022] [Indexed: 06/15/2023]
Abstract
The conformational preferences of Leu-enkephalin (Leu-Enk) were explored by the conformational search and density functional theory (DFT) calculations. By a combination of low-energy conformers of each residue, the initial structures of the neutral Leu-Enk were generated and optimized using the ECEPP3 force field in the gas phase. These structures were reoptimized at the HF/3-21G(d) and M06-2X levels of theory with 6-31G(d) and 6-31+G(d) basis functions. We finally located the 139 structures with the relative energy <10 kcal mol-1 in the gas phase, from which the structures of the corresponding zwitterionic Leu-Enk were generated and reoptimized at the M06-2X/6-31+G(d) level of theory using the implicit solvation model based on density (SMD) in water. The conformational preferences of Leu-Enk were analyzed using Gibbs free energies corrected by single-point energies calculated at the double-hybrid DSD-PBEP86-D3BJ/def2-TZVP level of theory in the gas phase and in water. The neutral Leu-Enk dominantly adopted a folded structure in the gas phase stabilized by three H-bonds with a βII'-bend-like motif at the Gly3-Phe4 sequence and a close contact between the side chains of Phe4 and Leu5. The zwitterionic Leu-Enk exhibited a folded structure in water stabilized by three H-bonds with double β-bends such as a βII' bend at the Gly2-Gly3 sequence and a βI bend at the Gly3-Phe4 sequence. The calculated ensemble-averaged distance between CGly2 α and CLeu5 α of the zwitterionic Leu-Enk in water is consistent with the value estimated from the simulated annealing using the distance constraints derived from nuclear Overhauser effect spectroscopy (NOESY) spectra in water. Interestingly, the preferred conformations of the neutral and zwitterionic Leu-Enk are new folded structures not predicted by earlier computational studies. According to the refined model of the zwitterionic Leu-Enk bound to δ-opioid receptor (δOR), there were favorable interactions of the terminal charged groups of Leu-Enk with the side chains of charged residues of δOR as well as a favorable CAryl···H interaction of the Phe4 residue of Leu-Enk with Trp284 of δOR. Hence, these favorable interactions would induce the folded structure of the zwitterionic Leu-Enk with double β-bends isolated in water into the "bioactive conformation" like an extended structure when binding to δOR.
Collapse
Affiliation(s)
- Hae Sook Park
- Department
of Nursing, Cheju Halla University, 38 Halladaehak-ro, Jeju, Jeju-do 63092, Republic of Korea
| | - Byung Jin Byun
- Drug
Discovery Center, JW Pharmaceutical Co.
Ltd., 2477 Nambusunhwan-ro, Seocho-gu, Seoul 06725, Republic
of Korea
| | - Young Kee Kang
- Department
of Chemistry, Chungbuk National University, 1 Chungdae-ro,
Seowon-gu, Cheongju, Chungbuk 28644, Republic of Korea
| |
Collapse
|
2
|
Capability of ganglioside GM1 in modulating interactions, structure, location and dynamics of peptides/proteins: biophysical approaches. Glycoconj J 2014; 31:435-47. [DOI: 10.1007/s10719-014-9554-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
3
|
Unruh JR, Kuczera K, Johnson CK. Conformational heterogeneity of a leucine enkephalin analogue in aqueous solution and sodium dodecyl sulfate micelles: comparison of time-resolved FRET and molecular dynamics simulations. J Phys Chem B 2009; 113:14381-92. [PMID: 19780516 PMCID: PMC2780025 DOI: 10.1021/jp903302k] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We have undertaken time-resolved Förster resonance energy transfer (FRET) and molecular dynamics simulations to analyze conformations and conformational heterogeneity of an analogue of leucine enkephalin in solution and in the presence of sodium dodecyl sulfate (SDS) micelles. Enkephalins are opioid pentapeptides that interact with opioid receptors in the central nervous system. We used time-correlated single-photon counting to detect energy transfer between the N-terminal tyrosine and a tryptophan residue substituted for phenylalanine at the 4 position. FRET from Tyr to Trp was measured over a temperature range from 5 to 55 degrees C in aqueous solution. By taking into account Tyr rotamer interconversion rates measured previously, we determined average distances between Tyr and Trp for the two populated rotameric conformations of Tyr. Molecular dynamics simulations (100 ns) support this analysis and indicate extensive conformational heterogeneity. The simulations also predict that the FRET orientational factor is correlated with the Tyr-Trp separation. Failure to account for the correlation between orientation and distance results in errors that appear to be largely offset in the leucine enkephalin analogue (YGGWL) by a weighting bias inherent in the R(-6) dependence of the energy-transfer rate. The Tyr lifetimes decrease upon titration of the peptides with SDS, indicating formation of compact conformations of the peptide in the micelle environment. This result is consistent with the conjecture that the lipid environment may induce formation of bioactive conformations of the peptide.
Collapse
Affiliation(s)
- Jay R. Unruh
- Department of Chemistry, University of Kansas, Lawrence, KS 66045
| | - Krzysztof Kuczera
- Department of Chemistry, University of Kansas, Lawrence, KS 66045
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66045
| | - Carey K. Johnson
- Department of Chemistry, University of Kansas, Lawrence, KS 66045
| |
Collapse
|
4
|
Gayen A, Mukhopadhyay C. Evidence for effect of GM1 on opioid peptide conformation: NMR study on leucine enkephalin in ganglioside-containing isotropic phospholipid bicelles. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2008; 24:5422-5432. [PMID: 18412380 DOI: 10.1021/la704056d] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
Enkephalins are endogenous neuropeptides that have opioid-like activities and compete with morphines for the receptor binding. The binding of these neuropeptides to membrane appears crucial since enkephalins interact with the nerve cell membranes to achieve bioactive conformations that fit onto multiple receptor sites (micro, delta, and kappa). Using NMR spectroscopy, we have determined the solution structure of the small opiate pentapeptide leucine enkephalin in the presence of isotropic phospholipid bicelles: phosphocholine bicelles (DMPC:CHAPS 1:4) and phosphocholine bicelles doped with ganglioside GM1 (DMPC:CHAPS:GM1 1:4:0.3). Bicelles containing GM1 were found to interact strongly with leucine enkephalin, whereas a somewhat weaker interaction was observed in the case of bicelles without GM1. Structure calculation from torsion angles, chemical shifts, and NOE-based distance constraints explored that the peptide could flexibly switch between several mu- and delta-selective conformations in both the bicelles though micro-selective conformations turned out to be geometrically preferred in each bicellar system. A detailed analysis of the structures presented supports the variance over the singly associated conformation of enkephalin in nerve cell membranes.
Collapse
Affiliation(s)
- Anindita Gayen
- Department of Chemistry, University of Calcutta, Kolkata 700 009, India
| | | |
Collapse
|
5
|
Fujitani N, Shimizu H, Matsubara T, Ohta T, Komata Y, Miura N, Sato T, Nishimura SI. Structural transition of a 15 amino acid residue peptide induced by GM1. Carbohydr Res 2007; 342:1895-903. [PMID: 17572397 DOI: 10.1016/j.carres.2007.05.026] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2007] [Revised: 05/11/2007] [Accepted: 05/15/2007] [Indexed: 11/15/2022]
Abstract
The ganglioside GM1-binding peptide, p3, with a sequence of VWRLLAPPFSNRLLP, displayed a clear structural alteration depending on the presence or absence of GM1 micelles. The three-dimensional structures of the p3 peptide in the free and GM1 bound states were analyzed using two-dimensional NMR spectroscopic experiments with distance-restrained simulated annealing calculations. The NMR experiments for the p3 peptide alone indicated that the peptide has two conformers derived from the exchange of cis and trans forms at Pro(7)-Pro(8). Further study with theoretical modeling revealed that the p3 peptide has a curb conformation without regular secondary structure. On the other hand, the NMR studies for the p3 peptide with the GM1 micelles elucidated a trans conformer and gave a structure stabilized by hydrophobic interactions of beta- and helical turns. Based on these structural investigations, tryptophan, a core residue of the hydrophobic cluster, might be an essential residue for the recognition of the GM1 saccharides. The dynamic transition of the p3 peptide may play an important role in the function of GM1 as a multiple receptor as in the traditional pathway of the infection by cholera toxin.
Collapse
Affiliation(s)
- Naoki Fujitani
- Division of Advanced Chemical Biology, Graduate School of Advanced Life Science, Frontier Research Center for Post-Genomic Science and Technology, Hokkaido University, Sapporo 011-0021, Japan
| | | | | | | | | | | | | | | |
Collapse
|
6
|
Begotka BA, Hunsader JL, Oparaeche C, Vincent JK, Morris KF. A pulsed field gradient NMR diffusion investigation of enkephalin peptide-sodium dodecyl sulfate micelle association. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2006; 44:586-93. [PMID: 16552718 DOI: 10.1002/mrc.1793] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/07/2023]
Abstract
Pulsed field gradient NMR (PFG-NMR) diffusion experiments were used to investigate the binding of leucine and methionine enkephalin peptides to anionic sodium dodecyl sulfate (SDS) micelles. The study was undertaken to investigate the mechanism of interaction between enkephalin peptides and SDS micelles and to determine if NMR-derived association constants, K(eq), can predict the elution order in electrokinetic chromatography (EKC). In EKC, peptides are separated on the basis of their interactions with micelles. The Leu-enkephalin peptide-micelle association constant increased from 130 +/- 8 to 1459 +/- 57 and 1744 +/- 64 M(-1), respectively, when an Arg or Lys was added to the C-terminus. The association constant of Leu-enkephalinamide was approximately equal to that of Leu-enkephalin-Arg. Substitution of Phe4 with a Trp or Gly2 with an Ala in the Leu-enkephalin peptides also increased the micelle binding affinity. These results confirm that the interaction of Leu-enkephalin peptides with SDS micelles is largely electrostatic and that the non-polar amino acid side chains interact with the hydrophobic micelle core. The peptide-micelle association constants for the cationic Met-enkephalin peptides were also greater than their zwitterionic counterparts. For example, the Met-enkephalin K(eq) value was 162 +/- 9 M(-1), while the association constants for Met-enkephalin-Arg, Met-enkephalin-Lys, and Met-enkephalinamide were, respectively, 674 +/- 31, 426 +/- 23, and 453 +/- 27 M(-1). In both Met-enkephalin and Met-enkephalinamide, replacing Gly2 with an Ala did not significantly increase the association constant. These results confirm that with the Met-enkephalin peptides, there was little or no interaction of the amino acid side chains with the micelle core. For both the Leu-enkephalin and Met-enkephalin peptides, the association constants were consistent with EKC results, in that the peptides with smaller K(eq) values were found to elute before those with larger association constants.
Collapse
|
7
|
Chatterjee C, Mukhopadhyay C. Interaction and structural study of kinin peptide bradykinin and ganglioside monosialylated 1 micelle. Biopolymers 2005; 78:197-205. [PMID: 15815977 DOI: 10.1002/bip.20278] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
Partitioning of small proteins and peptides from the aqueous to membrane phase is often coupled with folding. In this work we examine the binding and folding of the kinin peptide, bradykinin (BK), in the presence of the ganglioside monosialylated 1 (GM1) micelle. Using two-dimensional NMR techniques, we have shown that at low concentration, GM1 micelle is able to induce a turn conformation to BK. A pulsed-field gradient diffusion NMR study indicated that the peptide partitions into the GM1 micelle with a DeltaG(part) of -3.14 +/- 0.03 kcal/mol. A saturation transfer difference (STD) NMR study indicated that the binding is mostly through hydrophobic residues.
Collapse
Affiliation(s)
- Chiradip Chatterjee
- Department of Chemistry, University of Calcutta, 92, A. P. C. Road, Kolkata-700 009, India
| | | |
Collapse
|
8
|
Dastidar SG, Mukhopadhyay C. Anomalous behavior of water around sodium dodecyl sulphate micelles. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2004; 70:061901. [PMID: 15697396 DOI: 10.1103/physreve.70.061901] [Citation(s) in RCA: 18] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/26/2004] [Indexed: 05/24/2023]
Abstract
The dynamics, structural properties, and energetics of hydration water around a sodium dodecyl sulphate micelle have been investigated using molecular dynamics simulation. A clear revelation of the slow dynamics of the hydration water has been made by separate measurements of the rotational and translational properties. Calculated diffusion coefficients fall within the range of experimentally observed quantities. The water-micelle head group (MHG) hydrogen bond is more stable (by an amount approximately 7.0 kcal/mol) compared to the water-water hydrogen bond. The difference in stability of the water monomers forming different numbers of hydrogen bonds (n=0,1,2) with the MHG has clearly been shown from the analyses of their rotational relaxation, residence times, as well as the energy of interaction with different components of the system. The singly hydrogen-bonded water species is the most abundant and stable. The entropy plays the key role in controlling the relative abundance of the different species.
Collapse
Affiliation(s)
- Shubhra Ghosh Dastidar
- Department of Chemistry, University of Calcutta, 92. A.P.C. Road, Kolkata 700 009, India
| | | |
Collapse
|
9
|
Chatterjee C, Majumder B, Mukhopadhyay C. Pulsed-Field Gradient and Saturation Transfer Difference NMR Study of Enkephalins in the Ganglioside GM1 Micelle. J Phys Chem B 2004. [DOI: 10.1021/jp037553r] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Chiradip Chatterjee
- Department of Chemistry, University of Calcutta, 92, A. P. C. Road, Kolkata 700 009, India, and NMR Facility, Department of Biophysics, Bose Institute, P1/12, CIT Scheme 7M, Kolkata 700 054, India
| | - Barun Majumder
- Department of Chemistry, University of Calcutta, 92, A. P. C. Road, Kolkata 700 009, India, and NMR Facility, Department of Biophysics, Bose Institute, P1/12, CIT Scheme 7M, Kolkata 700 054, India
| | - Chaitali Mukhopadhyay
- Department of Chemistry, University of Calcutta, 92, A. P. C. Road, Kolkata 700 009, India, and NMR Facility, Department of Biophysics, Bose Institute, P1/12, CIT Scheme 7M, Kolkata 700 054, India
| |
Collapse
|
10
|
Chatterjee C, Mukhopadhyay C. Conformational alteration of bradykinin in presence of GM1 micelle. Biochem Biophys Res Commun 2004; 315:866-71. [PMID: 14985092 DOI: 10.1016/j.bbrc.2004.01.134] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2004] [Indexed: 10/26/2022]
Abstract
We report here the interaction of bradykinin with ganglioside GM1 by circular dichroism, steady-state fluorescence, and one-dimensional 1H NMR spectroscopy. Circular dichroism spectroscopy is indicative of a turn formation of bradykinin backbone in the presence of GM1 micelle. The fluorescence quenching efficiencies of iodide and acrylamide are substantially reduced, indicating a shielding of phenylalanine residue of bradykinin from aqueous environment. Significant line broadening of NMR resonances of bradykinin, suggestive of motional restriction, is observed.
Collapse
Affiliation(s)
- Chiradip Chatterjee
- Department of Chemistry, University of Calcutta 92, A.P.C. Road, Kolkata 700 009, India
| | | |
Collapse
|