Acquisition Of Functions Research Articles

Bayesian optimization for personalized dose-finding trials with combination therapies

Abstract Identification of optimal dose combinations in early-phase dose-finding trials is challenging, due to the trade-off between precisely estimating the many parameters required to flexibly model the possibly nonmonotonic dose-response surface, and the small sample sizes in early-phase trials. This difficulty is even more pertinent in the context of personalized dose-finding, where patient characteristics are used to identify tailored optimal dose combinations. To overcome these challenges, we propose the use of Bayesian optimization for finding optimal dose combinations in standard (one size fits all) and personalized multi-agent dose-finding trials. Bayesian optimization is a method for estimating the global optima of expensive-to-evaluate objective functions. The objective function is approximated by a surrogate model, commonly a Gaussian process, paired with a sequential design strategy to select the next point via an acquisition function. This work is motivated by an industry-sponsored problem, where the focus is on optimizing a dual-agent therapy in a setting featuring minimal toxicity. To compare the performance of the standard and personalized methods under this setting, simulation studies are performed for a variety of scenarios. Our study concludes that taking a personalized approach is highly beneficial in the presence of heterogeneity.

Diatoms exhibit dynamic chloroplast calcium signals in response to high light and oxidative stress

Abstract Diatoms are a group of silicified algae that play a major role in marine and freshwater ecosystems. Diatom chloroplasts were acquired by secondary endosymbiosis and exhibit important structural and functional differences from the primary plastids of land plants and green algae. Many functions of primary plastids, including photoacclimation and inorganic carbon acquisition, are regulated by calcium-dependent signalling processes. Calcium signalling has also been implicated in the photoprotective responses of diatoms; however, the nature of calcium elevations in diatom chloroplasts and their wider role in cell signalling remains unknown. Using genetically encoded calcium indicators, we find that the diatom Phaeodactylum tricornutum exhibits dynamic calcium elevations within the chloroplast stroma. Stromal calcium ([Ca2+]str) acts independently from the cytosol and is not elevated by stimuli that induce large cytosolic calcium ([Ca2+]cyt) elevations. In contrast, high light and exogenous hydrogen peroxide (H2O2) induce large, sustained [Ca2+]str elevations that are not replicated in the cytosol. Measurements using the fluorescent H2O2 sensor roGFP2-Oxidant Receptor Peroxidase 1 (Orp1) indicate that [Ca2+]str elevations induced by these stimuli correspond to the accumulation of H2O2 in the chloroplast. [Ca2+]str elevations were also induced by adding methyl viologen, which generates superoxide within the chloroplast, and by treatments that disrupt non-photochemical quenching (NPQ). The findings indicate that diatoms generate specific [Ca2+]str elevations in response to high light and oxidative stress that likely modulate the activity of calcium-sensitive components in photoprotection and other regulatory pathways.

A study of the role of data and model uncertainty in active learning

Uncertainty-based active learning strategies have demonstrated significant superiority in small data research of materials domain. This study explores the effects of model uncertainty and data uncertainty separately on the performance of active learning strategies, specifically focusing on the number of iterations required to identify the optimal samples. For model uncertainty, three kinds of acquisition functions are compared, including predicted value strategy (PV), ranking of predicted value strategy (PR) and expected improvement strategy (EI). Among these, the active learning model utilizing PR requires the fewest average iterations (1.75). For data uncertainty, we evaluate the iterations of active learning by Gaussian process models that incorporate the uncertainty of the observations and noise samples that takes account into the uncertainty of the input features respectively. The results indicate that the active learning iterations of the three strategies converge to similar at the optimal weighting when the uncertainty of the observations is considered in the model (EI for 1.75, PV for 1.21 and PR for 1.18). In contrast, incorporating noise samples into the augmented dataset after the original samples would severely deteriorate the efficiency of active learning recommendations. Our findings aim to offer guidance for exploring more favorable acquisition functions and methods for active learning strategies.

Study on the Function of Leptin Nutrient Acquisition and Energy Metabolism of Zebrafish (Danio rerio)

