Advancing clinical translation of cardiac biomechanics models: a comprehensive review, applications and future pathways
Cardiac mechanics models are developed to represent a high level of detail, including refined anatomies, accurate cell mechanics models, and platforms to link microscale physiology to whole-organ function. However, cardiac biomechanics models still have limited clinical translation. In this review, we provide a picture of cardiac mechanics models, focusing on their clinical translation. We review the main experimental and clinical data used in cardiac models, as well as the steps followed in the literature to generate anatomical meshes ready for simulations. We describe the main models in active and passive mechanics and the different lumped parameter models to represent the circulatory system. Lastly, we provide a summary of the state-of-the-art in terms of ventricular, atrial, and four-chamber cardiac biomechanics models. We discuss the steps that may facilitate clinical translation of the biomechanics models we describe. A well-established software to simulate cardiac biomechanics is lacking, with all available platforms involving different levels of documentation, learning curves, accessibility, and cost. Furthermore, there is no regulatory framework that clearly outlines the verification and validation requirements a model has to satisfy in order to be reliably used in applications. Finally, better integration with increasingly rich clinical and/or experimental datasets as well as machine learning techniques to reduce computational costs might increase model reliability at feasible resources. Cardiac biomechanics models provide excellent opportunities to be integrated into clinical workflows, but more refinement and careful validation against clinical data are needed to improve their credibility. In addition, in each context of use, model complexity must be balanced with the associated high computational cost of running these models.
Comparison of translational and rotational modes towards passive rheology of the cytoplasm of MCF-7 cells using optical tweezers
A colloidal particle placed inside the cell cytoplasm is enmeshed within a network of cytoskeletal fibres immersed in the cytosolic fluid. The translational mode is believed to yield different rheological parameters than the rotational mode, given that these modes stretch the fibers differently. We compare the parameters for Michigan Cancer Foundation-7 (MCF-7) cells in this manuscript and find that the results are well comparable to each other. At low values of 0 Hz viscosity, the rotational and translational viscoelasticity matches well. However, discrepancies appear at higher values which may indicate that the cytoskeletal modes involved in rotation and translation of the particle are getting invoked. We also show that the 0 Hz viscosity increases as the cell ages under the conditions of constant room temperature of 25C on the sample chamber.
Overview of Methods for Noise and Heat Reduction in MRI Gradient Coils
Magnetic resonance imaging (MRI) gradient coils produce acoustic noise due to coil conductor vibrations caused by large Lorentz forces. Accurate sound pressure levels and modeling of heating are essential for the assessment of gradient coil safety. This work reviews the state-of-the-art numerical methods used in accurate gradient coil modeling and prediction of sound pressure levels (SPLs) and temperature rise. We review several approaches proposed for noise level reduction of high-performance gradient coils, with a maximum noise reduction of 20 decibels (dB) demonstrated. An efficient gradient cooling technique is also presented.
Applications of Machine Learning to Improve the Clinical Viability of Compton Camera Based Range Verification in Proton Radiotherapy
We studied the application of a deep, fully connected Neural Network (NN) to process prompt gamma (PG) data measured by a Compton camera (CC) during the delivery of clinical proton radiotherapy beams. The network identifies 1) recorded "bad" PG events arising from background noise during the measurement, and 2) the correct ordering of PG interactions in the CC to help improve the fidelity of "good" data used for image reconstruction. PG emission from a tissue-equivalent target during irradiation with a 150 MeV proton beam delivered at clinical dose rates was measured with a prototype CC. Images were reconstructed from both the raw measured data and the measured data that was further processed with a neural network (NN) trained to identify "good" and "bad" PG events and predict the ordering of individual interactions within the good PG events. We determine if NN processing of the CC data could improve the reconstructed PG images to a level in which they could provide clinically useful information about the range and range shifts of the proton beams delivered at full clinical dose rates. Results showed that a deep, fully connected NN improved the achievable contrast to noise ratio (CNR) in our images by more than a factor of 8x. This allowed the path, range, and lateral width of the clinical proton beam within a tissue equivalent target to easily be identified from the PG images, even at the highest dose rates of a 150 MeV proton beam used for clinical treatments. On average, shifts in the beam range as small as 3 mm could be identified. However, when limited by the amount of PG data measured with our prototype CC during the delivery of a single proton pencil beam (~1 × 10 protons), the uncertainty in the reconstructed PG images limited the identification of range shift to ~5 mm. Substantial improvements in CC images were obtained during clinical beam delivery through NN pre-processing of the measured PG data. We believe this shows the potential of NNs to help improve and push CC-based PG imaging toward eventual clinical application for proton RT treatment delivery verification.
Minimally Invasive Image-Guided Gut Transport Function Measurement Probe
Diseases such as celiac disease, environmental enteric dysfunction, infectious gastroenteritis, type II diabetes and inflammatory bowel disease are associated with increased gut permeability. Dual sugar absorption tests, such as the lactulose to rhamnose ratio (L:R) test, are the current standard for measuring gut permeability. Although easy to administer in adults, the L:R test has a number of drawbacks. These include an inability to assess for spatial heterogeneity in gut permeability that may distinguish different disease severity or pathology, additional sample collection for immunoassays, and challenges in carrying out the test in certain populations such as infants and small children. Here, we demonstrate a minimally invasive probe for real-time localized gut permeability evaluation through gut potential difference (GPD) measurement.
Toward the Development of an On-Chip Acoustic Focusing Fluorescence Lifetime Flow Cytometer
Conventional flow cytometry is a valuable quantitative tool. Flow cytometers reveal physical and biochemical information from cells at a high throughput, which is quite valuable for many biomedical, biological, and diagnostic research fields. Flow cytometers range in complexity and typically provide multiparametric data for the user at rates of up to 50,000 cells measured per second. Cytometry systems are configured such that fluorescence or scattered light signals are collected per-cell, and the integrated optical signal at a given wavelength range indicates a particular cellular feature such as phenotype or morphology. When the timing of the optical signal is measured, the cytometry system becomes "time-resolved." Time-resolved flow cytometry (TRFC) instruments can detect fluorescence decay kinetics, and such measurements are consequential for Förster Resonance Energy Transfer (FRET) studies, multiplexing, and metabolic mapping, to name a few. TRFC systems capture fluorescence lifetimes at rates of thousands of cells per-second, however the approach is challenged at this throughput by terminal cellular velocities. High flow rates limit the total number of photons integrated per-cell, reducing the reliability of the average lifetime as a cytometric parameter. In this contribution, we examine an innovative approach to address this signal-to-noise issue. The technology merges time-resolved hardware with microfluidics and acoustics. We present an "acoustofluidic" time-resolved flow cytometer so that cellular velocities can be adjusted on the fly with a standing acoustic wave (SAW). Our work shows that acoustic control can be combined with time-resolved features to appropriately balance the throughput with the optical signals necessary for lifetime data.
Regulation of Actin Bundle Mechanics and Structure by Intracellular Environmental Factors
The mechanical and structural properties of actin cytoskeleton drive various cellular processes, including structural support of the plasma membrane and cellular motility. Actin monomers assemble into double-stranded helical filaments as well as higher-ordered structures such as bundles and networks. Cells incorporate macromolecular crowding, cation interactions, and actin-crosslinking proteins to regulate the organization of actin bundles. Although the roles of each of these factors in actin bundling have been well-known individually, how combined factors contribute to actin bundle assembly, organization, and mechanics is not fully understood. Here, we describe recent studies that have investigated the mechanisms of how intracellular environmental factors influence actin bundling. This review highlights the effects of macromolecular crowding, cation interactions, and actin-crosslinking proteins on actin bundle organization, structure, and mechanics. Understanding these mechanisms is important in determining actin biophysics and providing insights into cell physiology.
Systems Biology Modeling of the Complement System Under Immune Susceptible Pathogens
The complement system is assembled from a network of proteins that function to bring about the first line of defense of the body against invading pathogens. However, complement deficiencies or invasive pathogens can hijack complement to subsequently increase susceptibility of the body to infections. Moreover, invasive pathogens are increasingly becoming resistant to the currently available therapies. Hence, it is important to gain insights into the highly dynamic interaction between complement and invading microbes in the frontlines of immunity. Here, we developed a mathematical model of the complement system composed of 670 ordinary differential equations with 328 kinetic parameters, which describes all three complement pathways (alternative, classical, and lectin) and includes description of mannose-binding lectin, collectins, ficolins, factor H-related proteins, immunoglobulin M, and pentraxins. Additionally, we incorporate two pathogens: (type 1) complement susceptible pathogen and (type 2) located in either nasopharynx or bloodstream. In both cases, we generate time profiles of the pathogen surface occupied by complement components and the membrane attack complex (MAC). Our model shows both pathogen types in bloodstream are saturated by complement proteins, whereas MACs occupy <<1.0% of the pathogen surface. Conversely, the MAC production in nasopharynx occupies about 1.5-10% of the total surface, thus making nasal MAC levels at least about eight orders of magnitude higher. Altogether, we predict complement-imbalance, favoring overactivation, is associated with nasopharynx homeostasis. Conversely, orientating toward complement-balance may cause disruption to the nasopharynx homeostasis. Thus, for sporadic meningococcal disease, our model predicts rising nasal levels of complement regulators as early infection biomarkers.
Surface Texturing and Combinatorial Approaches to Improve Biocompatibility of Implanted Biomaterials
Biomaterial associated microbial infection and blood thrombosis are two of the barriers that inhibit the successful use of implantable medical devices in modern healthcare. Modification of surface topography is a promising approach to combat microbial infection and thrombosis without altering bulk material properties necessary for device function and without contributing to bacterial antibiotic resistance. Similarly, the use of other antimicrobial techniques such as grafting poly(ethylene glycol) (PEG) and nitric oxide (NO) release also improve the biocompatibility of biomaterials. In this review, we discuss the development of surface texturing techniques utilizing ordered submicron-size pillars for controlling bacterial adhesion and biofilm formation, and we present combinatorial approaches utilizing surface texturing in combination with poly(ethylene glycol) (PEG) grafting and NO release to improve the biocompatibility of biomaterials. The manuscript also discusses efforts towards understanding the molecular mechanisms of bacterial adhesion responses to the surface texturing and NO releasing biomaterials, focusing on experimental aspects of the approach.
Feasibility of Hepatic Fat Quantification Using Proton Density Fat Fraction by Multi-Echo Chemical-Shift-Encoded MRI at 7T
Fat fraction quantification and assessment of its distribution in the hepatic tissue become more important with the growing epidemic of obesity, and the increasing prevalence of diabetes mellitus type 2 and non-alcoholic fatty liver disease. At 3Tesla, the multi-echo, chemical-shift-encoded magnetic resonance imaging (CSE-MRI)-based acquisition allows the measurement of proton density fat-fraction (PDFF) even in clinical protocols. Further improvements in SNR can be achieved by the use of phased array coils and increased static magnetic field. The purpose of the study is to evaluate the feasibility of PDFF imaging using a multi-echo CSE-MRI technique at ultra-high magnetic field (7Tesla). Thirteen volunteers (M/F) with a broad range of age, body mass index, and hepatic PDFF were measured at 3 and 7T by multi-gradient-echo MRI and single-voxel spectroscopy MRS. All measurements were performed in breath-hold (exhalation); the MRI protocols were optimized for a short measurement time, thus minimizing motion-related problems. 7T data were processed off-line using Matlab® (MRI:multi-gradient-echo) and jMRUI (MRS), respectively. For quantitative validation of the PDFF results, a similar protocol was performed at 3T, including on-line data processing provided by the system manufacturer, and correlation analyses between 7 and 3T data were performed off-line. The multi-echo CSE-MRI measurements at 7T with a phased-array coil configuration and an optimal post-processing yielded liver volume coverage ranging from 30 to 90% for high- and low-BMI subjects, respectively. PDFFs ranged between 1 and 20%. We found significant correlations between 7T MRI and -MRS measurements (R ≅ 0.97; 0.005), and between MRI-PDFF at 7T and 3T fields (R ≅ 0.94; < 0.005) in the evaluated volumes. Based on the measurements and analyses performed, the multi-echo CSE-MRI method using a 32-channel coil at 7T showed its aptitude for MRI-based quantitation of PDFF in the investigated volumes. The results are the first step toward qMRI of the whole liver at 7T with further improvements in hardware.
Disentangling the effects of restriction and exchange with diffusion exchange spectroscopy
Diffusion exchange spectroscopy (DEXSY) is a multidimensional NMR technique that can reveal how water molecules exchange between compartments within heterogeneous media, such as biological tissue. Data from DEXSY experiments is typically processed using numerical inverse Laplace transforms (ILTs) to produce a diffusion-diffusion spectrum. A tacit assumption of this ILT approach is that the signal behavior is Gaussian - i.e., the spin echo intensity decays exponentially with the degree of diffusion weighting. The assumptions that underlie Gaussian signal behavior may be violated, however, depending on the gradient strength applied and the sample under study. We argue that non-Gaussian signal behavior due to restrictions is to be expected in the study of biological tissue using diffusion NMR. Further, we argue that this signal behavior can produce confounding features in the diffusion-diffusion spectra obtained from numerical ILTs of DEXSY data - entangling the effects of restriction and exchange. Specifically, restricted signal behavior can result in broadening of peaks and in the appearance of illusory exchanging compartments with distributed diffusivities, which pearl into multiple peaks if not highly regularized. We demonstrate these effects on simulated data. That said, we suggest the use of features in the signal acquisition domain that can be used to rapidly probe exchange without employing an ILT. We also propose a means to characterize the non-Gaussian signal behavior due to restrictions within a sample using DEXSY measurements with a near zero mixing time or storage interval. We propose a combined acquisition scheme to independently characterize restriction and exchange with various DEXSY measurements, which we term Restriction and Exchange from Equally-weighted Double and Single Diffusion Encodings (REEDS-DE). We test this method on neonatal mouse spinal cord - a sample consisting primarily of gray matter - using a low-field, static gradient NMR system. In sum, we highlight critical shortcomings of prevailing DEXSY analysis methods that conflate the effects of restriction and exchange, and suggest a viable experimental approach to disentangle them.
The role of antifreeze glycoprotein (AFGP) and polyvinyl alcohol/polyglycerol (X/Z-1000) as ice modulators during partial freezing of rat livers
The current liver organ shortage has pushed the field of transplantation to develop new methods to prolong the preservation time of livers from the current clinical standard of static cold storage. Our approach, termed partial freezing, aims to induce a thermodynamically stable frozen state at high subzero storage temperatures (-10°C to -15°C), while simultaneously maintaining a sufficient unfrozen fraction to limit ice-mediated injury.
The Accuracy of Cerenkov Photons Simulation in Geant4/Gate Depends on the Parameterization of Primary Electron Propagation
Energetic electrons traveling in a dispersive medium can produce Cerenkov radiation. Cerenkov photons' prompt emission, combined with their predominantly forward emission direction with respect to the parent electron, makes them extremely promising to improve radiation detector timing resolution. Triggering gamma detections based on Cerenkov photons to achieve superior timing resolution is challenging due to the low number of photons produced per interaction. Monte Carlo simulations are fundamental to understanding their behavior and optimizing their pathway to detection. Therefore, accurately modeling the electron propagation and Cerenkov photons emission is crucial for reliable simulation results. In this work, we investigated the physics characteristics of the primary electrons (velocity, energy) and those of all emitted Cerenkov photons (spatial and timing distributions) generated by 511 keV photoelectric interactions in a bismuth germanate crystal using simulations with Geant4/GATE. Geant4 uses a stepwise particle tracking approach, and users can limit the electron velocity change per step. Without limiting it (default Geant4 settings), an electron mean step length of ~250 μm was obtained, providing only macroscopic modeling of electron transport, with all Cerenkov photons emitted in the forward direction with respect to the incident gamma direction. Limiting the electron velocity change per step reduced the electron mean step length (~0.200 μm), leading to a microscopic approach to its transport which more accurately modeled the electron physical properties in BGO at 511 keV. The electron and Cerenkov photons rapidly lost directionality, affecting Cerenkov photons' transport and, ultimately, their detection. Results suggested that a deep understanding of low energy physics is crucial to perform accurate optical Monte Carlo simulations and ultimately use them in TOF PET detectors.
Motor-driven advection competes with crowding to drive spatiotemporally heterogeneous transport in cytoskeleton composites
The cytoskeleton-a composite network of biopolymers, molecular motors, and associated binding proteins-is a paradigmatic example of active matter. Particle transport through the cytoskeleton can range from anomalous and heterogeneous subdiffusion to superdiffusion and advection. Yet, recapitulating and understanding these properties-ubiquitous to the cytoskeleton and other out-of-equilibrium soft matter systems-remains challenging. Here, we combine light sheet microscopy with differential dynamic microscopy and single-particle tracking to elucidate anomalous and advective transport in actomyosin-microtubule composites. We show that particles exhibit multi-mode transport that transitions from pronounced subdiffusion to superdiffusion at tunable crossover timescales. Surprisingly, while higher actomyosin content increases the range of timescales over which transport is superdiffusive, it also markedly increases the degree of subdiffusion at short timescales and generally slows transport. Corresponding displacement distributions display unique combinations of non-Gaussianity, asymmetry, and non-zero modes, indicative of directed advection coupled with caged diffusion and hopping. At larger spatiotemporal scales, particles in active composites exhibit superdiffusive dynamics with scaling exponents that are robust to changing actomyosin fractions, in contrast to normal, yet faster, diffusion in networks without actomyosin. Our specific results shed important new light on the interplay between non-equilibrium processes, crowding and heterogeneity in active cytoskeletal systems. More generally, our approach is broadly applicable to active matter systems to elucidate transport and dynamics across scales.
Editorial: Innovations in MR hardware from ultra-low to ultra-high field
Editorial: Capturing Biological Complexity and Heterogeneity Using Multidimensional MRI
Toward a 3D physical model of diffusive polymer chains
Recent studies in polymer physics have created macro-scale analogs to solute microscopic polymer chains like DNA by inducing diffusive motion on a chain of beads. These bead chains have persistence lengths of O(10) links and undergo diffusive motion under random fluctuations like vibration. We present a bead chain model within a new stochastic forcing system: an air fluidizing bed of granular media. A chain of spherical 6 mm resin beads crimped onto silk thread are buffeted randomly by the multiphase flow of grains and low density rising air "bubbles". We "thermalize" bead chains of various lengths at different fluidizing airflow rates, while X-ray imaging captures a projection of the chains' dynamics within the media. With modern 3D printing techniques, we can better represent complex polymers by geometrically varying bead connections and their relative strength, e.g., mimicking the variable stiffness between adjacent nucleotide pairs of DNA. We also develop Discrete Element Method (DEM) simulations to study the 3D motion of the bead chain, where the bead chain is represented by simulated spherical particles connected by linear and angular spring-like bonds. In experiment, we find that the velocity distributions of the beads follow exponential distributions rather than the Gaussian distributions expected from polymers in solution. Through use of the DEM simulation, we find that this difference can likely be attributed to the distributions of the forces imparted onto the chain from the fluidized bed environment. We anticipate expanding this study in the future to explore a wide range of chain composition and confinement geometry, which will provide insights into the physics of large biopolymers.
Non-contact and label-free biomechanical imaging: Stimulated Brillouin microscopy and beyond
Brillouin microscopy based on spontaneous Brillouin scattering has emerged as a unique elastography technique because of its merit of non-contact, label-free, and high-resolution mechanical imaging of biological cell and tissue. Recently, several new optical modalities based on stimulated Brillouin scattering have been developed for biomechanical research. As the scattering efficiency of the stimulated process is much higher than its counterpart in the spontaneous process, stimulated Brillouin-based methods have the potential to significantly improve the speed and spectral resolution of existing Brillouin microscopy. Here, we review the ongoing technological advancements of three methods, including continuous wave stimulated Brillouin microscopy, impulsive stimulated Brillouin microscopy, and laser-induced picosecond ultrasonics. We describe the physical principle, the representative instrumentation, and biological application of each method. We further discuss the current limitations as well as the challenges for translating these methods into a visible biomedical instrument for biophysics and mechanobiology.
Critical brain wave dynamics of neuronal avalanches
Analytical expressions for scaling of brain wave spectra derived from the general non-linear wave Hamiltonian form show excellent agreement with experimental "neuronal avalanche" data. The theory of the weakly evanescent non-linear brain wave dynamics reveals the underlying collective processes hidden behind the phenomenological statistical description of the neuronal avalanches and connects together the whole range of brain activity states, from oscillatory wave-like modes, to neuronal avalanches, to incoherent spiking, showing that the neuronal avalanches are just the manifestation of the different non-linear side of wave processes abundant in cortical tissue. In a more broad way these results show that a system of wave modes interacting through all possible combinations of the third order non-linear terms described by a general wave Hamiltonian necessarily produces anharmonic wave modes with temporal and spatial scaling properties that follow scale free power laws. To the best of our knowledge this has never been reported in the physical literature and may be applicable to many physical systems that involve wave processes and not just to neuronal avalanches.
Active condensation of filaments under spatial confinement
Living systems exhibit self-organization, a phenomenon that enables organisms to perform functions essential for life. The interior of living cells is a crowded environment in which the self-assembly of cytoskeletal networks is spatially constrained by membranes and organelles. Cytoskeletal filaments undergo active condensation in the presence of crosslinking motor proteins. In past studies, confinement has been shown to alter the morphology of active condensates. Here, we perform simulations to explore systems of filaments and crosslinking motors in a variety of confining geometries. We simulate spatial confinement imposed by hard spherical, cylindrical, and planar boundaries. These systems exhibit non-equilibrium condensation behavior where crosslinking motors condense a fraction of the overall filament population, leading to coexistence of vapor and condensed states. We find that the confinement lengthscale modifies the dynamics and condensate morphology. With end-pausing crosslinking motors, filaments self-organize into half asters and fully-symmetric asters under spherical confinement, polarity-sorted bilayers and bottle-brush-like states under cylindrical confinement, and flattened asters under planar confinement. The number of crosslinking motors controls the size and shape of condensates, with flattened asters becoming hollow and ring-like for larger motor number. End pausing plays a key role affecting condensate morphology: systems with end-pausing motors evolve into aster-like condensates while those with non-end-pausing crosslinking motor proteins evolve into disordered clusters and polarity-sorted bundles.
A framework of computer vision-enhanced microfluidic approach for automated assessment of the transient sickling kinetics in sickle red blood cells
The occurrence of vaso-occlusive crisis greatly depends on the competition between the sickling delay time and the transit time of individual sickle cells, i.e., red blood cells (RBCs) from sickle cell disease (SCD) patients, while they are traversing the circulatory system. Many drugs for treating SCD work by inhibiting the polymerization of sickle hemoglobin (HbS), effectively delaying the sickling process in sickle cells (SS RBCs). Most previous studies on screening anti-sickling drugs, such as voxelotor, rely on testing of sickling characteristics, often conducted under prolonged deoxygenation for up to 1 hour. However, since the microcirculation of RBCs typically takes less than 1 minute, the results of these studies may be less accurate and less relevant for in vitro-in vivo correlation. In our current study, we introduce a computer vision-enhanced microfluidic framework designed to automatically capture the transient sickling kinetics of SS RBCs within a 1-min timeframe. Our study has successfully detected differences in the transient sickling kinetics between vehicle control and voxelotor-treated SS RBCs. This approach has the potential for broader applications in screening anti-sickling therapies.