Complete spatiotemporal quantification of cardiac motion in mice through enhanced acquisition and super-resolution reconstruction.

The quantification of cardiac motion using cardiac magnetic resonance imaging (CMR) has shown promise as an early-stage marker for cardiovascular diseases. Despite the growing popularity of CMR-based myocardial strain calculations, measures of complete spatiotemporal strains (i.e., three-dimensional strains over the cardiac cycle) remain elusive. Complete spatiotemporal strain calculations are primarily hampered by poor spatial resolution, with the rapid motion of the cardiac wall also challenging the reproducibility of such strains. We hypothesize that a super-resolution reconstruction (SRR) framework that leverages combined image acquisitions at multiple orientations will enhance the reproducibility of complete spatiotemporal strain estimation. Two sets of CMR acquisitions were obtained for five wild-type mice, combining short-axis scans with radial and orthogonal long-axis scans. Super-resolution reconstruction, integrated with tissue classification, was performed to generate full four-dimensional (4D) images. The resulting enhanced and full 4D images enabled complete quantification of the motion in terms of 4D myocardial strains. Additionally, the effects of SRR in improving accurate strain measurements were evaluated using an in-silico heart phantom. The SRR framework revealed near isotropic spatial resolution, high structural similarity, and minimal loss of contrast, which led to overall improvements in strain accuracy. In essence, a comprehensive methodology was generated to quantify complete and reproducible myocardial deformation, aiding in the much-needed standardization of complete spatiotemporal strain calculations.

Multiscale modeling of unfolding and bond dissociation of rubredoxin metalloprotein.

Mechanical properties of proteins that have a crucial effect on their operation. This study used a molecular dynamics simulation package to investigate rubredoxin unfolding on the atomic scale. Different simulation techniques were applied, and due to the dissociation of covalent/hydrogen bonds, this protein demonstrates several intermediate states in force-extension behavior. A conceptual model based on the cohesive finite element method was developed to consider the intermediate damages that occur during unfolding. This model is based on force-displacement curves derived from molecular dynamics results. The proposed conceptual model is designed to accurately identify bond rupture points and determine the associated forces. This is achieved by conducting a thorough comparison between molecular dynamics and cohesive finite element results. The utilization of a viscoelastic cohesive zone model allows for the consideration of loading rate effects. This rate-dependent model can be further developed and integrated into the multiscale modeling of large assemblies of metalloproteins, providing a comprehensive understanding of mechanical behavior while maintaining a reduced computational cost.

Decisive structural elements in water and ion permeation through mechanosensitive channels of large conductance: insights from molecular dynamics simulation.

In this paper, a series of equilibrium molecular dynamics simulations (EMD), steered molecular dynamics (SMD), and computational electrophysiology methods are carried out to explore water and ion permeation through mechanosensitive channels of large conductance (MscL). This research aims to identify the pore-lining side chains of the channel in different conformations of MscL homologs by analyzing the pore size. The distribution of permeating water dipole angles through the pore domains enclosed by VAL21 and GLU104 demonstrated that water molecules are oriented toward the charged oxygen headgroups of GLU104 from their hydrogen atoms to retain this interaction in a stabilized fashion. Although, this behavior was not perceived for VAL21. Numerical assessments of the secondary structure clarified that, during the ion permeation, in addition to the secondary structure alterations, the structure of Tb-MscL would also undergo significant conformational changes. It was elucidated that VAL21, GLU104, and water molecules accomplish a fundamental task in ion permeation. The mentioned residues hinder ion permeation so that the pulling SMD force is increased remarkably when the ions permeate through the domains enclosed by VAL21 and GLU102. The hydration level and potassium diffusivity in the hydrophobic gate of the transmembrane domain were promoted by applying the external electric field. Furthermore, the implementation of an external electric field altered the distribution pattern for potassium ions in the system while intensifying the accumulation of Cl- in the vicinity of ARG11 and ARG98.

Unraveling Flow Separation at the Water-Carbon Nanotube Interface: An Atomic-Scale Overview by Molecular Dynamics Simulation.

Flow separation near the fluid-solid surface has attracted attention for decades. It is critical to understand the behavior of separated flow adjacent to the solid walls to broaden its range of potential applications. Therefore, we conducted molecular dynamics investigations to consider water flow separation at the water-carbon nanotube (CNT) interface for different diameters of CNTs between 13 and 50 Å and different pressures of 0.1-1.254 GPa. Density heat maps indicated that water flow separation is observed for all CNTs under high pressures, and an empty space of water molecules or evacuation is formed behind the CNTs. It is shown that in CNTs with small diameters, (10, 10) and (20, 20), the structure of the first layer (FL) of water molecules or hydrated layer adjacent to the CNT wall is completely preserved, indicating that evacuation occurs from behind the CNTs. In (30, 30) and (40, 40) CNTs, flow separation occurred from the FL of water molecules near the solid surface, and the layered structure of water around CNTs is completely destroyed. Our findings of fluid-solid and fluid-fluid interaction energies suggested that the flow separation can be due to an attraction between the FL of water molecules and CNT and a repulsion between the water molecules in the hydrated layer and the outer layers. Moreover, analyzing the relationship between the CNT size and flow separation revealed that in the case of small CNTs, there are extra water molecules that contribute to the structural stability of the hydrated layer by strengthening the repulsive interaction in the liquid-liquid surface.