Leptin plays an indispensable role in energy homeostasis, and its involvement in metabolic activities has been extensively explored in fish. We generated mutant lines of leptina (−5 bp) and leptinb (+8 bp) in zebrafish using CRISPR/Cas9 technology to explore the metabolic characteristics of lepa and lepb mutant zebrafish in response to high glucose nutritional stress induced by high levels of carbohydrates. The results were as follows: the body weight and food intake of adult zebrafish of the two mutant species were increased; the visceral fat accumulation, whole-body crude lipid, and crude protein contents of lepb−/− were increased; and the visceral fat accumulation and crude lipid in lepa−/− zebrafish were decreased. The blood glucose levels of the two mutant zebrafish were increased, the mRNA expression levels of glycolytic genes pk and gck were decreased in the two mutant zebrafish, and there were differences between lepa−/− and lepb−/− zebrafish. The expressions of glycogen synthesis and decomposition genes were inhibited and promoted, respectively. The expression of adipose synthesis genes in the liver and muscle was stimulated in lepb−/− zebrafish but suppressed in lepa−/− zebrafish. Lipolysis and oxidation genes were also stimulated in lepa−/− zebrafish livers, while the livers of lepb−/− zebrafish were stimulated but muscle was inhibited. In conclusion, the results indicate that lepa plays a major role in glucose metabolism, which is conducive to promoting glucose utilization and lipogenesis, while lepb mainly promotes lipolysis and oxidation, regulates protein generation, and plays a minor role in glucose metabolism.

Study on non-destructive visualization identification of concrete internal cracks based on SH array ultrasound

Ultrasonic nondestructive testing (NDT) technology was one of the most frequently used and fastest-developing NDT techniques in detection field. This method had high detection sensitivity, more accurate defect localization, and was harmless to the human body. With the development of ultrasonic testing technology, modern ultrasonic imaging technology generally had automatic data acquisition and processing functions, which allowed for intuitive and clear recognition of internal defects by the human eyes. The overall purpose of this study was to visualize the cracks in concrete by synthetic aperture focusing technique (SAFT) and full-waveform inversion (FWI) algorithms on the basis of ultrasonic nondestructive testing technology. At present, there have been studies on the use of ultrasonic instruments to detect the internal defects of concrete. However, the accuracy and authenticity of this method have not been compared in previous studies. Therefore, in this experiment, four concrete specimens with different buried depths, angles, shapes, and thicknesses of cracks were designed. Based on SH ultrasonic wave NDT technology, two types of operations were performed on the collected ultrasonic wave data: one was the application of SAFT method, and the other was the first application of FWI method for detecting internal cracks in concrete. The study found that both algorithms were most accurate in determining the burial depth of cracks, with error results not exceeding 10 %. They also showed good recognition effects for the angle, shape, and thickness of the buried cracks. This proved that both identification methods had the capability to recognize internal cracks. The Introview software bundled with the ultrasonic instrument could identify cracks using SAFT directly and quickly. When FWI was used with MATLAB computation, the results were more precise and intuitive. The overall purpose of this study was to detect cracks in concrete by SAFT and FWI algorithms on the basis of ultrasonic nondestructive testing technology. Combining the two crack identification techniques could be of great significance for ensuring the long-term safety and stability of concrete structures, providing strong support for structural maintenance and reinforcement.

Expected coordinate improvement for high-dimensional Bayesian optimization

Bayesian optimization (BO) algorithm is very popular for solving low-dimensional expensive optimization problems. Extending Bayesian optimization to high dimension is a meaningful but challenging task. One of the major challenges is that it is difficult to find good infill solutions as the acquisition functions are also high-dimensional. In this work, we propose the expected coordinate improvement (ECI) criterion for high-dimensional Bayesian optimization. The proposed ECI criterion measures the potential improvement we can get by moving the current best solution along one coordinate. The proposed approach selects the coordinate with the highest ECI value to refine in each iteration and covers all the coordinates gradually by iterating over the coordinates. The greatest advantage of the proposed ECI-BO (expected coordinate improvement based Bayesian optimization) algorithm over the standard BO algorithm is that the infill selection problem of the proposed algorithm is always a one-dimensional problem thus can be easily solved. Numerical experiments show that the proposed algorithm can achieve significantly better results than the standard BO algorithm and competitive results when compared with five high-dimensional BOs and six surrogate-assisted evolutionary algorithms. This work provides a simple but efficient approach for high-dimensional Bayesian optimization. A Matlab implementation of our ECI-BO is available at https://github.com/zhandawei/Expected_Coordinate_Improvement.

Bayesian Optimization Using Simulation-Based Multiple Information Sources over Combinatorial Structures

