Relative Contributions of Solubility and Mobility to the Stability of Amorphous Solid Dispersions of Poorly Soluble Drugs: A Molecular Dynamics Simulation Study

Predicting Drug-Polymer Compatibility in Amorphous Solid Dispersions by MD Simulation: On the Trap of Solvation Free Energies

Amorphous solid dispersions (ASDs) are a prevalent method for increasing the bioavailability and apparent solubility of poorly soluble drugs. Consequently, extensive research, encompassing both experimental and computational approaches, has been dedicated to developing methods for assessing the key factors influencing their stability, notably drug-polymer interactions. A common computational approach to rank the compatibility of a drug with a set of solvents or polymers is to compare thermodynamic observables, such as solvation free energies at infinite dilution. However, the impact of the molecular weight of the polymer excipient on these interactions remains underexplored. This study delves into this impact through atomistic simulations of Indomethacin in PVP(-VA) and HPMC, and through simulations using a coarse-grained model, emphasizing its critical importance. First, we demonstrate that the molecular weight of the polymer plays a pivotal role in determining the solvation free energy of the drug, at times exerting a more significant influence than the specific chemical identity of the polymer. Additionally, our simulations suggest that higher molecular weight polymers lead to lower solvation free energies and, thus, suggest better compatibility with the drug. Yet, the lower free energy of solvation of the drug in longer polymers does not translate into a higher solubility. This work highlights the subtle role polymer molecular weight plays when measuring thermodynamic observables in amorphous solid dispersions, a role which must be considered when optimizing pharmaceutical formulations.

Quantitative Analysis of Physical Stability Mechanisms of Amorphous Solid Dispersions by Molecular Dynamic Simulation

Amorphous solid dispersions (ASDs) represent a promising strategy for enhancing the solubility of poorly soluble drugs. However, the mechanisms underlying the physical stability of ASDs remain insufficiently understood. This study aims to investigate these mechanisms and propose quantitative thresholds to predict the maximum stable drug loading using molecular dynamics simulations. Poly(vinylpyrrolidone) (PVP) and poly (vinylpyrrolidone-co-vinyl acetate) (PVPVA64) are selected as polymeric carriers, while naproxen and acetaminophen serve as model drugs, resulting in the formulation of 18 distinct ASDs across four types for comparison with experimental results. Our findings indicate that the molecular mobility of active pharmaceutical ingredients (APIs) is the primary determinant of solid dispersion stability. High polymer concentrations limit drug molecular mobility through spatial structural constraints and ASD viscosity. As drug loading increases, the polymer concentration reaches a critical threshold (C*), beyond which drug-rich regions form, leading to potential aggregation, rearrangement, and recrystallization of drug molecules into more energetically stable forms. Notably, both the interaction energy and diffusion coefficient show sharp fluctuations at the maximum stable drug loading, which can serve as predictive indicators for ASD stability. Additionally, a search strategy is used to identify potential pre-crystalline sites. By integrating kinetic, thermodynamic, and pre-crystalline analyses through molecular dynamics simulations, this study provides a foundation for more accurate predictions of ASD stability, significantly aiding future formulation development.

Molecular simulation for food protein-ligand interactions: A comprehensive review on principles, current applications, and emerging trends

In recent years, investigations on molecular interaction mechanisms between food proteins and ligands have attracted much interest. The interaction mechanisms can supply much useful information for many fields in the food industry, including nutrient delivery, food processing, auxiliary detection, and others. Molecular simulation has offered extraordinary insights into the interaction mechanisms. It can reflect binding conformation, interaction forces, binding affinity, key residues, and other information that physicochemical experiments cannot reveal in a fast and detailed manner. The simulation results have proven to be consistent with the results of physicochemical experiments. Molecular simulation holds great potential for future applications in the field of food protein-ligand interactions. This review elaborates on the principles of molecular docking and molecular dynamics simulation. Besides, their applications in food protein-ligand interactions are summarized. Furthermore, challenges, perspectives, and trends in molecular simulation of food protein-ligand interactions are proposed. Based on the results of molecular simulation, the mechanisms of interfacial behavior, enzyme-substrate binding, and structural changes during food processing can be reflected, and strategies for hazardous substance detection and food flavor adjustment can be generated. Moreover, molecular simulation can accelerate food development and reduce animal experiments. However, there are still several challenges to applying molecular simulation to food protein-ligand interaction research. The future trends will be a combination of international cooperation and data sharing, quantum mechanics/molecular mechanics, advanced computational techniques, and machine learning, which contribute to promoting food protein-ligand interaction simulation. Overall, the use of molecular simulation to study food protein-ligand interactions has a promising prospect.

