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.

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.

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.

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.

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.

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.

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.

Short-time expansion of one-dimensional Fokker-Planck equations with heterogeneous diffusion.

We formulate a short-time expansion for one-dimensional Fokker-Planck equations with spatially dependent diffusion coefficients, derived from stochastic processes with Gaussian white noise, for general values of the discretization parameter 0≤α≤1 of the stochastic integral. The kernel of the Fokker-Planck equation (the propagator) can be expressed as a product of a singular and a regular term. While the singular term can be given in closed form, the regular term can be computed from a Taylor expansion whose coefficients obey simple ordinary differential equations. We illustrate the application of our approach with examples taken from statistical physics and biophysics. Furthermore, we show how our formalism allows us to define a class of stochastic equations which can be treated exactly. The convergence of the expansion cannot be guaranteed independently from the discretization parameter α.

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.

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.

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.

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.

Theory and Monte Carlo simulations for the stretching of flexible and semiflexible single polymer chains under external fields.

Recent developments of microscopic mechanical experiments allow the manipulation of individual polymer molecules in two main ways: uniform stretching by external forces and non-uniform stretching by external fields. Many results can be thereby obtained for specific kinds of polymers and specific geometries. In this work, we describe the non-uniform stretching of a single, non-branched polymer molecule by an external field (e.g., fluid in uniform motion, or uniform electric field) by a universal physical framework, which leads to general conclusions on different types of polymers. We derive analytical results both for the freely-jointed chain and the worm-like chain models based on classical statistical mechanics. Moreover, we provide a Monte Carlo numerical analysis of the mechanical properties of flexible and semiflexible polymers anchored at one end. The simulations confirm the analytical achievements, and moreover allow to study the situations where the theory cannot provide explicit and useful results. In all cases, we evaluate the average conformation of the polymer and its fluctuation statistics as a function of the chain length, bending rigidity, and field strength.

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.

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.

Exceptional plasticity of silicon nanobridges.

The plasticity of covalently bonded materials is a subject at the forefront of materials science, bearing on a wide range of technological and fundamental aspects. However, covalent materials fracture in a brittle manner when the deformation exceeds just a few per cent. It is predicted that a macroscopically brittle material like silicon can show nanoscale plasticity. Here we report the exceptional plasticity observed in silicon nanocontacts ('nanobridges') at room temperature using a special experimental setup combining a transmission electron microscope and a microelectromechanical system. When accounting for surface diffusion, we succeeded in elongating the nanocontact into a wire-like structure, with a fivefold increase in volume, up to more than twenty times the original length. Such a large plasticity was caused by the stress-assisted diffusion and the sliding of the intergranular, amorphous-like material among the nanocrystals.

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.

Atomic-scale modeling of the interaction between short polypeptides and carbon surfaces.

We performed a comparative study of the adsorption of an in vitro selected peptide on two different carbon surfaces: a flat graphene and a curved (0,15) nanotube. The sequence was selected from recent experiments, as the one giving the highest carbon affinity for carbon nanotubes. Rigid docking of the molecule on the two surfaces by a genetic algorithm was followed by molecular dynamics simulations with empirical force fields (OPLS-AA) in water at finite temperature. The total free energies of folding and adhesion and the quality of surface binding were determined, based on a combination of solvation energy, formation of hydrogen bonds, and amount of the apolar (hydrophobic) contact surface between peptide and carbon surface. For both cases, we find a strong adhesion energy and large nonpolar contact surface. Isoleucines and tryptophans are the most strongly bound residues to the two carbon surfaces, the latter one largely dominating. In the case of the carbon nanotube, the peptide shows several competing stable structures, corresponding to different possible molecular foldings, and a propensity to enhance the intramolecular stability by forming new hydrogen bonds. In both systems, different arrangements of the histidine and tryptophan residues enable a better adaptation to the carbon surfaces. These findings suggest that the experimentally observed surface specificity of the peptide on nanotubes may depend on its capability to support multiple strongly bound configurations.