Hydrodynamic effects on β-amyloid (16-22) peptide aggregation

Fluid flow and amyloid transport and aggregation in the brain interstitial space

The driving mechanisms at the base of the clearance of biological wastes in the brain interstitial space (ISS) are still poorly understood and an actively debated subject. A complete comprehension of the processes that lead to the aggregation of amyloid proteins in such environment, hallmark of the onset and progression of Alzheimer's disease, is of crucial relevance. Here we employ combined computational fluid dynamics and molecular dynamics techniques to uncover the role of fluid flow and proteins transport in the brain ISS. Our work identifies diffusion as the principal mechanism for amyloid-β proteins clearance, whereas fluid advection may lead transport for larger molecular bodies, like amyloid-β aggregates or extracellular vesicles. We also clearly quantify the impact of large nascent prefibrils on the fluid flowing and shearing. Finally, we show that, even in the irregular brain interstitial space (ISS), hydrodynamic interactions enhance amyloid-β aggregation at all stages of the aggregation pathway. Our results are key to understand the role of fluid flow and solvent-solute interplay on therapeutics like antibodies acting in the brain ISS.

Stability and deformation of biomolecular condensates under the action of shear flow

Biomolecular condensates play a key role in cytoplasmic compartmentalization and cell functioning. Despite extensive research on the physico-chemical, thermodynamic, or crowding aspects of the formation and stabilization of the condensates, one less studied feature is the role of external perturbative fluid flow. In fact, in living cells, shear stress may arise from streaming or active transport processes. Here, we investigate how biomolecular condensates are deformed under different types of shear flows. We first model Couette flow perturbations via two-way coupling between the condensate dynamics and fluid flow by deploying Lattice Boltzmann Molecular Dynamics. We then show that a simplified approach where the shear flow acts as a static perturbation (one-way coupling) reproduces the main features of the condensate deformation and dynamics as a function of the shear rate. With this approach, which can be easily implemented in molecular dynamics simulations, we analyze the behavior of biomolecular condensates described through residue-based coarse-grained models, including intrinsically disordered proteins and protein/RNA mixtures. At lower shear rates, the fluid triggers the deformation of the condensate (spherical to oblated object), while at higher shear rates, it becomes extremely deformed (oblated or elongated object). At very high shear rates, the condensates are fragmented. We also compare how condensates of different sizes and composition respond to shear perturbation, and how their internal structure is altered by external flow. Finally, we consider the Poiseuille flow that realistically models the behavior in microfluidic devices in order to suggest potential experimental designs for investigating fluid perturbations in vitro.

Multistep molecular mechanisms of Aβ16-22 fibril formation revealed by lattice Monte Carlo simulations

As a model of self-assembly from disordered monomers to fibrils, the amyloid-β fragment Aβ16-22 was subject to past numerous experimental and computational studies. Because dynamics information between milliseconds and seconds cannot be assessed by both studies, we lack a full understanding of its oligomerization. Lattice simulations are particularly well suited to capture pathways to fibrils. In this study, we explored the aggregation of 10 Aβ16-22 peptides using 65 lattice Monte Carlo simulations, each simulation consisting of 3 × 109 steps. Based on a total of 24 and 41 simulations that converge and do not converge to the fibril state, respectively, we are able to reveal the diversity of the pathways leading to fibril structure and the conformational traps slowing down the fibril formation.

Evolution of large Aβ16-22 aggregates at atomic details and potential of mean force associated to peptide unbinding and fragmentation events

Atomic characterization of large nonfibrillar aggregates of amyloid polypeptides cannot be determined by experimental means. Starting from β-rich aggregates of Y and elongated topologies predicted by coarse-grained simulations and consisting of more than 100 Aβ16-22 peptides, we performed atomistic molecular dynamics (MD), replica exchange with solute scaling (REST2), and umbrella sampling simulations using the CHARMM36m force field in explicit solvent. Here, we explored the dynamics within 3 μs, the free energy landscape, and the potential of mean force associated with either the unbinding of one single peptide in different configurations within the aggregate or fragmentation events of a large number of peptides. Within the time scale of MD and REST2, we find that the aggregates experience slow global conformational plasticity, and remain essentially random coil though we observe slow beta-strand structuring with a dominance of antiparallel beta-sheets over parallel beta-sheets. Enhanced REST2 simulation is able to capture fragmentation events, and the free energy of fragmentation of a large block of peptides is found to be similar to the free energy associated with fibril depolymerization by one chain for longer Aβ sequences.

