Mechanical Stability and Unfolding Pathways of Parallel Tetrameric G-Quadruplexes Probed by Pulling Simulations

Guanine quadruplex (GQ) is a noncanonical nucleic acid structure formed by guanine-rich DNA and RNA sequences. Folding of GQs is a complex process, where several aspects remain elusive, despite being important for understanding structure formation and biological functions of GQs. Pulling experiments are a common tool for acquiring insights into the folding landscape of GQs. Herein, we applied a computational pulling strategy─steered molecular dynamics (SMD) simulations─in combination with standard molecular dynamics (MD) simulations to explore the unfolding landscapes of tetrameric parallel GQs. We identified anisotropic properties of elastic conformational changes, unfolding transitions, and GQ mechanical stabilities. Using a special set of structural parameters, we found that the vertical component of pulling force (perpendicular to the average G-quartet plane) plays a significant role in disrupting GQ structures and weakening their mechanical stabilities. We demonstrated that the magnitude of the vertical force component depends on the pulling anchor positions and the number of G-quartets. Typical unfolding transitions for tetrameric parallel GQs involve base unzipping, opening of the G-stem, strand slippage, and rotation to cross-like structures. The unzipping was detected as the first and dominant unfolding event, and it usually started at the 3'-end. Furthermore, results from both SMD and standard MD simulations indicate that partial spiral conformations serve as a transient ensemble during the (un)folding of GQs.

Molecular Mechanism of Interaction between DNA Aptamer and Receptor-Binding Domain of Severe Acute Respiratory Syndrome Coronavirus 2 Variants Revealed by Steered Molecular Dynamics Simulations

The ongoing SARS-CoV-2 pandemic has underscored the urgent need for versatile and rapidly deployable antiviral strategies. While vaccines have been pivotal in controlling the spread of the virus, the emergence of new variants continues to pose significant challenges to global health. Here, our study focuses on a novel approach to antiviral therapy using DNA aptamers, short oligonucleotides with high specificity and affinity for their targets, as potential inhibitors against the spike protein of SARS-CoV-2 variants Omicron and JN.1. Our research utilizes steered molecular dynamics (SMD) simulations to elucidate the binding mechanisms of a specifically designed DNA aptamer, AM032-4, to the receptor-binding domain (RBD) of the aforementioned variants. The simulations reveal detailed molecular insights into the aptamer-RBD interaction, demonstrating the aptamer's potential to maintain effective binding in the face of rapid viral evolution. Our work not only demonstrates the dynamic interaction between aptamer-RBD for possible antiviral therapy but also introduces a computational method to study aptamer-protein interactions.

Measuring mechanical properties with high-speed atomic force microscopy

High-speed atomic force microscopy (HS-AFM) is now a widely used technique to study the dynamics of single biomolecules and complex structures. In the past, it has mainly been used to capture surface topography as structural analysis, leading to important discoveries not attainable by other methods. Similar to conventional AFM, the scope of HS-AFM was recently expanded to encompass quantities beyond topography, such as the measurement of mechanical properties. This review delves into various methodologies for assessing mechanical properties, ranging from semi-quantitative approaches to precise force measurements and their corresponding sample responses. We will focus on the application to single proteins such as bridging integrator-1, ion channels such as Piezo1, complex structures such as microtubules and supramolecular fibers. In all these examples, the unique combination of quantifiable force application and high spatiotemporal resolution allows to unravel mechanisms that cannot be investigated by conventional means.

Magnetic tweezers in cell mechanics

The chapter provides an overview of the applications of magnetic tweezers in living cells. It discusses the advantages and disadvantages of magnetic tweezers technology with a focus on individual magnetic tweezers configurations, such as electromagnetic tweezers. Solutions to the disadvantages identified are also outlined. The specific role of magnetic tweezers in the field of mechanobiology, such as mechanosensitivity, mechano-allostery and mechanotransduction are also emphasized. The specific usage of magnetic tweezers in mechanically probing cells via specific cell surface receptors, such as mechanosensitive channels is discussed and why mechanical probing has revealed the opening and closing of the channels. Finally, the future direction of magnetic tweezers is presented.

Visualizing Molecular Dynamics by High-Speed Atomic Force Microscopy

Dynamic processes and structural changes of biological molecules are essential to life. While conventional atomic force microscopy (AFM) is able to visualize molecules and supramolecular assemblies at sub-nanometer resolution, it cannot capture dynamics because of its low imaging rate. The introduction of high-speed atomic force microscopy (HS-AFM) solved this problem by providing a large increase in imaging velocity. Using HS-AFM, one is able to visualize dynamic molecular events with high spatiotemporal resolution under near-to physiological conditions. This approach opened new windows as finally dynamics of biomolecules at sub-nanometer resolution could be studied. Here we describe the working principles and an operation protocol for HS-AFM imaging and characterization of biological samples in liquid.

Complexity of Guanine Quadruplex Unfolding Pathways Revealed by Atomistic Pulling Simulations

Guanine quadruplexes (GQs) are non-canonical nucleic acid structures involved in many biological processes. GQs formed in single-stranded regions often need to be unwound by cellular machinery, so their mechanochemical properties are important. Here, we performed steered molecular dynamics simulations of human telomeric GQs to study their unfolding. We examined four pulling regimes, including a very slow setup with pulling velocity and force load accessible to high-speed atomic force microscopy. We identified multiple factors affecting the unfolding mechanism, i.e.,: (i) the more the direction of force was perpendicular to the GQ channel axis (determined by GQ topology), the more the base unzipping mechanism happened, (ii) the more parallel the direction of force was, GQ opening and cross-like GQs were more likely to occur, (iii) strand slippage mechanism was possible for GQs with an all-anti pattern in a strand, and (iv) slower pulling velocity led to richer structural dynamics with sampling of more intermediates and partial refolding events. We also identified that a GQ may eventually unfold after a force drop under forces smaller than those that the GQ withstood before the drop. Finally, we found out that different unfolding intermediates could have very similar chain end-to-end distances, which reveals some limitations of structural interpretations of single-molecule spectroscopic data.

