Reorientation Motion and Preferential Interactions of a Peptide in Denaturants and Osmolyte

Solvent accessible surface area-assessed molecular basis of osmolyte-induced protein stability

In solvent-modulated protein folding, under certain physiological conditions, an equilibrium exists between the unfolded and folded states of the protein without any need to break or make a covalent bond. In this process, interactions between various protein groups (peptides) and solvent molecules are known to play a major role in determining the directionality of the chemical reaction. However, an understanding of the mechanism of action of the co(solvent) by a generic theoretical underpinning is lacking. In this study, a generic solvation model is developed based on statistical mechanics and the thermodynamic transfer free energy model by considering the microenvironment polarity of the interacting co(solvent)-protein system. According to this model, polarity and the fractional solvent-accessible surface areas contribute to the interaction energies. The present model includes various orientations of participating interactant solvent surfaces of suitable areas. As model systems, besides the backbone we consider naturally occurring amino acid residues solvated in ten different osmolytes, small organic compounds known to modulate protein stability. The present model is able to predict the correct trend of the osmolyte-peptide interactions ranging from stabilizing to destabilizing not only for the backbone but also for side chains. Our model predicts Asn, Gln, Asp, Glu, Arg and Pro to be highly stable in most of the protecting osmolytes while Ala, Val, Ile, Leu, Thr, Met, Lys, Phe, Trp and Tyr are predicted to be moderately stable, and Ser, Cys and Histidine are predicted to be least stable. However, in denaturing solvents, both backbone and side chain models show similar stabilities in urea and guanidine. One of the important aspects of this model is that it is essentially parameter-free and consistent with the electrostatics of the interaction partners that make this model suitable for estimating any solute-solvent interaction energies.

Effects of Proline on Internal Friction in Simulated Folding Dynamics of Several Alanine-Based α-Helical Peptides

We have studied in silico the effect of proline, a model cosolvent, on local and global friction coefficients in (un)folding of several typical alanine-based α-helical peptides. Local friction is related to dwell times of a single, ensemble-averaged hydrogen bond (HB) within each peptide. Global friction is related to energy dissipated in a series of configurational changes of each peptide experienced by increasing the number of HBs during folding. Both of these approaches are important in relation to future atomic force microscopic-based measurements of internal friction via force-clamp single-molecule force spectroscopy. Molecular dynamics (MD) simulations for six peptides, namely, ALA5, ALA8, ALA15, ALA21, (AAQAA)3, and H2N-GN(AAQAA)2G-COONH2, have been conducted at 2 and 5 M proline solutions in water. Using previously obtained MD data for these peptides in pure water as well as upgraded theoretical models, we obtained variations of local and global internal friction coefficients as a function of solution viscosity. The results showed the substantial role of proline in stabilizing the folded state and slowing the overall folding dynamics. Consequently, larger friction coefficients were obtained at larger viscosities. The local and global internal friction, i.e., respective, friction coefficients approximated to zero viscosity, was also obtained. The evolution of friction coefficients with viscosity was weakly dependent on the number of concurrent folding pathways but was rather dominated by a stabilizing effect of proline on the folded states. Obtained values of local and global internal friction showed qualitatively similar results and a clear dependency on the structure of the studied peptide.

Observing reorientation dynamics with Time-Resolved fluorescence and molecular dynamics in varying periodic boundary conditions

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 ¹L_b 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.

In silico studies of the human IAPP in the presence of osmolytes

The human islet amyloid polypeptide or amylin is secreted along with insulin by pancreatic islets. Under the drastic environmental conditions, amylin can aggregate to form amyloid fibrils. This amyloid plaque of hIAPP in the pancreatic cells is the cause of type II diabetes. Early stages of amylin aggregates are more cytotoxic than the matured fibrils. Here, we have used the all-atom molecular dynamic simulation to see the effect of water, TMAO, urea and urea/TMAO having ratio 2:1 of different concentrations on the amylin protein. Our study suggest that the amylin protein forms β-sheets in its monomeric form and may cause the aggregation of protein through the residue 13-17 and the C-terminal region. α-Helical content of protein increases with an increase in TMAO concentration by decreasing the SASA value of protein, increase in intramolecular hydrogen bonds and on making the short-range hydrophobic interactions. Electrostatic potential surfaces show that hydrophobic groups are buried and normalised configurational entropy of backbone, and side-chain atoms is lesser in the presence of TMAO, whereas opposite behaviour is obtained in the case of urea. Counteraction effect of TMAO using Kast model towards urea is also observed in ternary solution of urea/TMAO.

Hydration of Simple Model Peptides in Aqueous Osmolyte Solutions

The biology and chemistry of proteins and peptides are inextricably linked with water as the solvent. The reason for the high stability of some proteins or uncontrolled aggregation of others may be hidden in the properties of their hydration water. In this study, we investigated the effect of stabilizing osmolyte–TMAO (trimethylamine N-oxide) and destabilizing osmolyte–urea on hydration shells of two short peptides, NAGMA (N-acetyl-glycine-methylamide) and diglycine, by means of FTIR spectroscopy and molecular dynamics simulations. We isolated the spectroscopic share of water molecules that are simultaneously under the influence of peptide and osmolyte and determined the structural and energetic properties of these water molecules. Our experimental and computational results revealed that the changes in the structure of water around peptides, caused by the presence of stabilizing or destabilizing osmolyte, are significantly different for both NAGMA and diglycine. The main factor determining the influence of osmolytes on peptides is the structural-energetic similarity of their hydration spheres. We showed that the chosen peptides can serve as models for various fragments of the protein surface: NAGMA for the protein backbone and diglycine for the protein surface with polar side chains.

