Three Weaknesses for Three Perturbations: Comparing Protein Unfolding Under Shear, Force, and Thermal Stresses

Microfluidic Stress Device to Decouple the Synergistic Effect of Shear and Interfaces on Antibody Aggregation

Protein denaturation and aggregation resulting from the effects of interfacial stress, often enhanced by flow and shear stress, pose significant challenges in the production of therapeutic proteins and monoclonal antibodies. The influence of flow on protein stability is closely intertwined with interfacial effects. In this study, we have developed a microfluidic device capable of exposing low volume (< 320 µL) protein solutions to highly uniform shear. To disentangle the synergistic impact of flow and interfaces on protein aggregation, we fabricated two devices composed of different materials, namely poly(methyl methacrylate) (PMMA) and stainless steel. Upon application of shear, we observed formation of protein particles in the micron-size range. Notably, The number of particles generated in the steel devices was ∼ 3.5 fold lower than in the PMMA device, hinting at an interface-mediated effect. With increasing the protein concentration from 1 to 50 mg/mL we observed a saturation in the amount of aggregates, further confirming the key role of solid-liquid interfaces in inducing particle formation. Introduction of non-ionic surfactants prevented protein aggregation, even at the highest tested protein concentration and low surfactant concentrations of 0.05 mg/mL. Overall, our findings corroborate the synergistic impact of shear and interface effects on protein aggregation. The device developed in this study offers a small-scale platform for assessing the stability of antibody formulations throughout various stages of the development and manufacturing process.

Multiparameter screen optimizes immunoprecipitation

Immunoprecipitation (IP) coupled with mass spectrometry effectively maps protein-protein interactions when genome-wide, affinity-tagged cell collections are used. Such studies have recorded significant portions of the compositions of physiological protein complexes, providing draft 'interactomes'; yet many constituents of protein complexes still remain uncharted. This gap exists partly because high-throughput approaches cannot optimize each IP. A key challenge for IP optimization is stabilizing in vivo interactions during the transfer from cells to test tubes; failure to do so leads to the loss of genuine interactions during the IP and subsequent failure to detect. Our high-content screening method explores the relationship between in vitro chemical conditions and IP outcomes, enabling rapid empirical optimization of conditions for capturing target macromolecular assemblies.

Dynamics and Structures of Amyloid Aggregates under Fluid Flows

In this work, we investigate how fluid flows impact the aggregation mechanisms of Aβ40 proteins and Aβ16-22 peptides and mechanically perturb their (pre)fibrillar aggregates. We exploit the OPEP coarse-grained model for proteins and the Lattice Boltzmann Molecular Dynamics technique. We show that beyond a critical shear rate, amyloid aggregation speeds up in Couette flow because of the shorter collisions times between aggregates, following a transition from diffusion limited to advection dominated dynamics. We also characterize the mechanical deformation of (pre)fibrillar states due to the fluid flows (Couette and Poiseuille), confirming the capability of (pre)fibrils to form pathological loop-like structures as detected in experiments. Our findings can be of relevance for microfluidic applications and for understanding aggregation in the interstitial brain space.

Mechanical Stress Decreases the Amplitude of Twisting and Bending Fluctuations of Actin Filaments

A variety of biological roles of mechanical forces have been proposed in cell biology, such as cell signaling pathways for survival, development, growth, and differentiation. Mechanical forces alter the mechanical conditions within cells and their environment, which strongly influences the reorganization of the actin cytoskeleton. Single-molecule imaging studies of actin filaments have led to the hypothesis that the actin filament acts as a mechanosensor; e.g., increases in actin filament tension alter their conformation and affinity for regulatory proteins. However, our understanding of the molecular mechanisms underlying how tension modulates the mechanical behavior of a single actin filament is still incomplete. In this study, a direct measurement of the twisting and bending of a fluorescently labeled single actin filament under different tension levels by force application (0.8-3.4 pN) was performed using single-molecule fluorescence polarization (SMFP) microscopy. The results showed that the amplitude of twisting and bending fluctuations of a single actin filament decreased with increasing tension. Electron micrograph analysis of tensed filaments also revealed that the fluctuations in the crossover length of actin filaments decreased with increasing filament tension. Possible molecular mechanisms underlying these results involving the binding of actin-binding proteins, such as cofilin, to the filament are discussed.

Plant-based meat analogs and fat substitutes, structuring technology and protein digestion: A review

