hdANM: a new comprehensive dynamics model for protein hinges.

Hinge motions are essential for many protein functions, and their dynamics are important to understand underlying biological mechanisms. The ways that these motions are represented by various computational methods differ significantly. By focusing on a specific class of motion, we have developed a new hinge-domain anisotropic network model (hdANM) that is based on the prior identification of flexible hinges and rigid domains in the protein structure and the subsequent generation of global hinge motions. This yields a set of motions in which the relative translations and rotations of the rigid domains are modulated and controlled by the deformation of the flexible hinges, leading to a more restricted, specific view of these motions. hdANM is the first model, to our knowledge, that combines information about protein hinges and domains to model the characteristic hinge motions of a protein. The motions predicted with this new elastic network model provide important conceptual advantages for understanding the underlying biological mechanisms. As a matter of fact, the generated hinge movements are found to resemble the expected mechanisms required for the biological functions of diverse proteins. Another advantage of this model is that the domain-level coarse graining makes it significantly more computationally efficient, enabling the generation of hinge motions within even the largest molecular assemblies, such as those from cryo-electron microscopy. hdANM is also comprehensive as it can perform in the same way as the well-known protein dynamics models (anisotropic network model, rotations-translations of blocks, and nonlinear rigid block normal mode analysis), depending on the definition of flexible and rigid parts in the protein structure and on whether the motions are extrapolated in a linear or nonlinear fashion. Furthermore, our results indicate that hdANM produces more realistic motions as compared to the anisotropic network model. hdANM is an open-source software, freely available, and hosted on a user-friendly website.

Low-Frequency Harmonic Perturbations Drive Protein Conformational Changes.

Protein dynamics has been investigated since almost half a century, as it is believed to constitute the fundamental connection between structure and function. Elastic network models (ENMs) have been widely used to predict protein dynamics, flexibility and the biological mechanism, from which remarkable results have been found regarding the prediction of protein conformational changes. Starting from the knowledge of the reference structure only, these conformational changes have been usually predicted either by looking at the individual mode shapes of vibrations (i.e., by considering the free vibrations of the ENM) or by applying static perturbations to the protein network (i.e., by considering a linear response theory). In this paper, we put together the two previous approaches and evaluate the complete protein response under the application of dynamic perturbations. Harmonic forces with random directions are applied to the protein ENM, which are meant to simulate the single frequency-dependent components of the collisions of the surrounding particles, and the protein response is computed by solving the dynamic equations in the underdamped regime, where mass, viscous damping and elastic stiffness contributions are explicitly taken into account. The obtained motion is investigated both in the coordinate space and in the sub-space of principal components (PCs). The results show that the application of perturbations in the low-frequency range is able to drive the protein conformational change, leading to remarkably high values of direction similarity. Eventually, this suggests that protein conformational change might be triggered by external collisions and favored by the inherent low-frequency dynamics of the protein structure.

Structural compliance: A new metric for protein flexibility.

Proteins are the active players in performing essential molecular activities throughout biology, and their dynamics has been broadly demonstrated to relate to their mechanisms. The intrinsic fluctuations have often been used to represent their dynamics and then compared to the experimental B-factors. However, proteins do not move in a vacuum and their motions are modulated by solvent that can impose forces on the structure. In this paper, we introduce a new structural concept, which has been called the structural compliance, for the evaluation of the global and local deformability of the protein structure in response to intramolecular and solvent forces. Based on the application of pairwise pulling forces to a protein elastic network, this structural quantity has been computed and sometimes is even found to yield an improved correlation with the experimental B-factors, meaning that it may serve as a better metric for protein flexibility. The inverse of structural compliance, namely the structural stiffness, has also been defined, which shows a clear anticorrelation with the experimental data. Although the present applications are made to proteins, this approach can also be applied to other biomolecular structures such as RNA. This present study considers only elastic network models, but the approach could be applied further to conventional atomic molecular dynamics. Compliance is found to have a slightly better agreement with the experimental B-factors, perhaps reflecting its bias toward the effects of local perturbations, in contrast to mean square fluctuations. The code for calculating protein compliance and stiffness is freely accessible at https://jerniganlab.github.io/Software/PACKMAN/Tutorials/compliance.

Terahertz vibration modes in Na/K-ATPase.

Mechanical vibration in the Terahertz range is believed to be connected with protein functions. In this paper, we present the results of a normal-mode analysis (modal analysis) of a Na/K-ATPase all-atom model, focusing the attention on low-frequency vibration modes. The numerical model helps in the interpretation of experimental results previously obtained by the authors via Raman spectroscopy of Na/K-ATPase samples, where several unassigned peaks were found in the sub-500 cm-1 range. In particular, vibration modes corresponding to peaks at 27, 190 and 300 cm-1, found experimentally, are confirmed here numerically, together with some other modes at lower frequencies (wavenumbers) that were not possible to observe in the experimental test. All the aforementioned modes correspond to vibrations involving the protein ends, i.e. portions directly related to the operating mechanism of the sodium-potassium pump.