Analysis of stabilization mechanisms in β-lactoglobulin-based amorphous solid dispersions by experimental and computational approaches

Our previous work shows that β-lactoglobulin-stabilized amorphous solid dispersion (ASD) loaded with 70 % indomethacin remains stable for more than 12 months. The stability is probably due to hydrogen bond networks spread throughout the ASD, facilitated by the indomethacin which has both hydrogen donors and acceptors. To investigate the stabilization mechanisms further, here we tested five other drug molecules, including two without any hydrogen bond donors. A combination of experimental techniques (differential scanning calorimetry, X-ray power diffraction) and molecular dynamics simulations was used to find the maximum drug loadings for ASDs with furosemide, griseofulvin, ibuprofen, ketoconazole and rifaximin. This approach revealed the underlying stabilization factors and the capacity of computer simulations to predict ASD stability. We searched the ASD models for crystalline patterns, and analyzed diffusivity of the drug molecules and hydrogen bond formation. ASDs loaded with rifaximin and ketoconazole remained stable for at least 12 months, even at 90 % drug loading, whereas stable drug loadings for furosemide, griseofulvin and ibuprofen were at a maximum of 70, 50 and 40 %, respectively. Steric confinement and hydrogen bonding to the proteins were the most important stabilization mechanisms at low drug loadings (≤ 40 %). Inter-drug hydrogen bond networks (including those with induced donors), ionic interactions, and a high Tg of the drug molecule were additional factors stabilizing the ASDs at drug loading greater than 40 %.

Molecular Dynamics Simulations of Selected Amorphous Stilbenoids and Their Amorphous Solid Dispersions with Poly(Vinylpyrrolidone)

Amorphous solid dispersions (ASDs) are one of the promising strategies to improve the solubility and dissolution rate of poorly soluble compounds. In this study, Molecular Dynamics simulations were used to investigate the interactions between three selected stilbenoids with important biological activity (resveratrol, pinostilbene and pterostilbene) and poly(vinylpyrrolidone). The analysis of the pair distribution functions and hydrogen bond distributions reveals a significant weakening of the hydrogen bond network of the stilbenoids in ASDs compared to the pure (no polymer) amorphous systems. This is accompanied by an increase in the mobility of the stilbenoid molecules in the ASDs, both in the translational dynamics determined from the molecular mean square displacements, and in the molecular reorientations followed by analysing several torsional distributions.

Lipid-based solubilization technology via hot melt extrusion: promises and challenges

INTRODUCTION

Self-emulsifying drug delivery systems (SEDDS) are a promising strategy to improve the oral bioavailability of poorly water-soluble drugs (PWSD). The excipients of SEDDS enable permeation through the mucus and gastro-intestinal barrier, inhibiting efflux transporters (e.g. P-glycoprotein) of drugs. Poor drug loading capacity and formulation instability are the main setbacks of traditional SEDDS. The use of polymeric precipitation inhibitors was shown to create supersaturable SEDDS with increased drug payload, and their solidification can help to overcome the instability challenge. As an alternative to several existing SEDDS solidification technologies, hot melt extrusion (HME) holds the potential for lean and continuous manufacturing of supersaturable solid-SEDDS. Despite being ubiquitously applied in solid lipid and polymeric processing, HME has not yet been widely considered for the preparation of SEDDS.

AREAS COVERED

The review begins with the rationale why SEDDS as the preferred lipid-based delivery systems (LBDS) is suitable for the oral delivery of PWSD and discusses the common barriers to oral administration. The potential of LBDS to surmount them is discussed. SEDDS as the flagship of LBDS for PWSD is proposed with a special emphasis on solid-SEDDS. Finally, the opportunities and challenges of HME from the lipid-based excipient (LBE) processing and product performance standpoint are highlighted.

