Moiré-Tile Manipulation-Induced Friction Switch of Graphene on a Platinum Surface.

Friction control and technological advancement are intimately intertwined. Concomitantly, two-dimensional materials occupy a unique position for realizing quasi-frictionless contacts. However, the question arises of how to tune superlubric sliding. Drawing inspiration from twistronics, we propose to control superlubricity via moiré patterning. Friction force microscopy and molecular dynamics simulations unequivocally demonstrate a transition from a superlubric to dissipative sliding regime for different twist angles of graphene moirés on a Pt(111) surface triggered by the normal force. This follows from a novel mechanism at superlattice level where, beyond a critical load, moiré tiles are manipulated in a highly dissipative shear process connected to the twist angle. Importantly, the atomic detail of the dissipation associated with the moiré tile manipulation─i.e., enduring forced registry beyond a critical normal load─allows the bridging of disparate sliding regimes in a reversible manner, thus paving the road for a subtly intrinsic control of superlubricity.

A molecular view of DNA flexibility.

DNA dynamics can only be understood by taking into account its complex mechanical behavior at different length scales. At the micrometer level, the mechanical properties of single DNA molecules have been well-characterized by polymer models and are commonly quantified by a persistence length of 50 nm (~150 bp). However, at the base pair level (~3.4 Å), the dynamics of DNA involves complex molecular mechanisms that are still being deciphered. Here, we review recent single-molecule experiments and molecular dynamics simulations that are providing novel insights into DNA mechanics from such a molecular perspective. We first discuss recent findings on sequence-dependent DNA mechanical properties, including sequences that resist mechanical stress and sequences that can accommodate strong deformations. We then comment on the intricate effects of cytosine methylation and DNA mismatches on DNA mechanics. Finally, we review recently reported differences in the mechanical properties of DNA and double-stranded RNA, the other double-helical carrier of genetic information. A thorough examination of the recent single-molecule literature permits establishing a set of general 'rules' that reasonably explain the mechanics of nucleic acids at the base pair level. These simple rules offer an improved description of certain biological systems and might serve as valuable guidelines for future design of DNA and RNA nanostructures.

Flexible Superlubricity Unveiled in Sidewinding Motion of Individual Polymeric Chains.

A combination of low temperature atomic force microcopy and molecular dynamic simulations is used to demonstrate that soft designer molecules realize a sidewinding motion when dragged over a gold surface. Exploiting their longitudinal flexibility, pyrenylene chains are indeed able to lower diffusion energy barriers via on-surface directional locking and molecular strain. The resulting ultralow friction reaches values on the order of tens of pN reported so far only for rigid chains sliding on an incommensurate surface. Therefore, we demonstrate how molecular flexibility can be harnessed to realize complex nanomotion while retaining a superlubric character. This is in contrast with the paradigm guiding the design of most superlubric nanocontacts (mismatched rigid contacting surfaces).

Double-stranded RNA bending by AU-tract sequences.

Sequence-dependent structural deformations of the DNA double helix (dsDNA) have been extensively studied, where adenine tracts (A-tracts) provide a striking example for global bending in the molecule. However, in contrast to dsDNA, sequence-dependent structural features of dsRNA have received little attention. In this work, we demonstrate that the nucleotide sequence can induce a bend in a canonical Watson-Crick base-paired dsRNA helix. Using all-atom molecular dynamics simulations, we identified a sequence motif consisting of alternating adenines and uracils, or AU-tracts, that strongly bend the RNA double-helix. This finding was experimentally validated using atomic force microscopy imaging of dsRNA molecules designed to display macroscopic curvature via repetitions of phased AU-tract motifs. At the atomic level, this novel phenomenon originates from a localized compression of the dsRNA major groove and a large propeller twist at the position of the AU-tract. Moreover, the magnitude of the bending can be modulated by changing the length of the AU-tract. Altogether, our results demonstrate the possibility of modifying the dsRNA curvature by means of its nucleotide sequence, which may be exploited in the emerging field of RNA nanotechnology and might also constitute a natural mechanism for proteins to achieve recognition of specific dsRNA sequences.

Understanding the paradoxical mechanical response of in-phase A-tracts at different force regimes.