The Sacred Mountain of Varallo in Italy: seismic risk assessment by acoustic emission and structural numerical models.

We examine an application of Acoustic Emission (AE) technique for a probabilistic analysis in time and space of earthquakes, in order to preserve the valuable Italian Renaissance Architectural Complex named "The Sacred Mountain of Varallo." Among the forty-five chapels of the Renaissance Complex, the structure of the Chapel XVII is of particular concern due to its uncertain structural condition and due to the level of stress caused by the regional seismicity. Therefore, lifetime assessment, taking into account the evolution of damage phenomena, is necessary to preserve the reliability and safety of this masterpiece of cultural heritage. A continuous AE monitoring was performed to assess the structural behavior of the Chapel. During the monitoring period, a correlation between peaks of AE activity in the masonry of the "Sacred Mountain of Varallo" and regional seismicity was found. Although the two phenomena take place on very different scales, the AE in materials and the earthquakes in Earth's crust, belong to the same class of invariance. In addition, an accurate finite element model, performed with DIANA finite element code, is presented to describe the dynamic behavior of Chapel XVII structure, confirming visual and instrumental inspections of regional seismic effects.

Anisotropic linear elastic properties of fractal-like composites.

In this work, the anisotropic linear elastic properties of two-phase composite materials, made up of square inclusions embedded in a matrix, are investigated. The inclusions present a fractal hierarchical distribution and are supposed to have the same Poisson's ratio as the matrix but a different Young's modulus. The effective elastic moduli of the medium are computed at each fractal iteration by coupling a position-space renormalization-group technique with a finite element analysis. The study allows to obtain and generalize some fundamental properties of fractal composite materials.

Phenomenological approach to mechanical damage growth analysis.

The problem of characterizing damage evolution in a generic material is addressed with the aim of tracing it back to existing growth models in other fields of research. Based on energetic considerations, a system evolution equation is derived for a generic damage indicator describing a material system subjected to an increasing external stress. The latter is found to fit into the framework of a recently developed phenomenological universality (PUN) approach and, more specifically, the so-called U2 class. Analytical results are confirmed by numerical simulations based on a fiber-bundle model and statistically assigned local strengths at the microscale. The fits with numerical data prove, with an excellent degree of reliability, that the typical evolution of the damage indicator belongs to the aforementioned PUN class. Applications of this result are briefly discussed and suggested.

Multiscale stochastic simulations for tensile testing of nanotube-based macroscopic cables.

Thousands of multiscale stochastic simulations are carried out in order to perform the first in-silico tensile tests of carbon nanotube (CNT)-based macroscopic cables with varying length. The longest treated cable is the space-elevator megacable but more realistic shorter cables are also considered in this bottom-up investigation. Different sizes, shapes, and concentrations of defects are simulated, resulting in cable macrostrengths not larger than approximately 10 GPa, which is much smaller than the theoretical nanotube strength (approximately 100 GPa). No best-fit parameters are present in the multiscale simulations: the input at level 1 is directly estimated from nanotensile tests of CNTs, whereas its output is considered as the input for the level 2, and so on up to level 5, corresponding to the megacable. Thus, five hierarchical levels are used to span lengths from that of a single nanotube (approximately 100 nm) to that of the space-elevator megacable (approximately 100 Mm).

The dynamic evolution of the power exponent in a universal growth model of tumors.

We have previously reported that a universal growth law, as proposed by West and collaborators for all living organisms, appears to be able to describe also the growth of tumors in vivo after an initial exponential growth phase. In contrast to the assumption of a fixed power exponent p (assumed by West et al. to be equal to 3/4), we propose in this paper a dynamic evolution of p, using experimental data from the cancer literature. In analogy with results obtained by applying scaling laws to the study of fragmentation of solids, the dynamic behaviour of p is related to the evolution of the fractal topology of neoplastic vascular systems. Our model might be applied for diagnostic purposes to mark the emergence of an efficient neo-angiogenetic structure if the results of our in silico experiments are confirmed by clinical observations.

Are scaling laws on strength of solids related to mechanics or to geometry?

One of the largest controversial issues of the materials science community is the interpretation of scaling laws on material strength. In spite of the prevailing view, which considers mechanics as the real cause of such effects, here, we propose a different argument, purely based on geometry. Thus, as happened for relativity, geometry could again hold an unexpected and fundamental role.