Molecular modeling and simulation of Mycobacterium tuberculosis cell wall permeability

The low permeability of the mycobacterial cell wall is thought to contribute to the intrinsic drug resistance of mycobacteria. In this study, the permeability of the Mycobacterium tuberculosis cell wall is studied by computer simulation. Thirteen known drugs with diverse chemical structures were modeled as solutes undergoing transport across a model for the M. tuberculosis cell wall. The properties of the solute-membrane complexes were investigated by means of molecular dynamics simulation, especially the diffusion coefficients of the solute molecules inside the cell wall. The molecular shape of the solute was found to be an important factor for permeation through the M. tuberculosis cell wall. Predominant lateral diffusion within, as opposed to transverse diffusion across, the membrane/cell wall system was observed for some solutes. The extent of lateral diffusion relative to transverse diffusion of a solute within a biological cell membrane may be an important finding with respect to absorption distribution, metabolism, elimination, and toxicity properties of drug candidates. Molecular similarity measures among the solutes were computed, and the results suggest that compounds having high molecular similarity will display similar transport behavior in a common membrane/cell wall environment. In addition, the diffusion coefficients of the solute molecules across the M. tuberculosis cell wall model were compared to those across the monolayers of dipalmitoylphosphatidylethanolamine and dimyristoylphosphatidylcholine, are two common phospholipids in bacterial and animal membranes. The differences among these three groups of diffusion coefficients were observed and analyzed.     

Printable microfluidic systems using pressure sensitive adhesive material for biosensing devices

Background

In biosensors with a fluid analyte, the integration of a microfluidic system, which guides the analyte into the sensing area, is critical. Quicker and economical ways to build up microfluidic systems will make point of care diagnostics viable. Printing is a low-cost technology that is increasingly used in emerging organic and flexible electronics and biosensors. In this paper, we present printed fluidic systems on flexible substrates made with pressure sensitive adhesive materials.

Methods

Printable pressure sensitive adhesive materials have been used for making microfluidic systems. Flexible substrates have been used, and two types of adhesive materials, one thermally dried and another UV curable, have been tested. Top sealing layer was laminated directly on top of the printed microfluidic structure. Flow tests were done with deionized water.

Results

Flow tests with deionized water show that both adhesive materials are suitable for capillary flow driven fluidic devices. Flow test using water as dielectric material was also done successfully on a printed electrolyte gated organic field effect transistor with an integrated microfluidic system.

General significance

Due to its ease of process and low cost, printed microfluidic system is believed to find more applications in biosensing devices. This article is part of a Special Issue entitled Organic Bioelectronics—Novel Applications in Biomedicine.     

Theoretical analysis of molecular diffusion in pressure-driven laminar flow in microfluidic channels

The T-sensor is a microfluidic analytical device that operates at low Reynolds numbers to ensure entirely laminar flow. Diffusion of molecules between streams flowing side by side may be observed directly. The pressure-driven velocity profile in the duct-shaped device influences diffusive transport in ways that affect the use of the T-sensor to measure molecular properties. The primary effect is a position-dependent variation in the extent of diffusion that occurs due to the distribution of residence time among different fluid laminae. A more detailed characterization reveals that resultant secondary concentration gradients yield variations in the scaling behavior between diffusive displacement and elapsed time in different regions of the channel. In this study, the time-dependent evolution of analyte distribution has been quantified using a combination of one- and two-dimensional models. The results include an accurate portrayal of the shape of the interdiffusion region in a representative T-sensor assay, calculation of the diffusive scaling law across the width of the channel, and quantification of artifacts that occur when making diffusion coefficient measurements in the T-sensor.     

Finite Element Model of Oxygen Transport for the Design of Geometrically Complex Microfluidic Devices Used in Biological Studies

Red blood cells play a crucial role in the local regulation of oxygen supply in the microcirculation through the oxygen dependent release of ATP. Since red blood cells serve as an oxygen sensor for the circulatory system, the dynamics of ATP release determine the effectiveness of red blood cells to relate the oxygen levels to the vessels. Previous work has focused on the feasibility of developing a microfluidic system to measure the dynamics of ATP release. The objective was to determine if a steep oxygen gradient could be developed in the channel to cause a rapid decrease in hemoglobin oxygen saturation in order to measure the corresponding levels of ATP released from the red blood cells. In the present study, oxygen transport simulations were used to optimize the geometric design parameters for a similar system which is easier to fabricate. The system is composed of a microfluidic device stacked on top of a large, gas impermeable flow channel with a hole to allow gas exchange. The microfluidic device is fabricated using soft lithography in polydimethyl-siloxane, an oxygen permeable material. Our objective is twofold: (1) optimize the parameters of our system and (2) develop a method to assess the oxygen distribution in complex 3D microfluidic device geometries. 3D simulations of oxygen transport were performed to simulate oxygen distribution throughout the device. The simulations demonstrate that microfluidic device geometry plays a critical role in molecule exchange, for instance, changing the orientation of the short wide microfluidic channel results in a 97.17% increase in oxygen exchange. Since microfluidic devices have become a more prominent tool in biological studies, understanding the transport of oxygen and other biological molecules in microfluidic devices is critical for maintaining a physiologically relevant environment. We have also demonstrated a method to assess oxygen levels in geometrically complex microfluidic devices.     