High-yield ligation-free assembly of DNA constructs with nucleosome positioning sequence repeats for single-molecule manipulation assays

Force and torque spectroscopy have provided unprecedented insights into the mechanical properties, conformational transitions, and dynamics of DNA and DNA-protein complexes, notably nucleosomes. Reliable single-molecule manipulation measurements require, however, specific and stable attachment chemistries to tether the molecules of interest. Here, we present a functionalization strategy for DNA that enables high-yield production of constructs for torsionally constrained and very stable attachment. The method is based on two subsequent PCRs: first ∼380 bp long DNA strands are generated that contain multiple labels, which are used as "megaprimers" in a second PCR to generate ∼kbp long double-stranded DNA constructs with multiple labels at the respective ends. To achieve high-force stability, we use dibenzocyclooctyne-based click chemistry for covalent attachment to the surface and biotin-streptavidin coupling to the bead. The resulting tethers are torsionally constrained and extremely stable under load, with an average lifetime of 70 ± 3 h at 45 pN. The high yield of the approach enables nucleosome reconstitution by salt dialysis on the functionalized DNA, and we demonstrate proof-of-concept measurements on nucleosome assembly statistics and inner turn unwrapping under force. We anticipate that our approach will facilitate a range of studies of DNA interactions and nucleoprotein complexes under forces and torques.

Uncovering and engineering the mechanical properties of the adhesion GPCR ADGRG1 GAIN domain

Key cellular functions depend on the transduction of extracellular mechanical signals by specialized membrane receptors including adhesion G-protein coupled receptors (aGPCRs). While recently solved structures support aGPCR activation through shedding of the extracellular GAIN domain, the molecular mechanisms underpinning receptor mechanosensing remain poorly understood. When probed using single-molecule atomic force spectroscopy and molecular simulations, ADGRG1 GAIN dissociated from its tethered agonist at forces significantly higher than other reported signaling mechanoreceptors. Strong mechanical resistance was achieved through specific structural deformations and force propagation pathways under mechanical load. ADGRG1 GAIN variants computationally designed to lock the alpha and beta subdomains and rewire mechanically-induced structural deformations were found to modulate the GPS-Stachel rupture forces. Our study provides unprecedented insights into the molecular underpinnings of GAIN mechanical stability and paves the way for engineering mechanosensors, better understanding aGPCR function, and informing drug-discovery efforts targeting this important receptor class.

Recent advances in sensing the inter-biomolecular interactions at the nanoscale - A comprehensive review of AFM-based force spectroscopy

Biomolecular interactions underpin most processes inside the cell. Hence, a precise and quantitative understanding of molecular association and dissociation events is crucial, not only from a fundamental perspective, but also for the rational design of biomolecular platforms for state-of-the-art biomedical and industrial applications. In this context, atomic force microscopy (AFM) appears as an invaluable experimental technique, allowing the measurement of the mechanical strength of biomolecular complexes to provide a quantitative characterization of their interaction properties from a single molecule perspective. In the present review, the most recent methodological advances in this field are presented with special focus on bioconjugation, immobilization and AFM tip functionalization, dynamic force spectroscopy measurements, molecular recognition imaging and theoretical modeling. We expect this work to significantly aid in grasping the principles of AFM-based force spectroscopy (AFM-FS) technique and provide the necessary tools to acquaint the type of data that can be achieved from this type of experiments. Furthermore, a critical assessment is done with other nanotechnology techniques to better visualize the future prospects of AFM-FS.

Novel perspective for protein-drug interaction analysis: atomic force microscope

Proteins are major drug targets, and drug-target interaction identification and analysis are important factors for drug discovery. Atomic force microscopy (AFM) is a powerful tool making it possible to image proteins with nanometric resolution and probe intermolecular forces under physiological conditions. We review recent studies conducted in the field of target protein drug discovery using AFM-based analysis technology, including drug-driven changes in nanomechanical properties of protein morphology and interactions. Underlying mechanisms (including thermodynamic and kinetic parameters) of the drug-target interaction and drug-modulating protein-protein interaction (PPI) on the surfaces of models or living cells are discussed. Furthermore, challenges and the outlook for the field are likewise discussed. Overall, this insight into the mechanical properties of protein-drug interactions provides an unprecedented information framework for rational drug discovery in the pharmaceutical field.

Emulating Titin by a Multidomain DNA Structure

Titin, a giant protein containing multiple tandem domains, is essential in maintaining the superior mechanical performance of muscle. The consecutive and reversible unfolding and refolding of the domains are crucial for titin to serve as a modular spring. Since the discovery of the mechanical features of a single titin molecule, the exploration of biomimetic materials with titin-emulating modular structures has been an active field. However, it remains a challenge to prepare these modular polymers on a large scale due to the complex synthesis process. In this study, we propose modular DNA with multiple hairpins (MH-DNA) as the fundamental block for the bottom-up design of advanced materials. By analyzing the unfolding and refolding dynamics of modular hairpins by atomic force microscopy (AFM)-based single-molecule force spectroscopy (SMFS), we find that MH-DNA shows comparable stability to those of polyproteins like titin. The unique low hysteresis of modular hairpin makes it an ideal molecular spring with remarkable mechanical efficiency. On the basis of the well-established DNA synthesis techniques, we anticipate that MH-DNA can be used as a promising building block for advanced materials with a combination of superior structural stability, considerable extensibility, and high mechanical efficiency.

Unfolding and identification of membrane proteins in situ