Bayesian optimization due to its flexibility and sample efficiency has become a standard approach for simulation optimization. To reduce this problem, one can resort to cheaper surrogates of the objective function. Examples are ubiquitous, from protein engineering or material science to tuning machine learning algorithms, where one could use a subset of the full training set or even a smaller related dataset. Cheap information sources in the optimization scheme have been studied in the literature as the multi-fidelity optimization problem. Of course, cheaper sources may hold some promise toward tractability, but cheaper models offer an incomplete model inducing unknown bias and epistemic uncertainty. In this manuscript, we are concerned with the discrete case, where fx,wi is the value of the performance measure associated with the environmental condition wi and p(wi) represents the relevance of the condition wi (i.e., the probability of occurrence or the fraction of time this condition occurs). The main contribution of this paper is the proposal of a Gaussian-based framework, called augmented Gaussian process (AGP), based on sparsification, originally proposed for continuous functions and its generalization in this paper to stochastic optimization using different risk profiles for combinatorial optimization. The AGP leverages sample and cost-efficient Bayesian optimization (BO) of multiple information sources and supports a new acquisition function to select the new source–location pair considering the cost of the source and the (location-dependent) model discrepancy. An extensive set of computational results supports risk-aware optimization based on CVaR (conditional value-at-risk). Computational experiments confirm the actual performance of the MISO-AGP method and the hyperparameter optimization on benchmark functions and real-world problems.

Toward digitally twinning the process of creating machine controller digital twins – A G-code generation scenario

This paper proposes BatchEGO, an extension of Efficient Global Optimization (EGO), in an active learning context for creating and maintaining digital twins (DT). The motivation for this study is to address the challenges of obtaining high-quality yet cost and time effective data for developing accurate DTs of machine tool controllers that are robust to noise and physical system changes. BatchEGO is employed to generate manufacturing recipes (G-codes) recursively and iteratively, based on the spatial constraints of the physical machine, providing G-code based data for DTs. BatchEGO Interm is also proposed that extends BatchEGO by sampling additional intermediate points on the generated trajectory, to improve convergence and accuracy of DTs. To account for the stochastic nature of proposed methods, a total of 70 experiments based around repetitions of 4 correlation functions, 3 acquisition functions, and 4 sampling schemes are used to evaluate the performance of both methods. The results showcase impressive and faster convergence, in terms of sum of squared errors, for BatchEGO Interm with specific settings producing highly accurate and smooth G-codes. For instance, BatchEGO Interm with Matérn 3/2 correlation function and acquisition function of lower confidence bound can reach an error range of 0.0026 within 3-5 iterations. The work demonstrates the novelty and potential of such frameworks for enhancing data quality and quantity for DTs. The importance of G-code generation for DTs is also highlighted, as it can enable better control and optimization of machine tools in various industrial applications. The paper also discusses the different approaches for stitching the sampled points into G-code and compares their accuracy and efficiency.

Integrating small data and shape prior knowledge with gradient-enhanced Kriging through adaptive knowledge sampling

Integrating prior knowledge into surrogate models constitutes an advanced approach to tackling the challenges posed by small data. The existing gradient-enhanced Kriging (GEK) method can enhance accuracy by incorporating gradient points of shape prior knowledge. Nonetheless, our findings indicate that the strategy employed for setting gradient points significantly impacts accuracy, even when the number of gradient points remains constant. To enhance the efficient utilization of gradient information, we propose a novel Adaptive Knowledge Sampling (AKS) strategy with an innovative acquisition function that balances local exploitation and global exploration, effectively extracting gradient points from shape prior knowledge. The acquisition function consists of three parts: (1) The Local Information Filling Criterion (LIFC), designed to mitigate model error; (2) The Global Information Filling Criterion (GIFC), aimed at reducing model uncertainty in data-sparse areas; and (3) The Minimum Distance Constraint (MDC), implemented to prevent excessive clustering of gradient points. To strike a balance between model accuracy and computational cost, an adaptive stopping condition is introduced to regulate the number of gradient points sampled. The impact of hyperparameters in the AKS strategy was analyzed in detail, and the method’s performance was compared to existing approaches using six benchmark functions and a missile aerodynamic prediction case. The results demonstrate that AKS-GEK offers higher prediction accuracy and robustness when integrating small data with shape prior knowledge.

Design and Implementation of Image Acquisition and Signal Transmission Functions based on FPGA

