How force unfolding differs from chemical denaturation

Toward a Molecular Mechanism of Complementary RNA Duplexes Denaturation

RNA duplexes are relatively rare but play very important biological roles. As an end-product of template-based RNA replication, they also have key implications for hypothetical primitive forms of life. Unless they are specifically separated by enzymes, these duplexes denature upon a temperature increase. However, mechanistic and kinetic aspects of RNA (and DNA) duplex thermal denaturation remain unclear at the microscopic level. We propose an in silico strategy that probes the thermal denaturation of RNA duplexes and allows for an extensive conformational space exploration along a wide temperature range with atomistic precision. We show that this approach first accounts for the strong sequence and length dependence of the duplexes melting temperature, reproducing the trends seen in the experiments and predicted by nearest-neighbor models. The simulations are then instrumental at providing a molecular picture of the temperature-induced strand separation. The textbook canonical "all-or-nothing" two-state model, very much inspired by the protein folding mechanism, can be nuanced. We demonstrate that a temperature increase leads to significantly distorted but stable structures with extensive base-fraying at the extremities, and that the fully formed duplexes typically do not form around melting. The duplex separation therefore appears as much more gradual than commonly thought.

Effects of External Perturbations on Protein Systems: A Microscopic View

Protein folding can be viewed as the origami engineering of biology resulting from the long process of evolution. Even decades after its recognition, research efforts worldwide focus on demystifying molecular factors that underlie protein structure-function relationships; this is particularly relevant in the era of proteopathic disease. A complex co-occurrence of different physicochemical factors such as temperature, pressure, solvent, cosolvent, macromolecular crowding, confinement, and mutations that represent realistic biological environments are known to modulate the folding process and protein stability in unique ways. In the current review, we have contextually summarized the substantial efforts in unveiling individual effects of these perturbative factors, with major attention toward bottom-up approaches. Moreover, we briefly present some of the biotechnological applications of the insights derived from these studies over various applications including pharmaceuticals, biofuels, cryopreservation, and novel materials. Finally, we conclude by summarizing the challenges in studying the combined effects of multifactorial perturbations in protein folding and refer to complementary advances in experiment and computational techniques that lend insights to the emergent challenges.

Dithiol Based on l-Cysteine and Cysteamine as a Disulfide-Reducing Agent

We report the synthesis, chemical properties, and disulfide bond-reducing performance of a dithiol called NACMEAA, conceived as a hybrid of two biologically relevant thiols: cysteine and cysteamine. NACMEAA is conveniently prepared from inexpensive l-cystine in an efficient manner. As a nonvolatile, highly soluble, and neutral compound at physiological pH with the first thiol pK_a value of 8.0, NACMEAA is reactive and user-friendly. We also demonstrate that NACMEAA reduces disulfide bonds in GSSG and lysozyme.

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.

Measuring (biological) materials mechanics with atomic force microscopy. 3. Mechanical unfolding of biopolymers

Biopolymers, such as polynucleotides, polypeptides and polysaccharides, are macromolecules that direct most of the functions in living beings. Studying the mechanical unfolding of biopolymers provides important information about their molecular elasticity and mechanical stability, as well as their energy landscape, which is especially important in proteins, since their three‐dimensional structure is essential for their correct activity. In this primer, we present how to study the mechanical properties of proteins with atomic force microscopy and how to obtain information about their stability and energetic landscape. In particular, we discuss the preparation of polyprotein constructs suitable for AFM single molecule force spectroscopy (SMFS), describe the parameters used in our force‐extension SMFS experiments and the models and equations employed in the analysis of the data. As a practical example, we show the effect of the temperature on the unfolding force, the distance to the transition state, the unfolding rate at zero force, the height of the transition state barrier, and the spring constant of the protein for a construct containing nine repeats of the I27 domain from the muscle protein titin.

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.

Environmental factors modulating protein conformations and their role in protein aggregation diseases

The adverse physiological conditions have been long known to impact protein synthesis, folding and functionality. Major physiological factors such as the effect of pH, temperature, salt and pressure are extensively studied for their impact on protein structure and homeostasis. However, in the current scenario, the environmental risk factors (pollutants) have gained impetus in research because of their increasing concentrations in the environment and strong epidemiologic link with protein aggregation disorders. Here, we review the physiological and environmental risk factors for their impact on protein conformational changes, misfolding, aggregation, and associated pathological conditions, especially environmental risk factors associated pathologies.

Revealing the Underestimated Anticancer Effect of Azurin by Mechanical Unfolding

As a potential anticancer agent, azurin has attracted extensive attraction among chemists, physicists, and material scientists. Its structural and unfolding/folding information has been partially understood, but some detailed information, such as the difference in the unfolding processes between apo-azurin and holo-azurin, the mechanical stability, and the role of the copper cluster in its stability, has not been addressed adequately, especially at the single-molecule level. Here, we employed AFM-based single-molecule force spectroscopy to investigate the unfolding process of azurin in the apo and holo forms under an external force. The results indicated that the unfolding processes of apo-azurin and holo-azurin are different, and holo-azurin requires a stronger force to unfold than does apo-azurin. The copper cluster exhibited a more significant impact on the stability and the folding process of holo-azurin: the copper cluster was completely broken, and the copper ion left the unfolded azurin during the unfolding process of azurin. We suspected that the presence of the disulfide bond in azurin made the unfolding of the copper cluster different from that in pseudoazurin, which is also a type I copper protein like azurin. Rarely reported in previous studies, the mechanical strength of the Cu-N(His) bond of the copper cluster was obtained in this study, which is weaker than that of most metal-S(Cys) bonds but higher than that of the Fe-N(His) bond. Altogether, our results offer a possible new scenario for azurin to widely extend its anticancer activity.

