The fundamentals of protein folding: bringing together theory and experiment

Energy landscapes for proteins described by the UNRES coarse-grained potential

The self-assembly of proteins is encoded in the underlying potential energy surface (PES), from which we can predict structure, dynamics, and thermodynamic properties. However, the corresponding analysis becomes increasingly challenging with larger protein sizes, due to the computational time required, which grows significantly with the number of atoms. Coarse-grained models offer an attractive approach to reduce the computational cost. In this Feature Article, we describe our implementation of the UNited RESidue (UNRES) coarse-grained potential in the Cambridge energy landscapes software. We have applied this framework to explore the energy landscapes of four proteins that exhibit native states involving different secondary structures. Here we have tested the ability of the UNRES potential to represent the global energy landscape of proteins containing up to 100 amino acid residues. The resulting potential energy landscapes exhibit good agreement with experiment, with low-lying minima close to the PDB geometries and to results obtained using the all-atom AMBER force field. The new program interfaces will allow us to investigate larger biomolecules in future work, using the UNRES potential in combination with all the methodology available in the computational energy landscapes framework.

Modern Alchemical Free Energy Methods for Drug Discovery Explained

This Perspective provides a contextual explanation of the current state-of-the-art alchemical free energy methods and their role in drug discovery as well as highlights select emerging technologies. The narrative attempts to answer basic questions about what goes on "under the hood" in free energy simulations and provide general guidelines for how to run simulations and analyze the results. It is the hope that this work will provide a valuable introduction to students and scientists in the field.

Side Chain Geometry Determines the Fibrillation Propensity of a Minimal Two-Beads-per-Residue Peptide Model

The molecular mechanism of fibrillation is an important issue for understanding peptide aggregation. In our previous work, we demonstrated that the interchain attraction and intrachain bending stiffness control the aggregation kinetics and transient aggregate morphologies of a one-bead-per-residue implicit solvent peptide model. However, that model did not lead to fibrillation. In this work, we study the molecular origin of fibril formation using a two-beads-per-residue model, where one bead represents the backbone residue atoms and the other the side chain atoms. We show that the side chain geometry determines the fibrillation propensity that is further modulated by the modified terminal beads. This allows us to bring out the effects of side chain geometry and terminal capping on the fibrillation propensity. Our model does not assume a secondary structure and is, perhaps, the simplest bead-based chain model leading to fibrillation.

Antibody Interfaces Revealed Through Structural Mining

Antibodies are fundamental effectors of humoral immunity, and have become a highly successful class of therapeutics. There is increasing evidence that antibodies utilize transient homotypic interactions to enhance function, and elucidation of such interactions can provide insights into their biology and new opportunities for their optimization as drugs. Yet the transitory nature of weak interactions makes them difficult to investigate. Capitalizing on their rich structural data and high conservation, we have characterized all the ways that antibody fragment antigen-binding (Fab) regions interact crystallographically. This approach led to the discovery of previously unrealized interfaces between antibodies. While diverse interactions exist, β-sheet dimers and variable-constant elbow dimers are recurrent motifs. Disulfide engineering enabled interactions to be trapped and investigated structurally and functionally, providing experimental validation of the interfaces and illustrating their potential for optimization. This work provides first insight into previously undiscovered oligomeric interactions between antibodies, and enables new opportunities for their biotherapeutic optimization.

Insights into Fluctuations of Structure of Proteins: Significance of Intermediary States in Regulating Biological Functions

Proteins are indispensable to cellular communication and metabolism. The structure on which cells and tissues are developed is deciphered from proteins. To perform functions, proteins fold into a three-dimensional structural design, which is specific and fundamentally determined by their characteristic sequence of amino acids. Few of them have structural versatility, allowing them to adapt their shape to the task at hand. The intermediate states appear momentarily, while protein folds from denatured (D) ⇔ native (N), which plays significant roles in cellular functions. Prolific effort needs to be taken in characterizing these intermediate species if detected during the folding process. Protein folds into its native structure through definite pathways, which involve a limited number of transitory intermediates. Intermediates may be essential in protein folding pathways and assembly in some cases, as well as misfolding and aggregation folding pathways. These intermediate states help to understand the machinery of proper folding in proteins. In this review article, we highlight the various intermediate states observed and characterized so far under in vitro conditions. Moreover, the role and significance of intermediates in regulating the biological function of cells are discussed clearly.

Structural Lessons From the Mutant Proinsulin Syndrome

Insight into folding mechanisms of proinsulin has been provided by analysis of dominant diabetes-associated mutations in the human insulin gene (INS). Such mutations cause pancreatic β-cell dysfunction due to toxic misfolding of a mutant proinsulin and impairment in trans of wild-type insulin secretion. Anticipated by the "Akita" mouse (a classical model of monogenic diabetes mellitus; DM), this syndrome illustrates the paradigm endoreticulum (ER) stress leading to intracellular proteotoxicity. Diverse clinical mutations directly or indirectly perturb native disulfide pairing leading to protein misfolding and aberrant aggregation. Although most introduce or remove a cysteine (Cys; leading in either case to an unpaired thiol group), non-Cys-related mutations identify key determinants of folding efficiency. Studies of such mutations suggest that the hormone's evolution has been constrained not only by structure-function relationships, but also by the susceptibility of its single-chain precursor to impaired foldability. An intriguing hypothesis posits that INS overexpression in response to peripheral insulin resistance likewise leads to chronic ER stress and β-cell dysfunction in the natural history of non-syndromic Type 2 DM. Cryptic contributions of conserved residues to folding efficiency, as uncovered by rare genetic variants, define molecular links between biophysical principles and the emerging paradigm of Darwinian medicine: Biosynthesis of proinsulin at the edge of non-foldability provides a key determinant of "diabesity" as a pandemic disease of civilization.

