An In Silico Transwell Device for the Study of Drug Transport and Drug–Drug Interactions

Recent Advances in Drug Design and Delivery Across Biological Barriers using Computational Models

Abstract:

The systemic delivery of pharmacological substances generally exhibits several significant limitations associated with the bio-distribution of active drugs in the body. As per consequence, human body’s defense mechanisms become impediments to drug delivery. Various technologies to overcome these limitations have been evolved including computational approaches and advanced drug delivery. As the body of human has evolved to defend itself from hostile biological as well as chemical invaders, along with that these biological barriers such as ocular barriers, blood-brain barriers, intestinal and skin barriers also limit the passage of drugs across desired sites. Therefore, efficient delivery remains an utmost challenge for researchers and scientists. The present review focuses on the techniques to deliver the drugs with efficient therapeutic efficacy at the targeted sites. This review article considered the insights into main biological barriers along with the application of computational or numerical methods dealing with different barriers by determining the drug flow, temperature and various other parameters. It also summarizes the advanced implantable drug delivery system to circumvent the inherent resistance showed by these biological barriers and in turn to improve the drug delivery.

In silico drug absorption tract: An agent-based biomimetic model for human oral drug absorption

Background

An agent-based modeling approach has been suggested as an alternative to traditional, equation-based modeling methods for describing oral drug absorption. It enables researchers to gain a better understanding of the pharmacokinetic (PK) mechanisms of a drug. This project demonstrates that a biomimetic agent-based model can adequately describe the absorption and disposition kinetics both of midazolam and clonazepam.

Methods

An agent-based biomimetic model, in silico drug absorption tract (ISDAT), was built to mimic oral drug absorption in humans. The model consisted of distinct spaces, membranes, and metabolic enzymes, and it was altogether representative of human physiology relating to oral drug absorption. Simulated experiments were run with the model, and the results were compared to the referent data from clinical equivalence trials. Acceptable similarity was verified by pre-specified criteria, which included 1) qualitative visual matching between the clinical and simulated concentration-time profiles, 2) quantitative similarity indices, namely, weighted root mean squared error (RMSE), and weighted mean absolute percentage error (MAPE) and 3) descriptive similarity which requires less than 25% difference between key PK parameters calculated by the clinical and the simulated concentration-time profiles. The model and its parameters were iteratively refined until all similarity criteria were met. Furthermore, simulated PK experiments were conducted to predict bioavailability (F). For better visualization, a graphical user interface for the model was developed and a video is available in Supporting Information.

Results

Simulation results satisfied all three levels of similarity criteria for both drugs. The weighted RMSE was 0.51 and 0.92, and the weighted MAPE was 5.99% and 8.43% for midazolam and clonazepam, respectively. Calculated PK parameter values, including area under the curve (AUC), peak plasma drug concentration (C_max), time to reach C_max (T_max), terminal elimination rate constant (Kel), terminal elimination half life (T_1/2), apparent oral clearance (CL/F), and apparent volume of distribution (V/F), were reasonable compared to the referent values. The predicted absolute oral bioavailability (F) was 44% for midazolam (literature reported value, 31–72%) and 93% (literature reported value, ≥ 90%) for clonazepam.

Conclusion

The ISDAT met all the pre-specified similarity criteria for both midazolam and clonazepam, and demonstrated its ability to describe absorption kinetics of both drugs. Therefore, the validated ISDAT can be a promising platform for further research into the use of similar in silico models for drug absorption kinetics.

An Agent-Based Approach to Dynamically Represent the Pharmacokinetic Properties of Baicalein

Baicalein, a typical flavonoid presented in Scutellariae radix, exhibits a unique metabolic profile during first-pass metabolism: parallel glucuronidation and sulfation pathways, with possible substrate inhibition in both pathways. In this project, we aimed to construct an agent-based model to dynamically represent baicalein pharmacokinetics and to verify the substrate inhibition hypothesis. The model consisted of three 3D spaces and two membranes: apical space (S1), intracellular space (S2), basolateral space (S3), apical membrane (M1), and basolateral membrane (M2). In silico enzymes (UDP-glucuronosyltransferases (UGTs) and sulfotransferases (SULTs)) and binder components were placed in S2. The model was then executed to simulate one-pass metabolism experiments of baicalein. With the implementation of a two-site enzyme design, the simulated results captured the preset qualitative and quantitative features of the wet-lab observations. The feasible parameter set showed that substrate inhibition happened in both conjugation pathways of baicalein. The simulation results suggested that the sulfation pathway was dominant at low concentrations and that SULT was more inclined to substrate inhibition than UGT. Cross-model validation was satisfactory. Our findings were consistent with a previously reported catenary model. We conclude that the mechanisms represented by our model are plausible. Our novel modeling approach could dynamically represent the metabolic pathways of baicalein in a Caco-2 system.