A V(iii)-induced metallogel with solvent stimuli-responsive properties: structural proof-of-concept with MD simulations.

A new solvent stimuli-responsive metallogel (VGel) was synthesized through the introduction of vanadium ions into an adenine (Ade) and 1,3,5-benzene tricarboxylic acid (BTC) organogel, and its supramolecular self-assembly was investigated from a computational viewpoint. A relationship between the synthesized VGel integrity and the self-assembly of its components is demonstrated by a broad range of molecular dynamics (MD) simulations, an aspect that has not yet been explored for such a complex metallogel in particular. MD simulations and Voronoi tessellation assessments, both in agreement with experimental data, confirm the gel formation. Based on excellent water stability and the ethanol/methanol stimuli-responsive feature of the VGel an easy-to-use visualization assay for the detection of counterfeit liquor with a 6% (v/v) methanol limit of detection in 40% (v/v) ethanol is reported. These observations provide a cheap and technically simple method and are a step towards the immersible screening of similar molecules in methanol-spiked beverages.

A label-free graphene-based impedimetric biosensor for real-time tracing of the cytokine storm in blood serum; suitable for screening COVID-19 patients.

Concurrent with the pandemic announcement of SARS-CoV-2 infection by the WHO, a variety of reports were published confirming the cytokine storm as the most mortal effect of the virus on the infected patients. Hence, cytokine storm as an evidenced consequence in most of the COVID-19 patients could offer a promising opportunity to use blood as a disease progression marker. Here, we have developed a rapid electrochemical impedance spectroscopy (EIS) sensor for quantifying the overall immune activity of the patients. Since during the cytokine storm many types of cytokines are elevated in the blood, there is no need for specific detection of a single type of cytokine and the collective behavior is just measured without any electrode functionalization. The sensor includes a monolayer graphene on a copper substrate as the working electrode (WE) which is able to distinguish between the early and severe stage of the infected patients. The charge transfer resistance (R CT) in the moderate and severe cases varies about 65% and 138% compared to the normal groups, respectively and a specificity of 77% and sensitivity of 100% based on ELISA results were achieved. The outcomes demonstrate a significant correlation between the total mass of the three main hypercytokinemia associated cytokines including IL-6, TNF-α and IFN-γ in patients and the R CT values. As an extra application, the biosensor's capability for diagnosis of COVID-19 patients was tested and a sensitivity of 92% and specificity of 50% were obtained compared to the RT-PCR results.

Enhanced wettability of long narrow carbon nanotubes in a double-walled hetero-structure: unraveling the effects of a boron nitride nanotube as the exterior.

Studying the structure and dynamics of nano-confined water inside carbon nanotubes has consistently attracted the wide-spread interest of researchers. In the present work, molecular dynamics simulations indicated internal nonwetting behavior for the central region of the long and narrow single-wall carbon nanotube (5,5) (SWNT) and showed that continuous single-file water molecules are not formed through it. Unlike the SWNT, by adding boron nitride nanotubes (6,6) as an outer wall to the SWNT, a continuously long single-file water chain is formed through the double-walled carbon and boron nitride hetero-nanotube (DWHNT) and thorough internal wetting of the DWHNT is observed. The position and the number of water molecules, electrostatic potential heatmap of the nanotube's wall, free energy profile of nano-confined water, and number of hydrogen bonds between them confirmed the aforementioned results and complete internal wetting of the DWHNT. After using the boron nitride nanotube (6,6) as the outer wall, an homogeneous electrostatic potential distribution in the DWHNT and increase in the hydrophilic characteristics of the nano-channel wall are observed, bringing about gradual trapping of more water molecules through it. Finally, water molecules occupied the central region of the DWHNT and a thorough single-file water chain is formed inside the nano-channel. Water dipole orientation inside the DWHNT and their radial distribution function asserted the occurrence of the liquid-solid quasi-phase transition of single-file water molecules confined inside the long and narrow carbon nanotube (5,5) under ambient conditions.

Evaporation of Water on Suspended Graphene: Suppressing the Effect of Physically Heterogeneous Surfaces.