Protein Unfolding: Denaturant vs. Force

While protein refolding has been studied for over 50 years since the pioneering work of Christian Anfinsen, there have been a limited number of studies correlating results between chemical, thermal, and mechanical unfolding. The limited knowledge of the relationship between these processes makes it challenging to compare results between studies if different refolding methods were applied. Our current work compares the energetic barriers and folding rates derived from chemical, thermal, and mechanical experiments using an immunoglobulin-like domain from the muscle protein titin as a model system. This domain, I83, has high solubility and low stability relative to other Ig domains in titin, though its stability can be modulated by calcium. Our experiments demonstrated that the free energy of refolding was equivalent with all three techniques, but the refolding rates exhibited differences, with mechanical refolding having slightly faster rates. This suggests that results from equilibrium-based measurements can be compared directly but care should be given comparing refolding kinetics derived from refolding experiments that used different unfolding methods.

Protons Are Fast and Smart; Proteins Are Slow and Dumb: On the Relationship of Electrospray Ionization Charge States and Conformations

We present simple considerations of how differences in time scales of motions of protons, the lightest and fastest chemical moiety, and the much longer time scales associated with the dynamics of proteins, among the heaviest and slowest analytes, may allow many protein conformations from solution to be kinetically trapped during the process of electrospraying protein solutions into the gas phase. In solution, the quantum nature of protons leads them to change locations by tunneling, an instantaneous process; moreover, the Grotthuss mechanism suggests that these small particles can respond nearly instantaneously to the dynamic motions of proteins that occur on much longer time scales. A conformational change is accompanied by favorable or unfavorable variations in the free energy of the system, providing the impetus for solvent ↔ protein proton exchange. Thus, as thermal distributions of protein conformations interconvert, protonation states rapidly respond, as specific acidic and basic sites are exposed or protected. In the vacuum of the mass spectrometer, protons become immobilized in locations that are specific to the protein conformations from which they were incorporated. In this way, conformational states from solution are preserved upon electrospraying them into the gas phase. These ideas are consistent with the exquisite sensitivity of electrospray mass spectra to small changes of the local environment that alter protein structure in solution. We might remember this approximation for the protonation of proteins in solution with the colloquial expression-protons are fast and smart; proteins are slow and dumb.

Characterizing the Interplay between Polymer Solvation and Conformation

Conformational transitions of flexible molecules, especially those driven by hydrophobic effects, tend to be hindered by desolvation barriers. For such transitions, it is thus important to characterize and understand the interplay between solvation and conformation. Using specialized molecular simulations, here we perform such a characterization for a hydrophobic polymer solvated in water. We find that an external potential, which unfavorably perturbs the polymer hydration waters, can trigger a coil-to-globule or collapse transition, and that the relative stabilities of the collapsed and extended states can be quantified by the strength of the requisite potential. Our results also provide mechanistic insights into the collapse transition, highlighting that the bottleneck to polymer collapse is the formation of a sufficiently large cluster, and the collective dewetting of such a cluster. We also study the collapse of the hydrophobic polymer in octane, a nonpolar solvent, and interestingly, we find that the mechanistic details of the transition are qualitatively similar to that in water.

Mechanical stress induced protein precipitation method for drug target screening

The process of protein precipitation can be used to decipher the interaction of ligand and protein. For example, the classic Thermal Proteome Profiling (TPP) method uses heating as the driving force for protein precipitation, to discover the drug target protein. Under heating or other denature forces, the target protein that binds with the drug compound will be more resistant to precipitation than the free protein. Similar to thermal stress, mechanical stress can also induce protein precipitation. Upon mechanical stress, protein will gradually precipitate along with protein conformational changes, which can be exploited for the study of the ligand-protein interaction. Herein, we proposed a Mechanical Stress Induced Protein Precipitation (MSIPP) method for drug target deconvolution. Its streamlined workflow allows in situ sample preparation on the surface of microparticles, from protein precipitation to digestion. The mechanical stress was generated by vortexing the slurry of protein solution and microparticle materials. The mechanical stress induced protein precipitate was captured by the microparticles, which guarantees the MSIPP method to be scalable and user-friendly. The MSIPP method was successfully applied to four drug compounds, Methotrexate, Raltitrexed, SHP099, Geldanamycin and a pan-inhibitor of protein kinases, Staurosporine. Besides, DHFR was demonstrated to be a target of Raltitrexed, which has not been revealed by any other modification-free drug target discovery method yet. Thus, MSIPP is a complementary method to other drug target screening methods.

Does protein unfolding play a functional role in vivo?