Magnetic resonance imaging of a microvascular-interstitium model on a microfluidic device

We report a microvascular-interstitium model on microfluidic devices to study leakage of drugs from blood vessels under in vivo-like flow conditions. We employed magnetic resonance imaging to demonstrate the compatibility of the model for experimental animals and humans. We observed transport of two types of different molecular-weight contrast agents into the model interstitium. The ratio of the transport rates of agents agreed with the ratio calculated from diffusion coefficients of the agents. We expect that the model will be useful for the estimation and evaluation of leakage of many kinds of agents in vivo.     

Local macromolecule diffusion coefficients in structurally non-uniform bacterial biofilms using fluorescence recovery after photobleaching (FRAP)

Pure culture Pseudomonas putida biofilms were cultivated under controlled conditions to a desired overall biofilm thickness, then employed within classical half-cell diffusion chambers to estimate, from transient solute concentrations, the effective diffusion coefficient for several macromolecules of increasing molecular weight and molecular complexity. Results of traditional half-cell studies were found to be erroneous due to the existence of microscopic water channels or crevasses that perforate the polysaccharidic gel matrix of the biofilm, sometimes completely to the supporting substratum. Thus, half-cell devices measure a composite transfer coefficient that may overestimate the true, local flux of solutes in the biofilm polysaccharide gel matrix. An alternative analytical technique was refined to determine the local diffusion coefficients on a micro-scale to avoid the errors created by the biofilm architectural irregularities. This technique is based upon the Fluorescence Return After Photobleaching (FRAP), which allows image analysis observation of the transport of fluorescently labeled macromolecules as they migrate into a micro-scale photobleached zone. The technique can be computerized and allows one to map the local diffusion coefficients of various solute molecules at different horizontal planes and depths in a biofilm. These mappings also indirectly indicate the distribution of water channels in the biofilm, which was corroborated independently by direct microscopic observation of the settling of fluorescently-labeled latex spheres within the biofilm. Fluorescence return after photobleaching results indicate a significant reduction in the solute transport coefficients in biofilm polymer gel vs. the same value in water, with the reduction being dependent on solute molecule size and shape.     

Comparative modeling and analysis of microfluidic and conventional DNA microarrays

A theoretical analysis was developed to predict molecular hybridization rates for microarrays where samples flow through microfluidic channels and for conventional microarrays where samples remain stationary during hybridization. The theory was validated by using a multiplexed microfluidic microarray where eight samples were hybridized simultaneously against eight probes using 60-mer DNA strands. Mass transfer coefficients ranged over three orders of magnitude where either kinetic reaction rates or molecular diffusion rates controlled overall hybridization rates. Probes were printed using microfluidic channels and also conventional spotting techniques. Consistent with the theoretical model, the microfluidic microarray demonstrated the ability to print DNA probes in less than 1 min and to detect 10-pM target concentrations with hybridization times in less than 5 min.     

Crossover from picosecond collective to single particle dynamics defines the mechanism of lateral lipid diffusion

It has been widely accepted that the thermally excited motions of the molecules in a cell membrane is the prerequisite for a cell to carry its biological functions. On the other hand, the detailed mapping of the ultrafast picosecond single-molecule and the collective lipid dynamics in a cell membrane remains rather elusive. Here, we report all-atom molecular dynamics simulations of a 1,2-dipalmitoyl-sn-glycero-3-phosphocholine bilayer over a wide range of temperature. We elucidate a molecular mechanism underlying the lateral lipid diffusion in a cell membrane across the gel, rippled, and liquid phases using an analysis of the longitudinal and transverse current correlation spectra, the velocity auto-correlation functions, and the molecules mean square displacements. The molecular mechanism is based on the anomalous ultrafast vibrational properties of lipid molecules at the viscous-to-elastic crossover. The macroscopic lipid diffusion coefficients predicted by the proposed diffusion model are in a good agreement with experimentally observed values. Furthermore, we unveil the role of water confined at the water-lipid interface in triggering collective vibrations in a lipid bilayer.     

