Nucleosome Array Deformation in Chromatin is Sustained by Bending, Twisting and Kinking of Linker DNA.

Chromatin in the nucleus undergoes mechanical stresses from different sources during the various stages of cell life. Here a trinucleosome array is used as the minimal model to study the mechanical response to applied stress at the molecular level. By using large-scale, all-atom steered-molecular dynamics simulations, we show that the largest part of mechanical stress in compression is accommodated by the DNA linkers joining pairs of nucleosomes, which store the elastic energy accumulated by the applied force. Different mechanical instabilities (Euler bending, Brazier kinking, twist-bending) can deform the DNA canonical structure, as a function of the increasing force load. An important role of the histone tails in assisting the DNA deformation is highlighted. The overall response of the smallest chromatin fragment to compressive stress leaves the nucleosome assembly with a substantial plastic deformation and localised defects, which can have a potential impact on DNA transcription, downstream signaling pathways, the regulation of gene expression, and DNA repair.

Single-Cell Electrochemical Aptasensor Array.

Despite several demonstrations of electrochemical devices with limits of detection (LOD) of 1 cell/mL, the implementation of single-cell bioelectrochemical sensor arrays has remained elusive due to the challenges of scaling up. In this study, we show that the recently introduced nanopillar array technology combined with redox-labeled aptamers targeting epithelial cell adhesion molecule (EpCAM) is perfectly suited for such implementation. Combining nanopillar arrays with microwells determined for single cell trapping directly on the sensor surface, single target cells are successfully detected and analyzed. This first implementation of a single-cell electrochemical aptasensor array, based on Brownian-fluctuating redox species, opens new opportunities for large-scale implementation and statistical analysis of early cancer diagnosis and cancer therapy in clinical settings.

Infinite ergodicity in generalized geometric Brownian motions with nonlinear drift.

Geometric Brownian motion is an exemplary stochastic processes obeying multiplicative noise, with widespread applications in several fields, e.g., in finance, in physics, and biology. The definition of the process depends crucially on the interpretation of the stochastic integrals which involves the discretization parameter α with 0≤α≤1, giving rise to the well-known special cases α=0 (Itô), α=1/2 (Fisk-Stratonovich), and α=1 (Hänggi-Klimontovich or anti-Itô). In this paper we study the asymptotic limits of the probability distribution functions of geometric Brownian motion and some related generalizations. We establish the conditions for the existence of normalizable asymptotic distributions depending on the discretization parameter α. Using the infinite ergodicity approach, recently applied to stochastic processes with multiplicative noise by E. Barkai and collaborators, we show how meaningful asymptotic results can be formulated in a transparent way.

The out-of-field dose in radiation therapy induces delayed tumorigenesis by senescence evasion.

A rare but severe complication of curative-intent radiation therapy is the induction of second primary cancers. These cancers preferentially develop not inside the planning target volume (PTV) but around, over several centimeters, after a latency period of 1-40 years. We show here that normal human or mouse dermal fibroblasts submitted to the out-of-field dose scattering at the margin of a PTV receiving a mimicked patient's treatment do not die but enter in a long-lived senescent state resulting from the accumulation of unrepaired DNA single-strand breaks, in the almost absence of double-strand breaks. Importantly, a few of these senescent cells systematically and spontaneously escape from the cell cycle arrest after a while to generate daughter cells harboring mutations and invasive capacities. These findings highlight single-strand break-induced senescence as the mechanism of second primary cancer initiation, with clinically relevant spatiotemporal specificities. Senescence being pharmacologically targetable, they open the avenue for second primary cancer prevention.

DNA Aptamers Block the Receptor Binding Domain at the Spike Protein of SARS-CoV-2.