Evaporation of water nanodroplets on a hydrophilically adjusted graphene sheet was studied based on a molecular dynamics approach. Suspended graphene was used as a physically heterogeneous surface, and fixed graphene was considered as an ideally flat surface. State of the triple-phase contact line (TPCL) and shape evolution were addressed at four different temperatures on both substrates. Additionally, contact angle (CA) was studied during 3 and 22.5 ns simulations in both closed and opened conditions. The observed constant contact angle regime was predictable for the fixed graphene. However, it was not expected for the suspended system and was attributed to the oscillations of the substrate atoms. The size of the nanodroplet also affects the constant-contact-angle mode in both systems, when the number of water molecules decreases to less than 500. The oscillations created a surface on which physical heterogeneities were varying through time. Examination of the evaporation and condensation processes revealed higher rates for the fixed systems. Local mass fluxes were calculated to reveal the contribution of TPCL and meridian surface (MS) of the nanodroplet to evaporation and condensation. The obtained results indicate similar values for the mass flux ratio at the TPCL, which remains twice as large as the MS for both suspended and fixed graphene. The results confirm the assumption that a surface with varying heterogeneities can overwhelm the droplet and act as an ideally flat surface.

Determinative factors in inhibition of aquaporin by different pharmaceuticals: Atomic scale overview by molecular dynamics simulation.

The inhibition of water permeation through aquaporins by ligands of pharmaceutical compounds is considered as a method to control the cell lifetime. The inhibition of aquaporin 1 (AQP1) by bacopaside-I and torsemide, was explored and its atomistic nature was elucidated by molecular docking and molecular dynamics (MD) simulation collectively along with Poisson-Boltzmann surface area (PBSA) method. Docking results revealed that torsemide has a lower level of docking energy in comparison with bacopaside-I at the cytoplasmic side. Furthermore, the effect of steric constraints on water permeation was accentuated. Bacopaside-I inhibits the channel properly due to the strong interaction with the channel and larger spatial volume, whereas torsemide blocks the cytoplasmic side of the channel imperfectly. The most probable active sites of AQP1 for the formation of hydrogen bonds between the inhibitor and the channel were identified by numerical analysis of the bonds. Eventually, free energy assessments indicate that binding of both inhibitors is favorable in complex with AQP1, and van der Waals interaction has an important contribution in stabilizing the complexes.

Linear and nonlinear analyses of femoral fractures: Computational/experimental study.

This paper describes two new methods for computational fracture analysis of human femur using Quantitative Computed Tomography (QCT) voxel-based finite element (FE) simulation. The paper also reports comprehensive mechanical testing for validation of the methods and evaluation of the required material properties. The analyses and tests were carried out on 15 human femurs under 11 different stance-type loading orientations. Several classical forms of subcapital, transcervical, basicervical, and intertrochanteric fractures plus a specific type of subtrochanteric fracture were created and analyzed. A new procedure was developed for prediction of the strengths and the fracture initiation patterns using a FE-based linear scheme. The predicted and observed fracture patterns were in correspondence, and the FE predictions of the fracture loads were in very good agreement with the experimental results. Moreover, the crack initiation and growth behaviors of two subtrochanteric fractures were successfully simulated through a novel implementation of the cohesive zone model (CZM) within a nonlinear FE analysis scheme. The CZM parameters were obtained through a series of experimental tests on different types of specimens and determination of a variety of material properties for different anatomic regions and orientations. The presented results indicated that the locations and patterns of crack initiation, the sequences of crack growth on different paths, and the compatibility of growth increments agreed very well with the observed specifications. Also, very good agreements were achieved between the measured and simulated fracture loads.

Analysis of strength and failure pattern of human proximal femur using quantitative computed tomography (QCT)-based finite element method.

This paper presents a novel method for fast and reliable prediction of the failure strength of human proximal femur, using the quantitative computed tomography (QCT)-based linear finite element analysis (FEA). Ten fresh frozen human femora (age: 34±16) were QCT-scanned and the pertinent 3D voxel-based finite element models were constructed. A specially-designed holding frame was used to define and maintain a unique geometrical reference system for both FEA and in-vitro mechanical testing. The analyses and tests were carried out at 8 different loading orientations. A new scheme was developed for assortment of the element risk factor (defined as the ratio of the strain energy density to the yield strain energy for each element) and implemented for the prediction of the failure strength. The predicted and observed failure patterns were in correspondence, and the FEA predictions of the failure loads were in very good agreement with the experimental results (R2=0.86, slope=0.96, p<0.01). The average computational time was 5 min (on a regular desktop personal computer) for an average element number of 197,000. Noting that the run-time for a similar nonlinear model is about 8h, it was concluded that the proposed linear scheme is overwhelmingly efficient in terms of computational costs. Thus, it can efficiently be used to predict the femoral failure strength with the same accuracy of similar nonlinear models.