Computational Fluid Dynamics Simulations at Micro-Scale Stenosis for Microfluidic Thrombosis Model Characterization

Platelet aggregation plays a central role in pathological thrombosis, preventing healthy physiological blood flow within the circulatory system. For decades, it was believed that platelet aggregation was primarily driven by soluble agonists such as thrombin, adenosine diphosphate and thromboxane A2. However, recent experimental findings have unveiled an intriguing but complementary biomechanical mechanism—the shear rate gradients generated from flow disturbance occurring at sites of blood vessel narrowing, otherwise known as stenosis, may rapidly trigger platelet recruitment and subsequent aggregation. In our Nature Materials 2019 paper [1], we employed microfluidic devices which incorporated micro-scale stenoses to elucidate the molecular insights underlying the prothrombotic effect of blood flow disturbance. Nevertheless, the rheological mechanisms associated with this stenotic microfluidic device are poorly characterized. To this end, we developed a computational fluid dynamics (CFD) simulation approach to systematically analyze the hemodynamic influence of bulk flow mechanics and flow medium. Grid sensitivity studies were performed to ensure accurate and reliable results. Interestingly, the peak shear rate was significantly reduced with the device thickness, suggesting that fabrication of microfluidic devices should retain thicknesses greater than 50 µm to avoid unexpected hemodynamic aberration, despite thicker devices raising the cost of materials and processing time of photolithography. Overall, as many groups in the field have designed microfluidic devices to recapitulate the effect of shear rate gradients and investigate platelet aggregation, our numerical simulation study serves as a guideline for rigorous design and fabrication of microfluidic thrombosis models.     

A microfluidic system for investigation of extravascular transport and cellular uptake of drugs in tumors

Three-dimensional (3D) tumor models have been established in various microfluidic systems for drug delivery and resistance studies in vitro. However, one of the main drawbacks of these models is non-uniform distribution of cells, leaving regions with very low cell density within the 3D structures. As a result, molecular diffusion in the cell compartments is faster than that observed in solid tumors. To solve this problem, we developed a new technique for preparation of 3D tumor models in vitro. It was based on a microfluidic device containing three parallel channels separated by narrowly spaced posts. Tumor cells were loaded into the central channel at high density. To test the system, B16.F10 melanoma cells were perfusion-cultured overnight and the resulting 3D structure was characterized in terms of viability, density, and morphology of cells as well as transport properties of small fluorescent molecules. Immediately upon loading of tumor cells, the cell density was comparable to those observed in B16.F10 tumor tissues in vivo; and the viability of tumor cells was maintained through the overnight culture. The tumor model displayed low extracellular space and high resistance to diffusion of small molecules. For membrane-permeant molecules (e.g., Hoechst 33342), the rate of interstitial penetration was extremely slow, compared to membrane-impermeant molecules (e.g., sodium fluorescein). This versatile tumor model could be applied to in vitro studies of transport and cellular uptake of drugs and genes.     

Size-Dependent Diffusion of Dextrans in Excised Porcine Corneal Stroma

Delivery of therapeutic agents to the eye requires efficient transport through cellular and extracellular barriers. We evaluated the rate of diffusive transport in excised porcine corneal stroma using fluorescently labeled dextran molecules with hydrodynamic radii ranging from 1.3 to 34 nm. Fluorescence correlation spectroscopy (FCS) was used to measure diffusion coefficients of dextran molecules in the excised porcine corneal stroma. The preferential sensitivity of FCS to diffusion along two dimensions was used to differentially probe diffusion along the directions parallel to and perpendicular to the collagen lamellae of the corneal stroma. In order to develop an understanding of how size affects diffusion in cornea, diffusion coefficients in cornea were compared to diffusion coefficients measured in a simple buffer solution. Dextran molecules diffuse more slowly in cornea as compared to buffer solution. The reduction in diffusion coefficient is modest however (67% smaller), and is uniform over the range of sizes that we measured. This indicates that, for dextrans in the 1.3 to 34 nm range, the diffusion landscape of corneal stroma can be represented as a simple liquid with a viscosity approximately 1.5 times that of water. Diffusion coefficients measured parallel vs. perpendicular to the collagen lamellae were indistinguishable. This indicates that diffusion in the corneal stroma is not highly anisotropic. Our results support the notion that the corneal stroma is highly permeable and isotropic to transport of hydrophilic molecules and particles with hydrodynamic radii up to at least 34 nm.     

SPR biosensing coupled to a digital microfluidic microstreaming system