Single-molecule force spectroscopy (SMFS) uses the cantilever tip of an AFM to apply a force able to unfold a single protein. The obtained force-distance curve encodes the unfolding pathway, and from its analysis it is possible to characterize the folded domains. SMFS has been mostly used to study the unfolding of purified proteins, in solution or reconstituted in a lipid bilayer. Here, we describe a pipeline for analyzing membrane proteins based on SMFS, that involves the isolation of the plasma membrane of single cells and the harvesting of force-distance curves directly from it. We characterized and identified the embedded membrane proteins combining, within a Bayesian framework, the information of the shape of the obtained curves, with the information from Mass Spectrometry and proteomic databases. The pipeline was tested with purified/reconstituted proteins and applied to five cell types where we classified the unfolding of their most abundant membrane proteins. We validated our pipeline by overexpressing 4 constructs, and this allowed us to gather structural insights of the identified proteins, revealing variable elements in the loop regions. Our results set the basis for the investigation of the unfolding of membrane proteins in situ, and for performing proteomics from a membrane fragment.

Computed Three-Dimensional Atomic Force Microscopy Images of Biopolymers Using the Jarzynski Equality

Three-dimensional atomic force microscopy (3D-AFM) has resolved three-dimensional distributions of solvent molecules at solid-liquid interfaces at the subnanometer scale. This method is now being extended to the imaging of biopolymer assemblies such as chromosomes or proteins in cells, with the expectation of being able to resolve their three-dimensional structures. Here, we have developed a computational method to simulate 3D-AFM images of biopolymers by using the Jarzynski equality. It is found that some parts of the fiber structure of biopolymers are indeed resolved in the 3D-AFM image. The dependency of 3D-AFM images on the vertical scanning velocity is investigated, and optimum scanning velocities are found. It is also clarified that forces in nonequilibrium processes are measured in 3D-AFM measurements when the dynamics of polymers are slower than the scanning of the probe.

SARS-CoV-2 Delta Variant Decreases Nanobody Binding and ACE2 Blocking Effectivity

The Delta variant spreads more rapidly than previous variants of SARS-CoV-2. This variant comprises several mutations on the receptor-binding domain (RBD_Delta) of its spike glycoprotein, which binds to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. The RBD-PD interaction has been targeted by antibodies and nanobodies to prevent viral infection, but their effectiveness against the Delta variant remains unclear. Here, we investigated RBD_Delta-PD interactions in the presence and absence of nanobodies H11-H4, H11-D4, and Ty1 by performing 21.8 μs of all-atom molecular dynamics simulations. Unbiased simulations revealed that Delta variant mutations strengthen RBD binding to ACE2 by increasing the hydrophobic interactions and salt bridge formation, but weaken interactions with H11-H4, H11-D4, and Ty1. Among these nanobodies H11-H4 and H11-D4 bind RBD without overlapping ACE2. They were unable to dislocate ACE2 from RBD_Delta when bound side by side with ACE2 on RBD. Steered molecular dynamics simulations at comparable loading rates to high-speed atomic force microscopy (AFM) experiments estimated lower rupture forces of the nanobodies from RBD_Delta compared to ACE2. Our results suggest that existing nanobodies are less effective to inhibit RBD_Delta-PD interactions and a new generation of nanobodies is needed to neutralize the Delta variant.

Elastic versus brittle mechanical responses predicted for dimeric cadherin complexes

Cadherins are a superfamily of adhesion proteins involved in a variety of biological processes that include the formation of intercellular contacts, the maintenance of tissue integrity, and the development of neuronal circuits. These transmembrane proteins are characterized by ectodomains composed of a variable number of extracellular cadherin (EC) repeats that are similar but not identical in sequence and fold. E-cadherin, along with desmoglein and desmocollin proteins, are three classical-type cadherins that have slightly curved ectodomains and engage in homophilic and heterophilic interactions through an exchange of conserved tryptophan residues in their N-terminal EC1 repeat. In contrast, clustered protocadherins are straighter than classical cadherins and interact through an antiparallel homophilic binding interface that involves overlapped EC1 to EC4 repeats. Here we present molecular dynamics simulations that model the adhesive domains of these cadherins using available crystal structures, with systems encompassing up to 2.8 million atoms. Simulations of complete classical cadherin ectodomain dimers predict a two-phased elastic response to force in which these complexes first softly unbend and then stiffen to unbind without unfolding. Simulated α, β, and γ clustered protocadherin homodimers lack a two-phased elastic response, are brittle and stiffer than classical cadherins and exhibit complex unbinding pathways that in some cases involve transient intermediates. We propose that these distinct mechanical responses are important for function, with classical cadherin ectodomains acting as molecular shock absorbers and with stiffer clustered protocadherin ectodomains facilitating overlap that favors binding specificity over mechanical resilience. Overall, our simulations provide insights into the molecular mechanics of single cadherin dimers relevant in the formation of cellular junctions essential for tissue function.

Collective mechanical responses of cadherin-based adhesive junctions as predicted by simulations

Cadherin-based adherens junctions and desmosomes help stabilize cell-cell contacts with additional function in mechano-signaling, while clustered protocadherin junctions are responsible for directing neuronal circuits assembly. Structural models for adherens junctions formed by epithelial cadherin (CDH1) proteins indicate that their long, curved ectodomains arrange to form a periodic, two-dimensional lattice stabilized by tip-to-tip trans interactions (across junction) and lateral cis contacts. Less is known about the exact architecture of desmosomes, but desmoglein (DSG) and desmocollin (DSC) cadherin proteins are also thought to form ordered junctions. In contrast, clustered protocadherin (PCDH)-based cell-cell contacts in neuronal tissues are thought to be responsible for self-recognition and avoidance, and structural models for clustered PCDH junctions show a linear arrangement in which their long and straight ectodomains form antiparallel overlapped trans complexes. Here, we report all-atom molecular dynamics simulations testing the mechanics of minimalistic adhesive junctions formed by CDH1, DSG2 coupled to DSC1, and PCDHγB4, with systems encompassing up to 3.7 million atoms. Simulations generally predict a favored shearing pathway for the adherens junction model and a two-phased elastic response to tensile forces for the adhesive adherens junction and the desmosome models. Complexes within these junctions first unbend at low tensile force and then become stiff to unbind without unfolding. However, cis interactions in both the CDH1 and DSG2-DSC1 systems dictate varied mechanical responses of individual dimers within the junctions. Conversely, the clustered protocadherin PCDHγB4 junction lacks a distinct two-phased elastic response. Instead, applied tensile force strains trans interactions directly, as there is little unbending of monomers within the junction. Transient intermediates, influenced by new cis interactions, are observed after the main rupture event. We suggest that these collective, complex mechanical responses mediated by cis contacts facilitate distinct functions in robust cell-cell adhesion for classical cadherins and in self-avoidance signaling for clustered PCDHs.