Optimized OPEP Force Field for Simulation of Crowded Protein Solutions

Macromolecular crowding has profound effects on the mobility of proteins, with strong implications on the rates of intracellular processes. To describe the dynamics of crowded environments, detailed molecular models are needed, capturing the structures and interactions arising in the crowded system. In this work, we present OPEPv7, which is a coarse-grained force field at amino-acid resolution, suited for rigid-body simulations of the structure and dynamics of crowded solutions formed by globular proteins. Using the OPEP protein model as a starting point, we have refined the intermolecular interactions to match the experimentally observed dynamical slowdown caused by crowding. The resulting force field successfully reproduces the diffusion slowdown in homogeneous and heterogeneous protein solutions at different crowding conditions. Coupled with the lattice Boltzmann technique, it allows the study of dynamical phenomena in protein assemblies and opens the way for the in silico rheology of protein solutions.

An operative framework to model mucus clearance in silico by coupling cilia motion with the liquid environment

Mucociliary clearance is the first defense mechanism of the respiratory tract against inhaled particles. This mechanism is based on the collective beating motion of cilia at the surface of epithelial cells. Impaired clearance, either caused by malfunctioning or absent cilia, or mucus defects, is a symptom of many respiratory diseases. Here, by exploiting the lattice Boltzmann particle dynamics technique, we develop a model to simulate the dynamics of multiciliated cells in a two-layer fluid. First, we tuned our model to reproduce the characteristic length- and time-scales of the cilia beating. We then check for the emergence of the metachronal wave as a consequence of hydrodynamic mediated correlations between beating cilia. Finally, we tune the viscosity of the top fluid layer to simulate the mucus flow upon cilia beating, and evaluate the pushing efficiency of a carpet of cilia. With this work, we build a realistic framework that can be used to explore several important physiological aspects of mucociliary clearance.

Ingenuity in performing replica permutation: How to order the state labels for improving sampling efficiency

In the replica-permutation method, an advanced version of the replica-exchange method, all combinations of replicas and parameters are considered for parameter permutation, and a list of all the combinations is prepared. Here, we report that the temperature transition probability depends on how the list is created, especially in replica permutation with solute tempering (RPST). We found that the transition probabilities decrease at large replica indices when the combinations are sequentially assigned to the state labels as in the originally proposed list. To solve this problem, we propose to modify the list by randomly assigning the combinations to the state labels. We performed molecular dynamics simulations of amyloid-β(16-22) peptides using RPST with the "randomly assigned" list (RPST-RA) and RPST with the "sequentially assigned" list (RPST-SA). The results show the decreases in the transition probabilities in RPST-SA are eliminated, and the sampling efficiency is improved in RPST-RA.

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.

Pathway Dependence of the Formation and Development of Prefibrillar Aggregates in Insulin B Chain

Amyloid fibrils have been an important subject as they are involved in the development of many amyloidoses and neurodegenerative diseases. The formation of amyloid fibrils is typically initiated by nucleation, whereas its exact mechanisms are largely unknown. With this situation, we have previously identified prefibrillar aggregates in the formation of insulin B chain amyloid fibrils, which have provided an insight into the mechanisms of protein assembly involved in nucleation. Here, we have investigated the formation of insulin B chain amyloid fibrils under different pH conditions to better understand amyloid nucleation mediated by prefibrillar aggregates. The B chain showed strong propensity to form amyloid fibrils over a wide pH range, and prefibrillar aggregates were formed under all examined conditions. In particular, different structures of amyloid fibrils were found at pH 5.2 and pH 8.7, making it possible to compare different pathways. Detailed investigations at pH 5.2 in comparison with those at pH 8.7 have suggested that the evolution of protofibril-like aggregates is a common mechanism. In addition, different processes of evolution of the prefibrillar aggregates have also been identified, suggesting that the nucleation processes diversify depending on the polymorphism of amyloid fibrils.