A-tracts are A:T rich DNA sequences that exhibit unique structural and mechanical properties associated with several functions in vivo. The crystallographic structure of A-tracts has been well characterized. However, the mechanical properties of these sequences is controversial and their response to force remains unexplored. Here, we rationalize the mechanical properties of in-phase A-tracts present in the Caenorhabditis elegans genome over a wide range of external forces, using single-molecule experiments and theoretical polymer models. Atomic Force Microscopy imaging shows that A-tracts induce long-range (∼200 nm) bending, which originates from an intrinsically bent structure rather than from larger bending flexibility. These data are well described with a theoretical model based on the worm-like chain model that includes intrinsic bending. Magnetic tweezers experiments show that the mechanical response of A-tracts and arbitrary DNA sequences have a similar dependence with monovalent salt supporting that the observed A-tract bend is intrinsic to the sequence. Optical tweezers experiments reveal a high stretch modulus of the A-tract sequences in the enthalpic regime. Our work rationalizes the complex multiscale flexibility of A-tracts, providing a physical basis for the versatile character of these sequences inside the cell.

Sequential Bending and Twisting around C-C Single Bonds by Mechanical Lifting of a Pre-Adsorbed Polymer.

Bending and twisting around carbon-carbon single bonds are ubiquitous in natural and synthetic polymers. Force-induced changes were so far not measured at the single-monomer level, owing to limited ways to apply local forces. We quantified down to the submolecular level the mechanical response within individual poly-pyrenylene chains upon their detachment from a gold surface with an atomic force microscope at 5 K. Computer simulations based on a dedicated force field reproduce the experimental traces and reveal symmetry-broken bent and rotated conformations of the sliding physisorbed segment besides steric hindrance of the just lifted monomer. Our study also shows that the tip-molecule bond remains intact but remarkably soft and links force variations to complex but well-defined conformational changes.

Understanding the mechanical response of double-stranded DNA and RNA under constant stretching forces using all-atom molecular dynamics.

Multiple biological processes involve the stretching of nucleic acids (NAs). Stretching forces induce local changes in the molecule structure, inhibiting or promoting the binding of proteins, which ultimately affects their functionality. Understanding how a force induces changes in the structure of NAs at the atomic level is a challenge. Here, we use all-atom, microsecond-long molecular dynamics to simulate the structure of dsDNA and dsRNA subjected to stretching forces up to 20 pN. We determine all of the elastic constants of dsDNA and dsRNA and provide an explanation for three striking differences in the mechanical response of these two molecules: the threefold softer stretching constant obtained for dsRNA, the opposite twist-stretch coupling, and its nontrivial force dependence. The lower dsRNA stretching resistance is linked to its more open structure, whereas the opposite twist-stretch coupling of both molecules is due to the very different evolution of molecules' interstrand distance with the stretching force. A reduction of this distance leads to overwinding in dsDNA. In contrast, dsRNA is not able to reduce its interstrand distance and can only elongate by unwinding. Interstrand distance is directly correlated with the slide base-pair parameter and its different behavior in dsDNA and dsRNA traced down to changes in the sugar pucker angle of these NAs.

DNA Crookedness Regulates DNA Mechanical Properties at Short Length Scales.

Sequence-dependent DNA conformation and flexibility play a fundamental role in the specificity of DNA-protein interactions. Here we quantify the DNA crookedness: a sequence-dependent deformation of DNA that consists of periodic bends of the base pair centers chain. Using extensive 100 µs-long, all-atom molecular dynamics simulations, we found that DNA crookedness and its associated flexibility are bijective, which unveils a one-to-one relation between DNA structure and dynamics. This allowed us to build a predictive model to compute the stretch moduli of different DNA sequences from solely their structure. Sequences with very little crookedness show extremely high stretching stiffness and have been previously shown to form unstable nucleosomes and promote gene expression. Interestingly, the crookedness can be tailored by epigenetic modifications, known to affect gene expression. Our results rationalize the idea that the DNA sequence is not only a chemical code, but also a physical one that allows finely regulating its mechanical properties and, possibly, its 3D arrangement inside the cell.

Ab initio electronic structure calculations of entire blue copper azurins.