There is currently an increasing trend in the consumption of meat analogs and fat substitutes due to the health hazards by excessive consumption of meat. Simulating the texture and mouthfeel of meat through structured plant-derived polymers has become a popular processing method. In this review, the mechanical structuring technology of plant polymers for completely replacing real meat is mainly introduced in this review, which mainly focuses on the parameters and principles of mechanical equipment for the production of vegan meat. The difference in composition between plant meat and real meat is mainly reflected in the protein, and particular attention should be paid to the digestive characteristics of plant meat protein in the gastrointestinal tract. Therefore, the differences in the protein digestibility properties of meat analogs and real meat is discussed in this review, focusing primarily on protein digestibility and peptide/amino acid composition of mechanically structured vegan meats. In terms of fat substitutes for meat products, the types of plant polymer colloidal systems used for meat fat substitutes is comprehensively introduced, including emulsion, hydrogel and oleogel.

Self-Assembly of Amyloid-Beta (Aβ) Peptides from Solution to Near In Vivo Conditions

Understanding the atomistic resolution changes during the self-assembly of amyloid peptides or proteins is important to develop compounds or conditions to alter the aggregation pathways and suppress the toxicity and potentially aid in the development of drugs. However, the complexity of protein aggregation and the transient order/disorder of oligomers along the pathways to fibril are very challenging. In this Perspective, we discuss computational studies of amyloid polypeptides carried out under various conditions, including conditions closely mimicking in vivo and point out the challenges in obtaining physiologically relevant results, focusing mainly on the amyloid-beta Aβ peptides.

Aggregation of Lactoferrin Caused by Droplet Atomization Process via a Two-Fluid Nozzle: The Detrimental Effect of Air-Water Interfaces

Biological macromolecules, especially therapeutic proteins, are delicate and highly sensitive to denaturation from stresses encountered during the manufacture of dosage forms. Thin-film freeze-drying (TFFD) and spray freeze-drying (SFD) are two processes used to convert liquid forms of protein into dry powders. In the production of inhalable dry powders that contain proteins, these potential stressors fall into three categories based on their occurrence during the primary steps of the process: (1) droplet formation (e.g., the mechanism of droplet formation, including spray atomization), (2) freezing, and (3) frozen water removal (e.g., sublimation). This study compares the droplet formation mechanism used in TFFD and SFD by investigating the effects of spraying on the stability of proteins, using lactoferrin as a model. This study considers various perspectives on the denaturation (e.g., conformation) of lactoferrin after subjecting the protein solution to the atomization process using a pneumatic two-fluid nozzle (employed in SFD) or a low-shear drop application through the nozzle. The surface activity of lactoferrin was examined to explore the interfacial adsorption tendency, diffusion, and denaturation process. Subsequently, this study also investigates the secondary and tertiary structure of lactoferrin and the quantification of monomers, oligomers, and, ultimately, aggregates. The spraying process affected the tertiary structure more negatively than the tightly woven secondary structure, resulting in the peak position corresponding to the tryptophan (Trp) residues red-shifting by 1.5 nm. This conformational change can either (a) be reversed at low concentrations via relaxation or (b) proceed to form irreversible aggregates at higher concentrations. Interestingly, when the sample was allowed to progress into micrometer-sized aggregates, such a dramatic change was not detected using methods such as size-exclusion chromatography, polyacrylamide gel electrophoresis, and dynamic light scattering at 173°. A more complete understanding of the heterogeneous protein sample was achieved only through a combination of 173 and 13° backward and forward scattering, a combination of derived count rate measurements, and microflow imaging (MFI). After studying the impact of droplet formation mechanisms on aggregation tendency of lactoferrin, we further investigated two additional model proteins with different surface activity: bovine IgG (serving as a non surface-active negative reference), and β-galactosidase (another surface-active protein). The results corroborated the lactoferrin findings that spray-atomization-related stress-induced protein aggregation was much more pronounced for proteins that are surface active (lactoferrin and β-galactosidase), but it was minimal for non-surface-active protein (bovine IgG). Finally, compared to the low-shear dripping used in the TFFD process, lactoferrin underwent a relatively fast conformational change upon exposure to the high air-water interface of the two-fluid atomization nozzle used in the SFD process as compared to the low shear dripping used in the TFFD process. The interfacial-induced denaturation that occurred during spraying was governed primarily by the size of the atomized droplets, regardless of the duration of exposure to air. The percentage of denatured protein population and associated activity loss, in the case of β-galactosidase, was determined to range from 2 to 10% depending on the air-flow rate of the spraying process.

Molecular interpretation of single-molecule force spectroscopy experiments with computational approaches

