Experimental free energy surfaces reveal the mechanisms of maintenance of protein solubility

Plant protein solubility: A challenge or insurmountable obstacle

Currently, there is an increasing focus on comprehending the solubility of plant-based proteins, driven by the rising demand for animal-free food formulations. The solubility of proteins plays a crucial role in impacting other functional properties of proteins and food processing. Consequently, understanding protein solubility in a deeper sense may allow a better usage of plant proteins. Herein, we discussed the definition of protein solubility from both thermodynamic and colloidal perspectives. A range of factors affecting solubility of plant proteins are generalized, including intrinsic factors (amino acids composition, hydrophobicity), and extrinsic factors (pH, ionic strength, extraction and drying methods). Current methods to enhance solubility are outlined, including microwave, high intensity ultrasound, hydrostatic pressure, glycation, pH-shifting, enzymatic hydrolysis, enzymatic cross-linking, complexation and modulation of amino acids. We base the discussion on diverse modified methods of nitrogen solubility index available to determine and analyze protein solubility followed by addressing how other indigenous components affect the solubility of plant proteins. Some nonproteinaceous constituents in proteins such as carbohydrates and polyphenols may exert positive or negative impact on protein solubility. Appropriate protein extraction and modification methods that meet consumer and manufacturers requirements concerning nutritious and eco-friendly foods with lower cost should be investigated and further explored.

First 3-D structural evidence of a native-like intertwined dimer in the acylphosphatase family

Acylphosphatase (AcP, EC 3.6.1.7) is a small model protein conformed by a ferredoxin-like fold, profoundly studied to get insights into protein folding and aggregation processes. Numerous studies focused on the aggregation and/or amyloidogenic properties of AcPs suggest the importance of edge-β-strands in the process. In this work, we present the first crystallographic structure of Escherichia coli AcP (EcoAcP), showing notable differences with the only available NMR structure for this enzyme. EcoAcP is crystalised as an intertwined dimer formed by replacing a single C-terminal β-strand between two protomers, suggesting a flexible character of the C-terminal edge of EcoAcP. Despite numerous works where AcP from different sources have been used as a model system for protein aggregation, our domain-swapped EcoAcP structure is the first 3-D structural evidence of native-like aggregated species for any AcP reported to date, providing clues on molecular determinants unleashing aggregation.

Metal interactions of α-synuclein probed by NMR amide-proton exchange

The aberrant aggregation of α-synuclein (αS), a disordered protein primarily expressed in neuronal cells, is strongly associated with the underlying mechanisms of Parkinson's disease. It is now established that αS has a weak affinity for metal ions and that these interactions alter its conformational properties by generally promoting self-assembly into amyloids. Here, we characterised the nature of the conformational changes associated with metal binding by αS using nuclear magnetic resonance (NMR) to measure the exchange of the backbone amide protons at a residue specific resolution. We complemented these experiments with ¹⁵N relaxation and chemical shift perturbations to obtain a comprehensive map of the interaction between αS and divalent (Ca²⁺, Cu²⁺, Mn²⁺, and Zn²⁺) and monovalent (Cu⁺) metal ions. The data identified specific effects that the individual cations exert on the conformational properties of αS. In particular, binding to calcium and zinc generated a reduction of the protection factors in the C-terminal region of the protein, whereas both Cu(II) and Cu(I) did not alter the amide proton exchange along the αS sequence. Changes in the R₂/R₁ ratios from ¹⁵N relaxation experiments were, however, detected as a result of the interaction between αS and Cu⁺ or Zn²⁺, indicating that binding to these metals induces conformational perturbations in distinctive regions of the protein. Collectively our data suggest that multiple mechanisms of enhanced αS aggregation are associated with the binding of the analysed metals.

The role of structural dynamics in the thermal adaptation of hyperthermophilic enzymes

Proteins from hyperthermophilic organisms are evolutionary optimised to adopt functional structures and dynamics under conditions in which their mesophilic homologues are generally inactive or unfolded. Understanding the nature of such adaptation is of crucial interest to clarify the underlying mechanisms of biological activity in proteins. Here we measured NMR residual dipolar couplings of a hyperthermophilic acylphosphatase enzyme at 80°C and used these data to generate an accurate structural ensemble representative of its native state. The resulting energy landscape was compared to that obtained for a human homologue at 37°C, and additional NMR experiments were carried out to probe fast (¹⁵N relaxation) and slow (H/D exchange) backbone dynamics, collectively sampling fluctuations of the two proteins ranging from the nanosecond to the millisecond timescale. The results identified key differences in the strategies for protein-protein and protein-ligand interactions of the two enzymes at the respective physiological temperatures. These include the dynamical behaviour of a β-strand involved in the protection against aberrant protein aggregation and concerted motions of loops involved in substrate binding and catalysis. Taken together these results elucidate the structure-dynamics-function relationship associated with the strategies of thermal adaptation of protein molecules.

Quantitative Measurement of the Affinity of Toxic and Nontoxic Misfolded Protein Oligomers for Lipid Bilayers and of its Modulation by Lipid Composition and Trodusquemine

Many neurodegenerative diseases are associated with the self-assembly of peptides and proteins into fibrillar aggregates. Soluble misfolded oligomers formed during the aggregation process, or released by mature fibrils, play a relevant role in neurodegenerative processes through their interactions with neuronal membranes. However, the determinants of the cytotoxicity of these oligomers are still unclear. Here we used liposomes and toxic and nontoxic oligomers formed by the same protein to measure quantitatively the affinity of the two oligomeric species for lipid membranes. To this aim, we quantified the perturbation to the lipid membranes caused by the two oligomers by using the fluorescence quenching of two probes embedded in the polar and apolar regions of the lipid membranes and a well-defined protein-oligomer binding assay using fluorescently labeled oligomers to determine the Stern-Volmer and dissociation constants, respectively. With both approaches, we found that the toxic oligomers have a membrane affinity 20-25 times higher than that of nontoxic oligomers. Circular dichroism, intrinsic fluorescence, and FRET indicated that neither oligomer type changes its structure upon membrane interaction. Using liposomes enriched with trodusquemine, a potential small molecule drug known to penetrate lipid membranes and make them refractory to toxic oligomers, we found that the membrane affinity of the oligomers was remarkably lower. At protective concentrations of the small molecule, the binding of the oligomers to the lipid membranes was fully prevented. Furthermore, the affinity of the toxic oligomers for the lipid membranes was found to increase and slightly decrease with GM1 ganglioside and cholesterol content, respectively, indicating that physicochemical properties of lipid membranes modulate their affinity for misfolded oligomeric species.

