A model for intracellular trafficking of adenoviral vectors

New Horizons of Model Informed Drug Development in Rare Diseases Drug Development

Model-informed approaches provide a quantitative framework to integrate all available nonclinical and clinical data, thus furnishing a totality of evidence approach to drug development and regulatory evaluation. Maximizing the use of all available data and information about the drug enables a more robust characterization of the risk-benefit profile and reduces uncertainty in both technical and regulatory success. This offers the potential to transform rare diseases drug development, where conducting large well-controlled clinical trials is impractical and/or unethical due to a small patient population, a significant portion of which could be children. Additionally, the totality of evidence generated by model-informed approaches can provide confirmatory evidence for regulatory approval without the need for additional clinical data. In the article, applications of novel quantitative approaches such as quantitative systems pharmacology, disease progression modeling, artificial intelligence, machine learning, modeling of real-world data using model-based meta-analysis and strategies such as external control and patient-reported outcomes as well as clinical trial simulations to optimize trials and sample collection are discussed. Specific case studies of these modeling approaches in rare diseases are provided to showcase applications in drug development and regulatory review. Finally, perspectives are shared on the future state of these modeling approaches in rare diseases drug development along with challenges and opportunities for incorporating such tools in the rational development of drug products.

Mechanistic modeling of viral particle production

Viral systems such as wild-type viruses, viral vectors, and virus-like particles are essential components of modern biotechnology and medicine. Despite their importance, the commercial-scale production of viral systems remains highly inefficient for multiple reasons. Computational strategies are a promising avenue for improving process development, optimization, and control, but require a mathematical description of the system. This article reviews mechanistic modeling strategies for the production of viral particles, both at the cellular and bioreactor scales. In many cases, techniques and models from adjacent fields such as epidemiology and wild-type viral infection kinetics can be adapted to construct a suitable process model. These process models can then be employed for various purposes such as in-silico testing of novel process operating strategies and/or advanced process control.

Application of physiologically-based pharmacokinetic models for therapeutic proteins and other novel modalities

The past two decades have seen diversification of drug development pipelines and approvals from traditional small molecule therapies to alternative modalities including monoclonal antibodies, engineered proteins, antibody drug conjugates (ADCs), oligonucleotides and gene therapies. At the same time, physiologically-based pharmacokinetic (PBPK) models for small molecules have seen increased industry and regulatory acceptance.This review focusses on the current status of the application of PBPK models to these newer modalities and give a perspective on the successes, challenges and future directions of this field.There is greatest experience in the development of PBPK models for therapeutic proteins, and PBPK models for ADCs benefit from prior experience for both therapeutic proteins and small molecules. For other modalities, the application of PBPK models is in its infancy.Challenges are discussed and a common theme is lack of availability of physiological and experimental data to characterise systems and drug parameters to enable a priori prediction of pharmacokinetics. Furthermore, sufficient clinical data are required to build confidence in developed models.The PBPK modelling approach provides a quantitative framework for integrating knowledge and data from multiple sources and can be built on as more data becomes available.

Probing the Intracellular Delivery of Nanoparticles into Hard-to-Transfect Cells

Hard-to-transfect cells are cells that are known to present special difficulties in intracellular delivery of exogenous entities. However, the special transport behaviors underlying the special delivery problem in these cells have so far not been examined carefully. Here, we combine single-particle motion analysis, cell biology studies, and mathematical modeling to investigate nanoparticle transport in bone marrow-derived mesenchymal stem cells (BMSCs), a technologically important type of hard-to-transfect cells. Tat peptide-conjugated quantum dots (QDs-Tat) were used as the model nanoparticles. Two different yet complementary single-particle methods, namely, pair-correlation function and single-particle tracking, were conducted on the same cell samples and on the same viewing stage of a confocal microscope. Our results reveal significant differences in each individual step of transport of QDs-Tat in BMSCs vs a commonly used model cell line, HeLa cells. Single-particle motion analysis demonstrates that vesicle escape and cytoplasmic diffusion are dramatically more difficult in BMSCs than in HeLa cells. Cell biology studies show that BMSCs use different biological pathways for the cellular uptake, vesicular transport, and exocytosis of QDs-Tat than HeLa cells. A reaction-diffusion-advection model is employed to mathematically integrate the individual steps of cellular transport and can be used to predict and design nanoparticle delivery in BMSCs. This work provides dissective, quantitative, and mechanistic understandings of nanoparticle transport in BMSCs. The investigative methods described in this work can help to guide the tailored design of nanoparticle-based delivery in specific types and subtypes of hard-to-transfect cells.

Current progress and limitations of AAV mediated delivery of protein therapeutic genes and the importance of developing quantitative pharmacokinetic/pharmacodynamic (PK/PD) models