Abstract In recent years, video network monitoring systems have evolved rapidly, from the first generation of digital-analogue systems to the fully IP-based systems. Based on the FPGA core, this design focuses on image acquisition and information transmission. The FPGA expansion board contains several peripheral modules, including the OV7670 camera captures images, DRAM stores data, and FIFO handles cross-clock domain signal transmission. While the reference clock is 40ns, in this system, the clock period is 20ns, making two clock cycles equal to one reference clock cycle. Then, two counters are used, one is the reference clock and another is the line synchronization (cnt1) to count 800 clock cycles per line. The end of a line is indicated by end_cnt0, and FIFO data is stored in SDRAM. Configuring the VGA display driver allows image data to be transmitted and displayed on the VGA screen. The system features a simple structure, high resource utilization, and strong reliability.

Exploring intraspecific and interspecific variation of coral reef algae using a novel trait‐based framework

Abstract The development of trait‐based approaches has accelerated our understanding of how communities assemble, respond to environmental change and may best be managed in the Anthropocene. Understanding the magnitude and pattern of interspecific variability forms a critical underpinning of trait‐based approaches while exploring intraspecific variability can identify the potential of species to adapt to changing environmental drivers. Our work is motivated by the critical need for a novel conceptual framework for understanding the functional ecology of macroalgae, as the current paradigm is still mired in functional group models developed in the 1980s. Our objective was to quantify interspecific and intraspecific functional trait variability in three common and morphologically diverse species of tropical marine macroalgae by exploring traits relating to the ecological functions of resource acquisition, resistance to herbivory, and resistance to physical disturbance and the trade‐offs between them. We quantified intraspecific and interspecific variability of 11 functional traits for three common and morphologically diverse species of tropical macroalgae from five fringing reefs of Mo'orea, French Polynesia that were likely to capture a wide range of environmental variability. Differences in traits among species and sites were determined with PERMANOVA, visualised with non‐metric multidimensional scaling, and trade‐offs between pairs of traits explored with correlation. Finally, spatial patterns among select traits across all species were quantified. Species clustered together in distinctly different trait spaces driven by trade‐offs among suites of functional traits. Two of three species had considerable intraspecific variability, though this variability occurred at different scales, while one clustered tightly. Exploration of individual traits across species and sites revealed trade‐offs between two strategies for resource acquisition, growing tall and strong versus investing in large surface area. Synthesis. We captured novel patterns of interspecific and intraspecific variability for tropical marine macroalgae. We found fundamental differences in traits between species that may represent ecological strategies while considerable intraspecific variability demonstrates a wide range in abilities to respond to environmental drivers. Overall, our work provides novel insights into intra and interspecific trait variation that form an essential underpinning for using a trait‐based framework in a taxon that is increasingly dominant on tropical reefs.

Integrating adaptive learning with post hoc model explanation and symbolic regression to build interpretable surrogate models

Abstract We develop a materials informatics workflow to build an interpretable surrogate model for micromagnetic simulations. Our goal is to predict the energy barrier of a moving isolated skyrmion in rare-earth-free $$\hbox {Mn}_4$$ Mn 4 N. Our approach integrates adaptive learning with post hoc model explanation and symbolic regression methods. We discuss an unexplored acquisition function (information condensing active learning) within the adaptive learning loop and compare it with the known standard deviation function for efficient navigation of the search space. Model-agnostic post hoc explanation techniques then uncover trends learned by the trained model, which we then leverage to constrain the expressions used for symbolic regression. Graphical abstract

Bayesian optimization acquisition functions for accelerated search of cluster expansion convex hull of multi-component alloys

Atomistic simulations are crucial for predicting material properties and understanding phase stability, essential for materials selection and development. However, the high computational cost of density functional theory calculations challenges the design of materials with complex structures and composition. This study introduces new data acquisition strategies using Bayesian-Gaussian optimization that efficiently integrate the geometry of the convex hull to optimize the yield of batch experiments. We developed uncertainty-based acquisition functions to prioritize the computation tasks of configurations of multi-component alloys, enhancing our ability to identify the ground-state line. Our methods were validated across diverse materials systems including Co-Ni alloys, Zr-O compounds, Ni-Al-Cr ternary alloys, and a planar defect system in intermetallic (Ni1−x, Cox)3Al. Compared to traditional genetic algorithms, our strategies reduce training parameters and user interaction, cutting the number of experiments needed to accurately determine the ground-state line by over 30%. These approaches can be expanded to multi-component systems and integrated with cost functions to further optimize experimental designs.

Enhancing solar energy conversion efficiency: Thermophysical property predicting of MXene/Graphene hybrid nanofluids via bayesian-optimized artificial neural networks