High-resolution Single-molecule Magnetic Tweezers

Single-molecule magnetic tweezers deliver magnetic force and torque to single target molecules, permitting the study of dynamic changes in biomolecular structures and their interactions. Because the magnetic tweezer setups can generate magnetic fields that vary slowly over tens of millimeters-far larger than the nanometer scale of the single molecule events being observed-this technique can maintain essentially constant force levels during biochemical experiments while generating a biologically meaningful force on the order of 1-100 pN. When using bead-tether constructs to pull on single molecules, smaller magnetic beads and shorter submicrometer tethers improve dynamic response times and measurement precision. In addition, employing high-speed cameras, stronger light sources, and a graphics programming unit permits true high-resolution single-molecule magnetic tweezers that can track nanometer changes in target molecules on a millisecond or even submillisecond time scale. The unique force-clamping capacity of the magnetic tweezer technique provides a way to conduct measurements under near-equilibrium conditions and directly map the energy landscapes underlying various molecular phenomena. High-resolution single-molecule magnetic tweezers can thus be used to monitor crucial conformational changes in single-protein molecules, including those involved in mechanotransduction and protein folding. Expected final online publication date for the Annual Review of Biochemistry, Volume 91 is June 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

Nano-Precision Tweezers for Mechanosensitive Proteins and Beyond

Mechanical forces play pivotal roles in regulating cell shape, function, and fate. Key players that govern the mechanobiological interplay are the mechanosensitive proteins found on cell membranes and in cytoskeleton. Their unique nanomechanics can be interrogated using single-molecule tweezers, which can apply controlled forces to the proteins and simultaneously measure the ensuing structural changes. Breakthroughs in high-resolution tweezers have enabled the routine monitoring of nanometer-scale, millisecond dynamics as a function of force. Undoubtedly, the advancement of structural biology will be further fueled by integrating static atomic-resolution structures and their dynamic changes and interactions observed with the force application techniques. In this minireview, we will introduce the general principles of single-molecule tweezers and their recent applications to the studies of force-bearing proteins, including the synaptic proteins that need to be categorized as mechanosensitive in a broad sense. We anticipate that the impact of nano-precision approaches in mechanobiology research will continue to grow in the future.

Angle-dependent strength of a single chemical bond by stereographic force spectroscopy

A wealth of chemical bonds and polymers have been studied with single-molecule force spectroscopy, usually by applying a force perpendicular to the anchoring surface. However, the direction-dependence of the bond strength lacks fundamental understanding. Here we establish stereographic force spectroscopy to study the single-bond strength for various pulling angles. Surprisingly, we find that the apparent bond strength increases with increasing pulling angle relative to the anchoring surface normal, indicating a sturdy mechanical anisotropy of a chemical bond. This finding can be rationalized by a fixed pathway for the rupture of the bond, resulting in an effective projection of the applied pulling force onto a nearly fixed rupture direction. Our study is fundamental for the molecular understanding of the role of the direction of force application in molecular adhesion and friction. It is also a prerequisite for the nanoscale tailoring of the anisotropic strength of bottom-up designed materials.

Stereographic force spectroscopy reveals that a chemical bond ruptures along a fixed pathway such that the apparent bond strength strongly depends on the angle of force application.

Massively Parallelized Molecular Force Manipulation with On-Demand Thermal and Optical Control

In single-molecule force spectroscopy (SMFS), a tethered molecule is stretched using a specialized instrument to study how macromolecules extend under force. One problem in SMFS is the serial and slow nature of the measurements, performed one molecule at a time. To address this long-standing challenge, we report on the origami polymer force clamp (OPFC) which enables parallelized manipulation of the mechanical forces experienced by molecules without the need for dedicated SMFS instruments or surface tethering. The OPFC positions target molecules between a rigid nanoscale DNA origami beam and a responsive polymer particle that shrinks on demand. As a proof-of-concept, we record the steady state and time-resolved mechanical unfolding dynamics of DNA hairpins using the fluorescence signal from ensembles of molecules and confirm our conclusion using modeling.

Structure and Mechanical Stabilities of the Three-Way Junction Motifs in Prohead RNA

The core structure of phi29 prohead RNA (pRNA) is composed of three major helices organized into three-way junction pRNA (3WJ-pRNA) and has stout structural rigidity along the coaxial helices. Prohead RNAs of the other Bacillus subtilis bacteriophages such as GA1 and SF5 share similar secondary structure and function with phi29; whether these pRNAs have similar mechanical rigidity remains to be elucidated. In this study, we constructed the tertiary structures of GA1 and SF5 3WJ-pRNAs by comparative modeling. Both GA1 and SF5 3WJ-pRNAs adopt a similar structure, in which three helices are organized as the three-way junction and two of the three helices are stacked coaxially. Moreover, detailed structural features of GA1 and SF5 3WJ-pRNAs are also similar to those of phi29 3WJ-pRNA: all of the bases of the coaxial helices are paired, and all of the adenines in the junction region are paired, which eliminates the interference of A-minor tertiary interactions. Hence, the coaxial helices tightly join to each other, and the major groove between them is very narrow. Two Mg²⁺ ions can thus fit into this major groove and form double Mg clamps. A steered molecular dynamics simulation was used to study the mechanical properties of these 3WJ-pRNAs. Both GA1 and SF5 3WJ-pRNAs show strong resistance to applied force in the direction of their coaxial helices. Such mechanical stability can be attributed to the Mg clamps.

