Bourns College of Engineering

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

Heterogeneous networks, HetNets, are known to enhance the bandwidth efficiency and throughput of wireless networks by more effectively utilizing the network resources. However, the higher density of users and access points in HetNets introduces significant inter-user interference that needs to be mitigated through complex and sophisticated interference cancellation schemes. Moreover, due to significant channel attenuation and the presence of hardware impairments, e.g. phase noise and amplifier nonlinearities, the vast bandwidth in the millimeterwave band has not been fully utilized to date. In order to enable the development of multi-Gigabit per second wireless networks, we introduce a novel millimeter-wave HetNet paradigm, termed hybrid HetNet, which exploits the vast bandwidth and propagation characteristics in the 60 GHz and 70-80 GHz bands to reduce the impact of interference in HetNets. Simulation results are presented to illustrate the performance advantage of hybrid HetNets with respect to traditional networks. Next, two specific transceiver structures that enable hand-offs from the 60 GHz band, i.e. the V-band to the 70-80 GHz band, i.e. the E-band, and vice versa are proposed. Finally, the practical and regulatory challenges for establishing a hybrid HetNet are outlined.

The antimalarial drug-resistance conundrum which threatens to reverse the great strides taken to curb the malaria scourge warrants an urgent need to find novel chemical scaffolds to serve as templates for the development of new antimalarial drugs. Plants represent a viable alternative source for the discovery of unique potential antiplasmodial chemical scaffolds. To expedite the discovery of new antiplasmodial compounds from plants, the aim of this study was to use phylogenetic analysis to identify higher plant orders and families that can be rationally prioritised for antimalarial drug discovery. We queried the PubMed database for publications documenting antiplasmodial properties of natural compounds isolated from higher plants. Thereafter, we manually collated compounds reported along with plant species of origin and relevant pharmacological data. We systematically assigned antiplasmodial-associated plant species into recognised families and orders, and then computed the resistance index, selectivity index and physicochemical properties of the compounds from each taxonomic group. Correlating the generated phylogenetic trees and the biological data of each clade allowed for the identification of 3 ‘hot’ plant orders and families. The top 3 ranked plant orders were the (i) Caryophyllales, (ii) Buxales, and (iii) Chloranthales. The top 3 ranked plant families were the (i) Ancistrocladaceae, (ii) Simaroubaceae, and (iii) Buxaceae. The highly active natural compounds (IC50 ≤ 1 µM) isolated from these plant orders and families are structurally unique to the ‘legacy’ antimalarial drugs. Our study was able to identify the most prolific taxa at order and family rank that we propose be prioritised in the search for potent, safe and drug-like antimalarial molecules.

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

This review article provides an overview of the current status of velocity-selective arterial spin labeling (VSASL) perfusion MRI and is part of a wider effort arising from the International Society for Magnetic Resonance in Medicine (ISMRM) Perfusion Study Group. Since publication of the 2015 consensus paper on arterial spin labeling (ASL) for cerebral perfusion imaging, important advancements have been made in the field. The ASL community has, therefore, decided to provide an extended perspective on various aspects of technical development and application. Because VSASL has the potential to become a principal ASL method because of its unique advantages over traditional approaches, an in-depth discussion was warranted. VSASL labels blood based on its velocity and creates a magnetic bolus immediately proximal to the microvasculature within the imaging volume. VSASL is, therefore, insensitive to transit delay effects, in contrast to spatially selective pulsed and (pseudo-) continuous ASL approaches. Recent technical developments have improved the robustness and the labeling efficiency of VSASL, making it a potentially more favorable ASL approach in a wide range of applications where transit delay effects are of concern. In this review article, we (1) describe the concepts and theoretical basis of VSASL; (2) describe different variants of VSASL and their implementation; (3) provide recommended parameters and practices for clinical adoption; (4) describe challenges in developing and implementing VSASL; and (5) describe its current applications. As VSASL continues to undergo rapid development, the focus of this review is to summarize the fundamental concepts of VSASL, describe existing VSASL techniques and applications, and provide recommendations to help the clinical community adopt VSASL.

