Computational Prediction of Circular Dichroism Spectra and Quantification of Helicity Loss upon Peptide Adsorption on Silica

Abstract

Circular dichroism (CD) spectroscopy is one of the few experimental techniques sensitive to the structural changes that peptides undergo when they adsorb on inorganic material surfaces, a problem of deep significance in medicine, biotechnology, and materials science. Although the theoretical calculation of the CD spectrum of a molecule in a given conformation can be routinely performed, the inverse problem of extracting atomistic structural details from a measured spectrum is not uniquely determined. Especially complicated is the case of oligopeptides, whose folding/unfolding energy landscapes are often very broad and shallow. This means that the CD spectra measured for either dissolved or adsorbed peptides arise from a multitude of different structures, each present with a probability dictated by their relative free-energy variations, according to Boltzmann statistics. Here we present a modeling method based on replica exchange with solute tempering in combination with metadynamics, which allows us to predict both the helicity loss of a small peptide upon interaction with silica colloids in water and to compute the full CD spectra of the adsorbed and dissolved states, in quantitative agreement with experimental measurements. In our method, the CD ellipticity Θ for any given wavelength λ is calculated as an external collective variable by means of reweighting the biased trajectory obtained using the peptide-SiO2 surface distance and the structural helicity as two independent, internal collective variables. Our results also provide support for the often-employed hypothesis that the Θ intensity at λ = 222 nm is linearly correlated with the peptides' fractional helicity.

Quantitative first principles calculations of protein circular dichroism in the near-ultraviolet

“…The influence of conformational dynamics could also be explored using molecular dynamics simulations. Some individual cases have been studied 96 – 98 and we are planning a more comprehensive investigation.…”

Lowest matric potential in quartz: Metadynamics evidence

“…The free‐energy difference between the pure bulk water state and the desorbed state is the self‐hydration free energy of pure bulk water (Δ F 2 ). These two free energies are frequently used in chemistry and material science [e.g., Hermans et al ., ; Meißner et al ., ]. Thereafter, Δ F can be calculated with Δ F 1 and Δ F 2 : …”

What makes a good graphene-binding peptide? Adsorption of amino acids and peptides at aqueous graphene interfaces

“…Second, many materials-binding peptides are thought to be intrinsically disordered; the corresponding potential energy landscape describing the peptide conformational ensemble is anticipated to be complex. Therefore careful conformational sampling is an indispensable tool for ensuring meaningful results [17][18][19][20][21][22][23][24][25][26][27][28][29] , since insufficient sampling may give rise to misleading conclusions. The challenge of extensive conformational sampling is made more acute by the need to use an explicit description of liquid water, 30,31 since spatial structuring of interfacial solvent is thought to be a key determinant in peptide-materials binding.…”