While protein therapeutics are one of the most successful class of drug molecules, they are expensive and not suited for treating chronic disorders that require long-term dosing. Adeno-associated virus (AAV) mediated in vivo gene therapy represents a viable alternative, which can deliver the genes of protein therapeutics to produce long-term expression of proteins in target tissues. Ongoing clinical trials and recent regulatory approvals demonstrate great interest in these therapeutics, however, there is a lack of understanding regarding their cellular disposition, whole-body disposition, dose-exposure relationship, exposure-response relationship, and how product quality and immunogenicity affects these important properties. In addition, there is a lack of quantitative studies to support the development of pharmacokinetic-pharmacodynamic models, which can support the discovery, development, and clinical translation of this delivery system. In this review, we have provided a state-of-the-art overview of current progress and limitations related to AAV mediated delivery of protein therapeutic genes, along with our perspective on the steps that need to be taken to improve clinical translation of this therapeutic modality.

Adenovirus Entry: From Infection to Immunity

More than 80 different adenovirus (AdV) types infect humans through the respiratory, ocular, or gastrointestinal tracts. They cause acute clinical mani-festations or persist under humoral and cell-based immunity. Immuno-suppressed individuals are at risk of death from an AdV infection. Concepts about cell entry of AdV build on strong foundations from molecular and cellular biology-and increasingly physical virology. Here, we discuss how virions enter and deliver their genome into the nucleus of epithelial cells. This process breaks open the virion at distinct sites because the particle has nonisometric mechanical strength and reacts to specific host factors along the entry pathway. We further describe how macrophages and dendritic cells resist AdV infection yet enhance productive entry into polarized epithelial cells. A deep understanding of the viral mechanisms and cell biological and biophysical principles will continue to unravel how epithelial and antigen-presenting cells respond to AdVs and control inflammation and persistence in pathology and therapy.

Kinetic Modeling of Virus Growth in Cells

When a virus infects a host cell, it hijacks the biosynthetic capacity of the cell to produce virus progeny, a process that may take less than an hour or more than a week. The overall time required for a virus to reproduce depends collectively on the rates of multiple steps in the infection process, including initial binding of the virus particle to the surface of the cell, virus internalization and release of the viral genome within the cell, decoding of the genome to make viral proteins, replication of the genome, assembly of progeny virus particles, and release of these particles into the extracellular environment. For a large number of virus types, much has been learned about the molecular mechanisms and rates of the various steps. However, in only relatively few cases during the last 50 years has an attempt been made-using mathematical modeling-to account for how the different steps contribute to the overall timing and productivity of the infection cycle in a cell. Here we review the initial case studies, which include studies of the one-step growth behavior of viruses that infect bacteria (Qβ, T7, and M13), human immunodeficiency virus, influenza A virus, poliovirus, vesicular stomatitis virus, baculovirus, hepatitis B and C viruses, and herpes simplex virus. Further, we consider how such models enable one to explore how cellular resources are utilized and how antiviral strategies might be designed to resist escape. Finally, we highlight challenges and opportunities at the frontiers of cell-level modeling of virus infections.

How Computational Models Enable Mechanistic Insights into Virus Infection

An implicit aim in cellular infection biology is to understand the mechanisms how viruses, microbes, eukaryotic parasites, and fungi usurp the functions of host cells and cause disease. Mechanistic insight is a deep understanding of the biophysical and biochemical processes that give rise to an observable phenomenon. It is typically subject to falsification, that is, it is accessible to experimentation and empirical data acquisition. This is different from logic and mathematics, which are not empirical, but built on systems of inherently consistent axioms. Here, we argue that modeling and computer simulation, combined with mechanistic insights, yields unprecedented deep understanding of phenomena in biology and especially in virus infections by providing a way of showing sufficiency of a hypothetical mechanism. This ideally complements the necessity statements accessible to empirical falsification by additional positive evidence. We discuss how computational implementations of mathematical models can assist and enhance the quantitative measurements of infection dynamics of enveloped and non-enveloped viruses and thereby help generating causal insights into virus infection biology.

A blueprint for robust crosslinking of mobile species in biogels with weakly adhesive molecular anchors

Biopolymeric matrices can impede transport of nanoparticulates and pathogens by entropic or direct adhesive interactions, or by harnessing "third-party" molecular anchors to crosslink nanoparticulates to matrix constituents. The trapping potency of anchors is dictated by association rates and affinities to both nanoparticulates and matrix; the popular dogma is that long-lived, high-affinity bonds to both species facilitate optimal trapping. Here we present a contrasting paradigm combining experimental evidence (using IgG antibodies and Matrigel®), a theoretical framework (based on multiple timescale analysis), and computational modeling. Anchors that bind and unbind rapidly from matrix accumulate on nanoparticulates much more quickly than anchors that form high-affinity, long-lived bonds with matrix, leading to markedly greater trapping potency of multiple invading species without saturating matrix trapping capacity. Our results provide a blueprint for engineering molecular anchors with finely tuned affinities to effectively enhance the barrier properties of biogels against diverse nanoparticulate species.Biological polymeric matrices often use molecular anchors, such as antibodies, to trap nanoparticulates. Here, the authors find that anchor-matrix bonds that are weak and short-lived confer superior trapping potency, contrary to the prevailing belief that effective molecular anchors should form strong bonds to both the matrix and the nanoparticulates.

Mathematical modeling and quantitative analysis of HIV-1 Gag trafficking and polymerization