A high-speed single-mode fiber-based polarization-sensitive optical coherence tomography (PS OCT) system was developed. With a polarization modulator, Stokes parameters of reflected flight for four input polarization states are measured as a function of depth. A phase modulator in the reference arm of a Michelson interferometer permits independent control of the axial scan rate and carrier frequency. In vivo PS OCT images of human skin are presented, showing subsurface structures that are not discernible in conventional OCT images. A phase retardation image in tissue is calculated based on the reflected Stokes parameters of the four input polarization states.

Microorganisms, in the most hyperarid deserts around the world, inhabit the inside of rocks as a survival strategy. Water is essential for life, and the ability of a rock substrate to retain water is essential for its habitability. Here we report the mechanism by which gypsum rocks from the Atacama Desert, Chile, provide water for its colonizing microorganisms. We show that the microorganisms can extract water of crystallization (i.e., structurally ordered) from the rock, inducing a phase transformation from gypsum (CaSO₄·2H₂O) to anhydrite (CaSO₄). To investigate and validate the water extraction and phase transformation mechanisms found in the natural geological environment, we cultivated a cyanobacterium isolate on gypsum rock samples under controlled conditions. We found that the cyanobacteria attached onto high surface energy crystal planes ({011}) of gypsum samples generate a thin biofilm that induced mineral dissolution accompanied by water extraction. This process led to a phase transformation to an anhydrous calcium sulfate, anhydrite, which was formed via reprecipitation and subsequent attachment and alignment of nanocrystals. Results in this work not only shed light on how microorganisms can obtain water under severe xeric conditions but also provide insights into potential life in even more extreme environments, such as Mars, as well as offering strategies for advanced water storage methods.

Targeting effectual epitopes is essential for therapeutic antibodies to accomplish their desired biological functions. This study developed a competitive dual color fluorescence-activated cell sorting (FACS) to maturate a matrix metalloprotease 14 (MMP-14) inhibitory antibody. Epitope-specific screening was achieved by selection on MMP-14 during competition with N-terminal domain of tissue inhibitor of metalloproteinase-2 (TIMP-2) (nTIMP-2), a native inhibitor of MMP-14 binding strongly to its catalytic cleft. 3A2 variants with high potency, selectivity, and improved affinity and proteolytic stability were isolated from a random mutagenesis library. Binding kinetics indicated that the affinity improvements were mainly from slower dissociation rates. In vitro degradation tests suggested the isolated variants had half lives 6-11-fold longer than the wt. Inhibition kinetics suggested they were competitive inhibitors which showed excellent selectivity toward MMP-14 over highly homologous MMP-9. Alanine scanning revealed that they bound to the vicinity of MMP-14 catalytic cleft especially residues F204 and F260, suggesting that the desired epitope was maintained during maturation. When converted to immunoglobulin G, B3 showed 5.0 nM binding affinity and 6.5 nM inhibition potency with in vivo half-life of 4.6 days in mice. In addition to protease inhibitory antibodies, the competitive FACS described here can be applied for discovery and engineering biosimilars, and in general for other circumstances where epitope-specific modulation is needed.

Nanoparticles (NPs) dispersed within a conductive host are essential for a range of applications including electrochemical energy storage, catalysis, and energetic devices. However, manufacturing high quality NPs in an efficient manner remains a challenge, especially due to agglomeration during assembly processes. Here we report a rapid thermal shock method to in situ synthesize well-dispersed NPs on a conductive fiber matrix using metal precursor salts. The temperature of the carbon nanofibers (CNFs) coated with metal salts was ramped from room temperature to ∼2000 K in 5 ms, which corresponds to a rate of 400,000 K/s. Metal salts decompose rapidly at such high temperatures and nucleate into metallic nanoparticles during the rapid cooling step (cooling rate of ∼100,000 K/s). The high temperature duration plays a critical role in the size and distribution of the nanoparticles: the faster the process is, the smaller the nanoparticles are, and the narrower the size distribution is. We also demonstrated that the peak temperature of thermal shock can reach ∼3000 K, much higher than the decomposition temperature of many salts, which ensures the possibility of synthesizing various types of nanoparticles. This universal, in situ, high temperature thermal shock method offers considerable potential for the bulk synthesis of unagglomerated nanoparticles stabilized within a matrix.