Monitoring the surface tension by the pendant drop technique for detection of insulin fibrillogenesis

Monitoring the aggregation of amyloid-prone proteins is critical for understanding the mechanism of amyloid fibril formation. Insulin, when dissolved in low pH buffer, has a surface tension of 61-64 mN m^-1, as measured by the pendant drop technique. Formation of insulin amyloid fibrils resulted in the increase of the surface tension values up to 71.2-73.5 mN m^-1. The kinetics of fibril formation and fibril morphology were validated by ThT fluorescence and AFM, respectively. The results demonstrate that monitoring the surface tension by the pendant drop technique is a valuable tool for the detection of insulin amyloid aggregation.

Kinetics theories to understand the mechanism of aggregation of a protein and to design strategies for its inhibition

Protein aggregation phenomenon is closely related to the formation of amyloids which results in many neurodegenerative diseases like Alzheimer's, Parkinson's, Huntington's, and Amyotrophic Lateral Sclerosis. In order to prevent and treat these diseases, a clear understanding of the mechanism of misfolding and self-assembly of peptides and proteins is very crucial. The aggregation of a protein may involve various microscopic events. Multiple simulations utilizing the solutions of the master equation have given a better understanding of the kinetic profiles involved in the presence and absence of a particular microscopic event. This review focuses on understanding the contribution of these molecular events to protein aggregation based on the analysis of kinetic profiles of aggregation. We also discuss the effect of inhibitors, which target various species of aggregation pathways, on the kinetic profile of protein aggregation. At the end of this review, some strategies for the inhibition of aggregation that can be utilized by combining the chemical kinetics approach with thermodynamics are proposed.

Understanding molecular mechanisms of biologics drug delivery and stability from NMR spectroscopy

Protein therapeutics carry inherent limitations of membrane impermeability and structural instability, despite their predominant role in the modern pharmaceutical market. Effective formulations are needed to overcome physiological and physicochemical barriers, respectively, for improving bioavailability and stability. Knowledge of membrane affinity, cellular internalization, encapsulation, and release of drug-loaded carrier vehicles uncover the structural basis for designing and optimizing biopharmaceuticals with enhanced delivery efficiency and therapeutic efficacy. Understanding stabilizing and destabilizing interactions between protein drugs and formulation excipients provide fundamental mechanisms for ensuring the stability and quality of biological products. This article reviews the molecular studies of biologics using solution and solid-state NMR spectroscopy on structural attributes pivotal to drug delivery and stability. In-depth investigation of the structure-function relationship of drug delivery systems based on cell-penetrating peptides, lipid nanoparticles and polymeric colloidal, and biophysical and biochemical stability of peptide, protein, monoclonal antibody, and vaccine, as the integrative efforts on drug product design, will be elaborated.

A Review on Computational Approaches for Analyzing Hydrogen- Deuterium (H/D) Exchange of Proteins

Native state Hydrogen-Deuterium (H/D) exchange method has been used to study the structures and the unfolding pathways for quite a number of proteins. The H/D exchange method is generally monitored using nuclear magnetic resonance (NMR) spectroscopy and mass spectrometry (MS) techniques. NMR-assisted H/D exchange methods primarily monitor the residue level fluctuation of proteins, whereas MS-assisted H/D exchange methods analyze multifold ensemble conformations of proteins. In this connection, quite a large number of computational tools and algorithms have been developed for processing and analyzing huge amount of the H/D exchange data generated from these techniques. In this review, most of the freely available computational tools associated with the H/D exchange of proteins have been comprehensively reviewed and scopes to improve/ develop novel computational approaches for analyzing the H/D exchange data of proteins have also been brought into fore.

Cytosolic protein quality control machinery: Interactions of Hsp70 with a network of co-chaperones and substrates

The chaperone heat shock protein 70 (Hsp70) and its network of co-chaperones serve as a central hub of cellular protein quality control mechanisms. Domain organization in Hsp70 dictates ATPase activity, ATP dependent allosteric regulation, client/substrate binding and release, and interactions with co-chaperones. The protein quality control activities of Hsp70 are classified as foldase, holdase, and disaggregase activities. Co-chaperones directly assisting protein refolding included J domain proteins and nucleotide exchange factors. However, co-chaperones can also be grouped and explored based on which domain of Hsp70 they interact. Here we discuss how the network of cytosolic co-chaperones for Hsp70 contributes to the functions of Hsp70 while closely looking at their structural features. Comparison of domain organization and the structures of co-chaperones enables greater understanding of the interactions, mechanisms of action, and roles played in protein quality control.