EXPERT OPINION

HME can be a continuous, solvent-free, cost-effective, and scalable technology for manufacturing solid supersaturable SEDDS. Several critical formulations and process parameters in successfully preparing SEDDS via HME are identified.

Study on the Effect of Polymer Excipients on the Dispersibility, Interaction, Solubility, and Scavenging Reactive Oxygen Species of Myricetin Solid Dispersion: Experiment and Molecular Simulation

Although the preparation of amorphous solid dispersions can improve the solubility of crystalline drugs, there is still a lack of guidance on the micromechanism in the screening and evaluation of polymer excipients. In this study, a particular method of experimental characterization combined with molecular simulation was attempted on solubilization of myricetin (MYR) by solid dispersion. According to the analysis of the dispersibility and hydrogen-bond interaction, the effectiveness of the solid dispersion and the predicted sequence of poly(vinyl pyrrolidone) (PVP) > hypromellose (HPMC) > poly(ethylene glycol) (PEG) as the polymer excipient were verified. Through the dissolution, cell viability, and reactive oxygen species (ROS)-level detection, the reliability of simulation and micromechanism analysis was further confirmed. This work not only provided the theoretical guidance and screening basis for the miscibility of solid dispersions from the microscopic level but also served as a reference for the modification of new drugs.

Pharmaceutical profiling and molecular dynamics simulations reveal crystallization effects in amorphous formulations

Robust and reliable in vivo performance of medicines based on amorphous solid dispersions (ASDs) depend on maintenance of physical stability and efficient supersaturation. However, molecular drivers of these two kinetic processes are poorly understood. Here we used molecular dynamics (MD) simulations coupled with experimental assessments to explore supersaturation, nucleation, and crystal growth. The effect of drug loading on physical stability and supersaturation potential was highly drug specific. Storage under humid conditions influenced crystallization, but also resulted in morphological changes and particle fusion. This led to increased particle size, which significantly reduced dissolution rate. MD simulations identified the importance of nano-compartmentalization in the crystallization rate of the ASDs. Nucleation during storage did not inherently compromise the ASD. Rather, the poorer performance resulted from a combination of properties of the compound, nanostructures formed in the formulation, and crystallization.

Polyelectrolyte Matrices in the Modulation of Intermolecular Electrostatic Interactions for Amorphous Solid Dispersions: A Comprehensive Review

Polyelectrolyte polymers have been widely used in the pharmaceutical field as excipients to facilitate various drug delivery systems. Polyelectrolytes have been used to modulate the electrostatic environment and enhance favorable interactions between the drug and the polymer in amorphous solid dispersions (ASDs) prepared mainly by hot-melt extrusion. Polyelectrolytes have been used alone, or in combination with nonionic polymers as interpolyelectrolyte complexes, or after the addition of small molecular additives. They were found to enhance physical stability by favoring stabilizing intermolecular interactions, as well as to exert an antiplasticizing effect. Moreover, they not only enhance drug dissolution, but they have also been used for maintaining supersaturation, especially in the case of weakly basic drugs that tend to precipitate in the intestine. Additional uses include controlled and/or targeted drug release with enhanced physical stability and ease of preparation via novel continuous processes. Polyelectrolyte matrices, used along with scalable manufacturing methods in accordance with green chemistry principles, emerge as an attractive viable alternative for the preparation of ASDs with improved physical stability and biopharmaceutic performance.

Analytical and Computational Methods for the Determination of Drug-Polymer Solubility and Miscibility