Signatures of sluggish dynamics and local structural ordering during ice nucleation

We investigate the microscopic pathway of spontaneous crystallization in the ST2 model of water under deeply supercooled conditions via unbiased classical molecular dynamics simulations. After quenching below the liquid–liquid critical point, the ST2 model spontaneously separates into low-density liquid (LDL) and high-density liquid phases, respectively. The LDL phase, which is characterized by lower molecular mobility and enhanced structural order, fosters the formation of a sub-critical ice nucleus that, after a stabilization time, develops into the critical nucleus and grows. Polymorphic selection coincides with the development of the sub-critical nucleus and favors the formation of cubic (Ic) over hexagonal (Ih) ice. We rationalize polymorphic selection in terms of geometric arguments based on differences in the symmetry of second neighbor shells of ice Ic and Ih, which are posited to favor formation of the former. The rapidly growing critical nucleus absorbs both Ic and Ih crystallites dispersed in the liquid phase, a crystal with stacking faults. Our results are consistent with, and expand upon, recent observations of non-classical nucleation pathways in several systems.

Computational models for studying physical instabilities in high concentration biotherapeutic formulations

Computational prediction of the behavior of concentrated protein solutions is particularly advantageous in early development stages of biotherapeutics when material availability is limited and a large set of formulation conditions needs to be explored. This review provides an overview of the different computational paradigms that have been successfully used in modeling undesirable physical behaviors of protein solutions with a particular emphasis on high-concentration drug formulations. This includes models ranging from all-atom simulations, coarse-grained representations to macro-scale mathematical descriptions used to study physical instability phenomena of protein solutions such as aggregation, elevated viscosity, and phase separation. These models are compared and summarized in the context of the physical processes and their underlying assumptions and limitations. A detailed analysis is also given for identifying protein interaction processes that are explicitly or implicitly considered in the different modeling approaches and particularly their relations to various formulation parameters. Lastly, many of the shortcomings of existing computational models are discussed, providing perspectives and possible directions toward an efficient computational framework for designing effective protein formulations.

Replica permutation with solute tempering for molecular dynamics simulation and its application to the dimerization of amyloid-β fragments

We propose the replica permutation with solute tempering (RPST) by combining the replica-permutation method (RPM) and the replica exchange with solute tempering (REST). Temperature permutations are performed among more than two replicas in RPM, whereas temperature exchanges are performed between two replicas in the replica-exchange method (REM). The temperature transition in RPM occurs more efficiently than in REM. In REST, only the temperatures of the solute region, the solute temperatures, are exchanged to reduce the number of replicas compared to REM. Therefore, RPST is expected to be an improved method taking advantage of these methods. For comparison, we applied RPST, REST, RPM, and REM to two amyloid-β(16–22) peptides in explicit water. We calculated the transition ratio and the number of tunneling events in the temperature space and the number of dimerization events of amyloid-β(16–22) peptides. The results indicate that, in RPST, the number of replicas necessary for frequent random walks in the temperature and conformational spaces is reduced compared to the other three methods. In addition, we focused on the dimerization process of amyloid-β(16–22) peptides. The RPST simulation with a relatively small number of replicas shows that the two amyloid-β(16–22) peptides form the intermolecular antiparallel β-bridges due to the hydrophilic side-chain contact between Lys and Glu and hydrophobic side-chain contact between Leu, Val, and Phe, which stabilizes the dimer of the peptides.

Dynamics of Amyloid Formation from Simplified Representation to Atomistic Simulations

Amyloid fibril formation is an intrinsic property of short peptides, non-disease proteins, and proteins associated with neurodegenerative diseases. Aggregates of the Aβ and tau proteins, the α-synuclein protein, and the prion protein are observed in the brain of Alzheimer's, Parkinson's, and prion disease patients, respectively. Due to the transient short-range and long-range interactions of all species and their high aggregation propensities, the conformational ensemble of these devastating proteins, the exception being for the monomeric prion protein, remains elusive by standard structural biology methods in bulk solution and in lipid membranes. To overcome these limitations, an increasing number of simulations using different sampling methods and protein models have been performed. In this chapter, we first review our main contributions to the field of amyloid protein simulations aimed at understanding the early aggregation steps of short linear amyloid peptides, the conformational ensemble of the Aβ40/42 dimers in bulk solution, and the stability of Aβ aggregates in lipid membrane models. Then we focus on our studies on the interactions of amyloid peptides/inhibitors to prevent aggregation, and long amyloid sequences, including new results on a monomeric tau construct.

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.