Diabetes mellitus due to toxic misfolding of proinsulin variants

BACKGROUND

Dominant mutations in the human insulin gene (INS) lead to pancreatic β-cell dysfunction and diabetes mellitus (DM) due to toxic misfolding of a mutant proinsulin. Analogous to a classical mouse model of monogenic DM ("Akita"), this syndrome highlights the susceptibility of β-cells to endoreticulum (ER) stress due to protein misfolding and aberrant aggregation.

SCOPE OF REVIEW

Diverse clinical mutations directly or indirectly perturb native disulfide pairing. Whereas most introduce or remove a cysteine (Cys; leading in either case to an unpaired thiol group), non-Cys-related mutations identify key determinants of folding efficiency. Studies of such mutations suggest that the hormone's evolution has been constrained not only by structure-function relationships but also by the susceptibility of its single-chain precursor to impaired foldability. An intriguing hypothesis posits that INS overexpression in response to peripheral insulin resistance likewise leads to chronic ER stress and β-cell dysfunction in the natural history of nonsyndromic Type 2 DM.

MAJOR CONCLUSIONS

Cryptic contributions of conserved residues to folding efficiency, as uncovered by rare genetic variants, define molecular links between biophysical principles and the emerging paradigm of Darwinian medicine: Biosynthesis of proinsulin at the edge of nonfoldability provides a key determinant of "diabesity" as a pandemic disease of civilization.

Cryptic intermediates and metastable states of proteins as predicted by OneG computational method

Understanding structural excursions of proteins under folding conditions is crucial to map energy landscapes of proteins. In the present study, OneG computational tool has been used for analyzing possible existence of cryptic intermediates and metastable states of 26 proteins for which three prerequisite inputs of the OneG such as atomic coordinates of proteins, free energy of unfolding (ΔG_U) and free energy of exchange (ΔG_HX) determined in the absence of denaturant were available during the course of the study. The veraciousness of the tool on predicting the partially folded states of the proteins has been comprehensively described using experimental data available for 15 of the 26 proteins. Meanwhile, possible existence of partially structured states in the folding pathways of 11 other proteins has also been delineated as predicted by the OneG. In addition to mapping the folding pathways of proteins, the salient merits of the tool on systematically addressing the discrepancy between the ΔG_U and the ΔG_HX of the proteins have also been dealt.Communicated by Ramaswamy H. Sarma.

Identification of the Primary Factors Determining theSpecificity of Human VKORC1 Recognition by Thioredoxin-Fold Proteins

Redox (reduction-oxidation) reactions control many important biological processes in all organisms, both prokaryotes and eukaryotes. This reaction is usually accomplished by canonical disulphide-based pathways involving a donor enzyme that reduces the oxidised cysteine residues of a target protein, resulting in the cleavage of its disulphide bonds. Focusing on human vitamin K epoxide reductase (hVKORC1) as a target and on four redoxins (protein disulphide isomerase (PDI), endoplasmic reticulum oxidoreductase (ERp18), thioredoxin-related transmembrane protein 1 (Tmx1) and thioredoxin-related transmembrane protein 4 (Tmx4)) as the most probable reducers of VKORC1, a comparative in-silico analysis that concentrates on the similarity and divergence of redoxins in their sequence, secondary and tertiary structure, dynamics, intraprotein interactions and composition of the surface exposed to the target is provided. Similarly, hVKORC1 is analysed in its native state, where two pairs of cysteine residues are covalently linked, forming two disulphide bridges, as a target for Trx-fold proteins. Such analysis is used to derive the putative recognition/binding sites on each isolated protein, and PDI is suggested as the most probable hVKORC1 partner. By probing the alternative orientation of PDI with respect to hVKORC1, the functionally related noncovalent complex formed by hVKORC1 and PDI was found, which is proposed to be a first precursor to probe thiol-disulphide exchange reactions between PDI and hVKORC1.

Single-particle virology

The development of advanced experimental methodologies, such as optical tweezers, scanning-probe and super-resolved optical microscopies, has led to the evolution of single-molecule biophysics, a field of science that allows direct access to the mechanistic detail of biomolecular structure and function. The extension of single-molecule methods to the investigation of particles such as viruses permits unprecedented insights into the behavior of supramolecular assemblies. Here we address the scope of viral exploration at the level of individual particles. In an era of increased awareness towards virology, single-particle approaches are expected to facilitate the in-depth understanding, and hence combating, of viral diseases.

Improving double-ended transition state searches for soft-matter systems

Transitions between different stable configurations of biomolecules are important in understanding disease mechanisms, structure-function relations, and novel molecular-scale engineering. The corresponding pathways can be characterized efficiently using geometry optimization schemes based on double-ended transition state searches. An interpolation is first constructed between the known states and then refined, yielding a band that contains transition state candidates. Here, we analyze an example where various interpolation schemes lead to bands with a single step transition, but the correct pathway actually proceeds via an intervening, low-energy minimum. We compare a number of different interpolation schemes for this problem. We systematically alter the number of discrete images in the interpolations and the spring constants used in the optimization and test two schemes for adjusting the spring constants and image distribution, resulting in a total of 2760 different connection attempts. Our results confirm that optimized bands are not necessarily a good description of the transition pathways in themselves, and further refinement to actually converge transition states and establish their connectivity is required. We see an improvement in the optimized bands if we employ the adjustment of spring constants with doubly-nudged elastic band and a smaller improvement from the image redistribution. The example we consider is representative of numerous cases we have encountered in a wide variety of molecular and condensed matter systems.