We present a theoretical study of the blue-copper azurin extracted from Pseudomonas aeruginosa and several of its single amino acid mutants. For the first time, we consider the whole structure of this kind of protein rather than limiting our analysis to the copper complex only. This is accomplished by combining fully ab initio calculations based on density functional theory with atomic-scale molecular dynamics simulations. Beyond the main features arising from the copper complex, our study reveals the role played by the peripheral parts of the proteins. In particular, we find that oxygen atoms belonging to carboxyl groups which are distributed all over the protein contribute to electronic states near the HOMO. The contribution of the outer regions to the electronic structure of azurins had so far been overlooked. Our results highlight the need to investigate them thoroughly; this is especially important in prospect of understanding complex processes such as the electronic transport through metal-metalloprotein-metal junctions.

Bioengineering a Single-Protein Junction.

Bioelectronics moves toward designing nanoscale electronic platforms that allow in vivo determinations. Such devices require interfacing complex biomolecular moieties as the sensing units to an electronic platform for signal transduction. Inevitably, a systematic design goes through a bottom-up understanding of the structurally related electrical signatures of the biomolecular circuit, which will ultimately lead us to tailor its electrical properties. Toward this aim, we show here the first example of bioengineered charge transport in a single-protein electrical contact. The results reveal that a single point-site mutation at the docking hydrophobic patch of a Cu-azurin causes minor structural distortion of the protein blue Cu site and a dramatic change in the charge transport regime of the single-protein contact, which goes from the classical Cu-mediated two-step transport in this system to a direct coherent tunneling. Our extensive spectroscopic studies and molecular-dynamics simulations show that the proteins' folding structures are preserved in the single-protein junction. The DFT-computed frontier orbital of the relevant protein segments suggests that the Cu center participation in each protein variant accounts for the different observed charge transport behavior. This work is a direct evidence of charge transport control in a protein backbone through external mutagenesis and a unique nanoscale platform to study structurally related biological electron transfer.

Time Strengthening of Crystal Nanocontacts.

We demonstrate how an exponentially saturating increase of the contact area between a nanoasperity and a crystal surface, occurring on time scales governed by the Arrhenius equation, is consistent with measurements of the static friction and lateral contact stiffness on a model alkali-halide surface at different temperatures in ultrahigh vacuum. The "contact ageing" effect is attributed to atomic attrition and is eventually broken by thermally activated slip of the nanoasperity on the surface. The combination of the two effects also leads to regions of strengthening and weakening in the velocity dependence of the friction, which are well-reproduced by an extended version of the Prandtl-Tomlinson model.

Albumin (BSA) adsorption onto graphite stepped surfaces.

Nanomaterials are good candidates for the design of novel components with biomedical applications. For example, nano-patterned substrates may be used to immobilize protein molecules in order to integrate them in biosensing units. Here, we perform long MD simulations (up to 200 ns) using an explicit solvent and physiological ion concentrations to characterize the adsorption of bovine serum albumin (BSA) onto a nano-patterned graphite substrate. We have studied the effect of the orientation and step size on the protein adsorption and final conformation. Our results show that the protein is stable, with small changes in the protein secondary structure that are confined to the contact area and reveal the influence of nano-structuring on the spontaneous adsorption, protein-surface binding energies, and protein mobility. Although van der Waals (vdW) interactions play a dominant role, our simulations reveal the important role played by the hydrophobic lipid-binding sites of the BSA molecule in the adsorption process. The complex structure of these sites, that incorporate residues with different hydrophobic character, and their flexibility are crucial to understand the influence of the ion concentration and protein orientation in the different steps of the adsorption process. Our study provides useful information for the molecular engineering of components that require the immobilization of biomolecules and the preservation of their biological activity.

Albumin (BSA) Adsorption over Graphene in Aqueous Environment: Influence of Orientation, Adsorption Protocol, and Solvent Treatment.

We report 150 ns explicit solvent MD simulations of the adsorption on graphene of albumin (BSA) in two orientations and using two different adsorption protocols, i.e., free and forced adsorption. Our results show that free adsorption occurs with little structural rearrangements. Even taking adsorption to an extreme, by forcing it with a 5 nN downward force applied during the initial 20 ns, we show that along a particular orientation BSA is able to preserve the structural properties of the majority of its binding sites. Furthermore, in all the cases considered in this work, the ibuprofen binding site has shown a strong resilience to structural changes. Finally, we compare these results with implicit solvent simulations and find that the latter predicts an extreme protein unfolding upon adsorption. The origin of this discrepancy is attributed to a poor description of the water entropic forces at interfaces in the implicit solvent methods.