This article reports on a proof-of-concept system composed of a droplet based surface plasmon resonance (SPR) system coupled to a surface acoustic wave (SAW) microfluidic plateform. It is now well established that surface based binding analyses such as SPR are highly influenced by the transport of analyte to the sensing surface. Further, obtaining reliable equilibrium in flow cells to realize quantification studies is not straightforward. An original solution compared to generally used pressure driven flows is then proposed to favourably cope with these issues. Efficiency of SAW microstreaming coupled to SPR biosensing is considered, in order to improve the accuracy of kinetic parameter estimation in mass transport limited regime and to realize reliable quantification studies. First, the droplet based SPR technique and its advantages are presented. Then, the integration of the microstreaming on the system is discussed. Streptavidin binding is then monitored in static mode and under SAW streaming mode.     

A molecular machine biosensor: construction, predictive models and experimental studies

This paper describes the construction, operation and predictive modeling of a molecular machine, functioning as a high sensitivity biosensor. Embedded gramicidin A (gA) ionchannels in a self-assembled tethered lipid bilayer act as biological switches in response to target molecules and provide a signal amplification mechanism that results in high sensitivity molecular detection. The biosensor can be used as a rapid and sensitive point of care diagnostic device in different media such as human serum, plasma and whole blood without the need for pre and post processing steps required in an enzyme-linked immunosorbent assay. The electrical reader of the device provides the added advantage of objective measurement. Novel ideas in the construction of the molecular machine, including fabrication of biochip arrays, and experimental studies of its ability to detect analyte molecules over a wide range of concentrations are presented. Remarkably, despite the complexity of the device, it is shown that the response can be predicted by modeling the analyte fluid flow and surface chemical reactions. The derived predictive models for the sensing dynamics also facilitate determining important variables in the design of a molecular machine such as the ion channel lifetime and diffusion dynamics within the bilayer lipid membrane as well as the bio-molecular interaction rate constants.     

Fringe pattern photobleaching, a new method for the measurement of transport coefficients of biological macromolecules.

The conventional method of studying mass transport in membranes by spot photobleaching and then following the recovery of fluorescence has disadvantages. Among them, the need for a high density of fluorescent molecules, the measurement of the beam profile, and a knowledge of the photobleaching processes are of a crucial importance. The application of a planar fringe pattern of light both for the bleaching and the monitoring of the fluorescent molecules solves these three major difficulties. Brownian diffusion coefficients and flow velocities can be measured independently and are averaged over the whole fringe pattern volume. These transport coefficients are explored over the wide range of experimentally accessible distances (from an interfringe spacing 0.5-50 micron). The quantification of the mobile and immobile components is further simplified by scanning the fringe pattern and detecting only a modulated fluorescence recovery signal. The fringe pattern photobleaching method is particularly adapted to the measurements of diffusion coefficients and flow velocity of membrane components, as well as of cytoplasmic proteins. The theoretical results and the test experiments with fluorescent bovine serum albumin are described.     

Microfluidic delivery of small molecules into mammalian cells based on hydrodynamic focusing

Microfluidics-based cell assays offer high levels of automation and integration, and allow multiple assays to be run in parallel, based on reduced sample volumes. These characteristics make them attractive for studies associated with drug discovery. Controlled delivery of drug molecules or other exogenous materials into cells is a critical issue that needs to be addressed before microfluidics can serve as a viable platform for drug screening and studies. In this study, we report the application of hydrodynamic focusing for controlled delivery of small molecules into cells immobilized on the substrate of a microfluidic device. We delivered calcein AM which was permeant to the cell membrane into cells, and monitored its enzymatic conversion into fluorescent calcein during and after the delivery. Different ratios of the sample flow to the side flow were tested to determine how the conditions of hydrodynamic focusing affected the delivery. A 3D numerical model was developed to help understand the fluid flow, molecular diffusion due to hydrodynamic focusing in the microfluidic channel. The results from the simulation indicated that the calcein AM concentration on the outer surface of a cell was determined by the conditions of hydrodynamic focusing. By comparing the results from the simulation with those from the experiment, we found that the calcein AM concentration on the cell outer surface correlated very well with the amount of the molecules delivered into the cell. This suggests that hydrodynamic focusing provides an effective way for potentially quantitative delivery of exogenous molecules into cells at the single cell or subcellular level. We expect that our technique will pave the way to high-throughput drug screening and delivery on a microfluidic platform.     

Imaging Barriers to Diffusion by Pair Correlation Functions