Single molecule force-spectroscopy techniques have granted access to unprecedented molecular-scale details about biochemical and biological mechanisms. However, the interpretation of the experimental data is often challenging. Computational and simulation approaches (all-atom steered MD simulations in particular) are key to provide molecular details about the associated mechanisms, to help test different hypotheses and to predict experimental results. In this review, particular recent efforts directed towards the molecular interpretation of single-molecule force spectroscopy experiments on proteins and protein-related systems (often in close collaboration with experimental groups) will be presented. These results will be discussed in the broader context of the field, highlighting the recent achievements and the ongoing challenges for computational biophysicists and biochemists. In particular, I will focus on the input gained from molecular simulations approaches to rationalize the origin of the unfolded protein elasticity and the protein conformational behavior under force, to understand how force denaturation differs from chemical, thermal or shear unfolding, and to unravel the molecular details of unfolding events for a variety of systems. I will also discuss the use of models based on Langevin dynamics on a 1-D free-energy surface to understand the effect of protein segmentation on the work exerted by a force, or, at the other end of the spectrum of computational techniques, how quantum calculations can help to understand the reactivity of disulfide bridges exposed to force.

Assessing models of force-dependent unbinding rates via infrequent metadynamics

Protein–ligand interactions are crucial for a wide range of physiological processes. Many cellular functions result in these non-covalent “bonds” being mechanically strained, and this can be integral to proper cellular function. Broadly, two classes of force dependence have been observed—slip bonds, where the unbinding rate increases, and catch bonds, where the unbinding rate decreases. Despite much theoretical work, we cannot predict for which protein–ligand pairs, pulling coordinates, and forces a particular rate dependence will appear. Here, we assess the ability of MD simulations combined with enhanced sampling techniques to probe the force dependence of unbinding rates. We show that the infrequent metadynamics technique correctly produces both catch and slip bonding kinetics for model potentials. We then apply it to the well-studied case of a buckyball in a hydrophobic cavity, which appears to exhibit an ideal slip bond. Finally, we compute the force-dependent unbinding rate of biotin–streptavidin. Here, the complex nature of the unbinding process causes the infrequent metadynamics method to begin to break down due to the presence of unbinding intermediates, despite the use of a previously optimized sampling coordinate. Allowing for this limitation, a combination of kinetic and free energy computations predicts an overall slip bond for larger forces consistent with prior experimental results although there are substantial deviations at small forces that require further investigation. This work demonstrates the promise of predicting force-dependent unbinding rates using enhanced sampling MD techniques while also revealing the methodological barriers that must be overcome to tackle more complex targets in the future.

Recent Advances and Emerging Challenges in the Molecular Modeling of Mechanobiological Processes

Many biological processes result from the effect of mechanical forces on macromolecular structures and on their interactions. In particular, the cell shape, motion, and differentiation directly depend on mechanical stimuli from the extracellular matrix or from neighboring cells. The development of experimental techniques that can measure and characterize the tiny forces acting at the cellular scale and down to the single-molecule, biomolecular level has enabled access to unprecedented details about the involved mechanisms. However, because the experimental observables often do not provide a direct atomistic picture of the corresponding phenomena, particle-based simulations performed at various scales are instrumental in complementing these experiments and in providing a molecular interpretation. Here, we will review the recent key achievements in the field, and we will highlight and discuss the many technical challenges these simulations are facing, as well as suggest future directions for improvement.

Exposure of Von Willebrand Factor Cleavage Site in A1A2A3-Fragment under Extreme Hydrodynamic Shear

Von Willebrand Factor (vWf) is a giant multimeric extracellular blood plasma involved in hemostasis. In this work we present multi-scale simulations of its three-domains fragment A1A2A3. These three domains are essential for the functional regulation of vWf. Namely the A2 domain hosts the site where the protease ADAMTS13 cleavages the multimeric vWf allowing for its length control that prevents thrombotic conditions. The exposure of the cleavage site follows the elongation/unfolding of the domain that is caused by an increased shear stress in blood. By deploying Lattice Boltzmann molecular dynamics simulations based on the OPEP coarse-grained model for proteins, we investigated at molecular level the unfolding of the A2 domain under the action of a perturbing shear flow. We described the structural steps of this unfolding that mainly concerns the β-strand structures of the domain, and we compared the process occurring under shear with that produced by the action of a directional pulling force, a typical condition of single molecule experiments. We observe, that under the action of shear flow, the competition among the elongational and rotational components of the fluid field leads to a complex behaviour of the domain, where elongated structures can be followed by partially collapsed melted globule structures with a very different degree of exposure of the cleavage site. Our simulations pose the base for the development of a multi-scale in-silico description of vWf dynamics and functionality in physiological conditions, including high resolution details for molecular relevant events, e.g., the binding to platelets and collagen during coagulation or thrombosis.