Understanding the Role of Protein Glycation in the Amyloid Aggregation Process

Protein function and flexibility is directly related to the native distribution of its structural elements and any alteration in protein architecture leads to several abnormalities and accumulation of misfolded proteins. This phenomenon is associated with a range of increasingly common human disorders, including Alzheimer and Parkinson diseases, type II diabetes, and a number of systemic amyloidosis characterized by the accumulation of amyloid aggregates both in the extracellular space of tissues and as intracellular deposits. Post-translational modifications are known to have an active role in the in vivo amyloid aggregation as able to affect protein structure and dynamics. Among them, a key role seems to be played by non-enzymatic glycation, the most unwanted irreversible modification of the protein structure, which strongly affects long-living proteins throughout the body. This study provided an overview of the molecular effects induced by glycation on the amyloid aggregation process of several protein models associated with misfolding diseases. In particular, we analyzed the role of glycation on protein folding, kinetics of amyloid formation, and amyloid cytotoxicity in order to shed light on the role of this post-translational modification in the in vivo amyloid aggregation process.

Divergence Entropy-Based Evaluation of Hydrophobic Core in Aggressive and Resistant Forms of Transthyretin

The two forms of transthyretin differing slightly in the tertiary structure, despite the presence of five mutations, show radically different properties in terms of susceptibility to the amyloid transformation process. These two forms of transthyretin are the object of analysis. The search for the sources of these differences was carried out by means of a comparative analysis of the structure of these molecules in their native and early intermediate stage forms in the folding process. The criterion for assessing the degree of similarity and differences is the status of the hydrophobic core. The comparison of the level of arrangement of the hydrophobic core and its initial stages is possible thanks to the application of divergence entropy for the early intermediate stage and for the final forms. It was shown that the minimal differences observed in the structure of the hydrophobic core of the forms available in PDB, turned out to be significantly different in the early stage (ES) structure in folding process. The determined values of divergence entropy for both ES forms indicate the presence of the seed of hydrophobic core only in the form resistant to amyloid transformation. In the form of aggressively undergoing amyloid transformation, the structure lacking such a seed is revealed, being a stretched one with a high content of β-type structure. In the discussed case, the active presence of water in the structural transformation of proteins expressed in the fuzzy oil drop model (FOD) is of decisive importance for the generation of the final protein structure. It has been shown that the resistant form tends to generate a centric hydrophobic core with the possibility of creating a globular structure, i.e., a spherical micelle-like form. The aggressively transforming form reveals in the structure of its early intermediate, a tendency to form the ribbon-like micelle as observed in amyloid.

Structure-aware protein solubility prediction from sequence through graph convolutional network and predicted contact map

Protein solubility is significant in producing new soluble proteins that can reduce the cost of biocatalysts or therapeutic agents. Therefore, a computational model is highly desired to accurately predict protein solubility from the amino acid sequence. Many methods have been developed, but they are mostly based on the one-dimensional embedding of amino acids that is limited to catch spatially structural information. In this study, we have developed a new structure-aware method GraphSol to predict protein solubility by attentive graph convolutional network (GCN), where the protein topology attribute graph was constructed through predicted contact maps only from the sequence. GraphSol was shown to substantially outperform other sequence-based methods. The model was proven to be stable by consistent [Formula: see text] of 0.48 in both the cross-validation and independent test of the eSOL dataset. To our best knowledge, this is the first study to utilize the GCN for sequence-based protein solubility predictions. More importantly, this architecture could be easily extended to other protein prediction tasks requiring a raw protein sequence.

Soluble Prion Peptide 107-120 Protects Neuroblastoma SH-SY5Y Cells against Oligomers Associated with Alzheimer's Disease

Alzheimer's disease (AD) is the most prevalent form of dementia and soluble amyloid β (Aβ) oligomers are thought to play a critical role in AD pathogenesis. Cellular prion protein (PrPC) is a high-affinity receptor for Aβ oligomers and mediates some of their toxic effects. The N-terminal region of PrPC can interact with Aβ, particularly the region encompassing residues 95-110. In this study, we identified a soluble and unstructured prion-derived peptide (PrP107-120) that is external to this region of the sequence and was found to successfully reduce the mitochondrial impairment, intracellular ROS generation and cytosolic Ca2+ uptake induced by oligomeric Aβ42 ADDLs in neuroblastoma SH-SY5Y cells. PrP107-120 was also found to rescue SH-SY5Y cells from Aβ42 ADDL internalization. The peptide did not change the structure and aggregation pathway of Aβ42 ADDLs, did not show co-localization with Aβ42 ADDLs in the cells and showed a partial colocalization with the endogenous cellular PrPC. As a sequence region that is not involved in Aβ binding but in PrP self-recognition, the peptide was suggested to protect against the toxicity of Aβ42 oligomers by interfering with cellular PrPC and/or activating a signaling that protected the cells. These results strongly suggest that PrP107-120 has therapeutic potential for AD.

The Structure of Amyloid Versus the Structure of Globular Proteins