Gag, as the major structural protein of HIV-1, is necessary for the assembly of the HIV-1 sphere shell. An in-depth understanding of its trafficking and polymerization is important for gaining further insights into the mechanisms of HIV-1 replication and the design of antiviral drugs. We developed a mathematical model to simulate two biophysical processes, specifically Gag monomer and dimer transport in the cytoplasm and the polymerization of monomers to form a hexamer underneath the plasma membrane. Using experimental data, an optimization approach was utilized to identify the model parameters, and the identifiability and sensitivity of these parameters were then analyzed. Using our model, we analyzed the weight of the pathways involved in the polymerization reactions and concluded that the predominant pathways for the formation of a hexamer might be the polymerization of two monomers to form a dimer, the polymerization of a dimer and a monomer to form a trimer, and the polymerization of two trimers to form a hexamer. We then deduced that the dimer and trimer intermediates might be crucial in hexamer formation. We also explored four theoretical combined methods for Gag suppression, and hypothesized that the N-terminal glycine residue of the MA domain of Gag might be a promising drug target. This work serves as a guide for future theoretical and experimental efforts aiming to understand HIV-1 Gag trafficking and polymerization, and might help accelerate the efficiency of anti-AIDS drug design.

The human immunodeficiency virus (HIV-1) is a retrovirus that causes acquired immunodeficiency syndrome (AIDS), an infectious disease with high annual mortality. Gag protein is the major structural protein of HIV-1 and can self-assemble into the HIV-1 sphere shell. Therefore, an in-depth understanding of Gag protein trafficking and polymerization is important for gaining further insights into the mechanisms of HIV-1 replication and the design of antiviral drugs. Through mathematical modeling, optimization and quantitative analysis, we hypothesized the budding and release time of virus-like particles and revealed that the dimer and trimer intermediates might be crucial in hexamer formation. We also concluded that the predominant pathways in hexamer formation might be the polymerization of two monomers to form a dimer, the polymerization of a dimer and a monomer to form a trimer, and the polymerization of two trimers to form a hexamer. Our analysis also suggested that the N-terminal glycine residue of the MA domain of Gag might be a promising drug target. These results serve as a guide for further theoretical and experimental efforts aiming to understand HIV-1 Gag trafficking and polymerization and could aid anti-AIDS drug design.

Mathematical model for transport of DNA plasmids from the external medium up to the nucleus by electroporation

We propose a mathematical model for the transport of DNA plasmids from the extracellular matrix up to the cell nucleus. The model couples two phenomena: the electroporation process, describing the cell membrane permeabilization to plasmids and the intracellular transport enhanced by the presence of microtubules. Numerical simulations of cells with arbitrary geometry, in 2D and 3D, and a network of microtubules show numerically the importance of the microtubules and the electroporation on the effectiveness of the DNA transfection, as observed by previous biological data. The paper proposes efficient numerical tools for forthcoming optimized procedures of cell transfection.

Multiscale modeling of virus replication and spread

Replication and spread of human viruses is based on the simultaneous exploitation of many different host functions, bridging multiple scales in space and time. Mathematical modeling is essential to obtain a systems-level understanding of how human viruses manage to proceed through their life cycles. Here, we review corresponding advances for viral systems of large medical relevance, such as human immunodeficiency virus-1 (HIV-1) and hepatitis C virus (HCV). We will outline how the combination of mathematical models and experimental data has advanced our quantitative knowledge about various processes of these pathogens, and how novel quantitative approaches promise to fill remaining gaps.

A stochastic spatio-temporal (SST) model to study cell-to-cell variability in HIV-1 infection

Although HIV viremia in infected patients proceeds in a manner that may be accounted for by deterministic mathematical models, single virus-cell encounters following initial HIV exposure result in a variety of outcomes, only one of which results in a productive infection. The development of single molecule tracking techniques in living cells allows studies of intracellular transport of HIV, but it remains less clear what its impact may be on viral integration efficiency. Here, we present a stochastic intracellular mathematical model of HIV replication that incorporates microtubule transport of viral components. Using this model, we could study single round infections and observe how viruses entering cells reach one of three potential fates - degradation of the viral RNA genome, formation of LTR circles, or successful integration and establishment of a provirus. Our model predicts global trafficking properties, such as the probability and the mean time for a HIV viral particle to reach the nuclear pore. Interestingly, our model predicts that trafficking determines neither the probability or time of provirus establishment - instead, they are a function of vRNA degradation and reverse transcription reactions. Thus, our spatio-temporal model provides novel insights into the HIV infection process and may constitute a useful tool for the identification of promising drug targets.

Coarse-graining intermittent intracellular transport: Two- and three-dimensional models

Viruses and other cellular cargo that lack locomotion must rely on diffusion and cellular transport systems to navigate through a biological cell. Indeed, advances in single particle tracking have revealed that viral motion alternates between (a) diffusion in the cytoplasm and (b) active transport along microtubules. This intermittency makes quantitative analysis of trajectories difficult. Therefore, the purpose of this paper is to construct mathematical methods to approximate intermittent dynamics by effective stochastic differential equations. The coarse-graining method that we develop is more accurate than existing techniques and applicable to a wide range of intermittent transport models. In particular, we apply our method to two- and three-dimensional cell geometries (disk, sphere, and cylinder) and demonstrate its accuracy. In addition to these specific applications, we also explain our method in full generality for use on future intermittent models.