Accelerating Reaction Rates of Biomolecules by Using Shear Stress in Artificial Capillary Systems

Biomimetics is a design principle within chemistry, biology, and engineering, but chemistry biomimetic approaches have been generally limited to emulating nature's chemical toolkit while emulation of nature's physical toolkit has remained largely unexplored. To begin to explore this, we designed biophysically mimetic microfluidic reactors with characteristic length scales and shear stresses observed within capillaries. We modeled the effect of shear with molecular dynamics studies and showed that this induces specific normally buried residues to become solvent accessible. We then showed using kinetics experiments that rates of reaction of these specific residues in fact increase in a shear-dependent fashion. We applied our results in the creation of a new microfluidic approach for the multidimensional study of cysteine biomarkers. Finally, we used our approach to establish dissociation of the therapeutic antibody trastuzumab in a reducing environment. Our results have implications for the efficacy of existing therapeutic antibodies in blood plasma as well as suggesting in general that biophysically mimetic chemistry is exploited in biology and should be explored as a research area.

Molecular dynamics simulations of amyloid-β(16-22) peptide aggregation at air-water interfaces

Oligomers of amyloid-β (Aβ) peptides are known to be related to Alzheimer's disease, and their formation is accelerated at hydrophilic-hydrophobic interfaces, such as the cell membrane surface and air-water interface. Here, we report molecular dynamics simulations of aggregation of Aβ(16-22) peptides at air-water interfaces. First, 100 randomly distributed Aβ(16-22) peptides moved to the interface. The high concentration of peptides then accelerated their aggregation and formation of antiparallel β-sheets. Two layers of oligomers were observed near the interface. In the first layer from the interface, the oligomer with less β-bridges exposed the hydrophobic residues to the air. The second layer consisted of oligomers with more β-bridges that protruded into water. They are more soluble in water because the hydrophobic residues are covered by N- and C-terminal hydrophilic residues that are aligned well along the oligomer edge. These results indicate that amyloid protofibril formation mainly occurs in the second layer.

Polyphenol-Peptide Interactions in Mitigation of Alzheimer's Disease: Role of Biosurface-Induced Aggregation

Alzheimer's disease (AD) is the most common age-related neurodegenerative disorder, responsible for nearly two-thirds of all dementia cases. In this review, we report the potential AD treatment strategies focusing on natural polyphenol molecules (green chemistry) and more specifically on the inhibition of polyphenol-induced amyloid aggregation/disaggregation pathways: in bulk and on biosurfaces. We discuss how these pathways can potentially alter the structure at the early stages of AD, hence delaying the aggregation of amyloid-β (Aβ) and tau. We also discuss multidisciplinary approaches, combining experimental and modelling methods, that can better characterize the biochemical and biophysical interactions between proteins and phenolic ligands. In addition to the surface-induced aggregation, which can occur on surfaces where protein can interact with other proteins and polyphenols, we suggest a new concept referred as "confinement stability". Here, on the contrary, the adsorption of Aβ and tau on biosurfaces other than Aβ- and tau-fibrils, e.g., red blood cells, can lead to confinement stability that minimizes the aggregation of Aβ and tau. Overall, these mechanisms may participate directly or indirectly in mitigating neurodegenerative diseases, by preventing protein self-association, slowing down the aggregation processes, and delaying the progression of AD.

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.

Promotion and Inhibition of Amyloid-β Peptide Aggregation: Molecular Dynamics Studies