The issue of changing the structure of globular proteins into an amyloid form is in the focus of researchers' attention. Numerous experimental studies are carried out, and mathematical models to define the essence of amyloid transformation are sought. The present work focuses on the issue of the hydrophobic core structure in amyloids. The form of ordering the hydrophobic core in globular proteins is described by a 3D Gaussian distribution analog to the distribution of hydrophobicity in a spherical micelle. Amyloid fibril is a ribbon-like micelle made up of numerous individual chains, each representing a flat structure. The distribution of hydrophobicity within a single chain included in the fibril describes the 2D Gaussian distribution. Such a description expresses the location of polar residues on a circle with a center with a high level of hydrophobicity. The presence of this type of order in the amyloid forms available in Preotin Data Bank (PDB) (both in proto- and superfibrils) is demonstrated in the present work. In this system, it can be assumed that the amyloid transformation is a chain transition from 3D Gauss ordering to 2D Gauss ordering. This means changing the globular structure to a ribbon-like structure. This observation can provide a simple mathematical model for simulating the amyloid transformation of proteins.

The Status of Edge Strands in Ferredoxin-Like Fold

There is an opinion in professional literature that edge-strands in β-sheet are critical to the processes of amyloid transformation. Propagation of fibrillar forms mainly takes place on the basis of β-sheet type interactions. In many proteins, the edge strands represent only a partially matched form to the β-sheet. Therefore, the edge-strand takes slightly distorted forms. The assessment of the level of arrangement can be carried out based on studying the secondary structure as well as the structure of the hydrophobic core. For this purpose, a fuzzy oil drop model was used to determine the contribution of each fragment with a specific secondary structure to the construction of the system being the effect of a certain synergy, which results in the construction of a hydrophobic core. Studying the participation of β-sheets edge fragments in the hydrophobic core construction is the subject of the current analysis. Statuses of these edge fragments in β-sheets in ferredoxin-like folds are treated as factors that disturb the symmetry of the system.

An NMR based phosphodiesterase assay

We propose a phosphodiesterase assay based on 1D ¹H NMR to monitor the hydrolysis of cyclic nucleotides directly, without requiring tags or the addition of exogenous reagents. The method is suitable to measure phosphodiesterase K_M and k_cat parameters and to identify phosphodiesterase inhibitors.

Changes in hydrophobicity mainly promotes the aggregation tendency of ALS associated SOD1 mutants

Protein misfolding and aggregation due to mutations, are associated with fatal neurodegenerative disorders. The mutations in Cu/Zn superoxide dismutase (SOD1) causing its misfolding and aggregation are found linked to the motor neuron disorder, amyotrophic lateral sclerosis. Since the mutations are scattered throughout SOD1 structure, determining the exact molecular mechanism underlying the ALS pathology remains unresolved. In this study, we have investigated the major molecular factors that mainly contribute to SOD1 destabilization, intrinsic disorder, and misfolding using sequence and structural information. We have analysed 153 ALS causing SOD1 point mutants for aggregation tendency using four different aggregation prediction tools, viz., Aggrescan3D (A3D), CamSol, GAP and Zyggregator. Our results suggest that 74-79 mutants are susceptible to aggregation, due to distorted native interactions originated at the mutation site. Majority of the aggregation prone mutants are located in the buried regions of SOD1 molecule. Further, the mutations at the hydrophobic amino acids primarily promote the aggregation tendency of SOD1 protein through different destabilizing mechanisms including changes in hydrophobic free energy, loss of electrostatic interactions in the protein's surface and loss of hydrogen bonds that bridges the protein core and surface.

Microsecond sub-domain motions and the folding and misfolding of the mouse prion protein

Protein aggregation appears to originate from partially unfolded conformations that are sampled through stochastic fluctuations of the native protein. It has been a challenge to characterize these fluctuations, under native like conditions. Here, the conformational dynamics of the full-length (23-231) mouse prion protein were studied under native conditions, using photoinduced electron transfer coupled to fluorescence correlation spectroscopy (PET-FCS). The slowest fluctuations could be associated with the folding of the unfolded state to an intermediate state, by the use of microsecond mixing experiments. The two faster fluctuations observed by PET-FCS, could be attributed to fluctuations within the native state ensemble. The addition of salt, which is known to initiate the aggregation of the protein, resulted in an enhancement in the time scale of fluctuations in the core of the protein. The results indicate the importance of native state dynamics in initiating the aggregation of proteins.

Backbone NMR assignments of HypF-N under conditions generating toxic and non-toxic oligomers

The HypF protein is involved in the maturation and regulation of hydrogenases. The N-terminal domain of HypF (HypF-N) has served as a key model system to study the pathways of protein amyloid formation and the nature of the toxicity of pre-fibrilar protein oligomers. This domain can aggregate into two forms of oligomers having significantly different toxic effects when added to neuronal cultures. Here, NMR assignments of HypF-N backbone resonances are presented in its native state and under the conditions favouring the formation of toxic and non-toxic oligomers. The analyses of chemical shifts provide insights into the protein conformational state and the possible pathways leading to the formation of different types of oligomers.

Signature of Pareto optimization in the Escherichia coli proteome

Proteins have coevolved with cellular environments to improve or preserve their functions, maintaining at the same time the degree of hydrophobicity necessary to fold correctly and enough solubility to perform their biological roles. Here, we study the Escherichia coli proteome using a Pareto front analysis in the solubility-hydrophobicity space. The results indicate the existence of a Pareto optimal front, a triangle whose vertices correspond to archetypal proteins specialized in distinct tasks, such as regulatory processes, membrane transport, outer-membrane pore formation, catalysis, and binding. The vertices are further enriched with proteins that occupy different subcellular compartments, namely, cytoplasmic, inner membrane, outer membrane, and outer membrane bounded periplasmic space. The combination of various enriching features offers an interpretation of how bacteria use the physico-chemical properties of proteins, both to drive them into their final destination in the cell and to have their tasks accomplished.

The structural basis of nanobody unfolding reversibility and thermoresistance