Functional Modules of the Proteostasis Network

Cells invest in an extensive network of factors to maintain protein homeostasis (proteostasis) and prevent the accumulation of potentially toxic protein aggregates. This proteostasis network (PN) comprises the machineries for the biogenesis, folding, conformational maintenance, and degradation of proteins with molecular chaperones as central coordinators. Here, we review recent progress in understanding the modular architecture of the PN in mammalian cells and how it is modified during cell differentiation. We discuss the capacity and limitations of the PN in maintaining proteome integrity in the face of proteotoxic stresses, such as aggregate formation in neurodegenerative diseases. Finally, we outline various pharmacological interventions to ameliorate proteostasis imbalance.

Physical PEGylation to Prevent Insulin Fibrillation

Insulin is one of the most marketed therapeutic proteins worldwide. However, its formulation suffers from fibrillation, which affects the long-term storage limiting the development of novel devices for sustained delivery including portable infusion devices. We have investigated the effect of physical PEGylation on structural and colloidal stability of insulin by using 2 PEGylating agents terminating with polycyclic hydrophobic moieties, cholane and cholesterol: mPEG_5kDa-cholane and mPEG_5kDa-cholesterol, respectively. Microcalorimetric analyses showed that mPEG_5kDa-cholane and mPEG_5kDa-cholesterol efficiently bind insulin with binding constants (Ka) of 3.98 10⁴ and 1.14 10⁵ M^-1, respectively. At room temperature, the 2 PEGylating agents yielded comparable structural stabilization of α-helix conformation and decreased dimerization of insulin. However, melting studies showed that mPEG_5kDa-cholesterol has superior stabilizing effect of the protein conformation than mPEG_5kDa-cholane. Furthermore, the fibrillation study showed that at a 1:1 and 1:5 insulin/polymer molar ratios, mPEG_5kDa-cholesterol delays insulin fibrillation 40% and 26% more efficiently, respectively, as compared to mPEG_5kDa-cholane which was confirmed by transmission electron microscopy imaging. Insulin was released from the mPEG_5kDa-cholane and mPEG_5kDa-cholesterol assemblies with comparable kinetic profiles. The physical PEGylation has a beneficial effect on the stabilization and shielding of the insulin structure into the monomeric form, which is not prone to fibrillation and aggregation.

Chaperonin facilitates protein folding by avoiding initial polypeptide collapse

Chaperonins assist folding of many cellular proteins, including essential proteins for cell viability. However, it remains unclear how chaperonin-assisted folding is different from spontaneous folding. Chaperonin GroEL/GroES facilitates folding of denatured protein encapsulated in its central cage but the denatured protein often escapes from the cage to the outside during reaction. Here, we show evidence that the in-cage-folding and the escape occur diverging from the same intermediate complex in which polypeptide is tethered loosely to the cage and partly protrudes out of the cage. Furthermore, denatured proteins in the chaperonin cage are kept in more extended conformation than those initially formed in spontaneous folding. We propose that the formation of tethered intermediate of polypeptide is necessary to prevent polypeptide collapse at the expense of polypeptide escape. The tethering of polypeptide would allow freely mobile portions of tethered polypeptide to fold segmentally.

Dual effect of non-ionic detergent Triton X-100 on insulin amyloid formation

Atomic force microscopy, Thioflavin T (ThT) fluorescence assay, circular dichroism spectroscopy, differential scanning calorimetry, and molecular modeling techniques have been employed to investigate the amyloid aggregation of insulin in the presence of non-ionic detergent, Triton X-100 (TX-100). In contrast to recently described inhibition of lysozyme amyloid formation by non-ionic detergents (Siposova, 2017), the amyloid aggregation of insulin in the presence of sub-micellar TX-100 concentration exhibits two dissimilar phases. The first, inhibition phase, is observed at the protein to detergent molar ratio of 1:0.1 to 1:1. During this phase, the insulin amyloid fibril formation is inhibited by TX-100 up to ∼60%. The second, "morphological" phase, is observed at increasing detergent concentration, corresponding to protein:detergent molar ratio of ∼1:1 - 1:10. Under these conditions a significant increase of the steady-state ThT fluorescence intensities and a dramatically changed morphology of the insulin fibrils were observed. Increasing of the detergent concentration above the CMC led to complete inhibition of amyloidogenesis. Analysis of the experimental and molecular modeling results suggests an existence of up to six TX-100 binding sites within dimer of insulin with different binding energy. The physiological relevance of the results is discussed.

Energy Landscapes of Mini-Dumbbell DNA Octanucleotides