Molecular Paradigms for Biological Mechanosensing

Many proteins in living cells are subject to mechanical forces, which can be generated internally by molecular machines, or externally, e.g., by pressure gradients. In general, these forces fall in the piconewton range, which is similar in magnitude to forces experienced by a molecule due to thermal fluctuations. While we would naively expect such moderate forces to produce only minimal changes, a wide variety of "mechanosensing" proteins have evolved with functions that are responsive to forces in this regime. The goal of this article is to provide a physical chemistry perspective on protein-based molecular mechanosensing paradigms used in living systems, and how these paradigms can be explored using novel computational methods.

Two energy barriers and a transient intermediate state determine the unfolding and folding dynamics of cold shock protein

Cold shock protein (Csp) is a typical two-state folding model protein which has been widely studied by biochemistry and single molecule techniques. Recently two-state property of Csp was confirmed by atomic force microscopy (AFM) through direct pulling measurement, while several long-lifetime intermediate states were found by force-clamp AFM. We systematically studied force-dependent folding and unfolding dynamics of Csp using magnetic tweezers with intrinsic constant force capability. Here we report that Csp mostly folds and unfolds with a single step over force range from 5 pN to 50 pN, and the unfolding rates show different force sensitivities at forces below and above ~8 pN, which determines a free energy landscape with two barriers and a transient intermediate state between them along one transition pathway. Our results provide a new insight on protein folding mechanism of two-state proteins.

Transient Diffusive Interactions with a Protein Crowder Affect Aggregation Processes of Superoxide Dismutase 1 β-Barrel

Aggregate formation of superoxide dismutase 1 (SOD1) inside motor neurons is known as a major factor in onset of amyotrophic lateral sclerosis. The thermodynamic stability of the SOD1 β-barrel has been shown to decrease in crowded environments such as inside a cell, but it remains unclear how the thermodynamics of crowding-induced protein destabilization relate to SOD1 aggregation. Here we have examined the effects of a protein crowder, lysozyme, on fibril aggregate formation of the SOD1 β-barrel. We found that aggregate formation of SOD1 is decelerated even in mildly crowded solutions. Intriguingly, transient diffusive interactions with lysozyme do not significantly affect the static structure of the SOD1 β-barrel but stabilize an alternative excited "invisible" state. The net effect of crowding is to favor species off the aggregation pathway, thereby explaining the decelerated aggregation in the crowded environment. Our observations suggest that the intracellular environment may have a similar negative (inhibitory) effect on fibril formation of other amyloidogenic proteins in living cells. Deciphering how crowded intracellular environments affect aggregation and fibril formation of such disease-associated proteins will probably become central in understanding the exact role of aggregation in the etiology of these enigmatic diseases.

2D and 3D inkjet printing of biopharmaceuticals - A review of trends and future perspectives in research and manufacturing

There is an ongoing global shift in pharmaceutical business models from small molecule drugs to biologics. This increase in complexity is in response to advancements in our diagnoses and understanding of diseases. With the more targeted approach coupled with its inherently more costly development and manufacturing, 2D and 3D printing are being explored as suitable techniques to deliver more personalised and affordable routes to drug discovery and manufacturing. In this review, we explore first the business context underlying this shift to biopharmaceuticals and provide an update on the latest work exploring discovery and pharmaceutics. We then draw on multiple disciplines to help reveal the shared challenges facing researchers and firms aiming to develop biopharmaceuticals, specifically when using the most commonly explored manufacturing routes of drop-on-demand inkjet printing and pneumatic extrusion. This includes separating out how to consider mechanical and chemical influences during manufacturing, the role of the chosen hardware and the challenges of aqueous formulation based on similar challenges being faced by the printing industry. Together, this provides a review of existing work and guidance for researchers and industry to help with the de-risking and rapid development of future biopharmaceutical products.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

Intradimer forces and their implication for conformations of von Willebrand factor multimers