Surface hydration and preferential interaction directs the charged amino acids-induced changes in protein stability

In the present study, we investigate the interaction of amino acid osmolytes, Arg, Lys, Asp and Glu, and a denaturant, guanidinium chloride (Gdm) with proteins. To achieve this, molecular dynamics (MD) simulation of RNase A and α-lactalbumin was performed in the presence of three charged amino acids Arg, Lys, and Asp and the molecular mechanism of amino acid-induced (de)stabilization of the proteins was examined by combining with our earlier report on Glu. As Arg has the side chain similar to that of Gdm and destabilizes the proteins, MD simulation was carried out in the presence of Gdm as well. Radial distribution function and hydration fraction around the protein surface reveals that preferential hydration increases upon the addition of any of the cosolvent; however, the extent of increase is more in the presence of stabilizing cosolvents (stAAs: Lys, Asp and Glu) compared to destabilizing cosolvents (Arg and Gdm). Moreover, the preferential interaction of Arg and Gdm with the proteins is higher than that of stAAs. Residue-level interaction analysis suggests that stAAs preferably interacts with charged amino acids of the proteins whereas Arg and Gdm interactions could be found on almost all the surface exposed residues which might provide higher preferential interaction for these residues. From the results, we propose that the net outcome of preferential hydration versus preferential interaction of the amino acids might determine their effect on the stability of proteins.

Dynamic elements and kinetics: Most favorable conformations of peptides in solution with measurements and simulations

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% 3₁₀ 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 3₁₀ 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.

Crowding Stabilizes DMSO-Water Hydrogen-Bonding Interactions

Up to 40% of intracellular water is confined due to the dense packing of macromolecules, ions, and osmolytes. Despite the large body of work concerning the effect of additives on the biomolecular structure and stability, the role of crowding and heterogeneity in these interactions is not well understood. Here, infrared spectroscopy and molecular dynamics simulations are used to describe the mechanisms by which crowding modulates hydrogen bonding interactions between water and dimethyl sulfoxide (DMSO). Specifically, we use formamide and dimethylformamide (DMF) as molecular crowders and show that the S═O hydrogen bond populations in aqueous mixtures are increased by both amides. These additives increase the amount of water within the DMSO first solvation shell through two mechanisms: (a) directly stabilizing water-DMSO hydrogen bonds; (b) increasing water exposure by destabilizing DMSO-DMSO self-interactions. Further, we quantified the hydrogen bond enthalpies between the different components: DMSO-water (61 kJ/mol) > DMSO-formamide (32 kJ/mol) > water-water (23 kJ/mol) ≫ formamide-water (4.7 kJ/mol). Spectra of carbonyl stretching vibrations show that DMSO induces the dehydration of amides as a result of strong DMSO-water interactions, which has been suggested as the main mechanism of protein destabilization.

Effects of Trimethylamine-N-oxide on the Conformation of Peptides and its Implications for Proteins

To provide insights into the stabilizing mechanisms of trimethylamine-N-oxide (TMAO) on protein structures, we perform all-atom molecular dynamics simulations of peptides and the Trp-cage miniprotein. The effects of TMAO on the backbone and charged residues of peptides are found to stabilize compact conformations, whereas effects of TMAO on nonpolar residues lead to peptide swelling. This suggests competing mechanisms of TMAO on proteins, which accounts for hydrophobic swelling, backbone collapse, and stabilization of charge-charge interactions. These mechanisms are observed in Trp cage.

TMAO and urea in the hydration shell of the protein SNase

We performed all-atom MD simulations of the protein SNase in aqueous solution and in the presence of two major osmolytes, trimethylamine-N-oxide (TMAO) and urea, as cosolvents at various concentrations and compositions and at different pressures and temperatures.

Preferential Solvation of a Highly Medium Responsive Pentacyanoferrate(II) Complex in Binary Solvent Mixtures: Understanding the Role of Dielectric Enrichment and the Specificity of Solute-Solvent Interactions

In this work, the preferential solvation of an intensely solvatochromic ferrocyanide(II) dye involving a 4,4'-bipyridine-based ligand was examined in various binary solvent mixtures. Its solvatochromic behavior was rationalized in terms of specific and nonspecific solute-solvent interactions. An exceptional case of solvatochromic inversion was observed when going from alcohol/water to amide/water mixtures. These effects were quantified using Onsager's solvent polarity function. Furthermore, the sensitivity of the solvatochromism of the dye was determined using various solvatochromic parameters such as π* expressing the dipolarity/polarizability of solvents and α expressing the hydrogen-bond-donor acidity of solvents. This analysis was useful for the rationalization of the selective solvation phenomena occurring in the three types of alcohol/water and amide/water mixtures studied. Furthermore, two preferential solvation models were employed for the interpretation of the experimental spectral results in binary solvent mixtures, namely, the model of Suppan on dielectric enrichment [J. Chem. Soc. Faraday Trans. 1 1987, 83, 495-509] and the model of Bosch, Rosés, and co-workers [J. Chem. Soc., Perkin Trans. 2, 1995, 8, 1607-1615]. The first model successfully predicted the charge transfer energies of the dye in formamide/water and N-methylformamide/water mixtures, but in the case of MeOH/water mixtures, the prediction was less accurate because of the significant contribution of specific solute-solvent interactions in that case. The second model gave more insights for both specific solute-solvent as well as solvent-solvent interactions in the cybotactic region. The role of dielectric enrichment and specific interactions was discussed based on the findings.