Your browser doesn't support javascript.
loading
Show: 20 | 50 | 100
Results 1 - 20 de 167
Filter
1.
Article in English | MEDLINE | ID: mdl-38676536

ABSTRACT

A primary cilium, made of nine microtubule doublets enclosed in a cilium membrane, is a mechanosensing organelle that bends under an external mechanical load and sends an intracellular signal through transmembrane proteins activated by cilium bending. The nine microtubule doublets are the main load-bearing structural component, while the transmembrane proteins on the cilium membrane are the main sensing component. No distinction was made between these two components in all existing models, where the stress calculated from the structural component (nine microtubule doublets) was used to explain the sensing location, which may be totally misleading. For the first time, we developed a microstructure-based primary cilium model by considering these two components separately. First, we refined the analytical solution of bending an orthotropic cylindrical shell for individual microtubule, and obtained excellent agreement between finite element simulations and the theoretical predictions of a microtubule bending as a validation of the structural component in the model. Second, by integrating the cilium membrane with nine microtubule doublets and simulating the tip-anchored optical tweezer experiment on our computational model, we found that the microtubule doublets may twist significantly as the whole cilium bends. Third, besides being cilium-length-dependent, we found the mechanical properties of the cilium are also highly deformation-dependent. More important, we found that the cilium membrane near the base is not under pure in-plane tension or compression as previously thought, but has significant local bending stress. This challenges the traditional model of cilium mechanosensing, indicating that transmembrane proteins may be activated more by membrane curvature than membrane stretching. Finally, we incorporated imaging data of primary cilia into our microstructure-based cilium model, and found that comparing to the ideal model with uniform microtubule length, the imaging-informed model shows the nine microtubule doublets interact more evenly with the cilium membrane, and their contact locations can cause even higher bending curvature in the cilium membrane than near the base.

2.
Soft Matter ; 20(3): 599-608, 2024 Jan 17.
Article in English | MEDLINE | ID: mdl-38131477

ABSTRACT

We numerically investigate the hydrodynamics and membrane dynamics of a multicomponent vesicle in two strongly confined geometries. This serves as a simplified model for red blood cells undergoing large deformations while traversing narrow constrictions. We propose a new parameterization for the bending modulus that remains positive for all lipid phase parameter values. For a multicomponent vesicle passing through a stenosis, we establish connections between various properties: lipid phase coarsening, size and flow profile of the lubrication layers, excess pressure, and the tank-treading velocity of the membrane. For a multicomponent vesicle passing through a contracting channel, we find that the lipid always phase separates so that the vesicle is stiffer in the front as it passes through the constriction. For both cases of confinement we find that lipid coarsening is arrested under strong confinement, and resumes at a high rate upon relief from extreme confinement. The results may be useful for efficient sorting lipid domains using microfluidic flows by controlled release of vesicles passing through strong confinement.

3.
Drug Dev Res ; 84(8): 1724-1738, 2023 Dec.
Article in English | MEDLINE | ID: mdl-37756467

ABSTRACT

Signal transducer and activator of transcription 3 (STAT3) and STAT5 are the transcription factors that have been studied extensively in relevance to the development of cancers in humans. Suppression of either STAT3 or STAT5-mediated signaling events has been demonstrated to be effective in inducing cytotoxicity in cancer cells. Herein, new hybrids of triazolyl-indolo-quinoxaline are synthesized and examined for their effect on the activation of STAT3 and STAT5 pathways in gastric cancer (GC) cells. Among the newly synthesized compounds, 2,3-difluoro-6-((1-(3-fluorophenyl)-1H-1,2,3-triazol-5-yl)methyl)-6H-indolo[2,3-b]quinoxaline (DTI) displayed selective cytotoxicity against GC cells over their normal counterpart. Flow cytometric analysis, annexin-V-fluorescein isothiocyanate staining, terminal deoxynucleotidyl transferase dUTP nick end labeling assay, live and dead assay, and caspase activation experiments suggested DTI as a potent inducer of apoptosis. The mechanistic approach revealed that DTI imparts cytotoxicity via downregulating the phosphorylation of STAT3Y705 and STAT5Y694/699 . DTI significantly reduced the nuclear pool of STAT3/STAT5 and reduced the DNA interaction ability of STAT3/STAT5 as evidenced by immunofluorescence and electrophoretic mobility shift assay. Further investigation revealed that inhibitory effects towards STAT proteins were mediated through the suppression of upstream kinases such as JAK1, JAK2, and Src. Treatment of GC cells with pervanadate counteracted the DTI-driven STAT3/STAT5 inhibition suggesting the involvement of tyrosine phosphatase. Upon DTI exposure, there was a significant upregulation in the mRNA and protein expression of PTPεC, which is a negative regulator of the JAK-STAT pathway. Knockdown of PTPεC suppressed the DTI-induced STATs inhibition in GC cells. Taken together, triazolyl-indolo-quinoxaline is presented as a new inhibitor of the STAT3/STAT5 pathway in GC cells.