Mathematical modeling for intracellular transport and binding of HIV-1 Gag proteins

This paper presents a modeling study for the intracellular trafficking and trimerization of the HIV-1 Gag proteins. A set of differential equations including initial and boundary conditions is used to characterize the transport, diffusion, association and dissociation of Gag monomers and trimers for the time period from the initial production of Gag protein monomers to the initial appearance of immature HIV-1 virions near the cell membrane (the time duration Ta). The existence and stability of the steady-state solution of the initial boundary value problems provide a quantitative characterization of the tendency and equilibrium of Gag protein movement. The numerical simulation results further demonstrate Gag trimerization near the cell membrane. Our calculations of Ta are in good agreement with published experimental data. Sensitivity analysis of Ta to the model parameters indicates that the timing of the initial appearance of HIV-1 virions on the cell membrane is affected by the diffusion and transport processes. These results provide important information and insight into the Gag protein transport and binding and HIV-1 virion formation.

Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity

We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells.

Stochastic model-assisted development of efficient low-dose viral transduction in microfluidics

Adenoviruses are commonly used in vitro as gene transfer vectors in multiple applications. Nevertheless, issues such as low infection efficiency and toxicity effects on host cells have not been resolved yet. This work aims at developing a new versatile tool to enhance the expression of transduced genes while working at low viral doses in a sequential manner. We developed a microfluidic platform with automatically controlled sequential perfusion stages, which includes 10 independent channels. In addition, we built a stochastic mathematical model, accounting for the discrete nature of cells and viruses, to predict not only the percentage of infected cells, but also the associated infecting-virus distribution in the cell population. Microfluidic system and mathematical model were coupled to define an efficient experimental strategy. We used human foreskin fibroblasts, infected by replication-incompetent adenoviruses carrying EGFP gene, as the testing system. Cell characterization was performed through fluorescence microscopy, followed by image analysis. We explored the effect of different aspects: perfusion, multiplicity of infection, and temporal patterns of infection. We demonstrated feasibility of performing efficient viral transduction at low doses, by repeated pulses of cell-virus contact. This procedure also enhanced the exogenous gene expression in the sequential microfluidic infection system compared to a single infection at a higher, nontoxic, viral dose.

Non-Markovian polymer reaction kinetics

Describing the kinetics of polymer reactions, such as the formation of loops and hairpins in nucleic acids or polypeptides, is complicated by the structural dynamics of their chains. Although both intramolecular reactions, such as cyclization, and intermolecular reactions have been studied extensively, both experimentally and theoretically, there is to date no exact explicit analytical treatment of transport-limited polymer reaction kinetics, even in the case of the simplest (Rouse) model of monomers connected by linear springs. We introduce a new analytical approach to calculate the mean reaction time of polymer reactions that encompasses the non-Markovian dynamics of monomer motion. This requires that the conformational statistics of the polymer at the very instant of reaction be determined, which provides, as a by-product, new information on the reaction path. We show that the typical reactive conformation of the polymer is more extended than the equilibrium conformation, which leads to reaction times significantly shorter than predicted by the existing classical Markovian theory.

Influence of the nuclear membrane, active transport, and cell shape on the Hes1 and p53-Mdm2 pathways: insights from spatio-temporal modelling

There are many intracellular signalling pathways where the spatial distribution of the molecular species cannot be neglected. These pathways often contain negative feedback loops and can exhibit oscillatory dynamics in space and time. Two such pathways are those involving Hes1 and p53-Mdm2, both of which are implicated in cancer. In this paper we further develop the partial differential equation (PDE) models of Sturrock et al. (J. Theor. Biol., 273:15-31, 2011) which were used to study these dynamics. We extend these PDE models by including a nuclear membrane and active transport, assuming that proteins are convected in the cytoplasm towards the nucleus in order to model transport along microtubules. We also account for Mdm2 inhibition of p53 transcriptional activity. Through numerical simulations we find ranges of values for the model parameters such that sustained oscillatory dynamics occur, consistent with available experimental measurements. We also find that our model extensions act to broaden the parameter ranges that yield oscillations. Hence oscillatory behaviour is made more robust by the inclusion of both the nuclear membrane and active transport. In order to bridge the gap between in vivo and in silico experiments, we investigate more realistic cell geometries by using an imported image of a real cell as our computational domain. For the extended p53-Mdm2 model, we consider the effect of microtubule-disrupting drugs and proteasome inhibitor drugs, obtaining results that are in agreement with experimental studies.

Trafficking coordinate description of intracellular transport control of signaling networks

Many cellular networks rely on the regulated transport of their components to transduce extracellular information into precise intracellular signals. The dynamics of these networks is typically described in terms of compartmentalized chemical reactions. There are many important situations, however, in which the properties of the compartments change continuously in a way that cannot naturally be described by chemical reactions. Here, we develop an approach based on transport along a trafficking coordinate to precisely describe these processes and we apply it explicitly to the TGF-β signal transduction network, which plays a fundamental role in many diseases and cellular processes. The results of this newly introduced approach accurately capture the distinct TGF-β signaling dynamics of cells with and without cancerous backgrounds and provide an avenue to predict the effects of chemical perturbations in a way that closely recapitulates the observed cellular behavior.