The largest blood glycoprotein von Willebrand factor (VWF) responds to hydrodynamic stresses in the bloodstream with abrupt conformation changes, thus increasing its adhesivity to platelets and collagen. Arterial and microvascular hemostasis relies on mechanical and physicochemical properties of this macromolecule. Recently, it was discovered that the mechanical properties of VWF are controlled by multiple pH-dependent interactions with opposite trends within dimeric subunits. In this work, computer simulations reveal the effect of these intradimer forces on the conformation of VWF multimers in various hydrodynamic conditions. A coarse-grained computer model of VWF has been proposed and parameterized to give a good agreement with experimental data. The simulations suggest that strong attraction between VWF D4 domains increases the resistance to elongation under shear stress, whereas even intermediate attraction between VWF C domains contributes to VWF compaction in nonsheared fluid. It is hypothesized that the detailed subdimer dynamics of VWF concatamers may be one of the biophysical regulators of initial hemostasis and arterial thrombosis.

The Unfolding Journey of Superoxide Dismutase 1 Barrels under Crowding: Atomistic Simulations Shed Light on Intermediate States and Their Interactions with Crowders

The thermal stability of the superoxide dismutase 1 protein in a crowded solution is investigated by performing enhanced sampling molecular simulations. By complementing thermal unfolding experiments done close to physiological conditions (200 mg/mL), we provide evidence that the presence of the protein crowder bovine serum albumin in different packing states has only a minor, and essentially destabilizing, effect. The finding that quinary interactions counteract the pure stabilization contribution stemming from excluded volume is rationalized here by exploring the SOD1 unfolding mechanism in microscopic detail. In agreement with recent experiments, we unveil the importance of intermediate unfolded states as well as the correlation between protein conformations and local packing with the crowders. This link helps us to elucidate why certain SOD1 mutations involved in the ALS disease reverse the stability effect of the intracellular environment.

Aggregation of disease-related peptides

Protein misfolding and aggregation of amyloid proteins is the fundamental cause of more than 20 diseases. Molecular mechanisms of the self-assembly and the formation of the toxic aggregates are still elusive. Computer simulations have been intensively used to study the aggregation of amyloid peptides of various amino acid lengths related to neurodegenerative diseases. We review atomistic and coarse-grained simulations of short amyloid peptides aimed at determining their transient oligomeric structures and the early and late aggregation steps.

Thermal versus mechanical unfolding in a model protein

Force spectroscopy techniques are often used to learn about the free energy landscape of single biomolecules, typically by recovering free energy quantities that, extrapolated to zero force, are compared to those measured in bulk experiments. However, it is not always clear how the information obtained from a mechanically perturbed system can be related to the information obtained using other denaturants since tensioned molecules unfold and refold along a reaction coordinate imposed by the force, which is not likely to be meaningful in its absence. Here, we explore this dichotomy by investigating the unfolding landscape of a model protein, which is unfolded first mechanically through typical force spectroscopy-like protocols and next thermally. When unfolded by nonequilibrium force extension and constant force protocols, we recover a simple two-barrier landscape as the protein reaches the extended conformation through a metastable intermediate. Interestingly, folding-unfolding equilibrium simulations at low forces suggested a totally different scenario, where this metastable state plays little role in the unfolding mechanism, and the protein unfolds through two competing pathways [R. Tapia-Rojo et al., J. Chem. Phys. 141, 135102 (2014)]. Finally, we use Markov state models to describe the configurational space of the unperturbed protein close to the critical temperature. The thermal dynamics is well understood by a one-dimensional landscape along an appropriate reaction coordinate, however it is very different from the mechanical picture. In this sense, the results of our protein model for the mechanical and thermal descriptions provide incompatible views of the folding/unfolding landscape of the system, and the estimated quantities to zero force result are hard to interpret.

Modelling lipid systems in fluid with Lattice Boltzmann Molecular Dynamics simulations and hydrodynamics

In this work we present the coupling between Dry Martini, an efficient implicit solvent coarse-grained model for lipids, and the Lattice Boltzmann Molecular Dynamics (LBMD) simulation technique in order to include naturally hydrodynamic interactions in implicit solvent simulations of lipid systems. After validating the implementation of the model, we explored several systems where the action of a perturbing fluid plays an important role. Namely, we investigated the role of an external shear flow on the dynamics of a vesicle, the dynamics of substrate release under shear, and inquired the dynamics of proteins and substrates confined inside the core of a vesicle. Our methodology enables future exploration of a large variety of biological entities and processes involving lipid systems at the mesoscopic scale where hydrodynamics plays an essential role, e.g. by modulating the migration of proteins in the proximity of membranes, the dynamics of vesicle-based drug delivery systems, or, more generally, the behaviour of proteins in cellular compartments.