The antimalarial drug-resistance conundrum which threatens to reverse the great strides taken to curb the malaria scourge warrants an urgent need to find novel chemical scaffolds to serve as templates for the development of new antimalarial drugs. Plants represent a viable alternative source for the discovery of unique potential antiplasmodial chemical scaffolds. To expedite the discovery of new antiplasmodial compounds from plants, the aim of this study was to use phylogenetic analysis to identify higher plant orders and families that can be rationally prioritised for antimalarial drug discovery. We queried the PubMed database for publications documenting antiplasmodial properties of natural compounds isolated from higher plants. Thereafter, we manually collated compounds reported along with plant species of origin and relevant pharmacological data. We systematically assigned antiplasmodial-associated plant species into recognised families and orders, and then computed the resistance index, selectivity index and physicochemical properties of the compounds from each taxonomic group. Correlating the generated phylogenetic trees and the biological data of each clade allowed for the identification of 3 ‘hot’ plant orders and families. The top 3 ranked plant orders were the (i) Caryophyllales, (ii) Buxales, and (iii) Chloranthales. The top 3 ranked plant families were the (i) Ancistrocladaceae, (ii) Simaroubaceae, and (iii) Buxaceae. The highly active natural compounds (IC50 ≤ 1 µM) isolated from these plant orders and families are structurally unique to the ‘legacy’ antimalarial drugs. Our study was able to identify the most prolific taxa at order and family rank that we propose be prioritised in the search for potent, safe and drug-like antimalarial molecules.

MOTIVATION: Some economically important plants including wheat and cotton have more than two copies of each chromosome. With the decreasing cost and increasing read length of next-generation sequencing technologies, reconstructing the multiple haplotypes of a polyploid genome from its sequence reads becomes practical. However, the computational challenge in polyploid haplotyping is much greater than that in diploid haplotyping, and there are few related methods. RESULTS: This article models the polyploid haplotyping problem as an optimal poly-partition problem of the reads, called the Polyploid Balanced Optimal Partition model. For the reads sequenced from a k-ploid genome, the model tries to divide the reads into k groups such that the difference between the reads of the same group is minimized while the difference between the reads of different groups is maximized. When the genotype information is available, the model is extended to the Polyploid Balanced Optimal Partition with Genotype constraint problem. These models are all NP-hard. We propose two heuristic algorithms, H-PoP and H-PoPG, based on dynamic programming and a strategy of limiting the number of intermediate solutions at each iteration, to solve the two models, respectively. Extensive experimental results on simulated and real data show that our algorithms can solve the models effectively, and are much faster and more accurate than the recent state-of-the-art polyploid haplotyping algorithms. The experiments also show that our algorithms can deal with long reads and deep read coverage effectively and accurately. Furthermore, H-PoP might be applied to help determine the ploidy of an organism. AVAILABILITY AND IMPLEMENTATION: https://github.com/MinzhuXie/H-PoPG CONTACT: xieminzhu@hotmail.comSupplementary information: Supplementary data are available at Bioinformatics online.