Single-stranded DNA structures play a significant role in biological systems, in particular during replication, translation, and DNA repair. Tracts of simple repetitive DNA are associated with slipped-strand mispairing, which can lead to genetic diseases. Recent NMR studies of TTTA and CCTG repeats have shown that these sequences form mini-dumbbells (MDBs), leading to frameshift mutations. Here we explore the energy landscapes of (CCTG)₂ and (TTTA)₂, which are currently the smallest known molecules to form MDBs. While (CCTG)₂ MDBs are stable, (TTTA)₂ exhibits numerous other structures with lower energies. A key factor identified in the stabilization of MDB structures is the bonding strength between residues 1 and 4, and 5 and 8.

Insight into human insulin aggregation revisited using NMR derived translational diffusion parameters

The NMR derived translational diffusion coefficients were performed on unlabeled and uniformly labeled ¹³C,¹⁵N human insulin in water, both in neat, with zinc ions only, and in pharmaceutical formulation, containing only m-cresol as phenolic ligand, glycerol and zinc ions. The results show the dominant role of the pH parameter and the concentration on aggregation. The diffusion coefficient Dav was used for monitoring the overall average state of oligomeric ensemble in solution. The analysis of the experimental data of diffusion measurements, using the direct exponential curve resolution algorithm (DECRA) allows suggesting the two main components of the oligomeric ensemble. The 3D HSQC-iDOSY, (diffusion ordered HSQC) experiments performed on ¹³C, ¹⁵N-fully labeled insulin at the two pH values, 4 and 7.5, allow for the first time a more detailed experimental observation of individual components in the ensemble. The discussion involves earlier static and dynamic laser light scattering experiments and recent NMR derived translational diffusion results. The results bring new informations concerning the preparation of pharmaceutical formulation and in particular a role of Zn²⁺ ions. They also will enable better understanding and unifying the results of studies on insulin misfolding effects performed in solution by diverse physicochemical methods at different pH and concentration.

The protofilament architecture of a de novo designed coiled coil-based amyloidogenic peptide

Amyloid fibrils are polymers formed by proteins under specific conditions and in many cases they are related to pathogenesis, such as Parkinson's and Alzheimer's diseases. Their hallmark is the presence of a β-sheet structure. High resolution structural data on these systems as well as information gathered from multiple complementary analytical techniques is needed, from both a fundamental and a pharmaceutical perspective. Here, a previously reported de novo designed, pH-switchable coiled coil-based peptide that undergoes structural transitions resulting in fibril formation under physiological conditions has been exhaustively characterized by transmission electron microscopy (TEM), cryo-TEM, atomic force microscopy (AFM), wide-angle X-ray scattering (WAXS) and solid-state NMR (ssNMR). Overall, a unique 2-dimensional carpet-like assembly composed of large coexisiting ribbon-like, tubular and funnel-like structures with a clearly resolved protofilament substructure is observed. Whereas electron microscopy and scattering data point somewhat more to a hairpin model of β-fibrils, ssNMR data obtained from samples with selectively labelled peptides are in agreement with both, hairpin structures and linear arrangements.

Chaperones convert the energy from ATP into the nonequilibrium stabilization of native proteins

During and after protein translation, molecular chaperones require ATP hydrolysis to favor the native folding of their substrates and, under stress, to avoid aggregation and revert misfolding. Why do some chaperones need ATP, and what are the consequences of the energy contributed by the ATPase cycle? Here, we used biochemical assays and physical modeling to show that the bacterial chaperones GroEL (Hsp60) and DnaK (Hsp70) both use part of the energy from ATP hydrolysis to restore the native state of their substrates, even under denaturing conditions in which the native state is thermodynamically unstable. Consistently with thermodynamics, upon exhaustion of ATP, the metastable native chaperone products spontaneously revert to their equilibrium non-native states. In the presence of ATPase chaperones, some proteins may thus behave as open ATP-driven, nonequilibrium systems whose fate is only partially determined by equilibrium thermodynamics.

Finding the needle in the haystack: towards solving the protein-folding problem computationally

Prediction of protein tertiary structures from amino acid sequence and understanding the mechanisms of how proteins fold, collectively known as "the protein folding problem," has been a grand challenge in molecular biology for over half a century. Theories have been developed that provide us with an unprecedented understanding of protein folding mechanisms. However, computational simulation of protein folding is still difficult, and prediction of protein tertiary structure from amino acid sequence is an unsolved problem. Progress toward a satisfying solution has been slow due to challenges in sampling the vast conformational space and deriving sufficiently accurate energy functions. Nevertheless, several techniques and algorithms have been adopted to overcome these challenges, and the last two decades have seen exciting advances in enhanced sampling algorithms, computational power and tertiary structure prediction methodologies. This review aims at summarizing these computational techniques, specifically conformational sampling algorithms and energy approximations that have been frequently used to study protein-folding mechanisms or to de novo predict protein tertiary structures. We hope that this review can serve as an overview on how the protein-folding problem can be studied computationally and, in cases where experimental approaches are prohibitive, help the researcher choose the most relevant computational approach for the problem at hand. We conclude with a summary of current challenges faced and an outlook on potential future directions.

Solvent-Responsive Molecularly Imprinted Nanogels for Targeted Protein Analysis in MALDI-TOF Mass Spectrometry