Subject(s)
Signal Transduction , Stomach Neoplasms , Humans , STAT5 Transcription Factor/metabolism , STAT5 Transcription Factor/pharmacology , STAT3 Transcription Factor/metabolism , DNA-Binding Proteins/metabolism , Trans-Activators , Up-Regulation , Quinoxalines/pharmacology , Janus Kinases/metabolism , Janus Kinases/pharmacology , STAT Transcription Factors/metabolism , STAT Transcription Factors/pharmacology , Phosphorylation , Apoptosis
4.
bioRxiv ; 2023 Jul 15.
Article in English | MEDLINE | ID: mdl-37503231

ABSTRACT

A primary cilium, made of nine microtubule doublets enclosed in a cilium membrane, is a mechanosensing organelle that bends under an external mechanical load and sends an intracellular signal through transmembrane proteins activated by cilium bending. The nine microtubule doublets are the main load-bearing structural component, while the transmembrane proteins on the cilium membrane are the main sensing component. No distinction was made between these two components in all existing models, where the stress calculated from the structural component (nine microtubule doublets) was used to explain the sensing location, which may be totally misleading. For the first time, we developed a microstructure-based primary cilium model by considering these two components separately. First, we refined the analytical solution of bending an orthotropic cylindrical shell for individual microtubule, and obtained excellent agreement between finite element simulations and the theoretical predictions of a microtubule bending as a validation of the structural component in the model. Second, by integrating the cilium membrane with nine microtubule doublets, we found that the microtubule doublets may twist significantly as the whole cilium bends. Third, besides being cilium-length-dependent, we found the mechanical properties of the cilium are also highly deformation-dependent. More important, we found that the cilium membrane near the base is not under pure in-plane tension or compression as previously thought, but has significant local bending stress. This challenges the traditional model of cilium mechanosensing, indicating that transmembrane proteins may be activated more by membrane curvature than membrane stretching. Finally, we incorporated imaging data of primary cilia into our microstructure-based cilium model, and found that comparing to the ideal model with uniform microtubule length, the imaging-informed model shows the nine microtubule doublets interact more evenly with the cilium membrane, and their contact locations can cause even higher bending curvature in the cilium membrane than near the base. SIGNIFICANCE: Factors regulating the mechanical response of a primary cilium to fluid flow remain unclear. Modeling the microtubule doublet as a composite of two orthotropic shells and the ciliary axoneme as an elastic shell enclosing nine such microtubule doublets, we found that the length distribution of microtubule doublets (inferred from cryogenic electron tomography images) is the primary determining factor in the bending stiffness of primary cilia, rather than just the ciliary length. This implies ciliary-associated transmembrane proteins may be activated by membrane curvature changes rather than just membrane stretching. These insights challenge the traditional view of ciliary mechanosensation and expands our understanding of the different ways in which cells perceive and respond to mechanical stimuli.

5.
J Public Health (Oxf) ; 45(2): 529-534, 2023 06 14.
Article in English | MEDLINE | ID: mdl-37326349

ABSTRACT

BACKGROUND: Public Health registrars (SpRs) were an important component of the workforce that contributed to the COVID-19 response. This study explores their contribution and the impact the early stages of the pandemic had on their learning and training. METHODS: Data were collected from SpRs in the London and Kent, Surrey, Sussex training programme between July and September 2020 through a mixture of questionnaires and semi-structured interviews. A thematic analysis of interview transcripts was undertaken to identify themes. RESULTS: 35/128 SpRs responded to the survey and 11 were interviewed. SpRs were placed across a range of organizations and made a significant contribution to the COVID-19 response. Overall, SpRs learned important skills but working on the response may for some have impacted negatively on training. A number of facilitators and barriers to learning were identified. CONCLUSION: The study findings highlight the opportunities for learning created by the pandemic. However, changing projects and the desire of SpRs to contribute to the response meant the impacts on training were mixed. Future deployment of SpRs should consider the balance of responsibility and pace when delegating work, as well as the need to supervise effectively and support remote working to maintain good mental wellbeing.