Accurately predicting thermo-physical properties (TPPs) of MXene/graphene-based nanofluids is crucial for photovoltaic/thermal solar systems, driving focused research on developing precise TPP predictive models. This study presents optimized multi-layer perceptron neural network (MLPNN) models, leveraging Bayesian optimization to refine architectural and training hyperparameters, including hidden layers, neurons, activation functions, standardization, and regularization terms. A comparative analysis of Bayesian acquisition functions—the probability of improvement (POI), lower confidence bound (LCB), expected improvement (EI), expected improvement plus (EIP), expected improvement per second plus (EIPSP), and expected improvement per second (EIPS)—demonstrated that the POI-MLPNN achieves the most accurate results, as evidenced by the lowest MAPE of 1.0923 % and exceptional consistency with an R-value of 0.99811. The EI-MLPNN and EIP-MLPNN models recorded the same outputs. The EI/EIP-MLPNN (R = 0.99668) model excels in consistency over LCB-MLPNN (R = 0.99529) and EIPSP-MLPNN (R = 0.99667). The optimized models offer a reliable, cost-efficient alternate for experimental and computational TPP analyses. Leveraging insights from these models enables better control over nanofluid TPPs in solar systems, enhancing energy conversion efficiency.

BO-GRNN: Machine Learning for Underwater Acoustic Source Localization

This article proposes a Bayesian optimization-based generalized regression neural network (BO-GRNN) framework for underwater source localization in shallow water environment. Specially, by leveraging kernel regression analysis, BO-GRNN fits a conditional probability density function using training and test samples to infer the source location of unseen test samples. Furthermore, BO-GRNN models the objective function as a Gaussian process (GP) based on observed samples, utilizing an acquisition function to sample around the best objective function values and high variance regions. By evaluating these samples and updating the GP iteratively, the optimal hyperparameters are identified, ultimately improving the training efficiency and localization performance of the GRNN. Simulation experiments and error analysis validate this framework’s effectiveness. Further validation using SWellEx-96 experimental data demonstrates the superior localization performance of this framework, even with limited data snapshots. This framework estimates the range and depth of underwater sources within a 10 km range, achieving mean absolute errors of less than 0.15 km for range and 0.2 m for depth. Compared to GRNN, this framework demonstrates a 5 times improvement in training efficiency. These advantages make BO-GRNN an efficient, accurate, and robust solution for underwater source localization.

Research on key circuits of high‐precision servo channel loading card for wing fatigue testing

AbstractThe key circuits of servo channel loading card with high‐precision data acquisition and output functions for wing fatigue testing are studied and improved in this paper. To solve the problem of large offset voltage in some op‐amps, non‐inverting combined amplifier circuit (NCAC) and inverting combination amplifier circuit (ICAC) are proposed by utilizing precision op‐amps combined with wide voltage and high‐power output op‐amps. In the force signal acquisition section, a non‐inverting combined zeroing filter circuit (NCZFC) is proposed to improve the precision of signal acquisition. In the excitation source output section, a non‐inverting combined excitation output circuit (NCEOC) is proposed to improve the precision of the positive excitation output; an inverting combined excitation output circuit (ICEOC) is also proposed, which improves the precision of the negative excitation output and also realizes the inversion of the input and output signals. In the servo drive output section, a non‐inverting combined voltage and current conversion circuit is proposed to improve the precision of voltage and current outputs, which has two modes of non‐inverting combined voltage output (NCVO) and non‐inverting combined current output (NCIO) modes. Simulation experiments in Multisim and actual tests have been carried out on the above improved circuits to verify the stability and precision.

Toward global rooftop PV detection with Deep Active Learning

It is crucial to know the location of rooftop PV systems to monitor the regional progress toward sustainable societies and to ensure the integration of decentralized energy resources into the electricity grid. However, locations of PV are often unknown, which is why a large number of studies have proposed variants of Deep Learning to detect PV panels in remote sensing data using supervised Deep Learning. However, these methods are based on annotating datasets and therefore often require relabeling when fine-tuned or extended to a different region. Recent advances in Deep Active Learning offer the opportunity to significantly reduce the number of required annotated images by intelligently selecting the images to label next based on their informative value for the model. In this study, we compare different Deep Active Learning algorithms using a variety of datasets from different regions and compare different model training variants. In the simulations, the entropy-based acquisition function shows the highest performance with only 3% of the data needed in case-imbalanced data, while remaining simple to implement. We believe that Deep Active Learning provides an elegant solution to maintain high model accuracy while reducing annotation effort substantially. This facilitates the development of generalizable models for worldwide rooftop PV detection.

Targeted nano-energetic material exploration through active learning algorithm implementation