Molecular diffusion and transport are fundamental processes in physical, chemical, biochemical, and biological systems. However, current approaches to measure molecular transport in cells and tissues based on perturbation methods such as fluorescence recovery after photobleaching are invasive, fluctuation correlation methods are local, and single-particle tracking requires the observation of isolated particles for relatively long periods of time. We propose to detect molecular transport by measuring the time cross-correlation of fluctuations at a pair of locations in the sample. When the points are farther apart than two times the size of the point spread function, the maximum of the correlation is proportional to the average time a molecule takes to move from a specific location to another. We demonstrate the method by simulations, using beads in solution, and by measuring the diffusion of molecules in cellular membranes. The spatial pair cross-correlation method detects barriers to diffusion and heterogeneity of diffusion because the time of the correlation maximum is delayed in the presence of diffusion barriers. This noninvasive, sensitive technique follows the same molecule over a large area, thereby producing a map of molecular flow. It does not require isolated molecules, and thus many molecules can be labeled at the same time and within the point spread function.     

Apparent Diffusivity and Taylor Dispersion of Water and Solutes in Capillary Beds

A physical theory explaining the anisotropic dispersion of water and solutes in biological tissues is introduced based on the phenomena of Taylor dispersion, in which highly diffusive solutes cycle between flowing and stagnant regions in the tissue, enhancing dispersion in the direction of microvascular flow. An effective diffusion equation is derived, for which the coefficient of dispersion in the axial direction (direction of capillary orientation) depends on the molecular diffusion coefficient, tissue perfusion, and vessel density. This analysis provides a homogenization that represents three-dimensional transport in capillary beds as an effectively one-dimensional phenomenon. The derived dispersion equation may be used to simulate the transport of solutes in tissues, such as in pharmacokinetic modeling. In addition, the analysis provides a physically based hypothesis for explaining dispersion anisotropy observed in diffusion-weighted imaging (DWI) and diffusion-tensor magnetic resonance imaging (DTMRI) and suggests the means of obtaining quantitative functional information on capillary vessel density from measurements of dispersion coefficients. It is shown that a failure to account for flow-mediated dispersion in vascular tissues may lead to misinterpretations of imaging data and significant overestimates of directional bias in molecular diffusivity in biological tissues. Measurement of the ratio of axial to transverse diffusivity may be combined with an independent measurement of perfusion to provide an estimate of capillary vessel density in the tissue.     

Production and study of a high-temperature plasma in the central solenoid of the AMBAL-M device

Results are presented from experiments on the production and study of a hot dense plasma in the central solenoid of the AMBAL-M fully axisymmetric ambipolar magnetic confinement system. The hot plasma in the solenoid and end cell is produced by filling the system with a thermally insulated current-carrying plasma stream with developed low-frequency turbulence. The plasma stream is generated by a gas-discharge plasma source placed upstream from the magnetic mirror of the solenoid. As a result, an MHD-stabilized plasma with a length of 6 m, a diameter of 40 cm, a density of 2×10¹³ cm^?3, an ion energy of 250 eV, and an electron temperature of 60 eV is produced in the central solenoid. It is found that, in the quiescent decay phase, transverse plasma losses from the solenoid due to low-frequency oscillations and nonambipolar transport are rather small and comparable with the classical diffusion losses.     

Towards mimicking natural protein channels with aligned carbon nanotube membranes for active drug delivery

AimsCarbon nanotube (CNT) membranes offer an exciting opportunity to mimic natural protein channels due to 1) a mechanism of dramatically enhanced fluid flow 2) ability to place ‘gatekeeper’ chemistry at the entrance to pores 3) the ability for biochemical reactions to occur on gatekeeper molecules and 4) an ability to chemically functionalize each side of the membrane independently.Main methodsAligned CNT membranes were fabricated and CNT pore entrances modified with gatekeeper chemistry. Pressure driven fluid flow and diffusion experiments were performed to study the mechanisms of transport through CNTs.Key findingsThe transport mechanism through CNT membranes is primarily 1) ionic diffusion near bulk expectation 2) gas flow enhanced 1–2 orders of magnitude primarily due to specular reflection 3) fluid flow 4–5 orders of magnitude faster than conventional materials due to a nearly ideal slip-boundary interface. The transport can be modulated by ‘gatekeeper’ chemistry at the pore entrance using steric hindrance, electrostatic attraction/repulsion, or biochemical state. The conformation of charged tethered molecules can be modulated by applied bias setting the stage for programmable drug release devices.SignificanceThe membrane structure is mechanically far more robust than lipid bilayer films, allowing for large-scale chemical separations, delivery or sensing based on the principles of protein channels. The performance of protein channels is several orders of magnitude faster than conventional membrane materials. The fundamental requirements of mimicking protein channels are present in the CNT membrane system.