Molecular markers are increasingly being deployed to accelerate genetic gain in crop plants. The objective of this study was to assess the potential of a mid-density genotyping panel for molecular applications in cowpea breeding. A core set of 2,602 targeted diversity array technology (DArTag) single-nucleotide polymorphisms (SNPs) was designed from an existing 51,128 Cowpea iSelect Consortium Array. The panel's usefulness was assessed using 376 genotypes from different populations of known genetic backgrounds. The panel was informative, with over 78% of SNPs exceeding a minor allele frequency of 0.20. The panel decoded three stratifications in the constituted population, as was expected. Linkage disequilibrium (LD) decay was correctly depicted as slower in a biparental subset than in other populations. A known flower and seed coat color gene region was located on chromosome Vu07, suggesting that the mid-density panel may be used to hypothesize genomic regions underlying target traits in cowpea. Unexpected heterozygosity was detected in some lines and highly among F₁ progenies, divulging the panel's potential application in germplasm purity and hybridity verification. The study unveils the potential of an excellent genomic resource that can be tapped to enhance the development of improved cowpea cultivars.

Describing visual image contents by semantic concepts is an effective and straightforward way to facilitate various high level applications. Inferring semantic concepts from low-level pictorial feature analysis is challenging due to the semantic gap problem, while manually labeling concepts is unwise because of a large number of images in both online and offline collections. In this paper, we present a novel approach to automatically generate intermediate image descriptors by exploiting concept co-occurrence patterns in the pre-labeled training set that renders it possible to depict complex scene images semantically. Our work is motivated by the fact that multiple concepts that frequently co-occur across images form patterns which could provide contextual cues for individual concept inference. We discover the co-occurrence patterns as hierarchical communities by graph modularity maximization in a network with nodes and edges representing concepts and co-occurrence relationships separately. A random walk process working on the inferred concept probabilities with the discovered co-occurrence patterns is applied to acquire the refined concept signature representation. Through experiments in automatic image annotation and semantic image retrieval on several challenging datasets, we demonstrate the effectiveness of the proposed concept co-occurrence patterns as well as the concept signature representation in comparison with state-of-the-art approaches.

We propose a context guided belief propagation (BP) algorithm to perform high spatial resolution multispectral imagery (HSRMI) classification efficiently utilizing superpixel representation. One important characteristic of HSRMI is that different land cover objects possess a similar spectral property. This property is exploited to speed up the standard BP (SBP) in the classification process. Specifically, we leverage this property of HSRMI as context information to guide messages passing in SBP. Furthermore, the spectral and structural features extracted at the superpixel level are fed into a Markov random field framework to address the challenge of low interclass variation in HSRMI classification by minimizing the discrete energy through context guided BP (CBP). Experiments show that the proposed CBP is significantly faster than the SBP while retaining similar performance as compared with SBP. Compared to the baseline methods, higher classification accuracy is achieved by the proposed CBP when the context information is used with both spectral and structural features.

In plants, the robust maintenance of tissue structure is crucial to supporting its functionality. The multi-layered shoot apical meristem (SAM) of Arabidopsis, containing stem cells, is an approximately radially symmetric tissue whose shape and structure is maintained throughout the life of the plant. In this paper, a new biologically calibrated pseudo-three-dimensional (P3D) computational model of a longitudinal section of the SAM is developed. It includes anisotropic expansion and division of cells out of the cross-section plane, as well as representation of tension experienced by the SAM epidermis. Results from the experimentally calibrated P3D model provide new insights into maintenance of the structure of the SAM epidermal cell monolayer under tension and quantify dependence of epidermal and subepidermal cell anisotropy on the amount of tension. Moreover, the model simulations revealed that out-of-plane cell growth is important in offsetting cell crowding and regulating mechanical stresses experienced by tunica cells. Predictive model simulations show that tension-determined cell division plane orientation in the apical corpus may be regulating cell and tissue shape distributions needed for maintaining structure of the wild-type SAM. This suggests that cells' responses to local mechanical cues may serve as a mechanism to regulate cell- and tissue-scale patterning.

