Long-Range Intra-Protein Communication Can Be Transmitted by Correlated Side-Chain Fluctuations Alone

Allosteric regulation is a key component of cellular communication, but the way in which information is passed from one site to another within a folded protein is not often clear. While backbone motions have long been considered essential for long-range information conveyance, side-chain motions have rarely been considered. In this work, we demonstrate their potential utility using Monte Carlo sampling of side-chain torsional angles on a fixed backbone to quantify correlations amongst side-chain inter-rotameric motions. Results indicate that long-range correlations of side-chain fluctuations can arise independently from several different types of interactions: steric repulsions, implicit solvent interactions, or hydrogen bonding and salt-bridge interactions. These robust correlations persist across the entire protein (up to 60 Å in the case of calmodulin) and can propagate long-range changes in side-chain variability in response to single residue perturbations

Mechanism of Protein Kinetic Stabilization by Engineered Disulfide Crosslinks

The impact of disulfide bonds on protein stability goes beyond simple equilibrium thermodynamics effects associated with the conformational entropy of the unfolded state. Indeed, disulfide crosslinks may play a role in the prevention of dysfunctional association and strongly affect the rates of irreversible enzyme inactivation, highly relevant in biotechnological applications. While these kinetic-stability effects remain poorly understood, by analogy with proposed mechanisms for processes of protein aggregation and fibrillogenesis, we propose that they may be determined by the properties of sparsely-populated, partially-unfolded intermediates. Here we report the successful design, on the basis of high temperature molecular-dynamics simulations, of six thermodynamically and kinetically stabilized variants of phytase from Citrobacter braakii (a biotechnologically important enzyme) with one, two or three engineered disulfides. Activity measurements and 3D crystal structure determination demonstrate that the engineered crosslinks do not cause dramatic alterations in the native structure. The inactivation kinetics for all the variants displays a strongly non-Arrhenius temperature dependence, with the time-scale for the irreversible denaturation process reaching a minimum at a given temperature within the range of the denaturation transition. We show this striking feature to be a signature of a key role played by a partially unfolded, intermediate state/ensemble. Energetic and mutational analyses confirm that the intermediate is highly unfolded (akin to a proposed critical intermediate in the misfolding of the prion protein), a result that explains the observed kinetic stabilization. Our results provide a rationale for the kinetic-stability consequences of disulfide-crosslink engineering and an experimental methodology to arrive at energetic/structural descriptions of the sparsely populated and elusive intermediates that play key roles in irreversible protein denaturation.This work was supported by grants BIO2009-09562, CSD2009-00088 from the Spanish Ministry of Science and Innovation, and FEDER Funds (JMS-R)

Local Cooperativity in an Amyloidogenic State of Human Lysozyme Observed at Atomic Resolution

The partial unfolding of human lysozyme underlies its conversion from the soluble state into amyloid fibrils observed in a fatal hereditary form of systemic amyloidosis. To understand the molecular origins of the disease, it is critical to characterize the structural and physicochemical properties of the amyloidogenic states of the protein. Here we provide a high-resolution view of the unfolding process at low pH for three different lysozyme variants, the wild-type protein and the mutants I56T and I59T, which show variable stabilities and propensities to aggregate in vitro. Using a range of biophysical techniques that includes differential scanning calorimetry and nuclear magnetic resonance spectroscopy, we demonstrate that thermal unfolding under amyloidogenic solution conditions involves a cooperative loss of native tertiary structure, followed by progressive unfolding of a compact, molten globule-like denatured state ensemble as the temperature is increased. The width of the temperature window over which the denatured ensemble progressively unfolds correlates with the relative amyloidogenicity and stability of these variants, and the region of lysozyme that unfolds first maps to that which forms the core of the amyloid fibrils formed under similar conditions. Together, these results present a coherent picture at atomic resolution of the initial events underlying amyloid formation by a globular protein

Binding Free Energy Landscape of Domain-Peptide Interactions

Peptide recognition domains (PRDs) are ubiquitous protein domains which mediate large numbers of protein interactions in the cell. How these PRDs are able to recognize peptide sequences in a rapid and specific manner is incompletely understood. We explore the peptide binding process of PDZ domains, a large PRD family, from an equilibrium perspective using an all-atom Monte Carlo (MC) approach. Our focus is two different PDZ domains representing two major PDZ classes, I and II. For both domains, a binding free energy surface with a strong bias toward the native bound state is found. Moreover, both domains exhibit a binding process in which the peptides are mostly either bound at the PDZ binding pocket or else interact little with the domain surface. Consistent with this, various binding observables show a temperature dependence well described by a simple two-state model. We also find important differences in the details between the two domains. While both domains exhibit well-defined binding free energy barriers, the class I barrier is significantly weaker than the one for class II. To probe this issue further, we apply our method to a PDZ domain with dual specificity for class I and II peptides, and find an analogous difference in their binding free energy barriers. Lastly, we perform a large number of fixed-temperature MC kinetics trajectories under binding conditions. These trajectories reveal significantly slower binding dynamics for the class II domain relative to class I. Our combined results are consistent with a binding mechanism in which the peptide C terminal residue binds in an initial, rate-limiting step