DNA aptamers are versatile molecular species obtained by the folding of short single-stranded nucleotide sequences, with highly specific recognition capabilities against proteins. Here we test the ability of DNA aptamers to interact with the spike (S-)protein of the SARS-CoV-2 viral capsid. The S-protein, a trimer made up of several subdomains, develops the crucial function of recognizing the ACE2 receptors on the surface of human cells, and subsequent fusioning of the virus membrane with the host cell membrane. In order to achieve this, the S1 domain of one protomer switches between a closed conformation, in which the binding site is inaccessible to the cell receptors, and an open conformation, in which ACE2 can bind, thereby initiating the entry process of the viral genetic material in the host cell. Here we show, by means of state-of-the-art molecular simulations, that small DNA aptamers experimentally identified can recognize the S-protein of SARS-CoV-2, and characterize the details of the binding process. We find that their interaction with different subdomains of the S-protein can effectively block, or at least considerably slow down the opening process of the S1 domain, thereby significantly reducing the probability of virus-cell binding. We provide evidence that, as a consequence, binding of the human ACE2 receptor may be crucially affected under such conditions. Given the facility and low cost of fabrication of specific aptamers, the present findings could open the way to both an innovative viral screening technique with sub-nanomolar sensitivity, and to an effective and low impact curative strategy.

Agent-based model of multicellular tumor spheroid evolution including cell metabolism.

Computational models aiming at the spatio-temporal description of cancer evolution are a suitable framework for testing biological hypotheses from experimental data, and generating new ones. Building on our recent work (J. Theor. Biol. 389, 146 (2016)) we develop a 3D agent-based model, capable of tracking hundreds of thousands of interacting cells, over time scales ranging from seconds to years. Cell dynamics is driven by a Monte Carlo solver, incorporating partial differential equations to describe chemical pathways and the activation/repression of "genes", leading to the up- or down-regulation of specific cell markers. Each cell-agent of different kind (stem, cancer, stromal etc.) runs through its cycle, undergoes division, can exit to a dormant, senescent, necrotic state, or apoptosis, according to the inputs from its systemic network. The basic network at this stage describes glucose/oxygen/ATP cycling, and can be readily extended to cancer-cell specific markers. Eventual accumulation of chemical/radiation damage to each cell's DNA is described by a Markov chain of internal states, and by a damage-repair network, whose evolution is linked to the cell systemic network. Aimed at a direct comparison with experiments of tumorsphere growth from stem cells, the present model will allow to quantitatively study the role of transcription factors involved in the reprogramming and variable radio-resistance of simulated cancer-stem cells, evolving in a realistic computer simulation of a growing multicellular tumorsphere.

Mechanical evolution of DNA double-strand breaks in the nucleosome.

Double strand breaks (DSB) in the DNA backbone are the most lethal type of defect induced in the cell nucleus by chemical and radiation treatments of cancer. However, little is known about the outcomes of damage in nucleosomal DNA, and on its effects on damage repair. We performed microsecond-long molecular dynamics computer simulations of nucleosomes including a DSB at various sites, to characterize the early stages of the evolution of this DNA lesion. The damaged structures are studied by the essential dynamics of DNA and histones, and compared to the intact nucleosome, thus exposing key features of the interactions. All DSB configurations tend to remain compact, with only the terminal bases interacting with histone proteins. Umbrella sampling calculations show that broken DNA ends at the DSB must overcome a free-energy barrier to detach from the nucleosome core. Finally, by calculating the covariant mechanical stress, we demonstrate that the coupled bending and torsional stress can force the DSB free ends to open up straight, thus making it accessible to damage signalling proteins.

Stability of radiation-damaged DNA after multiple strand breaks.

Damage to the DNA backbone occurs from natural sources, and with exceedingly large density during radiotherapy, as typically used for cancer treatment. Here, we focus on the molecular-scale dynamics of the events immediately following the production of single- and double-strand breaks, since this early-stage evolution of the damage is crucial to determine the subsequent fate of the DNA fragment. While multiple cleavage of phosphodiester bonds is the first step, however the remaining hydrogen-bond and π-stacking interactions maintain a considerable DNA cohesion, and determine further defect evolution. We use all-atom molecular dynamics to simulate the force spectra and thermal stability of different single- and multiple-defect configurations, in a random 31 bp DNA sequence. Simulations reveal a complex dynamical behaviour of the defects, where collective bond-rearrangement phenomena dominate with respect to simple bond cleavage. Defects are stable against thermal disruption, unless very closely spaced. We establish the necessary conditions for the events ultimately leading to DNA fragmentation. Such findings impact the early stages of damage recognition and signalling by specialised proteins, also implying that the identification and counting of DSBs by different experimental methods is non-unique.

Thermal conductivity of glassy GeTe4 by first-principles molecular dynamics.