Unfolding and refolding of multidomain proteins under force have yet to be recognized as a major mechanism of function for proteins in vivo. In this review, we discuss the inherent properties of multidomain proteins under a force vector from a structural and functional perspective. We then characterize three main systems where multidomain proteins could play major roles through mechanical unfolding: muscular contraction, cellular mechanotransduction, and bacterial adhesion. We analyze how key multidomain proteins for each system can produce a gain-of-function from the perspective of a fine-tuned quantized response, a molecular battery, delivery of mechanical work through refolding, elasticity tuning, protection and exposure of cryptic sites, and binding-induced mechanical changes. Understanding how mechanical unfolding and refolding affect function will have important implications in designing mechano-active drugs against conditions such as muscular dystrophy, cancer, or novel antibiotics.

Tungsten Oxide Nanodots Exhibit Mild Interactions with WW and SH3 Modular Protein Domains

Tungsten oxide nanodot (WO_3-x ) is an active photothermal nanomaterial that has recently been discovered as a promising candidate for tumor theranostics and treatments. However, its potential cytotoxicity remains elusive and needs to be evaluated to assess its biosafety risks. Herein, we investigate the interactions between WO_3-x and two ubiquitous protein domains involved in protein-protein interactions, namely, WW and SH3 domains, using atomistic molecular dynamics simulations. Our results show that WO_3-x interacts only weakly with the key residues at the putative proline-rich motif (PRM) ligand-binding site of both domains. More importantly, our free energy landscape calculations reveal that the binding strength between WO_3-x and WW/SH3 is weaker than that of the native PRM ligand with WW/SH3, implying that WO_3-x has a limited inhibitory effect over PRM on both the WW and SH3 domains. These findings suggest that the cytotoxic effects of WO_3-x on the key modular protein domains could be very mild, which provides new insights for the future potential biomedical applications of this nanomaterial.

Differences in stability and calcium sensitivity of the Ig domains in titin's N2A region

Titin is a large filamentous protein that spans half a sarcomere, from Z-disk to M-line. The N2A region within the titin molecule exists between the proximal immunoglobulin (Ig) region and the PEVK region and protein-protein interactions involving this region are required for normal muscle function. The N2A region consists of four Ig domains (I80-I83) with a 105 amino acid linker region between I80 and I81 that has a helical nature. Using chemical stability measurements, we show that predicted differences between the adjacent Ig domains (I81-I83) correlate with experimentally determined differences in chemical stability and refolding kinetics. Our work further shows that I83 has the lowest ΔG_unfolding , which is increased in the presence of calcium (pCa 4.3), indicating that Ca²⁺ plays a role in stabilizing this immunoglobulin domain. The characteristics of N2A's three Ig domains provide insight into the stability of the binding sites for proteins that interact with the N2A region. This work also provides insights into how Ca²⁺ might influence binding events involving N2A.

Directional extraction and penetration of phosphorene nanosheets to cell membranes

Recently, phosphorene, a novel two-dimensional nanomaterial with a puckered surface morphology, was shown to exhibit cytotoxicity, but its underlying molecular mechanisms remain unknown. Herein, using large scale molecular dynamics simulations, we show that phosphorene nanosheets can penetrate into and extract large amounts of phospholipids from the cell membranes due to the strong dispersion interaction between phosphorene and lipid molecules, which would reduce cell viability. The extracted phospholipid molecules are aligned along the wrinkle direction of the phosphorene nanosheet because of its unique puckered structure. Our results also reveal that small phosphorene nanosheets penetrate into the cell membrane in a specific direction which is determined by the size and surface topography of phosphorene and the thickness of the membrane. These findings might shed light on understanding phosphorene's cytotoxicity and would be helpful for the future potential biomedical applications of phosphorene, such as biosensors and antibacterial agents.

Mechanical Unfolding of Spectrin Repeats Induces Water-Molecule Ordering

Mechanical processes are involved at many stages of the development of living cells, and often external forces applied to a biomolecule result in its unfolding. Although our knowledge of the unfolding mechanisms and the magnitude of the forces involved has evolved, the role that water molecules play in the mechanical unfolding of biomolecules has not yet been fully elucidated. To this end, we investigated with steered molecular dynamics simulations the mechanical unfolding of dystrophin's spectrin repeat 1 and related the changes in the protein's structure to the ordering of the surrounding water molecules. Our results indicate that upon mechanically induced unfolding of the protein, the solvent molecules become more ordered and increase their average number of hydrogen bonds. In addition, the unfolded structures originating from mechanical pulling expose an increasing amount of the hydrophobic residues to the solvent molecules, and the uncoiled regions adapt a convex surface with a small radius of curvature. As a result, the solvent molecules reorganize around the protein's small protrusions in structurally ordered waters that are characteristic of the so-called "small-molecule regime," which allows water to maintain a high hydrogen bond count at the expense of an increased structural order. We also determined that the response of water to structural changes in the protein is localized to the specific regions of the protein that undergo unfolding. These results indicate that water plays an important role in the mechanically induced unfolding of biomolecules. Our findings may prove relevant to the ever-growing interest in understanding macromolecular crowding in living cells and their effects on protein folding, and suggest that the hydration layer may be exploited as a means for short-range allosteric communication.

Differences in thermal structural changes and melting between mesophilic and thermophilic dihydrofolate reductase enzymes

The thermal resistance of two homolog enzymes is investigated, with an emphasis on their local stability and flexibility, and on the possible implications regarding their reactivity.

The nanomechanics of individual proteins

