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.

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.

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.

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).

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.

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.

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.

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.

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.

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.

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.

Sequence-dependent mechanical properties of double-stranded RNA.

The mechanical properties of double-stranded RNA (dsRNA) are involved in many of its biological functions and are relevant for future nanotechnology applications. DsRNA must tightly bend to fit inside viral capsids or deform upon the interaction with proteins that regulate gene silencing or the immune response against viral attacks. However, the question of how the nucleotide sequence affects the global mechanical properties of dsRNA has so far remained largely unexplored. Here, we have employed state-of-the-art atomistic molecular dynamics simulations to unveil the mechanical response of different RNA duplexes to an external force. Our results reveal that, similarly to dsDNA, the mechanical properties of dsRNA are highly sequence-dependent. However, we find that the nucleotide sequence affects in a strikingly different manner the stretching and twisting response of RNA and DNA duplexes under force. We find that the elastic response of dsRNA is dominated by the local high flexibility of pyrimidine-purine steps. Moreover, the flexibility of pyrimidine-purine steps is independent of the sequence context, and the global flexibility of the duplex reasonably scales with the number of this kind of base-pair dinucleotides. We conclude that disparities of the mechanical response of dinucleotides are responsible for the differences observed in the mechanical properties of RNA and DNA duplexes.

Tuning Structure and Dynamics of Blue Copper Azurin Junctions via Single Amino-Acid Mutations.

In the growing field of biomolecular electronics, blue-copper Azurin stands out as one of the most widely studied protein in single-molecule contacts. Interestingly, despite the paramount importance of the structure/dynamics of molecular contacts in their transport properties, these factors remain largely unexplored from the theoretical point of view in the context of single Azurin junctions. Here we address this issue using all-atom Molecular Dynamics (MD) of Pseudomonas Aeruginosa Azurin adsorbed to a Au(111) substrate. In particular, we focus on the structure and dynamics of the free/adsorbed protein and how these properties are altered upon single-point mutations. The results revealed that wild-type Azurin adsorbs on Au(111) along two well defined configurations: one tethered via cysteine groups and the other via the hydrophobic pocket surrounding the Cu 2 + . Surprisingly, our simulations revealed that single amino-acid mutations gave rise to a quenching of protein vibrations ultimately resulting in its overall stiffening. Given the role of amino-acid vibrations and reorientation in the dehydration process at the protein-water-substrate interface, we suggest that this might have an effect on the adsorption process of the mutant, giving rise to new adsorption configurations.

Mechanical Deformation and Electronic Structure of a Blue Copper Azurin in a Solid-State Junction.

Protein-based electronics is an emerging field which has attracted considerable attention over the past decade. Here, we present a theoretical study of the formation and electronic structure of a metal-protein-metal junction based on the blue-copper azurin from pseudomonas aeruginosa. We focus on the case in which the protein is adsorbed on a gold surface and is contacted, at the opposite side, to an STM (Scanning Tunneling Microscopy) tip by spontaneous attachment. This has been simulated through a combination of molecular dynamics and density functional theory. We find that the attachment to the tip induces structural changes in the protein which, however, do not affect the overall electronic properties of the protein. Indeed, only changes in certain residues are observed, whereas the electronic structure of the Cu-centered complex remains unaltered, as does the total density of states of the whole protein.

Assessing the Accuracy of Different Solvation Models To Describe Protein Adsorption.

In protein adsorption, the surrounding solvent has an important role in mediating protein-surface interactions. Therefore, it is of paramount importance that the solvent methods employed to model these kinds of processes are able to correctly capture the complex mechanisms occurring in the protein-water-surface interface. Here, we test the suitability of the two most popular implicit solvent methods based on the Generalized Born formalism to describe the adsorption process of the immunoglobulin G (IgG) on a hydrophobic graphene surface. Our results show that in both cases, IgG experiences an extreme and early (in less than 40 ns) unfolding as a result of the adsorption to the surface in contrast with previous experimental findings. A detailed energy decomposition analysis of explicit and implicit solvent simulations reveals that this discrepancy arises from the ill-characterization of two energy components in implicit solvent methods. These findings help to elucidate how implicit solvent models may be improved to accurately characterize the protein adsorption process.

Conformations and cryo-force spectroscopy of spray-deposited single-strand DNA on gold.

Cryo-electron microscopy can determine the structure of biological matter in vitrified liquids. However, structure alone is insufficient to understand the function of native and engineered biomolecules. So far, their mechanical properties have mainly been probed at room temperature using tens of pico-newton forces with a resolution limited by thermal fluctuations. Here we combine force spectroscopy and computer simulations in cryogenic conditions to quantify adhesion and intra-molecular properties of spray-deposited single-strand DNA oligomers on Au(111). Sub-nanometer resolution images reveal folding conformations confirmed by simulations. Lifting shows a decay of the measured stiffness with sharp dips every 0.2-0.3 nm associated with the sequential peeling and detachment of single nucleotides. A stiffness of 30-35 N m-1 per stretched repeat unit is deduced in the nano-newton range. This combined study suggests how to better control cryo-force spectroscopy of adsorbed heterogeneous (bio)polymer and to potentially enable single-base recognition in DNA strands only few nanometers long.

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.