A transient thermal regime is achieved in glassy GeTe4 by first-principles molecular dynamics following the recently proposed "approach-to-equilibrium" methodology. The temporal and spatial evolution of the temperature do comply with the time-dependent solution of the heat equation. We demonstrate that the time scales required to create the hot and the cold parts of the system and observe the resulting approach to equilibrium are accessible to first-principles molecular dynamics. Such a strategy provides the thermal conductivity from the characteristic decay time. We rationalize in detail the impact on the thermal conductivity of the initial temperature difference, the equilibration duration, and the main simulation features.

On the ability of molecular dynamics simulation and continuum electrostatics to treat interfacial water molecules in protein-protein complexes.

Interfacial waters are increasingly appreciated as playing a key role in protein-protein interactions. We report on a study of the prediction of interfacial water positions by both Molecular Dynamics and explicit solvent-continuum electrostatics based on the Dipolar Poisson-Boltzmann Langevin (DPBL) model, for three test cases: (i) the barnase/barstar complex (ii) the complex between the DNase domain of colicin E2 and its cognate Im2 immunity protein and (iii) the highly unusual anti-freeze protein Maxi which contains a large number of waters in its interior. We characterize the waters at the interface and in the core of the Maxi protein by the statistics of correctly predicted positions with respect to crystallographic water positions in the PDB files as well as the dynamic measures of diffusion constants and position lifetimes. Our approach provides a methodology for the evaluation of predicted interfacial water positions through an investigation of water-mediated inter-chain contacts. While our results show satisfactory behaviour for molecular dynamics simulation, they also highlight the need for improvement of continuum methods.

Tunneling mechanism and contact mechanics of colloidal nanoparticle assemblies.

Nanoparticle assemblies with thiol-terminated alkyl chains are studied by conducting atomic force microscopy (c-AFM) regarding their use as strain gauges for touch-sensitive panels. Current-force spectroscopy is used as a characterization tool complementary to the macroscopic setup since it allows a bias to be applied to a limited number of junctions, overcoming the Coulomb blockade energy and focusing on the contact electromechanics and the transport mechanism across the ligand. First, transition voltage spectroscopy is applied with varying force to target the underlying tunneling mechanism by observing whether the transition between the ohmic and exponential current-voltage behavior is force-dependent. Secondly, current-force spectroscopy in the ohmic range below the transition voltage is performed. The current-force behavior of the AFM probe in contact with a nanoparticle multilayer is associated with the spread of force and current within the nanoparticle lattice and at the level of adjacent particles by detailed contact mechanics treatment. The result is twofold: concerning the architecture of sensors, this work is a sample case of contact electromechanics at scales ranging from the device scale down to the individual ligand molecule. Regarding transport across the molecule, the vacuum tunneling mechanism is favored over the conduction by coherent molecular states, which is a decision-making aid for the choice of ligand in applications.

Real-time mechanical characterization of DNA degradation under therapeutic X-rays and its theoretical modeling.

The killing of tumor cells by ionizing radiation beams in cancer radiotherapy is currently based on a rather empirical understanding of the basic mechanisms and effectiveness of DNA damage by radiation. By contrast, the mechanical behaviour of DNA encompassing sequence sensitivity and elastic transitions to plastic responses is much better understood. A novel approach is proposed here based on a micromechanical Silicon Nanotweezers device. This instrument allows the detailed biomechanical characterization of a DNA bundle exposed to an ionizing radiation beam delivered here by a therapeutic linear particle accelerator (LINAC). The micromechanical device endures the harsh environment of radiation beams and still retains molecular-level detection accuracy. In this study, the first real-time observation of DNA damage by ionizing radiation is demonstrated. The DNA bundle degradation is detected by the micromechanical device as a reduction of the bundle stiffness, and a theoretical model provides an interpretation of the results. These first real-time observations pave the way for both fundamental and clinical studies of DNA degradation mechanisms under ionizing radiation for improved tumor treatment.

Stochastic mechanical degradation of multi-cracked fiber bundles with elastic and viscous interactions.