Subject(s)
COVID-19 , Pandemics , Humans , London/epidemiology , Public Health , COVID-19/epidemiology , Learning
6.
Soft Matter ; 19(4): 776-789, 2023 Jan 25.
Article in English | MEDLINE | ID: mdl-36625263

ABSTRACT

Contaminants and other agents are often present at the interface between two fluids, giving rise to rheological properties such as surface shear and dilatational viscosities. The dynamics of viscous drops with interfacial viscosities has attracted greater interest in recent years, due to the influence of surface rheology on deformation and the surrounding flows. We investigate the effects of shear and dilatational viscosities on the electro-deformation of a viscous drop using the Taylor-Melcher leaky dielectric model. We use a large deformation analysis to derive an ordinary differential equation for the drop shape. Our model elucidates the contributions of each force to the overall deformation of the drop and reveals a rich range of dynamic behaviors that show the effects of surface viscosities and their dependence on rheological and electrical properties of the system. We also examine the physical mechanisms underlying the observed behaviors by analyzing the surface dilatation and surface deformation.

7.
Soft Matter ; 18(25): 4786-4791, 2022 Jun 29.
Article in English | MEDLINE | ID: mdl-35708007

ABSTRACT

Nanoscale phenomena such as surface hydration and the molecular layering of liquids under strong nanoscale confinement play a critical role in liquid-mediated surface adhesion that is not accounted for by available models, which assume a uniform liquid density with or without considering surface forces and associated disjoining pressure effects. This work introduces an alternative theoretical description that via the potential of mean force (PMF) considers the strong spatial variation of the liquid number density under nanoscale confinement. This alternative description based on the PMF predicts a dual effect of surface hydration by producing: (i) strong spatial oscillations of the local liquid density and pressure and, more importantly, (ii) a configuration-dependent liquid-solid surface energy under nanoscale confinement. Theoretical analysis and molecular dynamics simulations for the case of an axisymmetric water bridge with nanoscale heights show that the latter hydration effect is critical for the accurate prediction of the surface energy and adhesion forces when a small volume of liquid is nanoscopically confined by two surfaces approaching contact.

8.
J Cell Biochem ; 123(7): 1222-1236, 2022 07.
Article in English | MEDLINE | ID: mdl-35621239

ABSTRACT

Epithelial-mesenchymal transition (EMT) is a key process, which can promote the transition of tumor cells into other organs by weakening the cell-cell junctions. Tumor cell invasion and metastasis arising because of EMT can determine the prognosis of cancer. EMT can be induced by several growth factors including transforming growth factor-ß (TGF-ß), which can exert their effects by affecting several cell-signaling pathways. Fangchinoline (FCN), a kind of bisbenzylisoquinoline, belongs to the family Menispermaceae. FCN can display substantial antitumor effects against various malignant cell lines but its possible impact on EMT has not been explored. We examined the potential impact of FCN in affecting the activation of EMT in human colon cancer cells. We evaluated the influence of FCN on EMT in colon cancer cells by using Western blot analysis and reverse transcription-polymerase chain reaction assays. The cellular invasion and migration were observed by Boyden chamber and wound healing assays. Thereafter, the effect of the drug on proliferation and invasion was also evaluated by real-time cell analysis. FCN suppressed the levels of TGF-ß-induced mesenchymal markers, such as fibronectin, vimentin, MMP-9, MMP-2, N-cadherin, Twist, and Snail. However, FCN markedly enhanced the expression of epithelial markers such as occludin and E-cadherin. These results imply that FCN can potentially inhibit tumor metastasis through abrogating EMT. In addition, FCN downregulated c-Met/PI3K/Akt/mTOR and Wnt/ß-catenin cell signaling pathways and mitigated tumor migration as well as invasion. Overall, our study suggests a potential novel role of FCN as an antimetastatic agent against human colon cancer cells.


Subject(s)
Benzylisoquinolines , Colonic Neoplasms , Benzylisoquinolines/pharmacology , Cell Line, Tumor , Cell Movement , Colonic Neoplasms/drug therapy , Epithelial-Mesenchymal Transition , Humans , Phosphatidylinositol 3-Kinases , Signal Transduction , Transforming Growth Factor beta/pharmacology
9.
Soft Matter ; 18(3): 554-565, 2022 Jan 19.
Article in English | MEDLINE | ID: mdl-34931640

ABSTRACT