Modeling of intracellular transport and compartmentation

The complexity and internal organization of mammalian cells as well as the regulation of intracellular transport processes has increasingly moved into the focus of investigation during the past two decades. Advanced staining and microscopy techniques help to shed light onto spatial cellular compartmentation and regulation, increasing the demand for improved modeling techniques. In this chapter, we summarize recent developments in the field of quantitative simulation approaches and frameworks for the description of intracellular transport processes. Special focus is therefore laid on compartmented and spatiotemporally resolved simulation approaches. The processes considered include free and facilitated diffusion of molecules, active transport via the microtubule and actin filament network, vesicle distribution, membrane transport, cell cycle-dependent cell growth and morphology variation, and protein production. Commercially and freely available simulation packages are summarized as well as model data exchange and harmonization issues.

Gene therapy: a pharmacokinetic/pharmacodynamic modelling overview

Since gene therapy started over 20 years ago, more than one-thousand clinical trials have been carried out. Nonviral vectors present interesting properties for their clinical application, but their efficiency in vivo is relatively low, and further improvements in these vectors are needed. Elucidating how nonviral vectors behave at the intracellular level is enlightening for vector improvement and optimization. Model-based approach is a powerful tool to understand and describe the different processes that gene transfer systems should overcome inside the body. Model-based approach allows for proposing and predicting the effect of parameter changes on the overall gene therapy response, as well as the known application of the pharmacokinetic/pharmacodynamic modelling in conventional therapies. The objective of this paper is to critically review the works in which the time-course of naked or formulated DNA have been quantitatively studied or modelled.

Optimal cytoplasmic transport in viral infections

For many viruses, the ability to infect eukaryotic cells depends on their transport through the cytoplasm and across the nuclear membrane of the host cell. During this journey, viral contents are biochemically processed into complexes capable of both nuclear penetration and genomic integration. We develop a stochastic model of viral entry that incorporates all relevant aspects of transport, including convection along microtubules, biochemical conversion, degradation, and nuclear entry. Analysis of the nuclear infection probabilities in terms of the transport velocity, degradation, and biochemical conversion rates shows how certain values of key parameters can maximize the nuclear entry probability of the viral material. The existence of such "optimal" infection scenarios depends on the details of the biochemical conversion process and implies potentially counterintuitive effects in viral infection, suggesting new avenues for antiviral treatment. Such optimal parameter values provide a plausible transport-based explanation of the action of restriction factors and of experimentally observed optimal capsid stability. Finally, we propose a new interpretation of how genetic mutations unrelated to the mechanism of drug action may nonetheless confer novel types of overall drug resistance.

Physical principles and models describing intracellular virus particle dynamics

Modeling in cellular biology benefits greatly from quantitative analysis that arise from the theory of diffusion and chemical reactions. Recent progress in single particle imaging enables the visualization of viral trajectories evolving in the cytoplasm. Biophysical models and mathematical analysis have been developed to unravel the complexity of single viral trajectories. We review here models of active motion of viruses along the cytoskeleton as well as their diffusion. We present resent efforts to estimate global trafficking properties, such as the probability and the mean time for a viral particle to reach a small nuclear pore. However, most signaling pathways involved in controlling viral motion remain undescribed and should be the goal of future modeling efforts.

Quantitative analysis of virus and plasmid trafficking in cells

Intracellular transport of DNA carriers is a fundamental step of gene delivery. By combining both theoretical and numerical approaches we study here single and several viruses and DNA particles trafficking in the cell cytoplasm to a small nuclear pore. We present a physical model to account for certain aspects of cellular organization, starting with the observation that a viral trajectory consists of epochs of pure diffusion and epochs of active transport along microtubules. We define a general degradation rate to describe the limitations of the delivery of plasmid or viral particles to a nuclear pore imposed by various types of direct and indirect hydrolysis activity inside the cytoplasm. By replacing the switching dynamics by a single steady state stochastic description, we obtain estimates for the probability and the mean time for the first one of many particles to go from the cell membrane to a small nuclear pore. Computational simulations confirm that our model can be used to analyze and interpret viral trajectories and estimate quantitatively the success of nuclear delivery.

Numerical modeling of molecular-motor-assisted transport of adenoviral vectors in a spherical cell

Viral gene delivery in a spherical cell is investigated numerically. The model of intracellular trafficking of adenoviruses is based on molecular-motor-assisted transport equations suggested by Smith and Simmons. These equations are presented in spherical coordinates and extended by accounting for the random component of motion of viral particles bound to filaments. This random component is associated with the stochastic nature of molecular motors responsible for locomotion of viral particles bound to filaments. The equations are solved numerically to simulate viral transport between the cell membrane and cell nucleus during initial stages of viral infection.

Mitochondrial imaging in dorsal root ganglion neurons following the application of inducible adenoviral vector expressing two fluorescent proteins