In the pharmaceutical industry, poorly water-soluble drugs require enabling technologies to increase apparent solubility in the biological environment. Amorphous solid dispersion (ASD) has emerged as an attractive strategy that has been used to market more than 20 oral pharmaceutical products. The amorphous form is inherently unstable and exhibits phase separation and crystallization during shelf life storage. Polymers stabilize the amorphous drug by antiplasticization, reducing molecular mobility, reducing chemical potential of drug, and increasing glass transition temperature in ASD. Here, drug-polymer miscibility is an important contributor to the physical stability of ASDs. The current Review discusses the basics of drug-polymer interactions with the major focus on the methods for the evaluation of solubility and miscibility of the drug in the polymer. Methods for the evaluation of drug-polymer solubility and miscibility have been classified as thermal, spectroscopic, microscopic, solid-liquid equilibrium-based, rheological, and computational methods. Thermal methods have been commonly used to determine the solubility of the drug in the polymer, while other methods provide qualitative information about drug-polymer miscibility. Despite advancements, the majority of these methods are still inadequate to provide the value of drug-polymer miscibility at room temperature. There is still a need for methods that can accurately determine drug-polymer miscibility at pharmaceutically relevant temperatures.

Molecular Simulation and Statistical Learning Methods toward Predicting Drug-Polymer Amorphous Solid Dispersion Miscibility, Stability, and Formulation Design

Amorphous solid dispersions (ASDs) have emerged as widespread formulations for drug delivery of poorly soluble active pharmaceutical ingredients (APIs). Predicting the API solubility with various carriers in the API-carrier mixture and the principal API-carrier non-bonding interactions are critical factors for rational drug development and formulation decisions. Experimental determination of these interactions, solubility, and dissolution mechanisms is time-consuming, costly, and reliant on trial and error. To that end, molecular modeling has been applied to simulate ASD properties and mechanisms. Quantum mechanical methods elucidate the strength of API-carrier non-bonding interactions, while molecular dynamics simulations model and predict ASD physical stability, solubility, and dissolution mechanisms. Statistical learning models have been recently applied to the prediction of a variety of drug formulation properties and show immense potential for continued application in the understanding and prediction of ASD solubility. Continued theoretical progress and computational applications will accelerate lead compound development before clinical trials. This article reviews in silico research for the rational formulation design of low-solubility drugs. Pertinent theoretical groundwork is presented, modeling applications and limitations are discussed, and the prospective clinical benefits of accelerated ASD formulation are envisioned.

Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development

Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Formulation of Curcumin-β-cyclodextrin-polyvinylpyrrolidone supramolecular inclusion complex: experimental, molecular docking, and preclinical anti-inflammatory assessment

This research planned to ameliorate an aqueous solubility and dissolution of Curcumin (CUR) by the formulation of inclusion complex with β-cyclodextrin (β-CD) and polyvinylpyrrolidone (PVP). The phase solubility study was performed to assess the solubility of CUR. The prepared CUR complex assessed for dissolution study, physicochemical evaluation, in-vitro antioxidant activity, molecular modeling, and anti-inflammatory assessment. The pivotal findings of phase-solubility studies demonstrate apparent stability constant (Kc) and complexation efficiency (CE) values for CUR-β-CD and CUR-β-CD-PVP complex was 175.4 M ^-1, 1.15% and 833.3.2 M ^-1 and 5.21%, respectively. The characterization results revealed amorphization of crystalline state (CUR) into amorphous state. The maximum drug release found with the ternary CUR complex (F7), i.e. 45.41 ± 3.78% in 6 h study. The chemical shift in the NMR supports that the aromatic ring of CUR is completely complexed inside the β-CD cavity. The antioxidant activity of pure CUR was found to be 58.02 ± 2.21% and CUR ternary complex (F7) showed significantly higher activity to 96.02 ± 2.46%. The in-vivo effect of CUR complex (F7) was also found significantly higher than that of pure CUR. The molecular modeling study depicted that PVP increased the stability of the ternary complex by forming the link between CUR and β-CD. Thus, the ternary inclusion complex of CUR-β-CD-PVP could contribute as an innovative outcome in the enhancement of solubility and in-vivo activity.

Formulation of amorphous ternary solid dispersions of dapagliflozin using PEG 6000 and Poloxamer 188: solid-state characterization, ex vivo study, and molecular simulation assessment