Sickle cell anemia (SCA) is a disease that affects red blood cells (RBCs). Healthy RBCs are highly deformable objects that under flow can penetrate blood capillaries smaller than their typical size. In SCA there is an impaired deformability of some cells, which are much stiffer and with a different shape than healthy cells, and thereby affect regular blood flow. It is known that blood from patients with SCA has a higher viscosity than normal blood. However, it is unclear how the rigidity of cells is related to the viscosity of blood, in part because SCA patients are often treated with transfusions of variable amounts of normal RBCs and only a fraction of cells will be stiff. Here, we report systematic experimental measurements of the viscosity of a suspension varying the fraction of rigid particles within a suspension of healthy cells. We also perform systematic numerical simulations of a similar mixed suspension of soft RBCs, rigid particles, and their hydrodynamic interactions. Our results show that there is a rheological signature within blood viscosity to clearly identify the fraction of rigidified cells among healthy deformable cells down to a 5% volume fraction of rigidified cells. Although aggregation of RBCs is known to affect blood rheology at low shear rates, and our simulations mimic this effect via an adhesion potential, we show that such adhesion, or aggregation, is unlikely to provide a physical rationalization for the viscosity increase observed in the experiments at moderate shear rates due to rigidified cells. Through numerical simulations, we also highlight that most of the viscosity increase of the suspension is due to the rigidity of the particles rather than their sickled or spherical shape. Our results are relevant to better characterize SCA, provide useful insights relevant to rheological consequences of blood transfusions, and, more generally, extend to the rheology of mixed suspensions having particles with different rigidities, as well as offering possibilities for developments in the field of soft material composites.


Subject(s)
Anemia, Sickle Cell , Blood Viscosity , Erythrocytes , Humans , Rheology , Viscosity
10.
Math Biosci Eng ; 18(3): 2849-2881, 2021 Mar 25.
Article in English | MEDLINE | ID: mdl-33892575

ABSTRACT

Active fluids consume fuel at the microscopic scale, converting this energy into forces that can drive macroscopic motions over scales far larger than their microscopic constituents. In some cases, the mechanisms that give rise to this phenomenon have been well characterized, and can explain experimentally observed behaviors in both bulk fluids and those confined in simple stationary geometries. More recently, active fluids have been encapsulated in viscous drops or elastic shells so as to interact with an outer environment or a deformable boundary. Such systems are not as well understood. In this work, we examine the behavior of droplets of an active nematic fluid. We study their linear stability about the isotropic equilibrium over a wide range of parameters, identifying regions in which different modes of instability dominate. Simulations of their full dynamics are used to identify their nonlinear behavior within each region. When a single mode dominates, the droplets behave simply: as rotors, swimmers, or extensors. When parameters are tuned so that multiple modes have nearly the same growth rate, a pantheon of modes appears, including zigzaggers, washing machines, wanderers, and pulsators.

11.
Math Biosci Eng ; 18(2): 1215-1237, 2021 01 15.
Article in English | MEDLINE | ID: mdl-33757184

ABSTRACT

Primary cilia are non-motile, solitary (one per cell) microtubule-based organelles that emerge from the mother centriole after cells have exited the mitotic cycle. Identified as a mechanosensing organelle that responds to both mechanical and chemical stimuli, the primary cilium provides a fertile ground for integrative investigations of mathematical modeling, numerical simulations, and experiments. Recent experimental findings revealed considerable complexity to the underlying mechanosensory mechanisms that transmit extracellular stimuli to intracellular signaling many of which include primary cilia. In this invited review, we provide a brief survey of experimental findings on primary cilia and how these results lead to various mathematical models of the mechanics of the primary cilium bent under an external forcing such as a fluid flow or a trap. Mathematical modeling of the primary cilium as a fluid-structure interaction problem highlights the importance of basal anchorage and the anisotropic moduli of the microtubules. As theoretical modeling and numerical simulations progress, along with improved state-of-the-art experiments on primary cilia, we hope that details of ciliary regulated mechano-chemical signaling dynamics in cellular physiology will be understood in the near future.


Subject(s)
Cilia , Microtubules , Biophysics , Mitosis , Models, Theoretical
12.
ChemMedChem ; 16(11): 1740-1743, 2021 06 07.
Article in English | MEDLINE | ID: mdl-33522135

ABSTRACT