Mitochondrial morphology and dynamics are known to vary considerably depending on the cell type and organism studied. The objective of this study was to assess the potential application of adenoviral-fluorescent protein constructs for long-term tracking of mitochondria in neurons. An adenoviral vector containing two fluorescent proteins, the enhanced green fluorescent protein (eGFP) targeted to the cytoplasm to highlight the neuronal processes, and the red fluorescent protein (RFP) directed to mitochondria under the control of an inducible promoter, facilitated an efficient and accurate method to study mitochondrial dynamics in long-term studies. Dorsal root ganglion neurons from rat embryos were cultured and infected. The infected neurons exhibited green fluorescence after 24h, while 16 h following induction with doxycycline, red fluorescence protein began to localize within mitochondria. The red fluorescent protein was transported into mitochondria at the cell body followed by distribution within processes. As the neurons aged, the expression of red fluorescent protein was confined to cytoplasmic vacuoles and not mitochondria. Further analysis suggested that the cytoplasmic vacuoles were likely of lysosomal origin. Taken together, the current study presents novel strategies to study the life history of cellular organelles such as mitochondria in long-term studies.

The method of separation of variables for solving equations describing molecular-motor-assisted transport of intracellular particles in a dendrite or axon

This paper presents an analytical solution of one-dimensional transient molecular-motor-assisted transport equations that describe transport of different organelles (such as transport vesicles loaded with a cargo of specific proteins) in a neuron's axon or dendrite. Large intracellular organelles are transported in the cytoplasm by a combined action of diffusion and motor-driven transport. In an axon, organelles are transported away from the neuron's body towards the axon's terminal by kinesin-family molecular motors running on tracks composed of microtubules (MTs); old and used components are carried back towards the neuron's body by dynein-family molecular motors. Using the method of separation of variables, a generalized Fourier series solution for this problem is obtained. The solution uses three different orthogonal sets of eigenfunctions to represent the concentration of free organelles transported by diffusion, MT-bound organelles transported away from the neuron's body, and MT-bound organelles transported towards the neuron's body. Binding/detachment kinetic processes between the organelles and the MT are specified by first-order rate constants; these lead to coupling between the three organelle concentrations.

Quantifying intermittent transport in cell cytoplasm

Active cellular transport is a fundamental mechanism for protein and vesicle delivery, cell cycle, and molecular degradation. Viruses can hijack the transport system and use it to reach the nucleus. Most transport processes consist of intermittent dynamics, where the motion of a particle, such as a virus, alternates between pure Brownian and directed movement along microtubules. In this Rapid Communication, we estimate the mean time for a particle to attach to a microtubule network. This computation leads to a coarse grained equation of the intermittent motion in radial and cylindrical geometries. Finally, by using the degradation activity inside the cytoplasm, we obtain refined asymptotic estimations for the probability and the mean time a virus reaches a small nuclear pore.

Microtubule-dependent distribution of mRNA in adult cardiocytes

Synthesis of myofibrillar proteins in the diffusion-restricted adult cardiocyte requires microtubule-based active transport of mRNAs as part of messenger ribonucleoprotein particles (mRNPs) to translation sites adjacent to nascent myofibrils. This is especially important for compensatory hypertrophy in response to hemodynamic overloading. The hypothesis tested here is that excessive microtubule decoration by microtubule-associated protein 4 (MAP4) after cardiac pressure overloading could disrupt mRNP transport and thus hypertrophic growth. MAP4-overexpressing and pressure-overload hypertrophied adult feline cardiocytes were infected with an adenovirus encoding zipcode-binding protein 1-enhanced yellow fluorescent protein fusion protein, which is incorporated into mRNPs, to allow imaging of these particles. Speed and distance of particle movement were measured via time-lapse microscopy. Microtubule depolymerization was used to study microtubule-based transport and distribution of mRNPs. Protein synthesis was assessed as radioautographic incorporation of [3H]phenylalanine. After microtubule depolymerization, mRNPs persist only perinuclearly and apparent mRNP production and protein synthesis decrease. Reestablishing microtubules restores mRNP production and transport as well as protein synthesis. MAP4 overdecoration of microtubules via adenovirus infection in vitro or following pressure overloading in vivo reduces the speed and average distance of mRNP movement. Thus cardiocyte microtubules are required for mRNP transport and structural protein synthesis, and MAP4 decoration of microtubules, whether directly imposed or accompanying pressure-overload hypertrophy, causes disruption of mRNP transport and protein synthesis. The dense, highly MAP4-decorated microtubule network seen in severe pressure-overload hypertrophy both may cause contractile dysfunction and, perhaps even more importantly, may prevent a fully compensatory growth response to hemodynamic overloading.

The contributions of microtubule stability and dynamic instability to adenovirus nuclear localization efficiency

