1
|
Kuczera K, Szoszkiewicz R, Jas GS. Microscopic effects of proline co-solvent on alanine homopeptide structure, solvation and helix folding dynamics. J Biomol Struct Dyn 2025:1-11. [PMID: 40351163 DOI: 10.1080/07391102.2025.2500681] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2025] [Accepted: 04/25/2025] [Indexed: 05/14/2025]
Abstract
We present a computational investigation to explore the influence of the protective osmolyte proline as a co-solvent on peptide structure and dynamics for a series of alanine-based peptides, (ALA)n of length n = 5, 8, 15, and 21 residues. Applying multi-microsecond molecular dynamics simulations in a 2 M proline solution, we evaluate peptide structure, solvation and helix folding dynamics and compare to behavior in pure water. Proline addition enhances helix content and significantly slows folding and unfolding times, correlating with a 1.9-fold increase in solvent viscosity. Notably, ALA15 helix content increases from 25% to 49% and relaxation time rises from 110 ns to 540 ns in proline relative to water. Microscopic solvation effects of proline include peptide compaction and dehydration, exclusion of proline from the backbone, formation of weak interactions with the ALA methyl sidechains, and strong interactions with water. The differences of these effects on the helix and coil states drive helix stabilization by proline. Low-dimensional kinetic modeling with Optimal Dimensionality Reduction predicts distinct folding mechanisms: shorter peptides (ALA5-ALA15) exhibit direct helix-coil transitions, and only the longest ALA21 follows a more complex folding pathway involving intermediates. Statistically, enhanced stability of hydrogen bonds in the peptide centers and strong correlation between transitions on neighboring residues are shared between water and proline solutions. However, there is a preference for helix initiation at the N-terminus under proline influence. Our analysis describes the molecular mechanisms of how proline modulates peptide behavior, offering atomistic insights into helix stabilization and folding mechanisms mediated by osmolytes.
Collapse
Affiliation(s)
- Krzysztof Kuczera
- Department of Chemistry, The University of Kansas, Lawrence, Kansas, USA
- Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas, USA
| | - Robert Szoszkiewicz
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Warsaw, Poland
| | - Gouri S Jas
- Department of Chemistry, The University of Kansas, Lawrence, Kansas, USA
- Department of Pharmaceutical Chemistry, The University of Kansas, Lawrence, Kansas, USA
- College of Pharmacy, University of Nebraska Medical Center, Omaha, Nebraska, USA
| |
Collapse
|
2
|
Jas GS, Childs EW, Middaugh CR, Kuczera K. Observing reorientation dynamics with Time-Resolved fluorescence and molecular dynamics in varying periodic boundary conditions. J Biomol Struct Dyn 2022; 40:10614-10628. [PMID: 34308794 DOI: 10.1080/07391102.2021.1947894] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
This work presents a combined study of time-resolved fluorescence spectroscopy and all-atom molecular dynamics simulation to investigate periodic boundary conditions' and water models' influence on the orientation dynamics and translational and rotational diffusion of peptides in solution. We have characterized the effects of solvent box size and water model choice on the dynamics of two peptide systems, NATA and WK5. Computationally, translational, and rotational diffusion and internal fluctuations are investigated through all-atom molecular dynamics simulation with two water models and different box sizes. These results are compared with time-resolved fluorescence anisotropy decay (FAD) measurements. The associated time constant and orientation dynamics from FAD measurement along the 1Lb axis provided baseline data to validate molecular dynamics simulation. The modeling results show that diffusion rates vary roughly in inverse proportion to water model viscosity, as one would expect. Corrections for finite box size are significant for translational diffusion and insignificant for rotational diffusion. This study also finds that internal dynamics described by autocorrelation functions and kinetic network models are relatively insensitive to both box size and water model properties. Our observation suggests that different peptide properties respond differently to a change in simulation conditions.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Gouri S Jas
- Department of Pharmaceutical Chemistry, The University of Kansas, Lawrence, KS, USA
| | - Ed W Childs
- Department of Surgery, Morehouse School of Medicine, Atlanta, GA, USA
| | - C Russell Middaugh
- Department of Pharmaceutical Chemistry, The University of Kansas, Lawrence, KS, USA
| | - Krzysztof Kuczera
- Department of Chemistry and Department of Molecular Biosciences, The University of Kansas, Lawrence, KS, USA
| |
Collapse
|
3
|
Jas GS, Vallejo-Calzada R, Johnson CK, Kuczera K. Dynamic elements and kinetics: Most favorable conformations of peptides in solution with measurements and simulations. J Chem Phys 2019; 151:225102. [PMID: 31837693 DOI: 10.1063/1.5131782] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
Small peptides in solution adopt a specific morphology as they function. It is of fundamental interest to examine the structural properties of these small biomolecules in solution and observe how they transition from one conformation to another and form functional structures. In this study, we have examined the structural properties of a simple dipeptide and a five-residue peptide with the application of far-UV circular dichroism (CD) spectroscopy as a function of temperature, fluorescence anisotropy, and all-atom molecular dynamics simulation. Analysis of the temperature dependent CD spectra shows that the simplest dipeptide N-acetyl-tryptophan-amide (NATA) adopts helical, beta sheet, and random coil conformations. At room temperature, NATA is found to have 5% alpha-helical, 37% beta sheet, and 58% random coil conformations. To our knowledge, this type of structural content in a simplest dipeptide has not been observed earlier. The pentapeptide (WK5) is found to have four major secondary structural elements with 8% 310 helix, 14% poly-L-proline II, 8% beta sheet, and 14% turns. A 56% unordered structural population is also present for WK5. The presence of a significant population of 310 helix in a simple pentapeptide is rarely observed. Fluorescence anisotropy decay (FAD) measurements yielded reorientation times of 45 ps for NATA and 120 ps for WK5. The fluorescence anisotropy decay measurements reveal the size differences between the two peptides, NATA and WK5, with possible contributions from differences in shape, interactions with the environment, and conformational dynamics. All-atom molecular dynamics simulations were used to model the structures and motions of these two systems in solution. The predicted structures sampled by both peptides qualitatively agree with the experimental findings. Kinetic modeling with optimal dimensionality reduction suggests that the slowest dynamic processes in the dipeptide involve sidechain transitions occurring on a 1 ns timescale. The kinetics in the pentapeptide monitors the formation of a distorted helical structure from an extended conformation on a timescale of 10 ns. Modeling of the fluorescence anisotropy decay is found to be in good agreement with the measured data and correlates with the main contributions of the measured reorientation times to individual conformers, which we define as dynamic elements. In NATA, the FAD can be well represented as a sum of contributions from representative conformers. This is not the case in WK5, where our analysis suggests the existence of coupling between conformational dynamics and global tumbling. The current study involving detailed experimental measurements and atomically detailed modeling reveals the existence of specific secondary structural elements and novel dynamical features even in the simplest peptide systems.
Collapse
Affiliation(s)
- Gouri S Jas
- Department of Pharmaceutical Chemistry, The University of Kansas, Lawrence, Kansas 66047, USA
| | | | - Carey K Johnson
- Department of Chemistry, The University of Kansas, Lawrence, Kansas 66045, USA
| | - Krzysztof Kuczera
- Department of Chemistry, The University of Kansas, Lawrence, Kansas 66045, USA
| |
Collapse
|
4
|
Molina-Bolívar JA, Galisteo-González F, Ruiz CC, Medina-O'Donnell M, Martínez A, Parra A. Maslinic acid conjugate with 7-amino-4-methylcoumarin as probe to monitor the temperature dependent conformational changes of human serum albumin by FRET. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2019; 214:161-169. [PMID: 30776717 DOI: 10.1016/j.saa.2019.02.014] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/14/2018] [Revised: 01/14/2019] [Accepted: 02/09/2019] [Indexed: 05/25/2023]
Abstract
Synthesis, characterization and spectroscopic investigation of maslinic acid labeled with fluorescent 7-amino-4-methylcoumarin is reported. It was found that the coumarin-maslinic derivative (MaCo) forms an excellent fluorescence resonance energy transfer (FRET) pair with the tryptophan (Trp) residue of human serum albumin (HSA). This feature allowed for monitoring HSA conformational alterations by measuring the distance between donor (Trp) and acceptor (MaCo) through Förster energy transfer mechanism. Displacement experiments confirmed that MaCo binds to subdomain IIA of HSA with independence of temperature. It was observed that, in the temperature range 35-45 °C, the fluorescence emission maximum of HSA-MaCo complex decreased, whereas in the range 45 °C-65 °C, an increment was detected. The concomitant change in the polarity of environment surrounding Trp was confirmed by red edge excitation shift experiments. Thermal denaturation of HSA was followed by time-resolved fluorescence spectroscopy. Average lifetime of Trp residue decreased with temperature due to the increment of solvent collisions and changes in the solvent exposure of Trp. To discriminate the importance of each effect, lifetime of N-Acetyl-L-tryptophanamide (NATA) at different temperatures was measured. Circular dichroism (CD) studies confirmed the loss of secondary structure of HSA with increasing temperature and showed a different trend in the conformational transformation below and above 45 °C, in agreement with steady-state and time-resolved fluorescence experiments.
Collapse
Affiliation(s)
- J A Molina-Bolívar
- Department of Applied Physics II, Engineering School, University of Málaga, 29071 Málaga, Spain.
| | | | - C Carnero Ruiz
- Department of Applied Physics II, Engineering School, University of Málaga, 29071 Málaga, Spain
| | - M Medina-O'Donnell
- Department of Organic Chemistry, University of Granada, 18071 Granada, Spain
| | - A Martínez
- Department of Organic Chemistry, University of Granada, 18071 Granada, Spain
| | - A Parra
- Department of Organic Chemistry, University of Granada, 18071 Granada, Spain
| |
Collapse
|
5
|
Biswas B, Muttathukattil AN, Reddy G, Singh PC. Contrasting Effects of Guanidinium Chloride and Urea on the Activity and Unfolding of Lysozyme. ACS OMEGA 2018; 3:14119-14126. [PMID: 31458105 PMCID: PMC6644995 DOI: 10.1021/acsomega.8b01911] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/06/2018] [Accepted: 10/12/2018] [Indexed: 05/03/2023]
Abstract
Cosolvents play an important role in regulating the stability and function of proteins present in the cell. We studied the role of cosolvents, urea and guanidinium chloride (GdmCl), which act as protein denaturants, in the catalytic activity and structural stability of the protein lysozyme using activity measurements, spectroscopy, and molecular dynamics simulations. We find that the activity of lysozyme increases on the addition of urea, whereas it decreases sharply on the addition of GdmCl. At low GdmCl concentrations ([GdmCl] < 4 M), the activity of lysozyme decreases, even though there is no significant perturbation in the structure of the lysozyme folded state. We find that this is due to the strong interaction of the Gdm+ ion with the residues Asp52 and Glu35, which are present in the lysozyme catalytic site. In contrast, urea interacts with Trp63 present in the loop region present near the active site of lysozyme, inducing minor conformational changes in lysozyme, which can increase the activity of lysozyme. At higher denaturant concentrations, experiments show that GdmCl completely denatures the protein, whereas the folded state is stable in the presence of urea. We further show that GdmCl denatures lysozyme with the disulfide bonds intact in the protein, whereas urea denatures the protein only when the disulfide bonds are broken using reducing agents.
Collapse
Affiliation(s)
- Biswajit Biswas
- Department
of Spectroscopy, Indian Association for
the Cultivation of Science, 2A & 2B Raja S.C. Mullick Road, Jadavpur, Kolkata 700032, India
| | - Aswathy N. Muttathukattil
- Solid
State and Structural Chemistry Unit, Indian
Institute of Science, Bengaluru 560012, Karnataka, India
| | - Govardhan Reddy
- Solid
State and Structural Chemistry Unit, Indian
Institute of Science, Bengaluru 560012, Karnataka, India
- E-mail: (G.R.)
| | - Prashant Chandra Singh
- Department
of Spectroscopy, Indian Association for
the Cultivation of Science, 2A & 2B Raja S.C. Mullick Road, Jadavpur, Kolkata 700032, India
- E-mail: (P.C.S.)
| |
Collapse
|