This paper presents an active learning-based method designed to guide experiments towards user-defined specific regions, termed ”regions of interest,” within vast and multi-dimensional thermite design spaces. Thermites composed of metallic reactant coupled to an inorganic oxidizer are non-explosive energetic materials which stay inert and stable until subjected to a sufficiently strong thermal stimulus, after which they undergo fast burning with release of high amount of chemical energy (up to 16 kJ.cm−3). They represent an interesting class of nano-engineered energetic material because of their high adiabatic flame temperature (> 2600 °C) and customizable combustion properties. We introduced a new acquisition function combining linearly two factors with the usual standard deviation of a Gaussian Process Regression algorithm. A first factor guides the sampling towards the defined zone of interest, and a second, is an incentive function that encourages the exploration of under-sampled regions of the design space. We found that our algorithm effectively provides up to several tens of nanothermites that achieve specific desired properties well-distributed within the design space, after only 200 samplings, whereas, Latin Hypercube Sampling procedure samples less than 10 points of interest. Our framework is tailored for typical discrete search spaces involving multiple measured physical properties and when only a small dataset is available, as it is the case in thermite materials. But more generally, this research represents a significant advancement for large-scale problems in energetic materials science and engineering, where predicting the effect of feature modification is desired but limited simulations or experiments can be afforded due to their high cost.

Development of a wireless multichannel miniature impedance measurement system and its application for bolt loosening detection

The piezoelectric impedance-based technique is always regarded as one of the most promising structural health monitoring and nondestructive evaluation methods. In recent years, impedance measurement chip AD5933 with the characteristics of high integration and cost-effectiveness makes it possible to address the huge and high-cost problems of the commercialized impedance measurement instrument during the process of structural health monitoring and defect identification. However, it still faces lots of challenges for the chips to be utilized in practical applications due to several limitations, such as short distance for data transmission, single measurement channel and artificial attendant requirement. In this paper, a wireless multichannel miniature impedance measurement system composed by the front-end measurement device and the remote measurement and control platform on the server, is firstly developed and presented with the functions of wireless data transmission, multi-channel acquisition and remote data post-processing. Subsequently, the design concept and composition of the system were introduced in detail. Then, a series of piezoelectric transducers related tests were conducted to validate its impedance measurement performance, especially when comparing with the ones measured by commercialized instruments. In addition, to verify its effectiveness and feasibility of the developed system for the structural damage detection, a bolt loosening detection experiment on the flange connection of a pipeline specimen was investigated for its damage localization and severity quantification. Finally, all the results demonstrated that the developed system provides a great possibility to be used as a convenient and portable impedance measurement tool for the civil structural health monitoring and damage identification in practical applications.

Constraint based Bayesian optimization of bioink precursor: a machine learning framework.

Current research practice for optimizing bioink involves exhaustive experimentation with multi-material composition for determining the printability, shape fidelity and biocompatibility. Predicting bioink properties can be beneficial to the research community but is a challenging task due to the non-Newtonian behavior in complex composition. Existing models such as Cross model become inadequate for predicting the viscosity for heterogeneous composition of bioinks. In this paper, we utilize a machine learning framework to accurately predict the viscosity of heterogeneous bioink compositions, aiming to enhance extrusion-based bioprinting techniques. Utilizing Bayesian optimization (BO), our strategy leverages a limited dataset to inform our model. This is a technique especially useful of the typically sparse data in this domain. Moreover, we have also developed a mask technique that can handle complex constraints, informed by domain expertise, to define the feasible parameter space for the components of the bioink and their interactions. Our proposed method is focused on predicting the intrinsic factor (e.g. viscosity) of the bioink precursor which is tied to the extrinsic property (e.g. cell viability) through the mask function. Through the optimization of the hyperparameter, we strike a balance between exploration of new possibilities and exploitation of known data, a balance crucial for refining our acquisition function. This function then guides the selection of subsequent sampling points within the defined viable space and the process continues until convergence is achieved, indicating that the model has sufficiently explored the parameter space and identified the optimal or near-optimal solutions. Employing this AI-guided BO framework, we have developed, tested, and validated a surrogate model for determining the viscosity of heterogeneous bioink compositions. This data-driven approach significantly reduces the experimental workload required to identify bioink compositions conducive to functional tissue growth. It not only streamlines the process of finding the optimal bioink compositions from a vast array of heterogeneous options but also offers a promising avenue for accelerating advancements in tissue engineering by minimizing the need for extensive experimental trials.