A strategy for creating potent and pan-genotypic stimulator of interferon genes (STING) agonists is described. Locking a bioactive U-shaped conformation of cyclic dinucleotides by introducing a transannular macrocyclic bridge between the nucleic acid bases leads to a topologically novel macrocycle-bridged STING agonist (MBSA). In addition to substantially enhanced potency, the newly designed MBSAs, exemplified by clinical candidate E7766, exhibit broad pan-genotypic activity in all major human STING variants. E7766 is shown to have potent antitumor activity with long lasting immune memory response in a mouse liver metastatic tumor model. Two complementary stereoselective synthetic routes to E7766 are also described.


Subject(s)
Antineoplastic Agents/pharmacology , Interferons/agonists , Macrocyclic Compounds/pharmacology , Animals , Antineoplastic Agents/chemical synthesis , Antineoplastic Agents/chemistry , Cell Line, Tumor , Cell Proliferation/drug effects , Drug Screening Assays, Antitumor , Humans , Macrocyclic Compounds/chemical synthesis , Macrocyclic Compounds/chemistry , Mice , Models, Molecular , Molecular Structure , Neoplasms, Experimental/drug therapy , Neoplasms, Experimental/pathology
13.
Biomech Model Mechanobiol ; 19(2): 445-460, 2020 Apr.
Article in English | MEDLINE | ID: mdl-31501964

ABSTRACT

The fluctuating position of an optically trapped cilium tip under untreated and Taxol-treated conditions was used to characterize mechanical properties of the cilium axoneme and its basal body by combining experimental, analytical, and computational tools. We provide, for the first time, evidence that the persistence length of a ciliary axoneme is length-dependent; longer cilia are stiffer than shorter cilia. We demonstrate that this apparent length dependence can be understood by a combination of modeling axonemal microtubules as anisotropic elastic shells and including actomyosin-driven stochastic basal body motion. Our results also demonstrate the possibility of using observable ciliary dynamics to probe interior cytoskeletal dynamics. It is hoped that our improved characterization of cilia will result in deeper understanding of the biological function of cellular flow sensing by this organelle.


Subject(s)
Cilia/metabolism , Animals , Dogs , Elastic Modulus , Image Processing, Computer-Assisted , Madin Darby Canine Kidney Cells , Microspheres , Models, Biological , Optical Tweezers , Paclitaxel/pharmacology , Stochastic Processes
14.
Phys Rev E ; 99(6-1): 063104, 2019 Jun.
Article in English | MEDLINE | ID: mdl-31330602

ABSTRACT

In this work we quantify the effects of surfactant transport on the deformation of a viscous drop under a DC electric field. We study how convective and diffusive transport of surfactants at drop surfaces influence the equilibrium and dynamic deformation of a leaky dielectric drop and a conducting drop. Focusing on the prolate drop shape (elongates along the electric field), we show the differences in equilibrium deformation and flow circulation between a leaky dielectric drop and a conducting drop. We quantify the drop electrodeformation via its dependence on the interior flow circulation and the dominant surfactant transport regime (characterized by the surface Péclet number Pe_{s}). For a leaky dielectric drop with dominant surfactant diffusion (Pe_{s}≪1), equator-to-pole (pole-to-equator) circulation yields smaller (larger) equilibrium deformation with increasing surfactant coverage, compared to a clean drop. However, when convection dominates (Pe_{s}≫1), the equilibrium drop deformation increases (decreases) with larger surfactant coverage for equator-to-pole (pole-to-equator) circulation. Larger equilibrium drop deformation is found for a leaky dielectric drop than a conducting drop when the interior flow is from equator to pole. For an interior flow from pole to equator, we identify cases where larger deformation is found for a conducting interior fluid. Finally, we study the effect of the surfactant transport on the dynamic evolution of drop shape. We found the drop undergoes an overshoot in the early deformation phase, before settling to its equilibrium shape-similar to the overshoot observed for unsteady Stokes flow.

15.
Clin Otolaryngol ; 43(4): 1171-1177, 2018 Aug.
Article in English | MEDLINE | ID: mdl-29718589
17.
Phys Rev E ; 93(5): 053114, 2016 May.
Article in English | MEDLINE | ID: mdl-27300985

ABSTRACT