Change in Allosteric Network Affects Binding Affinities of PDZ Domains: Analysis through Perturbation Response Scanning

The allosteric mechanism plays a key role in cellular functions of several PDZ domain proteins (PDZs) and is directly linked to pharmaceutical applications; however, it is a challenge to elaborate the nature and extent of these allosteric interactions. One solution to this problem is to explore the dynamics of PDZs, which may provide insights about how intramolecular communication occurs within a single domain. Here, we develop an advancement of perturbation response scanning (PRS) that couples elastic network models with linear response theory (LRT) to predict key residues in allosteric transitions of the two most studied PDZs (PSD-95 PDZ3 domain and hPTP1E PDZ2 domain). With PRS, we first identify the residues that give the highest mean square fluctuation response upon perturbing the binding sites. Strikingly, we observe that the residues with the highest mean square fluctuation response agree with experimentally determined residues involved in allosteric transitions. Second, we construct the allosteric pathways by linking the residues giving the same directional response upon perturbation of the binding sites. The predicted intramolecular communication pathways reveal that PSD-95 and hPTP1E have different pathways through the dynamic coupling of different residue pairs. Moreover, our analysis provides a molecular understanding of experimentally observed hidden allostery of PSD-95. We show that removing the distal third alpha helix from the binding site alters the allosteric pathway and decreases the binding affinity. Overall, these results indicate that (i) dynamics plays a key role in allosteric regulations of PDZs, (ii) the local changes in the residue interactions can lead to significant changes in the dynamics of allosteric regulations, and (iii) this might be the mechanism that each PDZ uses to tailor their binding specificities regulation

Population of nonnative States of lysozyme variants drives amyloid fibril formation.

The propensity of protein molecules to self-assemble into highly ordered, fibrillar aggregates lies at the heart of understanding many disorders ranging from Alzheimer's disease to systemic lysozyme amyloidosis. In this paper we use highly accurate kinetic measurements of amyloid fibril growth in combination with spectroscopic tools to quantify the effect of modifications in solution conditions and in the amino acid sequence of human lysozyme on its propensity to form amyloid fibrils under acidic conditions. We elucidate and quantify the correlation between the rate of amyloid growth and the population of nonnative states, and we show that changes in amyloidogenicity are almost entirely due to alterations in the stability of the native state, while other regions of the global free-energy surface remain largely unmodified. These results provide insight into the complex dynamics of a macromolecule on a multidimensional energy landscape and point the way for a better understanding of amyloid diseases

A Nanobody Binding to Non-amyloidogenic Regions of the Protein Human Lysozyme Enhances Partial Unfolding but Inhibits Amyloid Fibril Formation.

We report the effects of the interaction of two camelid antibody fragments, generally called nanobodies, namely cAb-HuL5 and a stabilized and more aggregation-resistant variant cAb-HuL5G obtained by protein engineering, on the properties of two amyloidogenic variants of human lysozyme, I56T and D67H, whose deposition in vital organs including the liver, kidney, and spleen is associated with a familial non-neuropathic systemic amyloidosis. Both NMR spectroscopy and X-ray crystallographic studies reveal that cAb-HuL5 binds to the α-domain, one of the two lobes of the native lysozyme structure. The binding of cAb-HuL5/cAb-HuL5G strongly inhibits fibril formation by the amyloidogenic variants; it does not, however, suppress the locally transient cooperative unfolding transitions, characteristic of these variants, in which the β-domain and the C-helix unfold and which represents key early intermediate species in the formation of amyloid fibrils. Therefore, unlike two other nanobodies previously described, cAb-HuL5/cAb-HuL5G does not inhibit fibril formation via the restoration of the global cooperativity of the native structure of the lysozyme variants to that characteristic of the wild-type protein. Instead, it inhibits a subsequent step in the assembly of the fibrils, involving the unfolding and structural reorganization of the α-domain. These results show that nanobodies can protect against the formation of pathogenic aggregates at different stages in the structural transition of a protein from the soluble native state into amyloid fibrils, illustrating their value as structural probes to study the molecular mechanisms of amyloid fibril formation. Combined with their amenability to protein engineering techniques to improve their stability and solubility, these findings support the suggestion that nanobodies can potentially be developed as therapeutics to combat protein misfolding diseases.Department of Applied Biology and Chemical Technolog