Nanobodies represent the variable binding domain of camelid heavy-chain antibodies and are employed in a rapidly growing range of applications in biotechnology and biomedicine. Their success is based on unique properties including their reported ability to reversibly refold after heat-induced denaturation. This view, however, is contrasted by studies which involve irreversibly aggregating nanobodies, asking for a quantitative analysis that clearly defines nanobody thermoresistance and reveals the determinants of unfolding reversibility and aggregation propensity. By characterizing nearly 70 nanobodies, we show that irreversible aggregation does occur upon heat denaturation for the large majority of binders, potentially affecting application-relevant parameters like stability and immunogenicity. However, by deriving aggregation propensities from apparent melting temperatures, we show that an optional disulfide bond suppresses nanobody aggregation. This effect is further enhanced by increasing the length of a complementarity determining loop which, although expected to destabilize, contributes to nanobody stability. The effect of such variations depends on environmental conditions, however. Nanobodies with two disulfide bonds, for example, are prone to lose their functionality in the cytosol. Our study suggests strategies to engineer nanobodies that exhibit optimal performance parameters and gives insights into general mechanisms which evolved to prevent protein aggregation.

Concurrent structural and biophysical traits link with immunoglobulin light chains amyloid propensity

Light chain amyloidosis (AL), the most common systemic amyloidosis, is caused by the overproduction and the aggregation of monoclonal immunoglobulin light chains (LC) in target organs. Due to genetic rearrangement and somatic hypermutation, virtually, each AL patient presents a different amyloidogenic LC. Because of such complexity, the fine molecular determinants of LC aggregation propensity and proteotoxicity are, to date, unclear; significantly, their decoding requires investigating large sets of cases. Aiming to achieve generalizable observations, we systematically characterised a pool of thirteen sequence-diverse full length LCs. Eight amyloidogenic LCs were selected as responsible for severe cardiac symptoms in patients; five non-amyloidogenic LCs were isolated from patients affected by multiple myeloma. Our comprehensive approach (consisting of spectroscopic techniques, limited proteolysis, and X-ray crystallography) shows that low fold stability and high protein dynamics correlate with amyloidogenic LCs, while hydrophobicity, structural rearrangements and nature of the LC dimeric association interface (as observed in seven crystal structures here presented) do not appear to play a significant role in defining amyloid propensity. Based on the structural and biophysical data, our results highlight shared properties driving LC amyloid propensity, and these data will be instrumental for the design of synthetic inhibitors of LC aggregation.

Comparing the energy landscapes for native folding and aggregation of PrP

Protein sequences are evolved to encode generally one folded structure, out of a nearly infinite array of possible folds. Underlying this code is a funneled free energy landscape that guides folding to the native conformation. Protein misfolding and aggregation are also a manifestation of free-energy landscapes. The detailed mechanisms of these processes are poorly understood, but often involve rare, transient species and a variety of different pathways. The inherent complexity of misfolding has hampered efforts to measure aggregation pathways and the underlying energy landscape, especially using traditional methods where ensemble averaging obscures important rare and transient events. We recently studied the misfolding and aggregation of prion protein by examining 2 monomers tethered in close proximity as a dimer, showing how the steps leading to the formation of a stable aggregated state can be resolved in the single-molecule limit and the underlying energy landscape thereby reconstructed. This approach allows a more quantitative comparison of native folding versus misfolding, including fundamental differences in the dynamics for misfolding. By identifying key steps and interactions leading to misfolding, it should help to identify potential drug targets. Here we describe the importance of characterizing free-energy landscapes for aggregation and the challenges involved in doing so, and we discuss how single-molecule studies can help test proposed structural models for PrP aggregates.

Direct Conversion of an Enzyme from Native-like to Amyloid-like Aggregates within Inclusion Bodies

The acylphosphatase from Sulfolobus solfataricus (Sso AcP) is a globular protein able to aggregate in vitro from a native-like conformational ensemble without the need for a transition across the major unfolding energy barrier. This process leads to the formation of assemblies in which the protein retains its native-like structure, which subsequently convert into amyloid-like aggregates. Here, we investigate the mechanism by which Sso AcP aggregates in vivo to form bacterial inclusion bodies after expression in E. coli. Shortly after the initiation of expression, Sso AcP is incorporated into inclusion bodies as a native-like protein, still exhibiting small but significant enzymatic activity. Additional experiments revealed that this overall process of aggregation is enhanced by the presence of the unfolded N-terminal region of the sequence and by destabilization of the globular segment of the protein. At later times, the Sso AcP molecules in the inclusion bodies lose their native-like properties and convert into β-sheet-rich amyloid-like structures, as indicated by their ability to bind thioflavin T and Congo red. These results show that the aggregation behavior of this protein is similar in vivo to that observed in vitro, and that, at least for a predominant part of the protein population, the transition from a native to an amyloid-like structure occurs within the aggregate state.

Water is an active matrix of life for cell and molecular biology

Szent-Győrgi called water the "matrix of life" and claimed that there was no life without it. This statement is true, as far as we know, on our planet, but it is not clear whether it must hold throughout the cosmos. To evaluate that question requires a close consideration of the many varied and subtle roles that water plays in living cells-a consideration that must be free of both an assumed essentialism that gives water an almost mystical life-giving agency and a traditional tendency to see it as a merely passive solvent. Water is a participant in the "life of the cell," and here I describe some of the features of that active agency. Water's value for molecular biology comes from both the structural and dynamic characteristics of its status as a complex, structured liquid as well as its nature as a polar, protic, and amphoteric reagent. Any discussion of water as life's matrix must, however, begin with an acknowledgment that our understanding of it as both a liquid and a solvent is still incomplete.

DNA intercalators as amyloid assembly modulators: mechanistic insights

DNA intercalators modulate amyloid assembly of proteins through specific hetero-aromatic interactions diverting them to form amorphous aggregates.

From the Evolution of Protein Sequences Able to Resist Self-Assembly to the Prediction of Aggregation Propensity