Binding Mechanism of Neutralizing Nanobodies Targeting SARS-CoV-2 Spike Glycoprotein

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) enters human cells upon binding of its spike (S) glycoproteins to ACE2 receptors. Several nanobodies neutralize SARS-CoV-2 infection by binding to the receptor-binding domain (RBD) of the S protein, but how their binding antagonizes S-ACE2 interactions is not well understood. Here, we identified interactions between the RBD and nanobodies H11-H4, H11-D4, and Ty1 by performing all-atom molecular dynamics simulations. H11-H4 and H11-D4 can bind to RBD without overlapping with ACE2. H11-H4, and to a lesser extent H11-D4, binding dislocates ACE2 from its binding site due to electrostatic repulsion. In comparison, Ty1 overlaps with ACE2 on RBD and has a similar binding strength to ACE2. Mutations in the Alpha variant of SARS-CoV-2 had a minor effect in RBD binding strengths of ACE2 and nanobodies, but reduced the ability of H11-H4 and H11-D4 to dislocate ACE2 from RBD. In comparison, the Beta variant weakened the RBD binding strengths of H11-H4 and H11-D4, which were less effective to dislocate ACE2 binding. Unexpectedly, mutations in Beta strengthened Ty1 binding to RBD, suggesting that this nanobody may be more effective to neutralize the Beta variant of SARS-CoV-2.

Interpretation of Single-Molecule Force Experiments on Proteins Using Normal Mode Analysis

Single-molecule force spectroscopy experiments allow protein folding and unfolding to be explored using mechanical force. Probably the most informative technique for interpreting the results of these experiments at the structural level makes use of steered molecular dynamics (MD) simulations, which can explicitly model the protein under load. Unfortunately, this technique is computationally expensive for many of the most interesting biological molecules. Here, we find that normal mode analysis (NMA), a significantly cheaper technique from a computational perspective, allows at least some of the insights provided by MD simulation to be gathered. We apply this technique to three non-homologous proteins that were previously studied by force spectroscopy: T4 lysozyme (T4L), Hsp70 and the glucocorticoid receptor domain (GCR). The NMA results for T4L and Hsp70 are compared with steered MD simulations conducted previously, and we find that we can recover the main results. For the GCR, which did not undergo MD simulation, our approach identifies substructures that correlate with experimentally identified unfolding intermediates. Overall, we find that NMA can make a valuable addition to the analysis toolkit for the structural analysis of single-molecule force experiments on proteins.

Single-molecule force spectroscopy reveals structural differences of heparan sulfate chains during binding to vitronectin

The syndecans represent an ongoing research field focused on their regulatory roles in normal and pathological conditions. The role of syndecans in cancer progression is well documented, implicating their importance in diagnosis and even proposing various potential cancer treatments. Thus, the characterization of the unbinding properties at the single-molecule level will appeal to their use as targets for therapeutics. In our study, syndecan-1 and syndecan-4 were measured during the interaction with the vitronectin HEP II binding site. Our findings show that syndecans are calcium ion dependent molecules that reveal distinct, unbinding properties indicating the alterations in the structure of heparan sulfate (HS) chains, possibly in the chain sequence or sulfation pattern. In this way, we suppose that HS chain affinity to extracellular matrix proteins may govern cancer invasion by altering the syndecans' ability to interact with cancer-related receptors present in the tumor microenvironment, thereby promoting the activation of various signaling cascades regulating tumor cell behavior.

Biophysical reviews top five: atomic force microscopy in biophysics

Since its invention in the late 1980s, atomic force microscopy (AFM), in which a nanometer-sized tip is used to physically interrogate the properties of a surface at high resolution, has brought about scientific revolutions in both surface science and biological physics. In response to a request from the journal, I have prepared a top-five list of scientific papers that I feel represent truly landmark developments in the use of AFM in the biophysics field. This selection is necessarily limited by number (just five) and subjective (my opinion) and I offer my apologies to those not appearing in this list.

Microbial production of megadalton titin yields fibers with advantageous mechanical properties

Manmade high-performance polymers are typically non-biodegradable and derived from petroleum feedstock through energy intensive processes involving toxic solvents and byproducts. While engineered microbes have been used for renewable production of many small molecules, direct microbial synthesis of high-performance polymeric materials remains a major challenge. Here we engineer microbial production of megadalton muscle titin polymers yielding high-performance fibers that not only recapture highly desirable properties of natural titin (i.e., high damping capacity and mechanical recovery) but also exhibit high strength, toughness, and damping energy - outperforming many synthetic and natural polymers. Structural analyses and molecular modeling suggest these properties derive from unique inter-chain crystallization of folded immunoglobulin-like domains that resists inter-chain slippage while permitting intra-chain unfolding. These fibers have potential applications in areas from biomedicine to textiles, and the developed approach, coupled with the structure-function insights, promises to accelerate further innovation in microbial production of high-performance materials.

Simulation of protein pulling dynamics on second time scale with boxed molecular dynamics

We demonstrate how recently developed Boxed Molecular Dynamics (BXD) and kinetics [D. V. Shalashilin et al., J. Chem. Phys. 137, 165102 (2012)] can provide a kinetic description of protein pulling experiments, allowing for a connection to be made between experiment and the atomistic protein structure. BXD theory applied to atomic force microscopy unfolding is similar in spirit to the kinetic two-state model [A. Noy and R. W. Friddle, Methods 60, 142 (2013)] but with some differences. First, BXD uses a large number of boxes, and therefore, it is not a two-state model. Second, BXD rate coefficients are obtained from atomistic molecular dynamics simulations. BXD can describe the dependence of the pulling force on pulling speed. Similar to Shalashilin et al. [J. Chem. Phys. 137, 165102 (2012)], we show that BXD is able to model the experiment at a very long time scale up to seconds, which is way out of reach for standard molecular dynamics.