Molecular imprinted poly(acrylamido)-derivative nanogels have shown their selectivity to bind the protein human serum transferrin (HTR) and also showed their capability for instantaneous solvent-induced modification upon the addition of acetonitrile. Integrated to matrix-assisted laser desorption/ionization time-of-flight mass analysis the HTR-imprinted solvent-responsive nanogels permitted the determination of HTR straight from serum and offered novel perspectives in targeted protein analysis.

Thermodynamics of Hydrophobic Amino Acids in Solution: A Combined Experimental-Computational Study

We present a joint experimental-computational study to quantitatively describe the thermodynamics of hydrophobic leucine amino acids in aqueous solution. X-ray scattering data were acquired at a series of solute and salt concentrations to effectively measure interleucine interactions, indicating that a major scattering peak is observed consistently at q = 0.83 Å^-1. Atomistic molecular dynamics simulations were then performed and compared with the scattering data, achieving high consistency at both small and wider scattering angles (q = 0-1.5 Å^-1). This experimental-computational consistence enables a first glimpse of the leucine-leucine interacting landscape, where two leucine molecules are aligned mostly in a parallel fashion, as opposed to antiparallel, but also allows us to derive effective leucine-leucine interactions in solution. Collectively, this combined approach of employing experimental scattering and molecular simulation enables quantitative characterization of effective intermolecular interactions of hydrophobic amino acids, critical for protein function and dynamics such as protein folding.

Protein Folding Structures: Formation of Folding Structures Based on Probability Theory

To the best of our knowledge, this is the first study that shows that the X-ray structures of proteins can be dissected into their continuous folding structure units. Each folding structure unit was designed such that both the terminal di- or tri-peptide sequences shared common sequences with the two adjacent folding structure units. To encode the folding structure information of proteins into their amino acid sequences, we proposed 44 kinds of folding elements, which covered all of the amino acids in the protein chains, and defined all folding structure units. The folding element was defined to mean a minimum structural piece, which covered the frame of the main chain of each amino acid in a protein chain. A folding structure unit of a local sequence could be fully characterized by the sequential combination of individual folding elements assigned to each amino acid. The folding structure information showed amino acid preferences in various positions in folding structure units. Folding structure formation proceeded on the basis of probability theory. Strikingly, relative formation ability analysis clearly indicated that we can decode the types and the chain length of folding structure units from the amino acid sequence of a protein.

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.

Under-folded proteins: Conformational ensembles and their roles in protein folding, function, and pathogenesis

For decades, protein function was intimately linked to the presence of a unique, aperiodic crystal‐like structure in a functional protein. The two only places for conformational ensembles of under‐folded (or partially folded) protein forms in this picture were either the end points of the protein denaturation processes or transiently populated folding intermediates. Recent years witnessed dramatic change in this perception and conformational ensembles, which the under‐folded proteins are, have moved from the shadow. Accumulated to date data suggest that a protein can exist in at least three global forms–functional and folded, functional and intrinsically disordered (nonfolded), and nonfunctional and misfolded/aggregated. Under‐folded protein states are crucial for each of these forms, serving as important folding intermediates of ordered proteins, or as functional states of intrinsically disordered proteins (IDPs) and IDP regions (IDPRs), or as pathology triggers of misfolded proteins. Based on these observations, conformational ensembles of under‐folded proteins can be classified as transient (folding and misfolding intermediates) and permanent (IDPs and stable misfolded proteins). Permanently under‐folded proteins can further be split into intentionally designed (IDPs and IDPRs) and unintentionally designed (misfolded proteins). Although intrinsic flexibility, dynamics, and pliability are crucial for all under‐folded proteins, the different categories of under‐foldedness are differently encoded in protein amino acid sequences. © 2013 Wiley Periodicals, Inc. Biopolymers 99: 870–887, 2013.

Ligand and interfacial dynamics in a homodimeric hemoglobin

The structural dynamics of dimeric hemoglobin (HbI) from Scapharca inaequivalvis in different ligand-binding states is studied from atomistic simulations on the μs time scale. The intermediates are between the fully ligand-bound (R) and ligand-free (T) states. Tertiary structural changes, such as rotation of the side chain of Phe97, breaking of the Lys96-heme salt bridge, and the Fe-Fe separation, are characterized and the water dynamics along the R-T transition is analyzed. All these properties for the intermediates are bracketed by those determined experimentally for the fully ligand-bound and ligand-free proteins, respectively. The dynamics of the two monomers is asymmetric on the 100 ns timescale. Several spontaneous rotations of the Phe97 side chain are observed which suggest a typical time scale of 50-100 ns for this process. Ligand migration pathways include regions between the B/G and C/G helices and, if observed, take place in the 100 ns time scale.

Human, vector and parasite Hsp90 proteins: A comparative bioinformatics analysis

•

Cytosolic Hsp90s are more conserved than those from mitochondrial and ER.

•

Cell environment plays a role in the overall physicochemical properties of Hsp90s.

•

Serine and tyrosine are favored phosphorylated residues of Hsp90s.

•

Mitochondrial and ER Hsp90s have motifs unique to specific organisms.