Aggregates of amyloid-β (Aβ) peptides are known to be related to Alzheimer’s disease. Their aggregation is enhanced at hydrophilic–hydrophobic interfaces, such as a cell membrane surface and air-water interface, and is inhibited by polyphenols, such as myricetin and rosmarinic acid. We review molecular dynamics (MD) simulation approaches of a full-length Aβ peptide, Aβ40, and Aβ(16–22) fragments in these environments. Since these peptides have both hydrophilic and hydrophobic amino acid residues, they tend to exist at the interfaces. The high concentration of the peptides accelerates the aggregation there. In addition, Aβ40 forms a β-hairpin structure, and this structure accelerates the aggregation. We also describe the inhibition mechanism of the Aβ(16–22) aggregation by polyphenols. The aggregation of Aβ(16–22) fragments is caused mainly by the electrostatic attraction between charged amino acid residues known as Lys16 and Glu22. Since polyphenols form hydrogen bonds between their hydroxy and carboxyl groups and these charged amino acid residues, they inhibit the aggregation.

Tracking Internal and Global Diffusive Dynamics During Protein Aggregation by High-Resolution Neutron Spectroscopy

Proteins can misfold and form either amorphous or organized aggregates with different morphologies and features. Aggregates of amyloid nature are pathological hallmarks in so-called protein conformational diseases, including Alzheimer's and Parkinson's. Evidence prevails that the transient early phases of the reaction determine the aggregate morphology and toxicity. As a consequence, real-time monitoring of protein aggregation is of utmost importance. Here, we employed time-resolved neutron backscattering spectroscopy to follow center-of-mass self-diffusion and nano- to picosecond internal dynamics of lysozyme during aggregation into a specific β-sheet rich superstructure, called particulates, formed at the isoelectric point of the protein. Particulate formation is found to be a one-step process, and protein internal dynamics, to remain unchanged during the entire aggregation process. The time-resolved neutron backscattering spectroscopy approach developed here, in combination with standard kinetics assays, provides a unifying framework in which dynamics and conformational transitions can be related to the different aggregation pathways.

Structures of the intrinsically disordered Aβ, tau and α-synuclein proteins in aqueous solution from computer simulations

Intrinsically disordered proteins (IDPs) play many biological roles in the human proteome ranging from vesicular transport, signal transduction to neurodegenerative diseases. The Aβ and tau proteins, and the α-synuclein protein, key players in Alzheimer's and Parkinson's diseases, respectively are fully disordered at the monomer level. The structural heterogeneity of the monomeric and oligomeric states and the high self-assembly propensity of these three IDPs have precluded experimental structural determination. Simulations have been used to determine the atomic structures of these IDPs. In this article, we review recent computer models to capture the equilibrium ensemble of Aβ, tau and α-synuclein proteins at different association steps in aqueous solution and present new results of the PEP-FOLD framework on α-synuclein monomer.

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.

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.

Recent Advances in Coarse-Grained Models for Biomolecules and Their Applications

Molecular dynamics simulations have emerged as a powerful tool to study biological systems at varied length and timescales. The conventional all-atom molecular dynamics simulations are being used by the wider scientific community in routine to capture the conformational dynamics and local motions. In addition, recent developments in coarse-grained models have opened the way to study the macromolecular complexes for time scales up to milliseconds. In this review, we have discussed the principle, applicability and recent development in coarse-grained models for biological systems. The potential of coarse-grained simulation has been reviewed through state-of-the-art examples of protein folding and structure prediction, self-assembly of complexes, membrane systems and carbohydrates fiber models. The multiscale simulation approaches have also been discussed in the context of their emerging role in unravelling hierarchical level information of biosystems. We conclude this review with the future scope of coarse-grained simulations as a constantly evolving tool to capture the dynamics of biosystems.

Multiscale Aggregation of the Amyloid Aβ16-22 Peptide: From Disordered Coagulation and Lateral Branching to Amorphous Prefibrils

In this work we investigate the multiscale dynamics of the aggregation process of an amyloid peptide, Aβ_16-22. By performing massive coarse-grained simulations at the quasi-atomistic resolution and including hydrodynamic effects, we followed the formation and growth of a large elongated aggregate and its slow structuring. The elongation proceeds via a two-step nucleation mechanism with disordered aggregates formed initially and subsequently fusing to elongate the amorphous prefibril. A variety of coagulation events coexist, including lateral growth. The latter mechanism, sustained by long-range hydrodynamics correlations, actually can create a large branched structure spanning a few tens of nanometers. Our findings confirm the experimental hypothesis of a critical contribution of lateral growth to the amyloid aggregation kinetics and the capability of our model to sample critical structures like prefibril hosting annular pores.

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