Protein nanomechanics in biological context

How proteins respond to pulling forces, or protein nanomechanics, is a key contributor to the form and function of biological systems. Indeed, the conventional view that proteins are able to diffuse in solution does not apply to the many polypeptides that are anchored to rigid supramolecular structures. These tethered proteins typically have important mechanical roles that enable cells to generate, sense, and transduce mechanical forces. To fully comprehend the interplay between mechanical forces and biology, we must understand how protein nanomechanics emerge in living matter. This endeavor is definitely challenging and only recently has it started to appear tractable. Here, I introduce the main in vitro single-molecule biophysics methods that have been instrumental to investigate protein nanomechanics over the last 2 decades. Then, I present the contemporary view on how mechanical force shapes the free energy of tethered proteins, as well as the effect of biological factors such as post-translational modifications and mutations. To illustrate the contribution of protein nanomechanics to biological function, I review current knowledge on the mechanobiology of selected muscle and cell adhesion proteins including titin, talin, and bacterial pilins. Finally, I discuss emerging methods to modulate protein nanomechanics in living matter, for instance by inducing specific mechanical loss-of-function (mLOF). By interrogating biological systems in a causative manner, these new tools can contribute to further place protein nanomechanics in a biological context.

Ultrastructure of influenza virus ribonucleoprotein complexes during viral RNA synthesis

The single-stranded, negative-sense, viral genomic RNA (vRNA) of influenza A virus is encapsidated by viral nucleoproteins (NPs) and an RNA polymerase to form a ribonucleoprotein complex (vRNP) with a helical, rod-shaped structure. The vRNP is responsible for transcription and replication of the vRNA. However, the vRNP conformation during RNA synthesis is not well understood. Here, using high-speed atomic force microscopy and cryo-electron microscopy, we investigated the native structure of influenza A vRNPs during RNA synthesis in vitro. Two distinct types of vRNPs were observed in association with newly synthesized RNAs: an intact, helical rod-shaped vRNP connected with a folded RNA and a deformed vRNP associated with a looped RNA. Interestingly, the looped RNA was a double-stranded RNA, which likely comprises a nascent RNA and the template RNA detached from NPs of the vRNP. These results suggest that while some vRNPs keep their helical structures during RNA synthesis, for the repeated cycle of RNA synthesis, others accidentally become structurally deformed, which likely results in failure to commence or continue RNA synthesis. Thus, our findings provide the ultrastructural feature of vRNPs during RNA synthesis.

Nakano et al. use high-speed AFM and cryo-EM to study influenza virus RNA synthesis in vitro. They observe the synthesis of two types of RNA: folded single stranded RNA associated to an intact viral ribonucleoprotein (vRNP) and looped, double stranded RNA associated with a partially rearranged vRNP structure.

Elucidating Protein Translocon Dynamics with Single-Molecule Precision

Translocons are protein assemblies that facilitate the targeting and transport of proteins into and across biological membranes. Our understanding of these systems has been advanced using genetics, biochemistry, and structural biology. Despite these classic advances, until recently we have still largely lacked a detailed understanding of how translocons recognize and facilitate protein translocation. With the advent and improvements of cryogenic electron microscopy (cryo-EM) single-particle analysis and single-molecule fluorescence microscopy, the details of how translocons function are finally emerging. Here, we introduce these methods and evaluate their importance in understanding translocon structure, function, and dynamics.

Critical Interactions Between the SARS-CoV-2 Spike Glycoprotein and the Human ACE2 Receptor

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infects human cells by binding its spike (S) glycoproteins to angiotensin-converting enzyme 2 (ACE2) receptors and causes the coronavirus disease 2019 (COVID-19). Therapeutic approaches to prevent SARS-CoV-2 infection are mostly focused on blocking S-ACE2 binding, but critical residues that stabilize this interaction are not well understood. By performing all-atom molecular dynamics (MD) simulations, we identified an extended network of salt bridges, hydrophobic and electrostatic interactions, and hydrogen bonds between the receptor-binding domain (RBD) of the S protein and ACE2. Mutagenesis of these residues on the RBD was not sufficient to destabilize binding but reduced the average work to unbind the S protein from ACE2. In particular, the hydrophobic end of RBD serves as the main anchor site and is the last to unbind from ACE2 under force. We propose that blocking the hydrophobic surface of RBD via neutralizing antibodies could prove to be an effective strategy to inhibit S-ACE2 interactions.

Structural dynamics of channels and transporters by high-speed atomic force microscopy

Channels and transporters are vital for transmembrane transport of ions and solutes, and also of larger compounds such as lipids and macromolecules. Therefore, they are crucial in many biological processes such as sensing, signal transduction, and the regulation of the distribution of molecules. Dysfunctions of these membrane proteins are associated to numerous diseases, and their interaction with drugs is critical in medicine. Understanding the behavior of channels and transporters requires structural and dynamic information to decipher the molecular mechanisms underlying their function. High-Speed Atomic Force Microscopy (HS-AFM) now allows the study of single transmembrane channels and transporters in action under physiological conditions, i.e., at ambient temperature and pressure, in physiological buffer and in a membrane, and in a most direct, label-free manner. In this chapter, we discuss the HS-AFM sample preparation, application, and data analysis protocols to study the structural and conformational dynamics of membrane-embedded channels and transporters.

Estimating the Directional Flexibility of Proteins from Equilibrium Thermal Fluctuations