Folding of polypeptide chains into biologically active entities is an astonishingly complex process, determined by the nature and the sequence of residues emerging from ribosomes. While it has been long believed that evolution has pressed genomes so that specific sequences could adopt unique, functional three-dimensional folds, it is now clear that complex protein machineries act as quality control system and supervise folding. Notwithstanding that, events such as erroneous folding, partial folding, or misfolding are frequent during the life of a cell or a whole organism, and they can escape controls. One of the possible outcomes of this misbehavior is cross-β aggregation, a super secondary structure which represents the hallmark of self-assembled, well organized, and extremely ordered structures termed amyloid fibrils. What if evolution would have not taken into account such possibilities? Twenty years of research point toward the idea that, in fact, evolution has constantly supervised the risk of errors and minimized their impact. In this review we tried to survey the major findings in the amyloid field, trying to describe what the real pitfalls of protein folding are-from an evolutionary perspective-and how sequence and structural features have evolved to balance the need for perfect, dynamic, functionally efficient structures, and the detrimental effects implicit in the dangerous process of folding. We will discuss how the knowledge obtained from these studies has been employed to produce computational methods able to assess, predict, and discriminate the aggregation properties of protein sequences.

Antibody Aggregation: Insights from Sequence and Structure

Monoclonal antibodies (mAbs) are the fastest-growing biological therapeutics with important applications ranging from cancers, autoimmunity diseases and metabolic disorders to emerging infectious diseases. Aggregation of mAbs continues to be a major problem in their developability. Antibody aggregation could be triggered by partial unfolding of its domains, leading to monomer-monomer association followed by nucleation and growth. Although the aggregation propensities of antibodies and antibody-based proteins can be affected by the external experimental conditions, they are strongly dependent on the intrinsic antibody properties as determined by their sequences and structures. In this review, we describe how the unfolding and aggregation susceptibilities of IgG could be related to their cognate sequences and structures. The impact of antibody domain structures on thermostability and aggregation propensities, and effective strategies to reduce aggregation are discussed. Finally, the aggregation of antibody-drug conjugates (ADCs) as related to their sequence/structure, linker payload, conjugation chemistry and drug-antibody ratio (DAR) is reviewed.

Rational design of mutations that change the aggregation rate of a protein while maintaining its native structure and stability

A wide range of human diseases is associated with mutations that, destabilizing proteins native state, promote their aggregation. However, the mechanisms leading from folded to aggregated states are still incompletely understood. To investigate these mechanisms, we used a combination of NMR spectroscopy and molecular dynamics simulations to compare the native state dynamics of Beta-2 microglobulin (β2m), whose aggregation is associated with dialysis-related amyloidosis, and its aggregation-resistant mutant W60G. Our results indicate that W60G low aggregation propensity can be explained, beyond its higher stability, by an increased average protection of the aggregation-prone residues at its surface. To validate these findings, we designed β2m variants that alter the aggregation-prone exposed surface of wild-type and W60G β2m modifying their aggregation propensity. These results allowed us to pinpoint the role of dynamics in β2m aggregation and to provide a new strategy to tune protein aggregation by modulating the exposure of aggregation-prone residues.

Accurate Determination of Conformational Transitions in Oligomeric Membrane Proteins

The structural dynamics governing collective motions in oligomeric membrane proteins play key roles in vital biomolecular processes at cellular membranes. In this study, we present a structural refinement approach that combines solid-state NMR experiments and molecular simulations to accurately describe concerted conformational transitions identifying the overall structural, dynamical, and topological states of oligomeric membrane proteins. The accuracy of the structural ensembles generated with this method is shown to reach the statistical error limit, and is further demonstrated by correctly reproducing orthogonal NMR data. We demonstrate the accuracy of this approach by characterising the pentameric state of phospholamban, a key player in the regulation of calcium uptake in the sarcoplasmic reticulum, and by probing its dynamical activation upon phosphorylation. Our results underline the importance of using an ensemble approach to characterise the conformational transitions that are often responsible for the biological function of oligomeric membrane protein states.

The inverted free energy landscape of an intrinsically disordered peptide by simulations and experiments

The free energy landscape theory has been very successful in rationalizing the folding behaviour of globular proteins, as this representation provides intuitive information on the number of states involved in the folding process, their populations and pathways of interconversion. We extend here this formalism to the case of the Aβ40 peptide, a 40-residue intrinsically disordered protein fragment associated with Alzheimer’s disease. By using an advanced sampling technique that enables free energy calculations to reach convergence also in the case of highly disordered states of proteins, we provide a precise structural characterization of the free energy landscape of this peptide. We find that such landscape has inverted features with respect to those typical of folded proteins. While the global free energy minimum consists of highly disordered structures, higher free energy regions correspond to a large variety of transiently structured conformations with secondary structure elements arranged in several different manners, and are not separated from each other by sizeable free energy barriers. From this peculiar structure of the free energy landscape we predict that this peptide should become more structured and not only more compact, with increasing temperatures, and we show that this is the case through a series of biophysical measurements.

AGGRESCAN3D (A3D): server for prediction of aggregation properties of protein structures

Protein aggregation underlies an increasing number of disorders and constitutes a major bottleneck in the development of therapeutic proteins. Our present understanding on the molecular determinants of protein aggregation has crystalized in a series of predictive algorithms to identify aggregation-prone sites. A majority of these methods rely only on sequence. Therefore, they find difficulties to predict the aggregation properties of folded globular proteins, where aggregation-prone sites are often not contiguous in sequence or buried inside the native structure. The AGGRESCAN3D (A3D) server overcomes these limitations by taking into account the protein structure and the experimental aggregation propensity scale from the well-established AGGRESCAN method. Using the A3D server, the identified aggregation-prone residues can be virtually mutated to design variants with increased solubility, or to test the impact of pathogenic mutations. Additionally, A3D server enables to take into account the dynamic fluctuations of protein structure in solution, which may influence aggregation propensity. This is possible in A3D Dynamic Mode that exploits the CABS-flex approach for the fast simulations of flexibility of globular proteins. The A3D server can be accessed at http://biocomp.chem.uw.edu.pl/A3D/.