Adenoviruses (Ads) utilize host cell microtubules to traverse the intracellular space and reach the nucleus in a highly efficient manner. Previous studies have shown that Ad infection promotes the formation of stable, posttranslationally modified microtubules by a RhoA-dependent mechanism. Ad infection also shifts key parameters of microtubule dynamic instability by a Rac1-dependent mechanism, resulting in microtubules with lower catastrophe frequencies, persistent growth phases, and a bias toward net growth compared to microtubules in uninfected cells. Until now it was unclear whether changes in RhoGTPase activity or microtubule dynamics had a direct impact on the efficiency of Ad microtubule-dependent nuclear localization. Here we have performed synchronous Ad infections and utilized confocal microscopy to analyze the individual contributions of RhoA activation, Rac1 activation, microtubule stability, dynamic behavior, and posttranslational modifications on Ad nuclear localization efficiency (NLE). We found that drug-induced suppression of microtubule dynamics impaired Ad NLE by disrupting the radial organization of the microtubule array. When the microtubule array was maintained, the suppression or enhancement of microtubule turnover did not significantly affect Ad NLE. Furthermore, RhoA activation or the formation of acetylated microtubules did not enhance Ad NLE. In contrast, active Rac1 was required for efficient Ad nuclear localization. Because Rac1 mediates persistent growth of microtubules to the lamellar regions of cells, we propose that Ad-induced activation of Rac1 enhances the ability of microtubules to "search and capture" incoming virus particles.

Viral strategies for intracellular trafficking: motors and microtubules

To overcome barriers to diffusion, many viruses utilize the microtubule-associated molecular motor cytoplasmic dynein 1 to drive transport towards the nucleus of a target cell. Cytoplasmic dynein 1 generates movement towards the minus end of microtubules located at the microtubule organizing centre (MTOC), a structure that is typically in close proximity to the nucleus. Physiological cargoes for cytoplasmic dynein include membranous organelles, protein complexes and aggregates of misfolded protein. In this review, we discuss the study of microtubule-based translocation of viruses and raise questions about the mechanisms for association with and then dissociation from cytoplasmic dynein with a goal of understanding whether viruses are seen by the intracellular trafficking machinery as functional protein complexes or misfolded protein aggregates.

Effect of inhibition of dynein function and microtubule-altering drugs on AAV2 transduction

Over the past decade, adeno-associated (AAV) virus has emerged as an important vector for gene therapy. As a result, understanding its basic biology, including intracellular trafficking, has become increasingly important. Here, we describe the effect of inhibiting dynein function or altering the state of microtubule polymerization on rAAV2 transduction. Overexpression of dynamitin, resulting in a functional inhibition of the minus-end-directed microtubule motor protein dynein, did not inhibit transduction. Equally, treatment of cells with nocodazole, or concentrations of vinblastine that result in the disruption of microtubules, had no significant effect on transduction. In contrast, high concentrations of Taxol and vinblastine, resulting in microtubule stabilization and the formation of tubulin paracrystals respectively, reduced rAAV2 transduction in a vector-dose-dependent manner. These results demonstrate that AAV2 can infect HeLa cells independently of dynein function or an intact microtubule network.

Cell-assisted assembly of colloidal crystallites

Many cells ingest foreign particles through a process known as phagocytosis. It now turns out that some cell types organize phagocytosed microparticles into crystalline arrays. Much like the classic crystallization of colloidal particles in a thermal bath, crystallization within the cell is driven by the spatial confinement of mutually repelling particles, in this case by the cell membrane. Cytoskeleton-driven motions exert a randomizing force, similar to but stronger than thermal forces; these motions anneal defects and purify the colloidal crystals within the cells. Bidisperse mixtures of microspheres phase separate within the cell, with the larger particles crystallizing around the nucleus and the smaller particles crystallizing around the perimeter of the large particle array. Mitochondria also participate in this kind of size segregation, which appears to be driven by membrane tension and curvature minimization. Beyond the curiosity of the phenomenon itself, cell-assisted colloidal assembly may prove useful as a new tool to study a variety of biophysical processes including cytoskeletal rearrangements, organelle-membrane interactions, the in vivo mechanics of microtubules, the cooperativity of molecular motors and intracellular traffic jams on cytoskeletal filaments.

Modeling of pattern regulation in melanophores

Melanosomes, pigment granules in melanophores, play a principal role in physiological color adaptation of fish and frog. Melanophores regulate melanosome trafficking on cytoskeletal filaments to generate a range of spatiotemporal patterns. Here, we present the first comprehensive model of spatiotemporal evolution of melanosome patterns. The model encompasses both physical and biochemical aspects of melanosome dynamics. It consists of (i) a kinetic description of biochemical reactions involved in intracellular signaling, (ii) a system of macroscopic reaction-diffusion-convection equations for melanosome concentration, and (iii) a set of constitutive relationships for coupling transport with the biochemical network. The model relates molecular-level regulatory actions to cell-level melanosome distribution, allowing unification of existing experimental observations and qualitative hypotheses into an integrated, consistent framework. The model reproduces salient features of melanosome patterns, both during transient and steady state. It gives useful insights into how cells coordinate motor-assisted transport to maintain and adapt spatial organization of intracellular organelles. In particular, we calculate the optimal transition paths from aggregation to dispersion in fish melanophores. The calculations suggest that fish melanophores optimally control intracellular signaling to maximize the efficiency of motor-assisted transport during dispersion.