Efficient Electron Hopping Transport through Azurin-Based Junctions.

We conducted a theoretical study of electron transport through junctions of the blue-copper azurin from Pseudomonas aeruginosa. We found that single-site hopping can lead to either higher or lower current values compared to fully coherent transport. This depends on the structural details of the junctions as well as the alignment of the protein orbitals. Moreover, we show how the asymmetry of the IV curves can be affected by the position of the tip in the junction and that, under specific conditions, such a hopping mechanism is consistent with a fairly low temperature dependence of the current. Finally, we show that increasing the number of hopping sites leads to higher hopping currents. Our findings, from fully quantum calculations, provide deep insight to help guide the interpretation of experimental IV curves on highly complex systems.

Accurate Quantum-Mechanically Derived Force-Fields through a Fragment-Based Approach: Balancing Specificity and Transferability in the Prediction of Self-Assembly in Soft Matter.

The wide range of time/length scales covered by self-assembly in soft matter makes molecular dynamics (MD) the ideal candidate for simulating such a supramolecular phenomenon at an atomistic level. However, the reliability of MD outcomes heavily relies on the accuracy of the adopted force-field (FF). The spontaneous re-ordering in liquid crystalline materials stands as a clear example of such collective self-assembling processes, driven by a subtle and delicate balance between supramolecular interactions and single-molecule flexibility. General-purpose transferable FFs often dramatically fail to reproduce such complex phenomena, for example, the error on the transition temperatures being larger than 100 K. Conversely, quantum-mechanically derived force-fields (QMD-FFs), specifically tailored for the target system, were recently shown (J. Phys. Chem. Lett.2022,13, 243) to allow for the required accuracy as they not only well reproduced transition temperatures but also yielded a quantitative agreement with the experiment on a wealth of structural, dynamic, and thermodynamic properties. The main drawback of this strategy stands in the computational burden connected to the numerous quantum mechanical (QM) calculations usually required for a target-specific parameterization, which has undoubtedly hampered the routine application of QMD-FFs. In this work, we propose a fragment-based strategy to extend the applicability of QMD-FFs, in which the amount of QM calculations is significantly reduced, being a single-molecule-optimized geometry and its Hessian matrix the only QM information required. To validate this route, a new FF is assembled for a large mesogen, exploiting the parameters obtained for two smaller liquid crystalline molecules, in this and previous work. Lengthy MD simulations are carried out with the new transferred QMD-FF, observing again a spontaneous re-orientation in the correct range of temperatures, with good agreement with the available experimental measures. The present results strongly suggest that a partial transfer of QMD-FF parameters can be invoked without a significant loss of accuracy, thus paving the way to exploit the method's intrinsic predictive capabilities in the simulation of novel soft materials.

Predicting Spontaneous Orientational Self-Assembly: In Silico Design of Materials with Quantum Mechanically Derived Force Fields.

De novo design of self-assembled materials hinges upon our ability to relate macroscopic properties to individual building blocks, thus characterizing in such supramolecular architectures a wide range of observables at varied time/length scales. This work demonstrates that quantum mechanical derived force fields (QMD-FFs) do satisfy this requisite and, most importantly, do so in a predictive manner. To this end, a specific FF, built solely based on the knowledge of the target molecular structure, is employed to reproduce the spontaneous transition to an ordered liquid crystal phase. The simulations deliver a multiscale portrait of such self-assembly processes, where conformational changes within the individual building blocks are intertwined with a 200 ns ensemble reorganization. The extensive characterization provided not only is in quantitative agreement with the experiment but also connects the time/length scales at which it was performed. Realizing QMD-FF predictive power and unmatched accuracy stands as an important leap forward for the bottom-up design of advanced supramolecular materials.

Practical Guide to Single-Protein AFM Nanomechanical Spectroscopy Mapping: Insights and Pitfalls As Unraveled by All-Atom MD Simulations on Immunoglobulin G.