Understanding the effect of alternative splicing in the folding and function of the second PDZ from protein tyrosine phosphatase-BL

PDZ domains are the most prominent biological structural domains involved in protein-protein interactions in the human cell. The second PDZ domain of the protein tyrosine phosphatase BL (PDZ2) interacts and binds the C-termini of the tumour suppressor protein APC and of the LIM domain-containing protein RIL. One isoform of PDZ2 (PDZ2as) involves an alternative spliced form that exhibits an insertion of 5 residues in a loop. PDZ2as abrogates binding to its partners, even if the insertion is directly located in its binding pocket. Here, we investigate the folding and function of PDZ2as, in comparison to the previously characterized PDZ2 domain. Data reveal that, whilst the thermodynamic stability of PDZ2as appears as nearly identical to that of PDZ2, the insertion of 5 amino acids induces formation of some weak transient non-native interactions in the folding transition state, as mirrored by a concomitant increase of both the folding and unfolding rate constants. From a functional perspective, we show that the decrease in affinity is caused by a pronounced decrease of the association rate constants (by nearly ten fold), with no effect on the microscopic dissociation rate constants. The results are briefly discussed in the context of previous work on PDZ domains.

Structural dynamics and conformational equilibria of SERCA regulatory proteins in membranes by solid-state NMR restrained simulations

Solid-state NMR spectroscopy is emerging as a powerful approach to determine structure, topology, and conformational dynamics of membrane proteins at the atomic level. Conformational dynamics are often inferred and quantified from the motional averaging of the NMR parameters. However, the nature of these motions is difficult to envision based only on spectroscopic data. Here, we utilized restrained molecular dynamics simulations to probe the structural dynamics, topology and conformational transitions of regulatory membrane proteins of the calcium ATPase SERCA, namely sarcolipin and phospholamban, in explicit lipid bilayers. Specifically, we employed oriented solid-state NMR data, such as dipolar couplings and chemical shift anisotropy measured in lipid bicelles, to refine the conformational ensemble of these proteins in lipid membranes. The samplings accurately reproduced the orientations of transmembrane helices and showed a significant degree of convergence with all of the NMR parameters. Unlike the unrestrained simulations, the resulting sarcolipin structures are in agreement with distances and angles for hydrogen bonds in ideal helices. In the case of phospholamban, the restrained ensemble sampled the conformational interconversion between T (helical) and R (unfolded) states for the cytoplasmic region that could not be observed using standard structural refinements with the same experimental data set. This study underscores the importance of implementing NMR data in molecular dynamics protocols to better describe the conformational landscapes of membrane proteins embedded in realistic lipid membranes.

Conformational recognition of an intrinsically disordered protein

There is a growing interest in understanding the properties of intrinsically disordered proteins (IDPs); however, the characterization of these states remains an open challenge. IDPs appear to have functional roles that diverge from those of folded proteins and revolve around their ability to act as hubs for protein-protein interactions. To gain a better understanding of the modes of binding of IDPs, we combined statistical mechanics, calorimetry, and NMR spectroscopy to investigate the recognition and binding of a fragment from the disordered protein Gab2 by the growth factor receptor-bound protein 2 (Grb2), a key interaction for normal cell signaling and cancer development. Structural ensemble refinement by NMR chemical shifts, thermodynamics measurements, and analysis of point mutations indicated that the population of preexisting bound conformations in the free-state ensemble of Gab2 is an essential determinant for recognition and binding by Grb2. A key role was found for transient polyproline II (PPII) structures and extended conformations. Our findings are likely to have very general implications for the biological behavior of IDPs in light of the evidence that a large fraction of these proteins possess a specific propensity to form PPII and to adopt conformations that are more extended than the typical random-coil states.

Thermodynamic characterization of the unfolding of the prion protein

The prion protein appears to be unusually susceptible to conformational change, and unlike nearly all other proteins, it can easily be made to convert to alternative misfolded conformations. To understand the basis of this structural plasticity, a detailed thermodynamic characterization of two variants of the mouse prion protein (moPrP), the full-length moPrP (23-231) and the structured C-terminal domain, moPrP (121-231), has been carried out. All thermodynamic parameters governing unfolding, including the changes in enthalpy, entropy, free energy, and heat capacity, were found to be identical for the two protein variants. The N-terminal domain remains unstructured and does not interact with the C-terminal domain in the full-length protein at pH 4. Moreover, the enthalpy and entropy of unfolding of moPrP (121-231) are similar in magnitude to values reported for other proteins of similar size. However, the protein has an unusually high native-state heat capacity, and consequently, the change in heat capacity upon unfolding is much lower than that expected for a protein of similar size. It appears, therefore, that the native state of the prion protein undergoes substantial fluctuations in enthalpy and hence, in structure.

Differential effects of glycation on protein aggregation and amyloid formation

Amyloids are a class of insoluble proteinaceous substances generally composed of linear un-branched fibrils that are formed from misfolded proteins. Conformational diseases such as Alzheimer's disease, transmissible spongiform encephalopathies, and familial amyloidosis are associated with the presence of amyloid aggregates in the affected tissues. The majority of the cases are sporadic, suggesting that several factors must contribute to the onset and progression of these disorders. Among them, in the past 10 years, non-enzymatic glycation of proteins has been reported to stimulate protein aggregation and amyloid deposition. In this review, we analyze the most recent advances in this field suggesting that the effects induced by glycation may not be generalized as strongly depending on the protein structure. Indeed, being a post-translational modification, glycation could differentially affects the aggregation process in promoting, accelerating and/or stabilizing on-pathway and off-pathway species.

Edge strand engineering prevents native-like aggregation in Sulfolobus solfataricus acylphosphatase

