Nonadditivity in the alpha-helix to coil transition

Processing helix–coil transition data: Account of chain length and solvent effects

Numerous nanobiotechnologies include manipulations of short polypeptide chains. The conformational properties of these polypeptides are studied in vitro by circular dichroism and time-resolved infrared spectroscopy. To find out the interaction parameters, the measured temperature dependence of normalized helicity degree needs to be further processed by fitting to a model. Using recent advances in the Hamiltonian formulation of the classical Zimm and Bragg model, we explicitly include chain length and solvent effects in the theoretical description. The expression for the helicity degree we suggest successfully fits the experimental data and provides hydrogen bonding energies and nucleation parameter values within the standards in the field.

Implicit water model within the Zimm-Bragg approach to analyze experimental data for heat and cold denaturation of proteins

Studies of biopolymer conformations essentially rely on theoretical models that are routinely used to process and analyze experimental data. While modern experiments allow study of single molecules in vivo, corresponding theories date back to the early 1950s and require an essential update to include the recent significant progress in the description of water. The Hamiltonian formulation of the Zimm-Bragg model we propose includes a simplified, yet explicit model of water-polypeptide interactions that transforms into the equivalent implicit description after performing the summation of solvent degrees of freedom in the partition function. Here we show that our model fits very well to the circular dichroism experimental data for both heat and cold denaturation and provides the energies of inter- and intra-molecular H-bonds, unavailable with other processing methods. The revealed delicate balance between these energies determines the conditions for the existence of cold denaturation and thus clarifies its absence in some proteins.

Ensemble properties of network rigidity reveal allosteric mechanisms

The distance constraint model (DCM) is a unique computational modeling paradigm that integrates mechanical and thermodynamic descriptions of macromolecular structure. That is, network rigidity calculations are used to account for nonadditivity within entropy components, thus restoring the utility of free-energy decomposition. The DCM outputs a large number of structural characterizations that collectively allow for quantified stability-flexibility relationships (QSFR) to be identified. In this review, we describe the theoretical underpinnings of the DCM and introduce several common QSFR metrics. Application of the DCM across protein families highlights the sensitivity within the set of protein structure residue-to-residue couplings. Further, we have developed a perturbation method to identify putative allosteric sites, where large changes in QSFR upon rigidification (mimicking ligand-binding) detect sites likely to invoke allosteric changes.

50 Years of Lifson-Roig Models: Application to Molecular Simulation Data

Simple helix-coil transition theories have been indispensable tools in the analysis of data reporting on the reversible folding of α-helical polypeptides. They provide a transferable means to not only characterize different systems but to also compare different techniques, viz., experimental probes monitoring helix-coil transitions in vitro or biomolecular force fields in silico. This article addresses several issues with the application of Lifson-Roig theory to helix-coil transition data. We use computer simulation to generate two sets of ensembles for the temperature-controlled, reversible folding of the 21-residue, alanine-rich FS peptide. Ensembles differ in the rigidity of backbone bond angles and are analyzed using two distinct descriptors of helicity. The analysis unmasks an underlying phase diagram that is surprisingly complex. The complexities give rise to fitted nucleation and propagation parameters that are difficult to interpret and that are inconsistent with the distribution of isolated residues in the α-helical basin. We show that enthalpies of helix formation are more robustly determined using van't Hoff analysis of simple measures of helicity rather than fitted propagation parameters. To overcome some of these issues, we design a simple variant of the Lifson-Roig model that recovers physical interpretability of the obtained parameters by allowing bundle formation to be described in simple fashion. The relevance of our results is discussed in relation to the applicability of Lifson-Roig models to both in silico and in vitro data.