The directional flexibility of proteins is an equilibrium molecular property which is accessible to both experiment and computation. Single molecule force spectroscopy (SMFS) experiments report effective directional spring constants to describe the collective anisotropic response of a protein structure to mechanical pulling forces applied along selected axes. On the other hand, computational methods have thus far employed either indirect force based nonequilibrium simulations or coarse-grained elastic network models (ENM) to predict protein directional spring constants. Here, we examine the ability of equilibrium atomistic Molecular Dynamics (MD) simulations to estimate the directional flexibility and mechanical anisotropy of proteins. MD-derived effective directional spring constants are found to correlate well with SMFS spring constants (ρ² = 0.97-0.99; Adj R² = 0.92-0.99) and unfolding forces (ρ² = 0.85-0.97; Adj R² = 0.63-0.91) for five different globular proteins. Specifically, the computed spring constants reproduce the mechanical anisotropy reported by SMFS along five different directions of green fluorescence protein (GFP) and six directions of the immunoglobulin-binding B1 domain of streptococcal protein G (GB1). Further, protein dynamics as captured in MD can be translated into spring constants which can distinguish the N-C directional flexibility of ubiquitin (Ub) from two structurally homologous small ubiquitin-like modifier (SUMO1 and SUMO2) isoforms. We apply our computational framework to study the mechanical anisotropy of Ub along the seven lysine-C-term directions which are functionally relevant. We show that Ub possesses two distinct flexibility scales along these directions which roughly differ by an order of magnitude. Further, our studies reveal that the mechanical anisotropy of Ub is modified in contrasting ways by the binding of two partner proteins (UBCH5A and UEV) which attach and recognize these biomolecular tag proteins. On the basis of equilibrium MD benchmarks for flexibility along 2485 bond vectors in Ub, we propose and validate a new covariance-propagation scheme to extract spring constants from ENM normal modes. We also critically examine the ability of ENM to predict directional flexibility of proteins and suggest modifications to improve these intuitive and scalable descriptions.

Free-energy changes of bacteriorhodopsin point mutants measured by single-molecule force spectroscopy

Single amino acid mutations provide quantitative insight into the energetics that underlie the dynamics and folding of membrane proteins. Chemical denaturation is the most widely used assay and yields the change in unfolding free energy (ΔΔG). It has been applied to >80 different residues of bacteriorhodopsin (bR), a model membrane protein. However, such experiments have several key limitations: 1) a nonnative lipid environment, 2) a denatured state with significant secondary structure, 3) error introduced by extrapolation to zero denaturant, and 4) the requirement of globally reversible refolding. We overcame these limitations by reversibly unfolding local regions of an individual protein with mechanical force using an atomic-force-microscope assay optimized for 2 μs time resolution and 1 pN force stability. In this assay, bR was unfolded from its native bilayer into a well-defined, stretched state. To measure ΔΔG, we introduced two alanine point mutations into an 8-amino-acid region at the C-terminal end of bR's G helix. For each, we reversibly unfolded and refolded this region hundreds of times while the rest of the protein remained folded. Our single-molecule-derived ΔΔG for mutant L223A (-2.3 ± 0.6 kcal/mol) quantitatively agreed with past chemical denaturation results while our ΔΔG for mutant V217A was 2.2-fold larger (-2.4 ± 0.6 kcal/mol). We attribute the latter result, in part, to contact between Val²¹⁷ and a natively bound squalene lipid, highlighting the contribution of membrane protein-lipid contacts not present in chemical denaturation assays. More generally, we established a platform for determining ΔΔG for a fully folded membrane protein embedded in its native bilayer.

Modulation of a protein-folding landscape revealed by AFM-based force spectroscopy notwithstanding instrumental limitations

Single-molecule force spectroscopy is a powerful tool for studying protein folding. Over the last decade, a key question has emerged: how are changes in intrinsic biomolecular dynamics altered by attachment to μm-scale force probes via flexible linkers? Here, we studied the folding/unfolding of α₃D using atomic force microscopy (AFM)-based force spectroscopy. α₃D offers an unusual opportunity as a prior single-molecule fluorescence resonance energy transfer (smFRET) study showed α₃D's configurational diffusion constant within the context of Kramers theory varies with pH. The resulting pH dependence provides a test for AFM-based force spectroscopy's ability to track intrinsic changes in protein folding dynamics. Experimentally, however, α₃D is challenging. It unfolds at low force (<15 pN) and exhibits fast-folding kinetics. We therefore used focused ion beam-modified cantilevers that combine exceptional force precision, stability, and temporal resolution to detect state occupancies as brief as 1 ms. Notably, equilibrium and nonequilibrium force spectroscopy data recapitulated the pH dependence measured using smFRET, despite differences in destabilization mechanism. We reconstructed a one-dimensional free-energy landscape from dynamic data via an inverse Weierstrass transform. At both neutral and low pH, the resulting constant-force landscapes showed minimal differences (∼0.2 to 0.5 k _B T) in transition state height. These landscapes were essentially equal to the predicted entropic barrier and symmetric. In contrast, force-dependent rates showed that the distance to the unfolding transition state increased as pH decreased and thereby contributed to the accelerated kinetics at low pH. More broadly, this precise characterization of a fast-folding, mechanically labile protein enables future AFM-based studies of subtle transitions in mechanoresponsive proteins.

Biological physics by high-speed atomic force microscopy

While many fields have contributed to biological physics, nanotechnology offers a new scale of observation. High-speed atomic force microscopy (HS-AFM) provides nanometre structural information and dynamics with subsecond resolution of biological systems. Moreover, HS-AFM allows us to measure piconewton forces within microseconds giving access to unexplored, fast biophysical processes. Thus, HS-AFM provides a tool to nourish biological physics through the observation of emergent physical phenomena in biological systems. In this review, we present an overview of the contribution of HS-AFM, both in imaging and force spectroscopy modes, to the field of biological physics. We focus on examples in which HS-AFM observations on membrane remodelling, molecular motors or the unfolding of proteins have stimulated the development of novel theories or the emergence of new concepts. We finally provide expected applications and developments of HS-AFM that we believe will continue contributing to our understanding of nature, by serving to the dialogue between biology and physics. This article is part of a discussion meeting issue 'Dynamic in situ microscopy relating structure and function'.