The mechanics of fiber bundles has been largely investigated in order to understand their complex failure modes. Under a mechanical load, the fibers fail progressively while the load is redistributed among the unbroken fibers. The classical fiber bundle model captures the most important features of this rupture process. On the other hand, the homogenization techniques are able to evaluate the stiffness degradation of bulk solids with a given population of cracks. However, these approaches are inadequate to determine the effective response of a degraded bundle where breaks are induced by non-mechanical actions. Here, we propose a method to analyze the behavior of a fiber bundle, undergoing a random distribution of breaks, by considering the intrinsic response of the fibers and the visco-elastic interactions among them. We obtain analytical solutions for simple configurations, while the most general cases are studied by Monte Carlo simulations. We find that the degradation of the effective bundle stiffness can be described by two scaling regimes: a first exponential regime for a low density of breaks, followed by a power-law regime at increasingly higher break density. For both regimes, we find analytical effective expressions described by specific scaling exponents.

Surface-induced optimal packing of two-dimensional molecular networks.

High-density packing in organic crystals is usually associated with an increase of the coordination between molecules. Such a concept is not necessarily extended to two-dimensional molecular networks self-assembled on a solid surface, for which we demonstrate the key role of the surface in inducing the optimal packing. By a combination of scanning tunneling microscopy experiments and multiscale computer simulations, we study the phase transition between two polymorphs. We find that, contrary to intuition, the structure with the lowest packing fraction corresponds to the highest molecular coordination number, due to the competition between surface and intermolecular forces. Having the lowest free energy, this structure spreads out as the most stable polymorph over a wide range of molecular concentrations.

Scaling shift in multicracked fiber bundles.

Bundles of fibers, wires, or filaments are ubiquitous structures in both natural and artificial materials. We investigate the bundle degradation induced by an external damaging action through a theoretical model describing an assembly of parallel fibers, progressively damaged by a random population of cracks. Fibers in our model interact by means of a lateral linear coupling, thus retaining structural integrity even after substantial damage. Monte Carlo simulations of the Young's modulus degradation for increasing crack density demonstrate a remarkable scaling shift between an exponential and a power-law regime. Analytical solutions of the model confirm this behavior, and provide a thorough understanding of the underlying physics.

Structure and mechanical characterization of DNA i-motif nanowires by molecular dynamics simulation.

We studied the structure and mechanical properties of DNA i-motif nanowires by means of molecular dynamics computer simulations. We built up to 230 nm-long nanowires, based on a repeated TC5 sequence from crystallographic data, fully relaxed and equilibrated in water. The unusual C⋅C(+) stacked structure, formed by four ssDNA strands arranged in an intercalated tetramer, is here fully characterized both statically and dynamically. By applying stretching, compression, and bending deformations with the steered molecular dynamics and umbrella sampling methods, we extract the apparent Young's and bending moduli of the nanowire, as well as estimates for the tensile strength and persistence length. According to our results, the i-motif nanowire shares similarities with structural proteins, as far as its tensile stiffness, but is closer to nucleic acids and flexible proteins, as far as its bending rigidity is concerned. Furthermore, thanks to its very thin cross section, the apparent tensile toughness is close to that of a metal. Besides their yet to be clarified biological significance, i-motif nanowires may qualify as interesting candidates for nanotechnology templates, due to such outstanding mechanical properties.

Silicon Nanotweezers with a microfluidic cavity for the real time characterization of DNA damage under therapeutic radiation beams.

We report the biomechanical characterization of λ-DNA bundle exposed to a therapeutic radiation beam by silicon Nanotweezers. The micromechanical device endures the harsh environment of radiation beams, and still retains molecular-level detection accuracy. The real-time DNA bundle degradation is observed in terms of biomechanical stiffness and viscosity reduction, both in air and in solution. These results pave the way for both fundamental and clinical studies of DNA degradation mechanisms under ionizing radiation for improved tumor treatment.

Elasticity of flexible and semiflexible polymers with extensible bonds in the Gibbs and Helmholtz ensembles.

Stretching experiments on single molecules of arbitrary length opened the way for studying the statistical mechanics of small systems. In many cases in which the thermodynamic limit is not satisfied, different macroscopic boundary conditions, corresponding to different statistical mechanics ensembles, yield different force-displacement curves. We formulate analytical expressions and develop Monte Carlo simulations to quantitatively evaluate the difference between the Helmholtz and the Gibbs ensembles for a wide range of polymer models of biological relevance. We consider generalizations of the freely jointed chain and of the worm-like chain models with extensible bonds. In all cases we show that the convergence to the thermodynamic limit upon increasing contour length is described by a suitable power law and a specific scaling exponent, characteristic of each model.