Modeling bioavailability to organs protected by biological barriers

Computational pharmacokinetic (PK) modeling gives access to drug concentration vs. time profiles in target organs and allows better interpretation of clinical observations of therapeutic or toxic effects. Physiologically-based PK (PBPK) models in particular, based on mechanistic descriptions of the body anatomy and physiology, may also help to extrapolate in vitro or animal data to human.

Once in the systemic circulation, a chemical has access to the microvasculature of every organ or tissue. However, its penetration in the brain, retina, thymus, spinal cord, testis, placenta,… may be limited or even fully prevented by dynamic physiological blood-tissue barriers. Those barriers are both physical (involving tight junctions between adjacent cells) and biochemical (involving metabolizing enzymes and transporters).

On those cases, correct mechanistic characterization of the passage (or not) of molecules through the barrier can be crucial for improved PBPK modeling and prediction.

In parallel, attempts to understand and quantitatively characterize the processes involved in drug penetration of physiological barriers have led to the development of several in vitro experimental models. Data from such assays are very useful to calibrate PBPK models.

We review here those in vitro and computational models, highlighting the challenges and perspectives for in vitro and computational models to better assess drug availability to target tissues.

A critical analysis of the interplay between cytochrome P450 3A and P-glycoprotein: recent insights from knockout and transgenic mice

CYP3A is one of the most important drug-metabolizing enzymes, determining the first-pass metabolism, oral bioavailability, and elimination of many drugs. It is also an important determinant of variable drug exposure and is involved in many drug-drug interactions. Recent studies with CYP3A knockout and transgenic mice have yielded a number of key insights that are important to consider during drug discovery and development. For instance, studies with tissue-specific CYP3A-transgenic mice have highlighted the importance of intestinal CYP3A-dependent metabolism. They also revealed that intestinal CYP3A plays an important role in the regulation of various drug-handling systems in the liver. Intestinal CYP3A activity can thus have far-reaching pharmacological effects. Besides CYP3A, the active drug efflux transporter P-glycoprotein also has a strong effect on the pharmacokinetics of numerous drugs. CYP3A and P-glycoprotein have an extensive overlap in their substrate spectrum. It has been hypothesized that for many drugs, the combined activity of CYP3A and P-glycoprotein makes for efficient intestinal first-pass metabolism of orally administered drugs as a result of a potentially synergistic collaboration. However, there is only limited in vitro and in vivo evidence for this hypothesis. There has also been some confusion in the field about what synergy actually means in this case. Our recent studies with Cyp3a/P-glycoprotein combination knockout mice have provided further insights into the CYP3A-P-glycoprotein interplay. We here present our view of the status of the synergy hypothesis and an attempt to clarify the existing confusion about synergy. We hope that this will facilitate further critical testing of the hypothesis and improve communication among researchers. Above all, the recent findings and insights into the interplay between CYP3A and P-glycoprotein may have implications for improving oral drug bioavailability and reducing adverse side effects.

Internet resources for agent-based modelling

The use of agent-based models (ABMs) is steadily increasing in all the disciplines including environmental chemistry and toxicology. This growth is mainly driven by their ability to address problems that conventional modelling techniques cannot, such as the change of scale or the emergence of unanticipated phenomena resulting from interactions between their constitutive goal-directed agents. After a brief introduction on the basic principles of agent-based modelling and the presentation of selected case studies, the main software resources available on the Internet are presented. An attempt is made to estimate the complexity of these tools versus their potentialities and flexibility.

Mechanistic insight from in silico pharmacokinetic experiments: roles of P-glycoprotein, Cyp3A4 enzymes, and microenvironments