This tutorial review provides an overview of the single protein force spectroscopy field, including the main techniques and the basic tools for analysing the data obtained from the single molecule experiments.

The molecular mechanism of robust macrophage immune responses induced by PEGylated molybdenum disulfide

Molybdenum disulfide (MoS2), a representative hexagonal transition metal dichalcogenide (TMD), has been extensively exploited in biomedical applications due to its unique physicochemical properties and biocompatibility. However, the lack of adequate data regarding how MoS2 activates immunological responses of macrophages remains a key concern for its risk assessment. Here, we employ a combined theoretical and experimental approach to investigate the interactions of MoS2 and PEGylated MoS2 (MoS2-PEG) with macrophages. We first perform molecular dynamics simulations to examine the atomic-detailed interactions of MoS2 and MoS2-PEG nanoflakes with a realistic model of the macrophage membrane. We show that a small MoS2 nanoflake (edge length of 2.86 nm) is capable of penetrating the macrophage membrane independent of its concentration. We also demonstrate that when initiated with a corner point-on configuration, the surface-bound PEG chains of MoS2-PEG hinder its membrane insertion process, leading to a prolonged passage through the membrane. Moreover, when placed in a face-on arrangement initially, the MoS2-PEG exhibits a lower binding free energy than pristine MoS2 after its adsorption on the membrane surface. The PEG chains can even insert and get buried in the outer leaflet of the membrane, providing additional contact for membrane adsorption. Our flow cytometric experiments then show that the responses of macrophages to either MoS2-PEG or MoS2 are significantly higher than that of the control (no nanomaterial stimulus), with MoS2-PEG eliciting stronger cytokine secretion than the pristine MoS2. The characteristics of slower/prolonged membrane penetration and stronger membrane adsorption of MoS2-PEG compared to pristine MoS2 explain why it triggers more sustained stimulation and higher cytokine secretion in macrophages as observed in our experiments. Our findings reveal the underlying molecular mechanism of how MoS2-PEG influences the immune responses and suggest its potential applications in nanomedicine involving immune stimulation.

Modeling and Structural Characterization of the Sweet Taste Receptor Heterodimer

Sweet taste receptor, a heterodimer belonging to the class C G-protein coupled receptor (GPCR) family and composed of the T1R2 and T1R3 subunits, is responsible for the perception of natural sugars, sweet proteins, various d-amino acids, as well as artificial sweeteners. Despite the critical importance of the sweet receptor not only in mediating gustation but also in its role in the food industry, the architecture of the T1R2-T1R3 complex and the mechanism by which extracellular stimuli induce conformational changes that are propagated to the intracellular milieu, i.e., the signal transduction pathway, remain largely unknown. Here, we constructed and characterized a full-length structural model of the T1R2-T1R3 receptor, including both the transmembrane (TM) and extracellular (EC) domains of the heterodimer, using comparative modeling and extensive all-atom molecular dynamics simulations. Several heterodimer interfaces were first examined for the TM domain, and conformational changes occurring at the intracellular side and associated with the receptor's activation were characterized. From the analysis on the simulated data, putative allosteric binding sites for ligands, ions, and cholesterol were proposed. Also, insights into the protein interface of the TM domain upon activation are provided. The EC domain of the heterodimer, including both the Venus flytrap and cysteine-rich domains, was also investigated. Several important intersubunit interactions located at regions responsible for the receptor's proper function were observed, which resemble those recently identified in other class C GPCR members. Integration of the results from the TM and EC domains facilitates the generation of a full-length T1R2-T1R3 receptor. These findings along with the full-length structural model of the T1R2-T1R3 receptor provide a structural framework that may assist in understanding the mechanistic details associated with the receptor activation process for the sweet T1R2-T1R3 receptor as well as other members of the same family.

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.

Different platinum crystal surfaces show very distinct protein denaturation capabilities

Different platinum (Pt) surfaces of nanocrystals usually exhibit significant distinctions with regard to various biological, physical, and chemical characteristics, such as bio-recognition, surface wetting, and catalytic activities. In this study, we report for the first time that two shape-controlled Pt nanocrystals with the most common low-index surfaces, Pt(100) and Pt(111), show very dissimilar protein denaturation capabilities based on all-atom molecular dynamics simulations employing the widely used model protein, villin headpiece (HP35). We demonstrate that HP35 is well preserved on the Pt(100) crystal surface, whereas it is severely disrupted on the Pt(111) crystal surface. This surprising difference originates from the distinct water behavior in the first solvation shell (FSS) of the two Pt crystal surfaces. Within the FSS of the Pt(100) crystal surface, water molecules form a very compact and stable monolayer through a highly uniform rhombic hydrogen-bond network. This water monolayer prefers the adsorption of acidic residues (such as Glu and Asp) and acts as a shield to prevent other residues from directly coming into contact with the metal surface. On the other hand, the hydrogen bond network in the water monolayer in the FSS of the Pt(111) crystal surface is very sparse and quite defective, which makes it more vulnerable to the penetration of various residues, particularly those with planar side chains such as Phe, Trp and Arg due to strong dispersion interactions, leading to subsequent protein unfolding. The binding free energy calculations for some key amino acids on the two different crystal surfaces further uncover the molecular origin behind their distinct protein denaturation capability. Our study reveals the vital importance of interfacial water in determining the structure of proteins when binding to different metal crystal surfaces. The discovered molecular mechanisms may be helpful for the future development of a bio-assisted programmable synthetic strategy of sophisticated Pt nanostructures for biomedical applications.