Three-dimensional variation in structural components or fiber alignments results in complex mechanical property distribution in tissues and biomaterials. In this paper, we use a physics-informed UNet-based neural network model (El-UNet) to discover the three-dimensional (3D) internal composition and space-dependent material properties of heterogeneous isotropic and transversely isotropic materials without a priori knowledge of the composition. We then show the capabilities of El-UNet by validating against data obtained from finite-element simulations of two soft tissues, namely, brain tissue and articular cartilage, under various loading conditions. We first simulated compressive loading of 3D brain tissue comprising of distinct white matter and gray matter mechanical properties undergoing small strains with isotropic linear elastic behavior, where El-UNet reached mean absolute relative errors under 1.5 % for elastic modulus and Poisson's ratio estimations across the 3D volume. We showed that the 3D solution achieved by El-UNet was superior to relative stiffness mapping by inverse of axial strain and two-dimensional plane stress/plane strain approximations. Additionally, we simulated a transversely isotropic articular cartilage with known fiber orientations undergoing compressive loading, and accurately estimated the spatial distribution of all five material parameters, with mean absolute relative errors under 5 %. Our work demonstrates the application of the computationally efficient physics-informed El-UNet in 3D elasticity imaging and provides methods for translation to experimental 3D characterization of soft tissues and other materials. The proposed El-UNet offers a powerful tool for both in vitro and ex vivo tissue analysis, with potential extensions to in vivo diagnostics. STATEMENT OF SIGNIFICANCE: Elasticity imaging is a technique that reconstructs mechanical properties of tissue using deformation and force measurements. Given the complexity of this reconstruction, most existing methods have mostly focused on 2D problems. Our work is the first implementation of physics-informed UNets to reconstruct three-dimensional material parameter distributions for isotropic and transversely isotropic linear elastic materials by having deformation and force measurements. We comprehensively validate our model using synthetic data generated using finite element models of biological tissues with high bio-fidelity-the brain and articular cartilage. Our method can be implemented in elasticity imaging scenarios for in vitro and ex vivo mechanical characterization of biomaterials and biological tissues, with potential extensions to in vivo diagnostics.

In contrast to inert systems, living biological systems have the advantage to adapt to their environment through growth and evolution. This transfiguration is evident during embryonic development, when the predisposed need to grow allows form to follow function. Alterations in the equilibrium state of biological systems breed disease and mutation in response to environmental triggers. The need to characterize the growth of biological systems to better understand these phenomena has motivated the continuum theory of growth and stimulated the development of computational tools in systems biology. Biological growth in development and disease is increasingly studied using the framework of morphoelasticity. Here, we demonstrate the potential for morphoelastic simulations through examples of volume, area, and length growth, inspired by tumor expansion, chronic bronchitis, brain development, intestine formation, plant shape, and myopia. We review the systems biology of living systems in light of biochemical and optical stimuli and classify different types of growth to facilitate the design of growth models for various biological systems within this generic framework. Exploring the systems biology of growth introduces a new venue to control and manipulate embryonic development, disease progression, and clinical intervention.

Controlling the interactions between macrophages and biomaterials is critical for modulating the response to implants. While it has long been thought that biomaterial surface chemistry regulates the immune response, recent studies have suggested that material geometry may in fact dominate. Our previous work demonstrated that elongation of macrophages regulates their polarization toward a pro-healing phenotype. In this work, we elucidate how surface topology might be leveraged to alter macrophage cell morphology and polarization state. Using a deep etch technique, we fabricated titanium surfaces containing micro- and nanopatterned grooves, which have been previously shown to promote cell elongation. Morphology, phenotypic markers, and cytokine secretion of murine bone marrow derived macrophages on different groove widths were analyzed. The results suggest that micro- and nanopatterned grooves influenced macrophage elongation, which peaked on substrates with 400-500 nm wide grooves. Surface grooves did not affect inflammatory activation but drove macrophages toward an anti-inflammatory, pro-healing phenotype. While secretion of TNF-alpha remained low in macrophages across all conditions, macrophages secreted significantly higher levels of anti-inflammatory cytokine, IL-10, on intermediate groove widths compared to cells on other Ti surfaces. Our findings highlight the potential of using surface topography to regulate macrophage function, and thus control the wound healing and tissue repair response to biomaterials.