The present study was designed to prepare dapagliflozin (DFG) loaded ternary solid dispersions (SDs) using the carrier blend polyethylene glycol 6000 (PEG 6000) and poloxamer 188 (PLX 188). The prepared DFG-SDs were evaluated for solubility study, physicochemical characterization and molecular simulation study. The prepared DFG-SDs showed significant higher solubility and dissolution vis-a-vis pure DFG and DFG physical mixture. The composition DFG:PEG:PLX (1:2.25:0.75 mM) showed the highest solubility (0.476 ± 0.016 mg/mL). The physicochemical characterization confirms the polymorphic transition of DFG from crystalline state to stable amorphous form. The prepared DFG-SDs showed a significantly higher dissolution (64.78 ± 2.34% to 78.41 ± 2.39%) than pure DFG (15.70 ± 3.54%). DFG-SD2 showed a significantly enhanced drug permeation (p<.05) (58.76 ± 4.65 µg/cm) as compared to pure DFG (14.97 ± 3.32 µg/cm). The molecular docking study result revealed a good hydrophobic interaction of DFG with the used carrier due to the lowest energy pose. The interaction occurs between the methylene bridges and the central hydrophobic chain of polyoxypropylene of the polymer. Therefore, DFG-SDs prepared by microwave irradiation method using hydrophilic carrier blend might be a promising strategy for improving the solubility and in vitro dissolution performance.

Developing HME-Based Drug Products Using Emerging Science: a Fast-Track Roadmap from Concept to Clinical Batch

This paper presents a rational workflow for developing enabling formulations, such as amorphous solid dispersions, via hot-melt extrusion in less than a year. First, our approach to an integrated product and process development framework is described, including state-of-the-art theoretical concepts, modeling, and experimental characterization described in the literature and developed by us. Next, lab-scale extruder setups are designed (processing conditions and screw design) based on a rational, model-based framework that takes into account the thermal load required, the mixing capabilities, and the thermo-mechanical degradation. The predicted optimal process setup can be validated quickly in the pilot plant. Lastly, a transfer of the process to any GMP-certified manufacturing site can be performed in silico for any extruder based on our validated computational framework. In summary, the proposed workflow massively reduces the risk in product and process development and shortens the drug-to-market time for enabling formulations.

Experimental, Thermodynamic, and Molecular Modeling Evaluation of Amorphous Simvastatin-Poly(vinylpyrrolidone) Solid Dispersions

A crucial step for the selection of proper amorphous solid dispersion (ASD) matrix carriers is the in-depth assessment of drug/polymer physicochemical properties. In this context, the present study extends the work of previously published attempts by evaluating the formation of simvastatin (SIM)-poly(vinylpyrrolidone) (PVP) ASDs with the aid of thermodynamic and molecular modeling. Specifically, the implementation of both Flory-Huggins lattice theory and molecular dynamics (MD) simulations was able to predict the miscibility between the two components (a finding that was experimentally verified via differential scanning calorimetry (DSC) and hot stage polarized microscopy), while a complete temperature-concentration phase-transition profile was constructed, leading to the identification of the thermodynamically metastable and unstable ASD zones. Furthermore, as in the case of previously published reports, the analysis of the ASDs via Fourier transform infrared spectroscopy did not clarify the type and extent of observed molecular interactions. Hence, in the present study, a computer-based MD simulation model was developed for the first time in order to gain an insight into the properties of the observed interactions. MD amorphous assemblies of SIM, PVP, and their mixtures were initially developed, and the calculated glass transition temperatures were in close agreement with experimentally obtained results, indicating that the developed models could be considered as realistic representations of the actual systems. Furthermore, molecular interactions evaluation via radial distribution function and radius of gyration analysis revealed that increasing SIM content results in a significant PVP chain shrinkage, which eventually leads to SIM-SIM amorphous intermolecular interactions, leading to the formation of amorphous drug zones. Finally, MD-based results were experimentally verified via DSC.

pH-sensitive loading/releasing of doxorubicin using single-walled carbon nanotube and multi-walled carbon nanotube: A molecular dynamics study

BACKGROUND AND OBJECTIVE

Doxorubicin is one of the drugs used to treat cancer, and many studies have been conducted to control its release. In this study, carbon nanotubes have been proposed as a doxorubicin carrier, and the effect of carboxyl functional group on the controlled release of doxorubicin has been studied.

METHODS

This study has been done by molecular dynamics simulation and was based on changing the pH as a mechanism controller.

RESULTS