Atomic force microscopy is an invaluable characterization tool in almost every biophysics laboratory. However, obtaining atomic/sub-nanometer resolution on single proteins has thus far remained elusive-a feat long achieved on hard substrates. In this regard, nanomechanical spectroscopy mapping may provide a viable approach to overcome this limitation. By complementing topography with mechanical properties measured locally, one may thus enhance spatial resolution at the single-protein level. In this work, we perform all-atom molecular dynamics simulations of the indentation process on a single immunoglobulin G (IgG) adsorbed on a graphene slab. Our simulations reveal three different stages as a function of strain: a noncontact regime-where the mechanical response is linked to the presence of the water environment- followed by an elastic response and a final plastic deformation regime. In the noncontact regime, we are able to identify hydrophobic/hydrophilic patches over the protein. This regime provides the most local mechanical information that allows one to discern different regions with similar height/topography and leads to the best spatial resolution. In the elastic regime, we conclude that the Young modulus is a well-defined property only within mechanically decoupled domains. This is caused by the fact that the elastic deformation is associated with a global reorganization of the domain. Differences in the mechanical response are large enough to clearly resolve domains within a single protein, such as the three subunits forming the IgG. Two events, unfolding or protein slipping, are observed in the plastic regime. Our simulations allow us to characterize these two processes and to provide a strategy to identify them in the force curves. Finally, we elaborate on possible challenges that could hamper the interpretation of such experiments/simulations and how to overcome them. All in all, our simulations provide a detailed picture of nanomechanical spectroscopy mapping on single proteins, showing its potential and the challenges that need to be overcome to unlock its full potential.

Automated Parameterization of Quantum Mechanically Derived Force Fields for Soft Materials and Complex Fluids: Development and Validation.

The reliability of molecular dynamics (MD) simulations in predicting macroscopic properties of complex fluids and soft materials, such as liquid crystals, colloidal suspensions, or polymers, relies on the accuracy of the adopted force field (FF). We present an automated protocol to derive specific and accurate FFs, fully based on ab initio quantum mechanical (QM) data. The integration of the Joyce and Picky procedures, recently proposed by our group to provide an accurate description of simple liquids, is here extended to larger molecules, capable of exhibiting more complex fluid phases. While the standard Joyce protocol is employed to parameterize the intramolecular FF term, a new automated procedure is here proposed to handle the computational cost of the QM calculations required for the parameterization of the intermolecular FF term. The latter is thus obtained by integrating the old Picky procedure with a fragmentation reconstruction method (FRM) that allows for a reliable, yet computationally feasible sampling of the intermolecular energy surface at the QM level. The whole FF parameterization protocol is tested on a benchmark liquid crystal, and the performances of the resulting quantum mechanically derived (QMD) FF were compared with those delivered by a general-purpose, transferable one, and by the third, "hybrid" FF, where only the bonded terms were refined against QM data. Lengthy atomistic MD simulations are carried out with each FF on extended 5CB systems in both isotropic and nematic phases, eventually validating the proposed protocol by comparing the resulting macroscopic properties with other computational models and with experiments. The QMD-FF yields the best performances, reproducing both phases in the correct range of temperatures and well describing their structure, dynamics, and thermodynamic properties, thus providing a clear protocol that may be explored to predict such properties on other complex fluids or soft materials.

Giant thermal expansion of a two-dimensional supramolecular network triggered by alkyl chain motion.

Thermal expansion, the response in shape, area or volume of a solid with heat, is usually large in molecular materials compared to their inorganic counterparts. Resulting from the intrinsic molecule flexibility, conformational changes or variable intermolecular interactions, the exact interplay between these mechanisms is however poorly understood down to the molecular level. Here, we investigate the structural variations of a two-dimensional supramolecular network on Au(111) consisting of shape persistent polyphenylene molecules equipped with peripheral dodecyl chains. By comparing high-resolution scanning probe microscopy and molecular dynamics simulations obtained at 5 and 300 K, we determine the thermal expansion coefficient of the assembly of 980 ± 110 × 10-6 K-1, twice larger than other molecular systems hitherto reported in the literature, and two orders of magnitude larger than conventional materials. This giant positive expansion originates from the increased mobility of the dodecyl chains with temperature that determine the intermolecular interactions and the network spacing.