Understanding intracellular transport processes pertinent to synthetic gene delivery via stochastic simulations and sensitivity analyses

A major challenge in synthetic gene delivery is to quantitatively predict the optimal design of polymer-based gene carriers (polyplexes). Here, we report a consistent, integrated, and fundamentally grounded computational methodology to address this challenge. This is achieved by accurately representing the spatio-temporal dynamics of intracellular structures and by describing the interactions between gene carriers and cellular components at a discrete, nanoscale level. This enables the applications of systems tools such as optimization and sensitivity analysis to search for the best combination of systems parameters. We validate the approach using DNA delivery by polyethylenimine as an example. We show that the cell topology (e.g., size, circularity, and dimensionality) strongly influences the spatiotemporal distribution of gene carriers, and consequently, their optimal intracellular pathways. The model shows that there exists an upper limit on polyplexes' intracellular delivery efficiency due to their inability to protect DNA until nuclear entry. The model predicts that even for optimally designed polyethylenimine vectors, only approximately 1% of total DNA is delivered to the nucleus. Based on comparison with gene delivery by viruses, the model suggests possible strategies to significantly improve transfection efficiencies of synthetic gene vectors.

Pharmacodynamics of adenovector distribution within the inner ear tissues of the mouse

Recent studies have demonstrated that delivery of genes to the inner ear can achieve a variety of effects ranging from support of auditory neuron survival to protection and restoration of hair cells, demonstrating the utility of vector based gene delivery. Translation of these findings to useful experimental systems or even clinical applications requires a detailed understanding of the pharmacokinetics of gene delivery in the inner ear. Ideal gene delivery systems will employ a well tolerated vector which efficiently transduces the appropriate target cells within a tissue, but spare non-target structures. Adenovectors based on serotype 5 (Ad 5) are commonly used vectors, are easy to construct and have a long track record of efficacious gene transfer in the inner ear. In this study we demonstrate that distribution of Ad5 vector occurs in a basal to apical gradient with rapid distribution of vector to the vestibule after delivery via a round window cochleostomy. Transduction of the vector and expression of the delivered transgene occurs by 10 min post vector delivery. At 24 h post delivery only 16% of vector that was initially detectable within the inner ear by quantitative PCR remained. Perilymph sampling was used to determine that vector concentrations in perilymph peaked at 30 min post delivery and then declined rapidly. Understanding these basic distribution patterns and parameters for delivery are important for the design of gene delivery vectors and vital for modeling dose responses to achieve safe efficacious delivery of a therapeutic agent.

Infection with replication-deficient adenovirus induces changes in the dynamic instability of host cell microtubules

Adenovirus translocation to the nucleus occurs through a well characterized minus end-directed transport along microtubules. Here, we show that the adenovirus infection process has a significant impact on the stability and dynamic behavior of host cell microtubules. Adenovirus-infected cells had elevated levels of acetylated and detyrosinated microtubules compared with uninfected cells. The accumulation of modified microtubules within adenovirus-infected cells required active RhoA. Adenovirus-induced changes in microtubule dynamics were characterized at the centrosome and at the cell periphery in living cells. Adenovirus infection resulted in a transient enhancement of centrosomal microtubule nucleation frequency. At the periphery of adenovirus-infected cells, the dynamic instability of microtubules plus ends shifted toward net growth, compared with the nearly balanced growth and shortening observed in uninfected cells. In infected cells, microtubules spent more time in growth, less time in shortening, and underwent catastrophes less frequently compared with those in uninfected cells. Drug-induced inhibition of Rac1 prevented most of these virus-induced shifts in microtubule dynamic instability. These results demonstrate that adenovirus infection induces a significant stabilizing effect on host cell microtubule dynamics, which involve, but are not limited to, the activation of the RhoGTPases RhoA and Rac1.

Theory of spatial patterns of intracellular organelles

Here we report on a generalized theory of spatial patterns of intracellular organelles, which are controlled by cells using cytoskeleton-based movements powered by molecular motors. The theory reveals that organelles exhibit one of the four distinct, stable patterns, namely aggregation, hyperdispersion, radial dispersion, and areal dispersion. Existence of specific patterns is determined by the contributions from three transport mechanisms, characterized by two Peclet numbers. The predicted patterns compare well with experimental data. This study provides a firm theoretical ground for classification of spatial patterns of organelles and understanding their regulation by cells.

A superhighway to virus infection

Microtubule-mediated transport of macromolecules and organelles (also known as "cargo") is essential for cells to function. Deficiencies in cytoplasmic transport are frequently associated with severe diseases and syndromes. Cytoplasmic transport also provides viruses with the means to reach their site of replication and is the route for newly assembled progeny to leave the infected cell. This parasitic relationship of viruses with the host cytoskeleton provides an excellent basis for cell biologists to unlock the secrets of cytoplasmic transport and unravel mechanisms of disease. Recent advances in live cell imaging and computational tracking of fluorescently labeled viruses are now revealing how complex the movements of single viruses are in infected cells. This review focuses on microtubule-based motility of viruses and highlights the mechanisms regulating cytoplasmic transport.