This work is intended to test the efficacy of this drug carrier for the release of doxorubicin. A comparison was also made between single-walled and double-walled carbon nanotubes to answer the question of which one can be a better carrier for doxorubicin. The study of DOXORUBICIN adsorption and release showed that the DOXORUBICIN adsorption on single-walled carbon nanotube and multi-walled carbon nanotube in neutral pH was stronger than it was in acidic pH, which could be due to the electrostatic interactions between the carboxyl group of nanotubes and DOXORUBICIN. Based on this and according to the investigation of hydrogen bonds, diffusion coefficients, and other results it was clear that the drug release in acidic pH was appropriate for body conditions. Since cancer tissues pH is acidic, this shows the suitability of carbon nanotube in drug delivery and DOXORUBICIN release in cancer tissues. In addition, it was shown that the blood pH (pH = 7) is suitable for DOXORUBICIN loading on the carbon nanotube and carbon nanotube-DOXORUBICIN linkage remained stable at this pH; accordingly, the carbon nanotube could deliver DOXORUBICIN in blood quite well and release it in cancerous tissues. This suggests the carbon nanotubes as a promising drug carrier in the cancer therapy which can be also investigated in experiments.

CONCLUSION

It was revealed that the bonds between multi-walled carbon nanotube and DOXORUBICIN was stronger and this complex had a slower release in the cancer tissues compared to the single-walled carbon nanotube; this can be regarded as an advantage over the single-walled carbon nanotube in the DOXORUBICIN delivery and release.

Predicting Density of Amorphous Solid Materials Using Molecular Dynamics Simulation

The true density of an amorphous solid is an important parameter for studying and modeling materials behavior. Experimental measurements of density using helium pycnometry are standard but may be prevented if the material is prone to rapid recrystallization, or preparation of gram quantities of reproducible pure component amorphous materials proves impossible. The density of an amorphous solid can be approximated by assuming it to be 95% of its respective crystallographic density; however, this can be inaccurate or impossible if the crystal structure is unknown. Molecular dynamic simulations were used to predict the density of 20 amorphous solid materials. The calculated density values for 10 amorphous solids were compared with densities that were experimentally determined using helium pycnometry. In these cases, the amorphous densities calculated using molecular dynamics had an average percent error of - 0.7% relative to the measured values, with a maximum error of - 3.48%. In contrast, comparisons of amorphous density approximated from crystallographic structures with pycnometrically measured values resulted in an average percent error of + 3.7%, with a maximum error of + 9.42%. These data suggest that the density of an amorphous solid can be accurately predicted using molecular dynamic simulations and allowed reliable calculation of density for the remaining 10 materials for which pycnometry could not be done.

A Tumor-on-a-Chip System with Bioprinted Blood and Lymphatic Vessel Pair

Current in vitro anti-tumor drug screening strategies are insufficiently portrayed lacking true perfusion and draining microcirculation systems, which may post significant limitation in reproducing the transport kinetics of cancer therapeutics explicitly. Herein, we report the fabrication of an improved tumor model consisting of bioprinted hollow blood vessel and lymphatic vessel pair, hosted in a three-dimensional (3D) tumor microenvironment-mimetic hydrogel matrix, termed as the tumor-on-a-chip with bioprinted blood and lymphatic vessel pair (TOC-BBL). The bioprinted blood vessel was perfusable channel with opening on both ends while the bioprinted lymphatic vessel was blinded on one end, both of which were embedded in a hydrogel tumor mass, with vessel permeability individually tunable through optimization of the composition of the bioinks. We demonstrated that systems with different combinations of these bioprinted blood/lymphatic vessels exhibited varying levels of diffusion profiles for biomolecules and anti-cancer drugs. Our TOC-BBL platform mimicking the natural pathway of drug-tumor interactions would have the drug introduced through the perfusable blood vessel, cross the vascular wall into the tumor tissue via diffusion, and eventually drained into the lymphatic vessel along with the carrier flow. Our results suggested that this unique in vitro tumor model containing the bioprinted blood/lymphatic vessel pair may have the capacity of simulating the complex transport mechanisms of certain pharmaceutical compounds inside the tumor microenvironment, potentially providing improved accuracy in future cancer drug screening.