Most of the existing numerical and theoretical investigations on the electrohydrodynamics of a viscous drop have focused on the creeping Stokes flow regime, where nonlinear inertia effects are neglected. In this work we study the inertia effects on the electrodeformation of a viscous drop under a DC electric field using a novel second-order immersed interface method. The inertia effects are quantified by the Ohnesorge number Oh, and the electric field is characterized by an electric capillary number Ca_{E}. Below the critical Ca_{E}, small to moderate electric field strength gives rise to steady equilibrium drop shapes. We found that, at a fixed Ca_{E}, inertia effects induce larger deformation for an oblate drop than a prolate drop, consistent with previous results in the literature. Moreover, our simulations results indicate that inertia effects on the equilibrium drop deformation are dictated by the direction of normal electric stress on the drop interface: Larger drop deformation is found when the normal electric stress points outward, and smaller drop deformation is found otherwise. To our knowledge, such inertia effects on the equilibrium drop deformation has not been reported in the literature. Above the critical Ca_{E}, no steady equilibrium drop deformation can be found, and often the drop breaks up into a number of daughter droplets. In particular, our Navier-Stokes simulations show that, for the parameters we use, (1) daughter droplets are larger in the presence of inertia, (2) the drop deformation evolves more rapidly compared to creeping flow, and (3) complex distribution of electric stresses for drops with inertia effects. Our results suggest that normal electric pressure may be a useful tool in predicting drop pinch-off in oblate deformations.

18.
Biol Open ; 4(12): 1733-8, 2015 Nov 24.
Article in English | MEDLINE | ID: mdl-26603473

ABSTRACT

Mechanosensation is crucial for cells to sense and respond to mechanical signals within their local environment. While adaptation allows a sensor to be conditioned by stimuli within the environment and enables its operation in a wide range of stimuli intensities, the mechanisms behind adaptation remain controversial in even the most extensively studied mechanosensor, bacterial mechanosensitive channels. Primary cilia are ubiquitous sensory organelles. They have emerged as mechanosensors across diverse tissues, including kidney, liver and the embryonic node, and deflect with mechanical stimuli. Here, we show that both mechanical and chemical stimuli can alter cilium stiffness. We found that exposure to flow stiffens the cilium, which deflects less in response to subsequent exposures to flow. We also found that through a process involving acetylation, the cell can biochemically regulate cilium stiffness. Finally, we show that this altered stiffness directly affects the responsiveness of the cell to mechanical signals. These results demonstrate a potential mechanism through which the cell can regulate its mechanosensing apparatus.

19.
Soft Matter ; 11(37): 7316-27, 2015 Oct 07.
Article in English | MEDLINE | ID: mdl-26264420

ABSTRACT

In this article we report on a study of the near-wall dynamics of suspended colloidal hard spheres over a broad range of volume fractions. We present a thorough comparison of experimental data with predictions based on a virial approximation and simulation results. We find that the virial approach describes the experimental data reasonably well up to a volume fraction of ϕ≈ 0.25 which provides us with a fast and non-costly tool for the analysis and prediction of evanescent wave DLS data. Based on this we propose a new method to assess the near-wall self-diffusion at elevated density. Here, we qualitatively confirm earlier results [Michailidou, et al., Phys. Rev. Lett., 2009, 102, 068302], which indicate that many-particle hydrodynamic interactions are diminished by the presence of the wall at increasing volume fractions as compared to bulk dynamics. Beyond this finding we show that this diminishment is different for the particle motion normal and parallel to the wall.


Subject(s)
Hydrodynamics , Models, Theoretical , Suspensions/chemistry
20.
Article in English | MEDLINE | ID: mdl-26172793

ABSTRACT

The coarse-grained molecular dynamics (MD) or Brownian dynamics (BD) simulation is a particle-based approach that has been applied to a wide range of biological problems that involve interactions with surrounding fluid molecules or the so-called hydrodynamic interactions (HIs). In this paper, an efficient algorithm is proposed to simulate the motion of a single DNA molecule in linear flows. The algorithm utilizes the integrating factor to cope with the effect of the linear flow of the surrounding fluid and applies the Metropolis method (MM) by Bou-Rabee, Donev, and Vanden-Eijnden [Multiscale Model. Simul. 12, 781 (2014)] to achieve more efficient BD simulation. Thus our method permits much larger time step size than previous methods while still maintaining the stability of the BD simulation, which is advantageous for long-time BD simulation. Our numerical results on λ-DNA agree very well with both experimental data and previous simulation results. Finally, when combined with fast algorithms such as the fast multipole method which has nearly optimal complexity in the total number of beads, the resulting method is parallelizable, scalable to large systems, and stable for large time step size, thus making the long-time large-scale BD simulation within practical reach. This will be useful for the study of membranes, long-chain molecules, and a large collection of molecules in the fluids.


Subject(s)
DNA/metabolism , Hydrodynamics , Molecular Dynamics Simulation , Algorithms , DNA/chemistry , Nucleic Acid Conformation
SELECTION OF CITATIONS
SEARCH DETAIL
...