The perturbation of a protein conformation by a physiological fluid flow is crucial in various biological processes including blood clotting and bacterial adhesion to human tissues. Investigating such mechanisms by computer simulations is thus of great interest, but it requires development of ad hoc strategies to mimic the complex hydrodynamic interactions acting on the protein from the surrounding flow. In this study, we apply the Lattice Boltzmann Molecular Dynamics (LBMD) technique built on the implicit solvent coarse-grained model for protein Optimized Potential for Efficient peptide structure Prediction (OPEP) and a mesoscopic representation of the fluid solvent, to simulate the unfolding of a small globular cold-shock protein in shear flow and to compare it to the unfolding mechanisms caused either by mechanical or thermal perturbations. We show that each perturbation probes a specific weakness of the protein and causes the disruption of the native fold along different unfolding pathways. Notably, the shear flow and the thermal unfolding exhibit very similar pathways, while because of the directionality of the perturbation, the unfolding under force is quite different. For force and thermal disruption of the native state, the coarse-grained simulations are compared to all-atom simulations in explicit solvent, showing an excellent agreement in the explored unfolding mechanisms. These findings encourage the use of LBMD based on the OPEP model to investigate how a flow can affect the function of larger proteins, for example, in catch-bond systems.

Advances in coarse-grained modeling of macromolecular complexes

Recent progress in coarse-grained (CG) molecular modeling and simulation has facilitated an influx of computational studies on biological macromolecules and their complexes. Given the large separation of length-scales and time-scales that dictate macromolecular biophysics, CG modeling and simulation are well-suited to bridge the microscopic and mesoscopic or macroscopic details observed from all-atom molecular simulations and experiments, respectively. In this review, we first summarize recent innovations in the development of CG models, which broadly include structure-based, knowledge-based, and dynamics-based approaches. We then discuss recent applications of different classes of CG models to explore various macromolecular complexes. Finally, we conclude with an outlook for the future in this ever-growing field of biomolecular modeling.

On the Applicability of Force Fields To Study the Aggregation of Amyloidogenic Peptides Using Molecular Dynamics Simulations

Molecular dynamics simulations play an essential role in understanding biomolecular processes such as protein aggregation at temporal and spatial resolutions which are not attainable by experimental methods. For a correct modeling of protein aggregation, force fields must accurately represent molecular interactions. Here, we study the effect of five different force fields on the oligomer formation of Alzheimer's Aβ_16-22 peptide and two of its mutants: Aβ_16-22(F19V,F20V), which does not form fibrils, and Aβ_16-22(F19L) which forms fibrils faster than the wild type. We observe that while oligomer formation kinetics depends strongly on the force field, structural properties, such as the most relevant protein-protein contacts, are similar between them. The oligomer formation kinetics obtained with different force fields differ more from each other than the kinetics between aggregating and nonaggregating peptides simulated with a single force field. We discuss the difficulties in comparing atomistic simulations of amyloid oligomer formation with experimental observables.

Etersalate prevents the formations of 6Aβ16-22 oligomer: An in silico study

Oligomerization of amyloid beta (Aβ) peptides has been considered as the crucially causative agent in the development of Alzheimer's disease. Etersalate, a nonsteroidal anti-inflammatory oral drug (United State Food and Drug Administration—Unique Ingredient Identifier: 653GN04T2G) was previously suggested to bind well to proto-fibrils of Aβ peptides in silico. Here, the effect of etersalate on the oligomerization of soluble Aβ_16–22 hexamer (6Aβ_16–22) were extensively investigated using temperature replica exchange molecular dynamics (REMD) simulations over ~16.8 μs in total for 48 replicas (350 ns per replica). The results reveal that etersalate can enter the inner space or bind on the surface of 6Aβ_16–22 conformations, which destabilizes the hexamer. Etersalate was predicted to able to cross the blood brain barrier using prediction of absorption, distribution, metabolism, and excretion—toxicity (preADMET) tools. Overall, although the investigation was performed with the low concentration of trial inhibitor, the obtained results indicate that etersalate is a potential drug candidate for AD through inhibiting formation of Aβ oligomers with the average binding free energy of -11.7 kcal/mol.

