scholarly journals Analysis of Two-State Folding Using Parabolic Approximation IV: Non-Arrhenius Kinetics of FBP28 WW Part-II

2016 ◽  
Author(s):  
Robert S Sade

A model which treats the denatured and the native conformers as being confined to harmonic Gibbs energy wells has been used to rationalize the physical basis for the non-Arrhenius behaviour of spontaneously-folding fixed two-state systems. It is shown that at constant pressure and solvent conditions: (i) the rate constant for folding will be a maximum when the heat released upon formation of net molecular interactions is exactly compensated by the heat absorbed to desolvate net polar and non-polar solvent accessible surface area (SASA), as the denatured conformers driven by thermal noise bury their SASA and diffuse on the Gibbs energy surface to reach the activated state; (ii) the rate constant for unfolding will be a minimum when the heat absorbed by the native conformers to break various net backbone and sidechain interactions is exactly compensated by the heat of hydration released due to the net increase in SASA, as the native conformers unravel to reach the activated state; (iii) the activation entropy for folding will be zero, and the Gibbs barrier to folding will be a minimum, when the decrease in the backbone and the sidechain mobility is exactly compensated by the increase in entropy due to solvent-release, as the denatured conformers bury their SASA to reach the activated state; (iv) the activation entropy for unfolding will be zero, and the Gibbs barrier to unfolding will be a maximum when the increase in the backbone and sidechain mobility is exactly compensated by the negentropy of solvent-capture on the protein surface, as the native conformers unravel to reach the activated state; (v) while cold denaturation is driven by solvent effects, heat denaturation is primarily due to chain effects; (vi) the speed-limit for the folding is ultimately due to conformational searching; and(vii) Levinthal's paradox may have little basis if the entropy of solvent-release that accompanies protein folding is taken into consideration.

2016 ◽  
Author(s):  
Robert S Sade

A model which treats the denatured and native conformers of spontaneously-folding fixed two-state systems as being confined to harmonic Gibbs energy-wells has been developed. Within the assumptions of this model the Gibbs energy functions of the denatured (DSE) and the native state (NSE) ensembles are described by parabolas, with the mean length of the reaction coordinate (RC) being given by the temperature-invariant denaturant m value. Consequently, the ensemble-averaged position of the transition state ensemble (TSE) along the RC, and the ensemble-averaged Gibbs energy of the TSE are determined by the intersection of the DSE and the NSE-parabolas. The equations derived enable equilibrium stability and the rate constants to be rationalized in terms of the mean and the variance of the Gaussian distribution of the solvent accessible surface area of the conformers in the DSE and the NSE. The implications of this model for protein folding are discussed.


2016 ◽  
Author(s):  
Robert S Sade

A model which treats the denatured and the native conformers as being confined to harmonic Gibbs energy wells has been used to analyse the non-Arrhenius behaviour of spontaneously-folding fixed two-state systems. The results demonstrate that when pressure and solvent are constant: (i) a two-state system is physically defined only for a finite temperature range; (ii) irrespective of the primary sequence, the 3-dimensional structure of the native conformer, the residual structure in the denatured state, and the magnitude of the folding and unfolding rate constants, the equilibrium stability of a two-state system is a maximum when its denatured conformers bury the least amount of solvent accessible surface area (SASA) to reach the activated state; (iii) the Gibbs barriers to folding and unfolding are not always due to the incomplete compensation of the activation enthalpies and entropies; (iv) the difference in heat capacity between the reaction-states is due to both the size of the solvent-shell and the non-covalent interactions; (v) the position of the transition state ensemble along the reaction coordinate (RC) depends on the choice of the RC; and (vi) the atomic structure of the transiently populated reaction-states cannot be inferred from perturbation-induced changes in their energetics.


2019 ◽  
Author(s):  
David Wright ◽  
Fouad Husseini ◽  
Shunzhou Wan ◽  
Christophe Meyer ◽  
Herman Van Vlijmen ◽  
...  

<div>Here, we evaluate the performance of our range of ensemble simulation based binding free energy calculation protocols, called ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) for use in fragment based drug design scenarios. ESMACS is designed to generate reproducible binding affinity predictions from the widely used molecular mechanics Poisson-Boltzmann surface area (MMPBSA) approach. We study ligands designed to target two binding pockets in the lactate dehydogenase A target protein, which vary in size, charge and binding mode. When comparing to experimental results, we obtain excellent statistical rankings across this highly diverse set of ligands. In addition, we investigate three approaches to account for entropic contributions not captured by standard MMPBSA calculations: (1) normal mode analysis, (2) weighted solvent accessible surface area (WSAS) and (3) variational entropy. </div>


2017 ◽  
Vol 95 (9) ◽  
pp. 991-998 ◽  
Author(s):  
Prabal K. Maiti

Using fully atomistic molecular dynamics simulation that are several hundred nanoseconds long, we demonstrate the pH-controlled sponge action of PAMAM dendrimer. We show how at varying pH levels, the PAMAM dendrimer acts as a wet sponge; at neutral or low pH levels, the dendrimer expands noticeably and the interior of the dendrimer opens up to host several hundreds to thousands of water molecules depending on the generation number. Increasing the pH (i.e., going from low pH to high pH) leads to the collapse of the dendrimer size, thereby expelling the inner water, which mimics the ‘sponge’ action. As the dendrimer size swells up at a neutral pH or low pH due to the electrostatic repulsion between the primary and tertiary amines that are protonated at this pH, there is dramatic increase in the available solvent accessible surface area (SASA), as well as solvent accessible volume (SAV).


2021 ◽  
Author(s):  
Cecylia Severin Lupala ◽  
Yongjin Ye ◽  
Hong Chen ◽  
Xiaodong Su ◽  
Haiguang Liu

The spreading of SARS-CoV-2 virus resulted the COVID-19 pandemic, which has caused more than 5 millions of death globally. Several major variants of SARS-CoV-2 have emerged and placed challenges in controlling the infections. The recently emerged Omicron variant raised serious concerns about reducing efficacy of antibodies or vaccines, due to its vast mutations. We modelled the complex structure of human ACE2 protein and the receptor binding domain of Omicron variant, then conducted atomistic molecular dynamics simulations to study the binding interactions. The analysis shows that the Omicron variant RBD binds more strongly to the human ACE2 protein than the original strain. The mutation at the ACE2-RBD interface enhanced the tight binding by increasing hydrogen bonding interaction and enlarging buried solvent accessible surface area.


2019 ◽  
Vol 48 (8) ◽  
pp. 773-779 ◽  
Author(s):  
Anatoliy Dragan ◽  
Peter Privalov ◽  
Colyn Crane-Robinson

Abstract The heat capacity change, ΔCp, accompanying the folding/unfolding of macromolecules reflects their changing state of hydration. Thermal denaturation of the DNA duplex is characterized by an increase in ΔCp but of much lower magnitude than observed for proteins. To understand this difference, the changes in solvent accessible surface area (ΔASA) have been determined for unfolding the B-form DNA duplex into disordered single strands. These showed that the polar component represents ~ 55% of the total increase in ASA, in contrast to globular proteins of similar molecular weight for which the polar component is only about 1/3rd of the total. As the exposure of polar surface results in a decrease of ΔCp, this explains the much reduced heat capacity increase observed for DNA and emphasizes the enhanced role of polar interactions in maintaining duplex structure. Appreciation of a non-zero ΔCp for DNA has important consequences for the calculation of duplex melting temperatures (Tm). A modified approach to Tm prediction is required and comparison is made of current methods with an alternative protocol.


Sign in / Sign up

Export Citation Format

Share Document