Surface Inhomogeneity of Graphene Oxide Influences Dissociation of Aβ16-21 Peptide Assembly

Abnormal peptide assembly and aggregation is associated with an array of neurodegenerative diseases including Alzheimer's disease (AD). A detailed understanding of how nanostructured materials such as oxidized graphene perturb the peptide assembly and subsequently induce fibril dissociation may open new directions for the development of potential AD treatments. Here, we investigate the impact of surface inhomogeneity of graphene oxide (GO) on the assembly of amyloid-beta Aβ_16-21 peptides on GO surfaces with different degrees of oxidation using molecular dynamics simulations. Interestingly, nonuniform GO nanosheets (in terms of oxidation sites) have a much stronger perturbation effect on the structure of Aβ_16-21 assembly. The Aβ peptides exhibit a remarkable tendency in binding to the scattered interfaces between unoxidized and oxidized regions, which induces the dissociation of Aβ amyloid fibril. These findings should deepen our understanding of surface-induced peptide dissociation and stimulate discovery of alternative AD treatments.

Lanosterol Disrupts the Aggregation of Amyloid-β Peptides

Lanosterol, an amphipathic molecule, was discovered only very recently to effectively hinder the aggregation of lens proteins and dissolve the extremely stable fibrillar aggregates in cataracts. Here, we combined computational and experimental approaches to study how lanosterol disrupts the aggregation of another important peptide, amyloid-β (Aβ) peptide, associated with the Alzheimer's Disease (AD). Molecular dynamics simulations using the core amyloidogenic segment (KLVFFA) of Aβ peptide revealed that lanosterol exhibits at least two types of inhibition mechanism on the self-assembly of Aβ peptides. First, lanosterol entangles with peptides and forms a hydrophobic core with residues Phe-19 and Phe-20 in particular. Second, it interferes with the steric zipper interaction at the β-sheet-β-sheet interface. These simulation data suggest that lanosterol induces the unfolding of the Aβ peptide and the separation of the β-sheet layers. This predicted inhibition effect of lanosterol was then confirmed by an in vitro ThT fluorescence assay and AFM imaging. The cell toxicity assay also showed that the treatment of lanosterol indeed mitigates the cytotoxicity of the Aβ peptide in PC-12 cells. Moreover, lanosterol shows a stronger suppression effect on Aβ peptides' aggregation than cholesterol because of its higher hydrophobicity. This result establishes a foundation for the development of lanosterol-based potential therapies for AD and other protein conformational diseases.

Parameterization of Molybdenum Disulfide Interacting with Water Using the Free Energy Perturbation Method

Water contact angles (WCA) are often used to parametrize force field parameters of novel 2D nanomaterials, such as molybdenum disulfide (MoS₂), which has emerged as a promising nanomaterial in many biomedical applications due to its unique and impressive properties. However, there is a wide range of water-MoS₂ contact angles in the literature depending on the aging process on the surface of a MoS₂ nanosheet and/or substrate material. In this study, we revisit and optimize existing parameters for the basal plane of MoS₂ with two popular water models, TIP3P and SPC/E, using the wide range of WCAs from various experiments. We develop and deploy the free energy perturbation method for parametrizing MoS₂ with experimentally determined WCAs for both fresh and aged surfaces. Energy decomposition analysis on the simulation trajectories reveals that MoS₂-water interaction is dominated by van der Waals interaction, which mainly comes from the top layer of MoS₂. We conclude that to describe both fresh and aged MoS₂ surfaces it is convenient to only adjust the Lennard-Jones parameter ε_S (the depth of the potential well of a sulfur atom), which displays a surprisingly linear correlation with WCAs.

Robust Antibacterial Activity of Tungsten Oxide (WO3-x) Nanodots

Antibacterial agents are an important tool in the prevention of bacterial infections. Inorganic materials are attractive due to their high stability under a variety of conditions compared to organic antibacterial agents. Herein tungsten oxide nanodots (WO_3-x), synthesized by a simple one-pot synthetic approach, were found to exhibit strong antibacterial capabilities. The analyses with colony-forming units (CFU) showed an excellent antibacterial activity of WO_3-x against both Gram-negative Escherichia coli (E. coli) and Gram-positive Staphylococcus aureus (S. aureus) strains. The scanning electron microscopy (SEM) and transmission electron microscopy (TEM) images revealed clear damages to the bacterial cell membranes, which was further confirmed by molecular dynamics simulations. Additionally, exposure to simulated sunlight was found to further increase the germicidal activity of WO_3-x nanodots, a 30 min exposure to sunlight combined with 50 μg/mL WO_3-x nanodots showed a 70% decrease in E. coli viability compared to without exposure. Electron spin resonance spectroscopy (ESR) was used to elucidate the underlying mechanism of this photocatalytic activity through the generation of hydroxyl radical species. The cell counting kit-8 (CCK-8) and the live/dead assay were further employed to evaluate the cytotoxicity of WO_3-x nanodots on eukaryotic cells, which demonstrated their general biocompatibility. In summary, our results suggest WO_3-x nanodots have considerable potential in antibacterial applications, while also being biocompatible at large.