Looking at the Disordered Proteins through the Computational Microscope

Intrinsically disordered proteins (IDPs) have attracted wide interest over the past decade due to their surprising prevalence in the proteome and versatile roles in cell physiology and pathology. A large selection of IDPs has been identified as potential targets for therapeutic intervention. Characterizing the structure-function relationship of disordered proteins is therefore an essential but daunting task, as these proteins can adapt transient structure, necessitating a new paradigm for connecting structural disorder to function. Molecular simulation has emerged as a natural complement to experiments for atomic-level characterizations and mechanistic investigations of this intriguing class of proteins. The diverse range of length and time scales involved in IDP function requires performing simulations at multiple levels of resolution. In this Outlook, we focus on summarizing available simulation methods, along with a few interesting example applications. We also provide an outlook on how these simulation methods can be further improved in order to provide a more accurate description of IDP structure, binding, and assembly.

Impact of hydrodynamic interactions on protein folding rates depends on temperature

We investigated the impact of hydrodynamic interactions (HI) on protein folding using a coarse-grained model. The extent of the impact of hydrodynamic interactions, whether it accelerates, retards, or has no effect on protein folding, has been controversial. Together with a theoretical framework of the energy landscape theory (ELT) for protein folding that describes the dynamics of the collective motion with a single reaction coordinate across a folding barrier, we compared the kinetic effects of HI on the folding rates of two protein models that use a chain of single beads with distinctive topologies: a 64-residue α/β chymotrypsin inhibitor 2 (CI2) protein, and a 57-residue β-barrel α-spectrin Src-homology 3 domain (SH3) protein. When comparing the protein folding kinetics simulated with Brownian dynamics in the presence of HI to that in the absence of HI, we find that the effect of HI on protein folding appears to have a "crossover" behavior about the folding temperature. This means that at a temperature greater than the folding temperature, the enhanced friction from the hydrodynamic solvents between the beads in an unfolded configuration results in lowered folding rate; conversely, at a temperature lower than the folding temperature, HI accelerates folding by the backflow of solvent toward the folded configuration of a protein. Additionally, the extent of acceleration depends on the topology of a protein: for a protein like CI2, where its folding nucleus is rather diffuse in a transition state, HI channels the formation of contacts by favoring a major folding pathway in a complex free energy landscape, thus accelerating folding. For a protein like SH3, where its folding nucleus is already specific and less diffuse, HI matters less at a temperature lower than the folding temperature. Our findings provide further theoretical insight to protein folding kinetic experiments and simulations.

Molecular Mechanism of Protein Unfolding under Shear: A Lattice Boltzmann Molecular Dynamics Study

Proteins are marginally stable soft-matter entities that can be disrupted using a variety of perturbative stresses, including thermal, chemical, or mechanical ones. Fluid under extreme flow conditions is a possible route to probe the weakness of biomolecules and collect information on the molecular cohesive interactions that secure their stability. Moreover, in many cases, physiological flow triggers the functional response of specialized proteins as occurring in blood coagulation or cell adhesion. We deploy the Lattice Boltzmann molecular dynamics technique based on the coarse-grained model for protein OPEP to study the mechanism of protein unfolding under Couette flow. Our simulations provide a clear view of how structural elements of the proteins are affected by shear, and for the simple study case, the β-hairpin, we exploited the analogy to pulling experiments to quantify the mechanical forces acting on the protein under shear.

Pathways of Amyloid-β Aggregation Depend on Oligomer Shape