Saquinavir exhibits paradoxical transport across modified Caco-2 cell monolayers (doi: 10.1124/jpet.103.056390) expressing P-glycoprotein and Cyp3A4. The data implicate complicated intracellular transport mechanisms. Drawing on recent discrete event modeling and simulation advances, we built an in silico analog of the confluent, asymmetric cell monolayer used in the cited work. We call it in silico experimental Caco-2 (cell monolayer) culture (ISECC). Concrete, working, hypothesized spatial mechanisms were implemented. Validation was achieved when in silico experimental results met similarity measure (SM) expectations that targeted 16 wet-lab experimental conditions. Initial mechanistic hypotheses turned out to be necessary parts of a more complicated explanation. We progressed through four stages of an iterative refinement and validation protocol that enabled and facilitated discovery of plausible, new mechanistic details. The process exercised abductive reasoning, a primary means of scientific knowledge creation and creative cognition. The ISECC that survived the most stringent SM challenge produced transport data that was statistically indistinguishable from referent wet-lab observations. It required a 7:1 ratio of apical transporters to metabolizing enzymes, a 97% reduction of efflux activity by an inhibitor, a biased distribution of metabolizing enzymes, heterogeneous intracellular spaces, and restrictions on intracellular drug movement. Experimenting on synthetic analogs such as ISECC provides a former unavailable means of discovering new mechanistic details and testing their plausibility. The approach thus provides a powerful new expansion of the scientific method: an independent, scientific means to challenge, explore, better understand, and improve any inductive mechanism and, importantly, the assumptions on which it rests.

At the biological modeling and simulation frontier

We provide a rationale for and describe examples of synthetic modeling and simulation (M&S) of biological systems. We explain how synthetic methods are distinct from familiar inductive methods. Synthetic M&S is a means to better understand the mechanisms that generate normal and disease-related phenomena observed in research, and how compounds of interest interact with them to alter phenomena. An objective is to build better, working hypotheses of plausible mechanisms. A synthetic model is an extant hypothesis: execution produces an observable mechanism and phenomena. Mobile objects representing compounds carry information enabling components to distinguish between them and react accordingly when different compounds are studied simultaneously. We argue that the familiar inductive approaches contribute to the general inefficiencies being experienced by pharmaceutical R&D, and that use of synthetic approaches accelerates and improves R&D decision-making and thus the drug development process. A reason is that synthetic models encourage and facilitate abductive scientific reasoning, a primary means of knowledge creation and creative cognition. When synthetic models are executed, we observe different aspects of knowledge in action from different perspectives. These models can be tuned to reflect differences in experimental conditions and individuals, making translational research more concrete while moving us closer to personalized medicine.

Computational strategies unravel and trace how liver disease changes hepatic drug disposition

Liver disease changes the disposition properties of drugs, complicating drug therapy management. We present normal and "diseased" versions of an abstract, agent-oriented In Silico Livers (ISLs), and validate their mechanisms against disposition data from perfused normal and diseased rat livers. Dynamic tracing features enabled spatiotemporal tracing of differences in dispositional events for diltiazem and sucrose across five levels, including interactions with representations of lobular microarchitectural features, cells, and intracellular factors that sequester and metabolize. Differences in attributes map to measures of histopathology. We measured disease-causing differences in local, intralobular ISL effects, obtaining until now unavailable views of how and where hepatic drug disposition may differ in normal and diseased rat livers from diltiazem's perspective. Exploration of disposition in less and more advanced stages of disease is feasible. The approach and technology represent an important step toward unraveling the complex changes from normal to disease states and their influences on drug disposition.

Mechanistic simulations explain paradoxical saquinavir metabolism during in vitro vectorial transport study

We were confronted by an unexpected observation: in a Transwell study of confluent, Cyp3A4 and P-gp-expressing Caco-2 cells, higher intracellular saquinavir levels, yet less first metabolite (M7) formation were observed following apical dosing, compared to basal dosing. To test two seemingly plausible hypothesized explanations, we constructed an in silico working analogue using the synthetic method. Neither mechanism alone was sufficient, but when combined and tuned within the analogue, the results generated were a semi-quantitative match to the experimental data. After 60 cycles, more of the simulated dose was present within analogue cells as parent drug after apical dosing. Furthermore, less M7 was present after apical dose. The paradox disappeared by having simulated drugs equilibrate among separate intracellular zones. Building, studying, and exploring mechanistic explanations for complex wet-lab phenomena using the new methods improved insight into the referent system, while providing a straightforward, scientific means of testing the plausibility of mechanistic hypotheses.