Role of Disulfide Bonds and Topological Frustration in the Kinetic Partitioning of Lysozyme Folding Pathways

Disulfide bonds in proteins can strongly influence the folding pathways by constraining the conformational space. Lysozyme has four disulfide bonds and is widely studied for its antibacterial properties. Experiments on lysozyme infer that the protein folds through a fast and a slow pathway. However, the reasons for the kinetic partitioning in the folding pathways are not completely clear. Using a coarse-grained protein model and simulations, we show that two out of the four disulfide bonds, which are present in the α-domain of lysozyme, are responsible for the slow folding pathway. In this pathway, a kinetically trapped intermediate state, which is close to the native state, is populated. In this state, the orientations of α-helices present in the α-domain are misaligned relative to each other. The protein in this state has to partially unfold by breaking down the interhelical contacts between the misaligned helices to fold to the native state. However, the topological constraints due to the two disulfide bonds present in the α-domain make the protein less flexible, and it is trapped in this conformation for hundreds of milliseconds. On disabling these disulfide bonds, we find that the kinetically trapped intermediate state and the slow folding pathway disappear. Simulations mimicking the folding of protein without disulfide bonds under oxidative conditions show that the native disulfide bonds are formed as the protein folds, indicating that folding guides the formation of disulfide bonds. The sequence of formation of the disulfide bonds is Cys64-Cys80 → Cys76-Cys94 → Cys30-Cys115 → Cys6-Cys127. Any disulfide bond that forms before its precursor in the sequence has to break and follow the sequence for the protein to fold. These results show that lysozyme also serves as a very good model system to probe the role of disulfide bonds and topological frustration in protein folding. The predictions from the simulations can be verified by single-molecule fluorescence resonance energy transfer or single-molecule pulling experiments, which can probe heterogeneity in the folding pathways.

Mechanobiology: protein refolding under force

The application of direct force to a protein enables to probe wide regions of its energy surface through conformational transitions as unfolding, extending, recoiling, collapsing, and refolding. While unfolding under force typically displayed a two-state behavior, refolding under force, from highly extended unfolded states, displayed a more complex behavior. The first recording of protein refolding at a force quench step displayed an initial rapid elastic recoil, followed by a plateau phase at some extension, concluding with a collapse to a final state, at which refolding occurred. These findings stirred a lively discussion, which led to further experimental and theoretical investigation of this behavior. It was demonstrated that the polymeric chain of the unfolded protein is required to fully collapse to a globular conformation for the maturation of native structure. This behavior was modeled using one-dimensional free energy landscape over the end-to-end length reaction coordinate, the collective measured variable. However, at low forces, conformational space is not well captured by such models, and using two-dimensional energy surfaces provides further insight into the dynamics of this process. This work reviews the main concepts of protein refolding under constant force, which is essential for understanding how mechanotransducing proteins operate in vivo.

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.

X-ray crystallographic analysis of time-dependent binding of guanidine hydrochloride to HEWL: First steps during protein unfolding

Time-dependent binding of guanidine hydrochloride (GuHCl) to hen egg-white lysozyme (HEWL), and effects of this binding on the protein structure have been investigated by solving X-ray structures of crystalline complexes. The complexes have been prepared by soaking, for different periods of time, native lysozyme crystals in solutions containing 2.5M GuHCl. In the refined structures, the number of water molecules in the protein's first solvent shell has progressively decreased from 152 to 115, showing protein's preference for guanidinium over water. Guanidinium ions preferentially hydrogen bond with the backbone carbonyl oxygen atoms. In their van der Waals interactions, they do not show any preference for apolar residues. Guanidinium ions have replaced water molecules that form cages around exposed hydrophobic residues. Guanidinium binding has decreased the average length of water-water hydrogen bond by 0.1Å. The hydrogen bonds between main chain atoms have been weakened by GuHCl, and this may be the reason for increased potency of GuHCl compared to urea. Guanidinium binding destabilizes the β-domain by causing loss of hydrogen bonds involving Asn 59 side chain. Interestingly, this loss is almost identical to that observed in structures of amyloidogenic variants of human lysozyme. Compounds preventing this loss could be anti-amyloidogenic.

Modeling protein folding in vivo