β-proteins are constantly threatened by the risk of aggregation because β-sheets are inherently structured for edge-to-edge interactions. To avoid native-like aggregation, evolution has resulted in a set of strategies that prevent intermolecular β-interactions. Acylphosphatase from Sulfolobus solfataricus (Sso AcP) represents a suitable model for the study of such a process. Under conditions promoting aggregation, Sso AcP acquires a native-like conformational state whereby an unstructured N-terminal segment interacts with the edge β-strand B4 of an adjacent Sso AcP molecule. Because B4 is poorly protected against aggregation, this interaction triggers the aggregation cascade without the need for unfolding. Recently, three single Sso AcP mutants (V84D, Y86E and V84P) were designed to engineer additional protection against aggregation in B4 and were observed to successfully impair native-like aggregation in all three variants at the expense of a lower stability. To understand the structural basis of the reduced aggregation propensity and lower stability, the crystal structures of the Sso AcP variants were determined in the present study. Structural analysis reveals that the V84D and Y86E mutations exert protection by the insertion of an edge negative charge. A conformationally less regular B4 underlies protection against aggregation in the V84P mutant. The thermodynamic basis of instability is discussed. Moreover, kinetic experiments indicate that aggregation of the three mutants is not native-like and is independent of the interaction between B4 and the unstructured N-terminal segment. The reported data rationalize previous evidence regarding Sso AcP native-like aggregation and provide a basis for the design of aggregation-free proteins.

DATABASE

The atomic coordinates and related experimental data for the Sso AcP mutants V84P, V84D, ΔN11 Y86E have been deposited in the Protein Data Bank under accession numbers 4OJ3, 4OJG and 4OJH, respectively.

STRUCTURED DIGITAL ABSTRACT

• Sso AcP and Sso AcP bind by fluorescence technology (View interaction).

Propensity to form amyloid fibrils is encoded as excitations in the free energy landscape of monomeric proteins

Protein aggregation, linked to many of diseases, is initiated when monomers access rogue conformations that are poised to form amyloid fibrils. We show, using simulations of src SH3 domain, that mechanical force enhances the population of the aggregation-prone (N(⁎)) states, which are rarely populated under force free native conditions but are encoded in the spectrum of native fluctuations. The folding phase diagrams of SH3 as a function of denaturant concentration ([C]), mechanical force (f), and temperature exhibit an apparent two-state behavior, without revealing the presence of the elusive N(⁎) states. Interestingly, the phase boundaries separating the folded and unfolded states at all [C] and f fall on a master curve, which can be quantitatively described using an analogy to superconductors in a magnetic field. The free energy profiles as a function of the molecular extension (R), which are accessible in pulling experiments, (R), reveal the presence of a native-like N(⁎) with a disordered solvent-exposed amino-terminal β-strand. The structure of the N(⁎) state is identical with that found in Fyn SH3 by NMR dispersion experiments. We show that the timescale for fibril formation can be estimated from the population of the N(⁎) state, determined by the free energy gap separating the native structure and the N(⁎) state, a finding that can be used to assess fibril forming tendencies of proteins. The structures of the N(⁎) state are used to show that oligomer formation and likely route to fibrils occur by a domain-swap mechanism in SH3 domain.

Role of loops connecting secondary structure elements in the stabilization of proteins isolated from thermophilic organisms

It has been recently discovered that the connection of secondary structure elements (ββ-unit, βα- and αβ-units) in proteins follows quite stringent principles regarding the chirality and the orientation of the structural units (Koga et al., Nature 2012;491:222-227). By exploiting these rules, a number of protein scaffolds endowed with a remarkable thermal stability have been designed (Koga et al., Nature 2012;491:222-227). By using structural databases of proteins isolated from either mesophilic or thermophilic organisms, we here investigate the influence of supersecondary associations on the thermal stability of natural proteins. Our results suggest that β-hairpins of proteins from thermophilic organisms are very frequently characterized by shortenings of the loops. Interestingly, this shortening leads to states that display a very strong preference for the most common connectivity of the strands observed in native protein hairpins. The abundance of selective states in these proteins suggests that they may achieve a high stability by adopting a strategy aimed to reduce the possible conformations of the unfolded ensemble. In this scenario, our data indicate that the shortening is effective if it increases the adherence to these rules. We also show that this mechanism may operate in the stabilization of well-known protein folds (thioredoxin and RNase A). These findings suggest that future investigations aimed at defining mechanism of protein stabilization should also consider these effects.

Long-range electron transfer with myoglobin immobilized at Au/mixed-SAM junctions: mechanistic impact of the strong protein confinement

Horse muscle myoglobin (Mb) was tightly immobilized at Au-deposited ~15-Å-thick mixed-type (1:1) alkanethiol SAMs, HS-(CH₂)₁₁-COOH/HS-(CH₂)₁₁-OH, and placed in contact with buffered H₂O or D₂O solutions. Fast-scan cyclic voltammetry (CV) and a Marcus-equation-based analysis were applied to determine unimolecular standard rate constants and reorganization free energies for electron transfer (ET), under variable-temperature (15-55 °C) and -pressure (0.01-150 MPa) conditions. The CV signal was surprisingly stable and reproducible even after multiple temperature and pressure cycles. The data analysis revealed the following values: standard rate constant, 33 s⁻¹ (25 °C, 0.01 MPa, H₂O); reorganization free energy, 0.5 ± 0.1 eV (throughout); activation enthalpy, 12 ± 3 kJ mol⁻¹; activation volume, -3.1 ± 0.2 cm³ mol⁻¹; and pH-dependent solvent kinetic isotope effect (k(H)⁰/k(D)⁰), 0.7-1.4. Furthermore, the values for the rate constant and reorganization free energy are very similar to those previously found for cytochrome c electrostatically immobilized at the monocomponent Au/HS-(CH₂)₁₁-COOH junction. In vivo, Mb apparently forms a natural electrostatic complex with cytochrome b₅ (cyt-b₅) through the "dynamic" (loose) docking pattern, allowing for a slow ET that is intrinsically coupled to the water's removal from the "defective" heme iron (altogether shaping the biological repair mechanism for Mb's "met" form). In contrary, our experiments rather mimic the case of a "simple" (tight) docking of the redesigned (mutant) Mb with cyt-b₅ (Nocek et al. J. Am. Chem. Soc. 2010, 132, 6165-6175). According to our analysis, in this configuration, Mb's distal pocket (linked to the "ligand channel") seems to be arrested within the restricted configuration, allowing the rate-determining reversible ET process to be coupled only to the inner-sphere reorganization (minimal elongation/shortening of an Fe-OH₂ bond) rather than the pronounced detachment (rebinding) of water and, hence, to be much faster.