The treatment of protozoan parasitic diseases is challenging, and thus identification and analysis of new drug targets is important. Parasites survive within host organisms, and some need intermediate hosts to complete their life cycle. Changing host environment puts stress on parasites, and often adaptation is accompanied by the expression of large amounts of heat shock proteins (Hsps). Among Hsps, Hsp90 proteins play an important role in stress environments. Yet, there has been little computational research on Hsp90 proteins to analyze them comparatively as potential parasitic drug targets. Here, an attempt was made to gain detailed insights into the differences between host, vector and parasitic Hsp90 proteins by large-scale bioinformatics analysis. A total of 104 Hsp90 sequences were divided into three groups based on their cellular localizations; namely cytosolic, mitochondrial and endoplasmic reticulum (ER). Further, the parasitic proteins were divided according to the type of parasite (protozoa, helminth and ectoparasite). Primary sequence analysis, phylogenetic tree calculations, motif analysis and physicochemical properties of Hsp90 proteins suggested that despite the overall structural conservation of these proteins, parasitic Hsp90 proteins have unique features which differentiate them from human ones, thus encouraging the idea that protozoan Hsp90 proteins should be further analyzed as potential drug targets.

How do chaperonins fold protein?

Protein folding is a biological process that is essential for the proper functioning of proteins in all living organisms. In cells, many proteins require the assistance of molecular chaperones for their folding. Chaperonins belong to a class of molecular chaperones that have been extensively studied. However, the mechanism by which a chaperonin mediates the folding of proteins is still controversial. Denatured proteins are folded in the closed chaperonin cage, leading to the assumption that denatured proteins are completely encapsulated inside the chaperonin cage. In contrast to the assumption, we recently found that denatured protein interacts with hydrophobic residues at the subunit interfaces of the chaperonin, and partially protrude out of the cage. In this review, we will explain our recent results and introduce our model for the mechanism by which chaperonins accelerate protein folding, in view of recent findings.

Probing conformational stability and dynamics of erythroid and nonerythroid spectrin: effects of urea and guanidine hydrochloride

We have studied the conformational stability of the two homologous membrane skeletal proteins, the erythroid and non-erythroid spectrins, in their dimeric and tetrameric forms respectively during unfolding in the presence of urea and guanidine hydrochloride (GuHCl). Fluorescence and circular dichroism (CD) spectroscopy have been used to study the changes of intrinsic tryptophan fluorescence, anisotropy, far UV-CD and extrinsic fluorescence of bound 1-anilinonapthalene-8-sulfonic acid (ANS). Chemical unfolding of both proteins were reversible and could be described as a two state transition. The folded erythroid spectrin and non-erythroid spectrin were directly converted to unfolded monomer without formation of any intermediate. Fluorescence quenching, anisotropy, ANS binding and dynamic light scattering data suggest that in presence of low concentrations of the denaturants (up-to 1M) hydrogen bonding network and van der Waals interaction play a role inducing changes in quaternary as well as tertiary structures without complete dissociation of the subunits. This is the first report of two large worm like, multi-domain proteins obeying twofold rule which is commonly found in small globular proteins. The free energy of stabilization (ΔGuH20) for the dimeric spectrin has been 20 kcal/mol lesser than the tetrameric from.

Technical decision making with higher order structure data: impact of a formulation change on the higher order structure and stability of a mAb

Changes in formulation may be required during the development of protein therapeutics. Some of the changes may alter the protein higher order structure (HOS). In this note, we show how the change from a trehalose-based formulation to an arginine-based formulation concomitantly impacted the tertiary structure and the thermal stability of a mAb (mAb1). The secondary structure was not disrupted by the formulation change. The destabilization of the tertiary structure did not affect the long-term stability or the bioactivity of mAb1. This indicates that loss of conformational stability was likely compensated by improvements in the colloidal stability of mAb1 in the arginine-based formulation. The formulation-induced changes in HOS were reversible as proven by measurements after dilution in a common buffer (phosphate-buffered saline). For aggregation driven by assembly of aggregates (colloidally limited), small changes in conformational structure and stability as measured by HOS methods may not necessarily be predictive of long-term stability.

Quantifying the Sources of Kinetic Frustration in Folding Simulations of Small Proteins

Experiments and atomistic simulations of polypeptides have revealed structural intermediates that promote or inhibit conformational transitions to the native state during folding. We invoke a concept of “kinetic frustration” to quantify the prevalence and impact of these behaviors on folding rates within a large set of atomistic simulation data for 10 fast-folding proteins, where each protein’s conformational space is represented as a Markov state model of conformational transitions. Our graph theoretic approach addresses what conformational features correlate with folding inhibition and therefore permits comparison among features within a single protein network and also more generally between proteins. Nonnative contacts and nonnative secondary structure formation can thus be quantitatively implicated in inhibiting folding for several of the tested peptides.

REMD and umbrella sampling simulations to probe the energy barrier of the folding pathways of engrailed homeodomain