A half century of studying protein folding in vitro and modeling it in silico has not provided us with a reliable computational method to predict the native conformations of proteins de novo, let alone identify the intermediates on their folding pathways. In this Opinion article, we suggest that the reason for this impasse is the over-reliance on current physical models of protein folding that are based on the assumption that proteins are able to fold spontaneously without assistance. These models arose from studies conducted in vitro on a biased sample of smaller, easier-to-isolate proteins, whose native structures appear to be thermodynamically stable. Meanwhile, the vast empirical data on the majority of larger proteins suggests that once these proteins are completely denatured in vitro, they cannot fold into native conformations without assistance. Moreover, they tend to lose their native conformations spontaneously and irreversibly in vitro, and therefore such conformations must be metastable. We propose a model of protein folding that is based on the notion that the folding of all proteins in the cell is mediated by the actions of the "protein folding machine" that includes the ribosome, various chaperones, and other components involved in co-translational or post-translational formation, maintenance and repair of protein native conformations in vivo. The most important and universal component of the protein folding machine consists of the ribosome in complex with the welcoming committee chaperones. The concerted actions of molecular machinery in the ribosome peptidyl transferase center, in the exit tunnel, and at the surface of the ribosome result in the application of mechanical and other forces to the nascent peptide, reducing its conformational entropy and possibly creating strain in the peptide backbone. The resulting high-energy conformation of the nascent peptide allows it to fold very fast and to overcome high kinetic barriers along the folding pathway. The early folding intermediates in vivo are stabilized by interactions with the ribosome and welcoming committee chaperones and would not be able to exist in vitro in the absence of such cellular components. In vitro experiments that unfold proteins by heat or chemical treatment produce denaturation ensembles that are very different from folding intermediates in vivo and therefore have very limited use in reconstructing the in vivo folding pathways. We conclude that computational modeling of protein folding should deemphasize the notion of unassisted thermodynamically controlled folding, and should focus instead on the step-by-step reverse engineering of the folding process as it actually occurs in vivo.

REVIEWERS

This article was reviewed by Eugene Koonin and Frank Eisenhaber.

Dynamic coarse-graining fills the gap between atomistic simulations and experimental investigations of mechanical unfolding

We present a dynamic coarse-graining technique that allows one to simulate the mechanical unfolding of biomolecules or molecular complexes on experimentally relevant time scales. It is based on Markov state models (MSMs), which we construct from molecular dynamics simulations using the pulling coordinate as an order parameter. We obtain a sequence of MSMs as a function of the discretized pulling coordinate, and the pulling process is modeled by switching among the MSMs according to the protocol applied to unfold the complex. This way we cover seven orders of magnitude in pulling speed. In the region of rapid pulling, we additionally perform steered molecular dynamics simulations and find excellent agreement between the results of the fully atomistic and the dynamically coarse-grained simulations. Our technique allows the determination of the rates of mechanical unfolding in a dynamical range from approximately 10^-8/ns to 1/ns thus reaching experimentally accessible time regimes without abandoning atomistic resolution.

Crowding Stabilizes DMSO-Water Hydrogen-Bonding Interactions

Up to 40% of intracellular water is confined due to the dense packing of macromolecules, ions, and osmolytes. Despite the large body of work concerning the effect of additives on the biomolecular structure and stability, the role of crowding and heterogeneity in these interactions is not well understood. Here, infrared spectroscopy and molecular dynamics simulations are used to describe the mechanisms by which crowding modulates hydrogen bonding interactions between water and dimethyl sulfoxide (DMSO). Specifically, we use formamide and dimethylformamide (DMF) as molecular crowders and show that the S═O hydrogen bond populations in aqueous mixtures are increased by both amides. These additives increase the amount of water within the DMSO first solvation shell through two mechanisms: (a) directly stabilizing water-DMSO hydrogen bonds; (b) increasing water exposure by destabilizing DMSO-DMSO self-interactions. Further, we quantified the hydrogen bond enthalpies between the different components: DMSO-water (61 kJ/mol) > DMSO-formamide (32 kJ/mol) > water-water (23 kJ/mol) ≫ formamide-water (4.7 kJ/mol). Spectra of carbonyl stretching vibrations show that DMSO induces the dehydration of amides as a result of strong DMSO-water interactions, which has been suggested as the main mechanism of protein destabilization.

Rare Dissipative Transitions Punctuate the Initiation of Chemical Denaturation in Proteins

Protein unfolding dynamics are bound by their degree of entropy production, a quantity that relates the amount of heat dissipated by a nonequilibrium process to a system's forward and time-reversed trajectories. We here explore the statistics of heat dissipation that emerge in protein molecules subjected to a chemical denaturant. Coupling large molecular dynamics datasets and Markov state models with the theory of entropy production, we demonstrate that dissipative processes can be rigorously characterized over the course of the urea-induced unfolding of the protein chymotrypsin inhibitor 2. By enumerating full entropy production probability distributions as a function of time, we first illustrate that distinct passive and dissipative regimes are present in the denaturation dynamics. Within the dissipative dynamical region, we next find that chymotrypsin inhibitor 2 is strongly driven into unfolded states in which the protein's hydrophobic core has been penetrated by urea molecules and disintegrated. Detailed analyses reveal that urea's interruption of key hydrophobic contacts between core residues causes many of the protein's native structural features to dissolve.

The force-dependent mechanism of DnaK-mediated mechanical folding

It is well established that chaperones modulate the protein folding free-energy landscape. However, the molecular determinants underlying chaperone-mediated mechanical folding remain largely elusive, primarily because the force-extended unfolded conformation fundamentally differs from that characterized in biochemistry experiments. We use single-molecule force-clamp spectroscopy, combined with molecular dynamics simulations, to study the effect that the Hsp70 system has on the mechanical folding of three mechanically stiff model proteins. Our results demonstrate that, when working independently, DnaJ (Hsp40) and DnaK (Hsp70) work as holdases, blocking refolding by binding to distinct substrate conformations. Whereas DnaK binds to molten globule-like forms, DnaJ recognizes a cryptic sequence in the extended state in an unanticipated force-dependent manner. By contrast, the synergetic coupling of the Hsp70 system exhibits a marked foldase behavior. Our results offer unprecedented molecular and kinetic insights into the mechanisms by which mechanical force finely regulates chaperone binding, directly affecting protein elasticity.

