Core-shell structures in single flexible-semiflexible block copolymers: Finding the free energy minimum for the folding transition

We investigate the folding transition of a single diblock copolymer consisting of a semiflexible and a flexible block. We obtain a {\it Saturn-shaped} core-shell conformation in the folded state, in which the flexible block forms a core and the semiflexible block wraps around it. We demonstrate two distinctive features of the core-shell structures: (i) The kinetics of the folding transition in the copolymer are significantly more efficient than those of a semiflexible homopolymer. (ii) The core-shell structure does not depend on the transition pathway

Random copolymer adsorption: Morita approximation compared to exact numerical simulations

We study the adsorption of ideal random lattice copolymers with correlations in the sequences on homogeneous substrates with two different methods: An analytical solution of the problem based on the constrained annealed approximation introduced by Morita in 1964 and the generating functional (GF) technique, and direct numerical simulations of lattice chains averaged over many realizations of random sequences. Both methods allow to calculate the free energy and different conformational characteristics of the adsorbed chain. The comparison of the results for random copolymers with different degree of correlations and different types of nonadsorbing monomers (neutral or repelling from the surface) shows not only qualitative but a very good quantitative agreement, especially in the cases of Bernoullian and quasi-alternating random sequences.Comment: 19 pages, 9 figure

Molecular Dynamics of "Fuzzy" Transcriptional Activator-Coactivator Interactions

Transcriptional activation domains (ADs) are generally thought to be intrinsically unstructured, but capable of adopting limited secondary structure upon interaction with a coactivator surface. The indeterminate nature of this interface made it hitherto difficult to study structure/function relationships of such contacts. Here we used atomistic accelerated molecular dynamics (aMD) simulations to study the conformational changes of the GCN4 AD and variants thereof, either free in solution, or bound to the GAL11 coactivator surface. We show that the AD-coactivator interactions are highly dynamic while obeying distinct rules. The data provide insights into the constant and variable aspects of orientation of ADs relative to the coactivator, changes in secondary structure and energetic contributions stabilizing the various conformers at different time points. We also demonstrate that a prediction of α-helical propensity correlates directly with the experimentally measured transactivation potential of a large set of mutagenized ADs. The link between α-helical propensity and the stimulatory activity of ADs has fundamental practical and theoretical implications concerning the recruitment of ADs to coactivators

Transition States in Protein Folding Kinetics: The Structural Interpretation of Phi-values

Phi-values are experimental measures of the effects of mutations on the folding kinetics of a protein. A central question is which structural information Phi-values contain about the transition state of folding. Traditionally, a Phi-value is interpreted as the 'nativeness' of a mutated residue in the transition state. However, this interpretation is often problematic because it assumes a linear relation between the nativeness of the residue and its free-energy contribution. We present here a better structural interpretation of Phi-values for mutations within a given helix. Our interpretation is based on a simple physical model that distinguishes between secondary and tertiary free-energy contributions of helical residues. From a linear fit of our model to the experimental data, we obtain two structural parameters: the extent of helix formation in the transition state, and the nativeness of tertiary interactions in the transition state. We apply our model to all proteins with well-characterized helices for which more than 10 Phi-values are available: protein A, CI2, and protein L. The model captures nonclassical Phi-values 1 in these helices, and explains how different mutations at a given site can lead to different Phi-values.Comment: 26 pages, 7 figures, 5 table

Thermal denaturation of fluctuating finite DNA chains: the role of bending rigidity in bubble nucleation

Statistical DNA models available in the literature are often effective models where the base-pair state only (unbroken or broken) is considered. Because of a decrease by a factor of 30 of the effective bending rigidity of a sequence of broken bonds, or bubble, compared to the double stranded state, the inclusion of the molecular conformational degrees of freedom in a more general mesoscopic model is needed. In this paper we do so by presenting a 1D Ising model, which describes the internal base pair states, coupled to a discrete worm like chain model describing the chain configurations [J. Palmeri, M. Manghi, and N. Destainville, Phys. Rev. Lett. 99, 088103 (2007)]. This coupled model is exactly solved using a transfer matrix technique that presents an analogy with the path integral treatment of a quantum two-state diatomic molecule. When the chain fluctuations are integrated out, the denaturation transition temperature and width emerge naturally as an explicit function of the model parameters of a well defined Hamiltonian, revealing that the transition is driven by the difference in bending (entropy dominated) free energy between bubble and double-stranded segments. The calculated melting curve (fraction of open base pairs) is in good agreement with the experimental melting profile of polydA-polydT. The predicted variation of the mean-square-radius as a function of temperature leads to a coherent novel explanation for the experimentally observed thermal viscosity transition. Finally, the influence of the DNA strand length is studied in detail, underlining the importance of finite size effects, even for DNA made of several thousand base pairs.Comment: Latex, 28 pages pdf, 9 figure

Kinetic pathway analysis of an α-helix in two protonation states: Direct observation and optimal dimensionality reduction

This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article appeared in J. Chem. Phys. 150, 074902 (2019); https://doi.org/10.1063/1.5082192 and may be found at https://aip.scitation.org/doi/full/10.1063/1.5082192Thermodynamically stable conformers of secondary structural elements make a stable tertiary/quaternary structure that performs its proper biological function efficiently. Formation mechanisms of secondary and tertiary/quaternary structural elements from the primary structure are driven by the kinetic properties of the respective systems. Here we have carried out thermodynamic and kinetic characterization of an alpha helical heteropeptide in two protonation states, created with the addition and removal of a proton involving a single histidine residue in the primary structure. Applying far-UV circular dichroism spectroscopy, the alpha helix is observed to be significantly more stable in the deprotonated state. Nanosecond laser temperature jump spectroscopy monitoring time-resolved tryptophan fluorescence on the protonated conformer is carried out to measure the kinetics of this system. The measured relaxation rates at a final temperature between 296K and 314 K generated a faster component of 20 ns–11 ns and a slower component of 314 ns–198 ns. Atomically detailed characterization of the helix-coil kinetic pathways is performed based on all-atom molecular dynamics trajectories of the two conformers. Application of clustering and kinetic coarse-graining with optimum dimensionality reduction produced description of the trajectories in terms of kinetic models with two to five states. These models include aggregate states corresponding to helix, coil, and intermediates. The “coil” state involves the largest number of conformations, consistent with the expected high entropy of this structural ensemble. The “helix” aggregate states are found to be mixed with the full helix and partially folded forms. The experimentally observed higher helix stability in the deprotonated form of the alpha helical heteropeptide is reflected in the nature of the “helix” aggregate state arising from the kinetic model. In the protonated form, the “coil” state exhibits the lowest free energy and longest lifetime, while in the deprotonated form, it is the “helix” that is found to be most stable. Overall, the coarse grained models suggest that the protonation of a single histidine residue in the primary structure induces significant changes in the free energy landscape and kinetic network of the studied helix-forming heteropeptide