Combining High-Pressure Perturbation with NMR Spectroscopy for a Structural and Dynamical Characterization of Protein Folding Pathways

High-hydrostatic pressure is an alternative perturbation method that can be used to destabilize globular proteins. Generally perfectly reversible, pressure exerts local effects on regions or domains of a protein containing internal voids, contrary to heat or chemical denaturant that destabilize protein structures uniformly. When combined with NMR spectroscopy, high pressure (HP) allows one to monitor at a residue-level resolution the structural transitions occurring upon unfolding and to determine the kinetic properties of the process. The use of HP-NMR has long been hampered by technical difficulties. Owing to the recent development of commercially available high-pressure sample cells, HP-NMR experiments can now be routinely performed. This review summarizes recent advances of HP-NMR techniques for the characterization at a quasi-atomic resolution of the protein folding energy landscape.

Atomic Force Microscopy-Based Force Spectroscopy and Multiparametric Imaging of Biomolecular and Cellular Systems

During the last three decades, a series of key technological improvements turned atomic force microscopy (AFM) into a nanoscopic laboratory to directly observe and chemically characterize molecular and cell biological systems under physiological conditions. Here, we review key technological improvements that have established AFM as an analytical tool to observe and quantify native biological systems from the micro- to the nanoscale. Native biological systems include living tissues, cells, and cellular components such as single or complexed proteins, nucleic acids, lipids, or sugars. We showcase the procedures to customize nanoscopic chemical laboratories by functionalizing AFM tips and outline the advantages and limitations in applying different AFM modes to chemically image, sense, and manipulate biosystems at (sub)nanometer spatial and millisecond temporal resolution. We further discuss theoretical approaches to extract the kinetic and thermodynamic parameters of specific biomolecular interactions detected by AFM for single bonds and extend the discussion to multiple bonds. Finally, we highlight the potential of combining AFM with optical microscopy and spectroscopy to address the full complexity of biological systems and to tackle fundamental challenges in life sciences.

Structural determinants of protocadherin-15 mechanics and function in hearing and balance perception

When sound vibrations reach the inner ear, fine protein filaments called “tip links” stretch and open cochlear hair-cell mechanosensitive channels that trigger sensory perception. Similarly, vestibular hair cells use tip links to sense mechanical stimuli produced by head motions. Tip links are formed by cadherin-23 and protocadherin-15, two large proteins involved in hearing loss and balance disorders. Here we present multiple structures, models, and simulations that depict the lower end of the tip link, including the complete protocadherin-15 ectodomain. These models show an essential connection between cadherin-23 and protocadherin-15 with dual molecular “handshakes” and various protein sites that are mutated in inherited deafness. The simulations also reveal how the tip link responds to force to mediate hearing and balance sensing.

The vertebrate inner ear, responsible for hearing and balance, is able to sense minute mechanical stimuli originating from an extraordinarily broad range of sound frequencies and intensities or from head movements. Integral to these processes is the tip-link protein complex, which conveys force to open the inner-ear transduction channels that mediate sensory perception. Protocadherin-15 and cadherin-23, two atypically large cadherins with 11 and 27 extracellular cadherin (EC) repeats, are involved in deafness and balance disorders and assemble as parallel homodimers that interact to form the tip link. Here we report the X-ray crystal structure of a protocadherin-15 + cadherin-23 heterotetrameric complex at 2.9-Å resolution, depicting a parallel homodimer of protocadherin-15 EC1-3 molecules forming an antiparallel complex with two cadherin-23 EC1-2 molecules. In addition, we report structures for 10 protocadherin-15 fragments used to build complete high-resolution models of the monomeric protocadherin-15 ectodomain. Molecular dynamics simulations and validated crystal contacts are used to propose models for the complete extracellular protocadherin-15 parallel homodimer and the tip-link bond. Steered molecular dynamics simulations of these models suggest conditions in which a structurally diverse and multimodal protocadherin-15 ectodomain can act as a stiff or soft gating spring. These results reveal the structural determinants of tip-link–mediated inner-ear sensory perception and elucidate protocadherin-15’s structural and adhesive properties relevant in disease.

Lipid bilayers: Phase behavior and nanomechanics

Lipid membranes are involved in many physiological processes like recognition, signaling, fusion or remodeling of the cell membrane or some of its internal compartments. Within the cell, they are the ultimate barrier, while maintaining the fluidity or flexibility required for a myriad of processes, including membrane protein assembly. The physical properties of in vitro model membranes as model cell membranes have been extensively studied with a variety of techniques, from classical thermodynamics to advanced modern microscopies. Here we review the nanomechanics of solid-supported lipid membranes with a focus in their phase behavior. Relevant information obtained by quartz crystal microbalance with dissipation monitoring (QCM-D) and atomic force microscopy (AFM) as complementary techniques in the nano/mesoscale interface is presented. Membrane morphological and mechanical characterization will be discussed in the framework of its phase behavior, phase transitions and coexistence, in simple and complex models, and upon the presence of cholesterol.

Rate-dependent force-extension models for single-molecule force spectroscopy experiments

Single-molecule force spectroscopy techniques allow for the measurement of several static and dynamic features of macromolecules of biological origin. In particular, atomic force microscopy, used with a variable pulling rate, provides valuable information on the folding/unfolding dynamics of proteins. We propose here two different models able to describe the out-of-equilibrium statistical mechanics of a chain composed of bistable units. These latter represent the protein domains, which can be either folded or unfolded. Both models are based on the Langevin approach and their implementation allows for investigating the effect of the pulling rate and of the device intrinsic elasticity on the chain unfolding response. The theoretical results (both analytical and numerical) have been compared with experimental data concerning the unfolding of the titin and filamin proteins, eventually obtaining a good agreement over a large range of the pulling rates.