Detecting Interactions between Nanomaterials and Cell Membranes by Synthetic Nanopores

Engineered nanomaterials have been increasingly utilized in industry for various consumer products, environmental treatments, energy storage, and biomedical applications. Meanwhile, it has been established that certain nanomaterials can be toxic to biological cells from extensive experimental and theoretical studies. Despite that the exact molecular mechanisms of this nanomaterial toxicity are still not well understood, it is ubiquitous that their interactions with cell membranes, through either endocytosis or penetration (and thus potential lysis), act as the first step toward the inflammation or even the death of a cell. To facilitate the study of nanomaterial-membrane interactions, here we demonstrate a nanopore-based single-molecule approach that can be applied to monitor a specific nanomaterial-membrane interaction in real time. Combined with molecular dynamics and experimental approaches, we show how an ionic current can be used to detect membrane damage by a graphene nanosheet and illustrate the underlying molecular mechanism. More generally, we expect that measured transmembrane ionic currents (both DC and AC) can signify many particle-induced membrane modifications, such as hole formation, particle adsorption, and protein insertion.

Mechanics of Protein Adaptation to High Temperatures

Inspired by Somero's corresponding state principle that relates protein enhanced thermal stability with mechanical rigidity, we deployed state of the art computational techniques (based on atomistic steered molecular dynamics and Hamiltonian-replica exchange simulations) to study the in silico realization of mechanical and thermal unfolding of two homologous Csp proteins that have evolved to thrive in different thermal environments. By complementing recent single-molecule experiments, we unambiguously show that, for these homologues whose structures are very similar, the increased thermal resistance of the thermophilic variant is not associated with an increased mechanical stability. Our approach provides microscopic insights that are otherwise inaccessible to experimental techniques, and explains why the protein weak spots for thermal and mechanical denaturation are distinct.

Studying the mechanical responses of proteins using magnetic tweezers

The mechanical stability of proteins has been extensively studied using AFM as a single-molecule force spectroscopy method. While this has led to many important results, these studies have been mainly limited to fast unfolding at a high-force regime due to the rapid mechanical drift in most AFM stretching experiments. Therefore, there is a gap between the knowledge obtained at a high-force regime and the mechanical properties of proteins at a lower force regime which is often more physiologically relevant. Recent studies have demonstrated that this gap can be addressed by stretching single protein molecules using magnetic tweezers, due to the excellent mechanical stability this technology offers. Here we review magnetic tweezers technology and its current application in studies of the force-dependent stability and interactions of proteins.

An activity transition from NADH dehydrogenase to NADH oxidase during protein denaturation

A decrease in the specific activity of an enzyme is commonly observed when the enzyme is inappropriately handled or is stored over an extended period. Here, we reported a functional transition of an FMN-bound diaphorase (FMN-DI) that happened during the long-term storage process. It was found that FMN-DI did not simply lose its β-nicotinamide adenine diphosphate (NADH) dehydrogenase activity after a long-time storage, but obtained a new enzyme activity of NADH oxidase. Further mechanistic studies suggested that the alteration of the binding strength of an FMN cofactor with a DI protein could be responsible for this functional switch of the enzyme.

The Role of Controlled Surface Topography and Chemistry on Mouse Embryonic Stem Cell Attachment, Growth and Self-Renewal

The success of stem cell therapies relies heavily on our ability to control their fate in vitro during expansion to ensure an appropriate supply. The biophysical properties of the cell culture environment have been recognised as a potent stimuli influencing cellular behaviour. In this work we used advanced plasma-based techniques to generate model culture substrates with controlled nanotopographical features of 16 nm, 38 nm and 68 nm in magnitude, and three differently tailored surface chemical functionalities. The effect of these two surface properties on the adhesion, spreading, and self-renewal of mouse embryonic stem cells (mESCs) were assessed. The results demonstrated that physical and chemical cues influenced the behaviour of these stem cells in in vitro culture in different ways. The size of the nanotopographical features impacted on the cell adhesion, spreading and proliferation, while the chemistry influenced the cell self-renewal and differentiation.

Exploring the Denatured State Ensemble by Single-Molecule Chemo-Mechanical Unfolding: The Effect of Force, Temperature, and Urea

While it is widely appreciated that the denatured state of a protein is a heterogeneous conformational ensemble, there is still debate over how this ensemble changes with environmental conditions. Here, we use single-molecule chemo-mechanical unfolding, which combines force and urea using the optical tweezers, together with traditional protein unfolding studies to explore how perturbants commonly used to unfold proteins (urea, force, and temperature) affect the denatured-state ensemble. We compare the urea m-values, which report on the change in solvent accessible surface area for unfolding, to probe the denatured state as a function of force, temperature, and urea. We find that while the urea- and force-induced denatured states expose similar amounts of surface area, the denatured state at high temperature and low urea concentration is more compact. To disentangle these two effects, we use destabilizing mutations that shift the T_m and C_m. We find that the compaction of the denatured state is related to changing temperature as the different variants of acyl-coenzyme A binding protein have similar m-values when they are at the same temperature but different urea concentration. These results have important implications for protein folding and stability under different environmental conditions.