Misfolding of amyloidogenic proteins and their interactions with membranes

In this paper, we discuss amyloidogenic proteins, their misfolding, resulting structures, and interactions with membranes, which lead to membrane damage and subsequent cell death. Many of these proteins are implicated in serious illnesses such as Alzheimer’s disease and Parkinson’s disease. Misfolding of amyloidogenic proteins leads to the formation of polymorphic oligomers and fibrils. Oligomeric aggregates are widely thought to be the toxic species, however, fibrils also play a role in membrane damage. We focus on the structure of these aggregates and their interactions with model membranes. Study of interactions of amlyoidogenic proteins with model and natural membranes has shown the importance of the lipid bilayer in protein misfolding and aggregation and has led to the development of several models for membrane permeabilization by the resulting amyloid aggregates. We discuss several of these models: formation of structured pores by misfolded amyloidogenic proteins, extraction of lipids, interactions with receptors in biological membranes, and membrane destabilization by amyloid aggregates perhaps analogous to that caused by antimicrobial peptides.

Delicate balance between functionally required flexibility and aggregation risk in a β-rich protein

Susceptibility to aggregation is general to proteins because of the potential for intermolecular interactions between hydrophobic stretches in their amino acid sequences. Protein aggregation has been implicated in several catastrophic diseases, yet we still lack in-depth understanding about how proteins are channeled to this state. Using a predominantly β-sheet protein whose folding has been explored in detail, cellular retinoic acid-binding protein 1 (CRABP1), as a model, we have tackled the challenge of understanding the links between a protein's natural tendency to fold, 'breathe', and function with its propensity to misfold and aggregate. We identified near-native dynamic species that lead to aggregation and found that inherent structural fluctuations in the native protein, resulting in opening of the ligand-entry portal, expose hydrophobic residues on the most vulnerable aggregation-prone sequences in CRABP1. CRABP1 and related intracellullar lipid-binding proteins have not been reported to aggregate inside cells, and we speculate that the cellular concentration of their open, aggregation-prone conformations is sufficient for ligand binding but below the critical concentration for aggregation. Our finding provides an example of how nature fine-tunes a delicate balance between protein function, conformational variability, and aggregation vulnerability and implies that with the evolutionary requirement for proteins to fold and function, aggregation becomes an unavoidable but controllable risk.

Substrate dynamics in enzyme action: rotations of monosaccharide subunits in the binding groove are essential for pectin methylesterase processivity

The dynamical behavior of biomacromolecules is a fundamental property regulating a large number of biological processes. Protein dynamics have been widely shown to play a role in enzyme catalysis; however, the interplay between substrate dynamics and enzymatic activity is less understood. We report insights into the role of dynamics of substrates in the enzymatic activity of PME from Erwinia chrysanthemi, a processive enzyme that catalyzes the hydrolysis of methylester groups from the galacturonic acid residues of homogalacturonan chains, the major component of pectin. Extensive molecular dynamics simulations of this PME in complex with decameric homogalacturonan chains possessing different degrees and patterns of methylesterification show how the carbohydrate substitution pattern governs the dynamics of the substrate in the enzyme's binding cleft, such that substrate dynamics represent a key prerequisite for the PME biological activity. The analyses reveal that correlated rotations around glycosidic bonds of monosaccharide subunits at and immediately adjacent to the active site are a necessary step to ensure substrate processing. Moreover, only substrates with the optimal methylesterification pattern attain the correct dynamical behavior to facilitate processive catalysis. This investigation is one of the few reported examples of a process where the dynamics of a substrate are vitally important.

Bacterial cell division regulation by Ser/Thr kinases: a structural perspective

Recent genetic, biochemical and structural studies have established that eukaryotic-like Ser/Thr protein-kinases are critical mediators of developmental changes and host pathogen interactions in bacteria. Although with lower abundance compared to their homologues from eukaryotes, Ser/Thr protein-kinases are widespread in gram-positive bacteria. These data underline a key role of reversible Ser/Thr phosphorylation in bacterial physiology and virulence. Numerous studies have revealed how phosphorylation/dephosphorylation of Ser/Thr protein-kinases governs cell division and cell wall biosynthesis and that Ser/Thr protein kinases are responsible for distinct phenotypes, dependent on different environmental signals. In this review we discuss the current understandings of Ser/Thr protein-kinases functional processes based on structural data.

The amyloid-cell membrane system. The interplay between the biophysical features of oligomers/fibrils and cell membrane defines amyloid toxicity

Amyloid cytotoxicity, structure and polymorphisms are themes of increasing importance. Present knowledge considers any peptide/protein able to undergo misfolding and aggregation generating intrinsically cytotoxic amyloids. It also describes growth and structure of amyloid fibrils and their possible disassembly, whereas reduced information is available on oligomer structure. Recent research has highlighted the importance of the environmental conditions as determinants of the amyloid polymorphisms and cytotoxicity. Another body of evidence describes chemical or biological surfaces as key sites of protein misfolding and aggregation or of interaction with amyloids and the resulting biochemical modifications inducing cell functional/viability impairment. In particular, the membrane lipid composition appears to modulate cell response to toxic amyloids, thus contributing to explain the variable vulnerability to the same amyloids of different cell types. Finally, a recent view describes amyloid toxicity as an emerging property dependent on a complex interplay between the biophysical features of early aggregates and the interacting cell membranes taken as a whole system.