Proteins fold by diverse pathways which depend on the energy barriers involved in reaching different intermediates. There has been a lot of development in the theoretical aspects of protein folding, from force-field to simulation techniques. One such simulation approach is replica exchange molecular dynamics simulation (REMD), which provides an efficient conformational sampling method to understand the events involved in protein folding. In this study, an attempt is made to explore the folding funnel of engrailed homeodomain protein (EnHD) using REMD simulations. EnHD is a 54 residue long helix bundle protein which has a folding time of about 15 μs. The protein was represented using the Amber United atom model in order to reduce the system size which helped to speed up the simulation. Individual replicas were simulated for 1.4-2 μs making cumulative time of more than 100 μs of REMD simulations. Free energy analysis was carried out to understand the folding behavior of EnHD protein. Effects of temperature range and exchange frequency in REMD simulations have been explored. In addition to this, multiple umbrella sampling (US) simulations of a total of 320 ns were also carried out, followed by weighted histogram analysis method (WHAM) to investigate the energy barriers involved during the folding of various intermediates. US studies were also carried on mutational variants of EnHD protein to see effect of the mutations on the folding pathway of the protein. The use of US technique may be helpful for predicting fast folding mutants or protein engineering. The combination of REMD with US may help in understanding the energetics between multiple pathways of fast folding proteins and their mutant counterparts.

Atomic force microscopy to study molecular mechanisms of amyloid fibril formation and toxicity in Alzheimer's disease

Alzheimer's disease (AD) is a devastating neurodegenerative disease characterized by dementia and memory loss for which no cure or effective prevention is currently available. Neurodegeneration in AD is linked to formation of amyloid plaques found in brain tissues of Alzheimer's patients during post-mortem examination. Amyloid plaques are composed of amyloid fibrils and small oligomers - insoluble protein aggregates. Although amyloid plaques are found on the neuronal cell surfaces, the mechanism of amyloid toxicity is still not well understood. Currently, it is believed that the cytotoxicity is a result of the nonspecific interaction of small soluble amyloid oligomers (rather than longer fibrils) with the plasma membrane. In recent years, nanotechnology has contributed significantly to understanding the structure and function of lipid membranes and to the study of the molecular mechanisms of membrane-associated diseases. We review the current state of research, including applications of the latest nanotechnology approaches, on the interaction of lipid membranes with the amyloid-β (Aβ) peptide in relation to amyloid toxicity. We discuss the interactions of Aβ with model lipid membranes with a focus to demonstrate that composition, charge and phase of the lipid membrane, as well as lipid domains and rafts, affect the binding of Aβ to the membrane and contribute to toxicity. Understanding the role of the lipid membrane in AD at the nanoscale and molecular level will contribute to the understanding of the molecular mechanism of amyloid toxicity and may aid into the development of novel preventive strategies to combat AD.

Dissecting the contributions of β-hairpin tyrosine pairs to the folding and stability of long-lived human γD-crystallins

Ultraviolet-radiation-induced damage to and aggregation of human lens crystallin proteins are thought to be a significant pathway to age-related cataract. The aromatic residues within the duplicated Greek key domains of γ- and β-crystallins are the main ultraviolet absorbers and are susceptible to direct and indirect ultraviolet damage. The previous site-directed mutagenesis studies have revealed a striking difference for two highly conserved homologous β-hairpin Tyr pairs, at the N-terminal domain (N-td) and C-terminal domain (C-td), respectively, in their contribution to the overall stability of HγD-Crys, but why they behave so differently still remains a mystery. In this paper, we systematically investigated the underlying molecular mechanism and detailed contributions of these two Tyr pairs with large scale molecular dynamics simulations. A series of different tyrosine-to-alanine pair(s) substitutions were performed in either the N-td, the C-td, or both. Our results suggest that the Y45A/Y50A pair substitution in the N-td mainly affects the stability of the N-td itself, while the Y133A/Y138A pair substitution in the C-td leads to a more cooperative unfolding of both N-td and C-td. The stability of motif 2 in the N-td is mainly determined by the interdomain interface, while motif 1 in the N-td or motifs 3 and 4 in the C-td are mainly stabilized by the intradomain hydrophobic core. The damage to any tyrosine pair(s) can directly introduce some apparent water leakage to the hydrophobic core at the interface, which in turn causes a serious loss in the stability of the N-td. However, for the C-td substitutions, it may further impair the stable "sandwich-like" Y133-R167-Y138 cluster (through cation-π interactions) in the wild-type, thus causing the loop regions near the residue A138 to undergo large fluctuations, which in turn results in the intrusion of water into the hydrophobic core of the C-td and induces the C-td to lose its stability. These findings help resolve the "mystery" on why these two Tyr pairs display such a striking difference in their contributions to the overall protein stability despite their highly homologous nature.

Modelling amyloid fibril formation kinetics: mechanisms of nucleation and growth

Amyloid and amyloid-like fibrils are self-assembling protein nanostructures, of interest for their robust material properties and inherent biological compatibility as well as their putative role in a number of debilitating mammalian disorders. Understanding fibril formation is essential to the development of strategies to control, manipulate or prevent fibril growth. As such, this area of research has attracted significant attention over the last half century. This review describes a number of different models that have been formulated to describe the kinetics of fibril assembly. We describe the macroscopic implications of mechanisms in which secondary processes such as secondary nucleation, fragmentation or branching dominate the assembly pathway, compared to mechanisms dominated by the influence of primary nucleation. We further describe how experimental data can be analysed with respect to the predictions of kinetic models.