One of the main research topics related to Alzheimer's disease is the aggregation of the amyloid-β peptide, which was shown to follow different pathways for the two major alloforms of the peptide, Aβ40 and the more toxic Aβ42. Experimental studies emphasized that oligomers of specific sizes appear in the early aggregation process in different quantities and might be the key toxic agents for each of the two alloforms. We use transition networks derived from all-atom molecular dynamics simulations to show that the oligomers leading to the type of oligomer distributions observed in experiments originate from compact conformations. Extended oligomers, on the other hand, contribute more to the production of larger aggregates thus driving the aggregation process. We further demonstrate that differences in the aggregation pathways of the two Aβ alloforms occur as early as during the dimer stage. The higher solvent-exposure of hydrophobic residues in Aβ42 oligomers contributes to the different aggregation pathways of both alloforms and also to the increased cytotoxicity of Aβ42.

HyRes: a coarse-grained model for multi-scale enhanced sampling of disordered protein conformations

Efficient coarse-grained (CG) models can be coupled with atomistic force fields to accelerate the sampling of atomistic energy landscapes in the multi-scale enhanced sampling (MSES) framework. This approach may be particularly suitable for generating atomistic conformational ensembles of intrinsically disordered proteins (IDPs). While MSES is relatively robust to inherent CG artifacts, achieving optimal sampling efficiency requires CG modeling to generate the local and long-range fluctuations that are largely consistent with those at the atomistic level. Here, we describe a new hybrid resolution CG model (HyRes) for MSES simulations of disordered protein states, which is specifically designed to provide semi-quantitative secondary structure propensities together with a qualitative description of long-range nonspecific interactions. The HyRes model contains an atomistic description of the backbone with intermediate resolution side chains. The secondary structure propensities are tuned by adjusting the backbone hydrogen-bonding strength and the ϕ/ψ torsion profile. The sizes and covalent geometries of the side chains are parameterized to reproduce distributions derived from atomistic simulations. Lennard-Jones parameters for sidechain beads are assigned to reproduce statistical potentials derived from the protein structural database, and then globally parameterized with nonspecific electrostatic interactions to reproduce the free energy profiles of pair wise interactions and the key conformational properties of model peptides. Application of HyRes to MSES simulations of small IDPs suggests that it is capable of driving faster structural transitions at the atomistic level and increasing the convergence rate compared to the Cα-only Gō-like models previously utilized. With further optimization, we believe that the new CG model could greatly improve the efficiency of MSES simulations of the larger and more complex IDPs frequently involved in cellular signalling and regulation.

Multi-scale simulations of biological systems using the OPEP coarse-grained model

Biomolecules are complex machines that are optimized by evolution to properly fulfill or contribute to a variety of biochemical tasks in the cellular environment. Computer simulations based on quantum mechanics and atomistic force fields have been proven to be a powerful microscope for obtaining valuable insights into many biological, physical, and chemical processes. Many interesting phenomena involve, however, a time scale and a number of degrees of freedom, notably if crowding is considered, that cannot be explored at an atomistic resolution. To bridge the gap between reality and simulation, many different advanced computational techniques and coarse-grained (CG) models have been developed. Here, we report some applications of the CG OPEP protein model to amyloid fibril formation, the response of catch-bond proteins to two types of fluid flow, and interactive simulations to fold peptides with well-defined 3D structures or with intrinsic disorder.

Simulations of disordered proteins and systems with conformational heterogeneity

Intrinsically disordered proteins (IDPs) and protein regions can facilitate a wide variety of complex physiological processes such as binding, signaling, and formation of membraneless organelles. They can however also play pathological roles by aggregating into cytotoxic oligomers and fibrils. Characterizing the structure and function of disordered proteins is an onerous task, primarily because these proteins adopt transient structures, which are difficult to capture in experiments. Simulations have emerged as a powerful tool for interpreting and augmenting experimental measurements of IDPs. In this review we focus on computer simulations of disordered protein structures, functions, assemblies, and emerging questions that, taken together, give an overview of the field as it exists today.

Multiscale simulation of molecular processes in cellular environments

We describe the recent advances in studying biological systems via multiscale simulations. Our scheme is based on a coarse-grained representation of the macromolecules and a mesoscopic description of the solvent. The dual technique handles particles, the aqueous solvent and their mutual exchange of forces resulting in a stable and accurate methodology allowing biosystems of unprecedented size to be simulated.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.