Topology-driven spatial organization of ring polymers under confinement.

Entropic repulsion between DNA ring polymers under confinement is a key mechanism governing the spatial segregation of bacterial DNA before cell division. Here we establish that "internal" loops within a modified-ring polymer architecture enhance entropic repulsion between two overlapping polymers confined in a cylinder. Interestingly, they also induce entropy-driven spatial organization of polymer segments as seen in vivo. Here we design polymers of different architectures in our simulations by introducing a minimal number of cross-links between specific monomers along the ring-polymer contour. The cross-links are likely induced by various bridging proteins inside living cells. We investigate the segregation of two polymers with modified topologies confined in a cylinder, which initially had spatially overlapping configurations. This helps us to identify the architectures that lead to higher success rates of segregation. We also establish the mechanism that leads to localization of specific polymer segments. We use the blob model to provide a theoretical understanding of why certain architectures lead to enhanced entropic repulsive forces between the polymers. Lastly, we establish a correspondence between the organizational patterns of the chromosome of the C.crescentus bacterium and our results for a specifically designed polymer architecture. However, the principles outlined here pertaining to the organization of polymeric segments are applicable to both synthetic and biological polymers.

Polymer architecture orchestrates the segregation and spatial organization of replicating E. coli chromosomes in slow growth.

The mechanism of chromosome segregation and organization in the bacterial cell cycle of E. coli is one of the least understood aspects in its life cycle. The E. coli chromosome is often modelled as a bead spring ring polymer. We introduce cross-links in the DNA-ring polymer, resulting in the formation of loops within each replicating bacterial chromosome. We use simulations to show that the chosen polymer-topology ensures its self-organization along the cell long-axis, such that various chromosomal loci get spatially localized as seen in vivo. The localization of loci arises due to entropic repulsion between polymer loops within each daughter DNA confined in a cylinder. The cellular addresses of the loci in our model are in fair agreement with those seen in experiments as given in J. A. Cass et al., Biophys. J., 2016, 110, 2597-2609. We also show that the adoption of such modified polymer architectures by the daughter DNAs leads to an enhanced propensity of their spatial segregation. Secondly, we match other experimentally reported results, including observation of the cohesion time and the ter-transition. Additionally, the contact map generated from our simulations reproduces the macro-domain like organization as seen in the experimentally obtained Hi-C map. Lastly, we have also proposed a plausible reconciliation of the 'Train Track' and the 'Replication Factory' models which provide conflicting descriptions of the spatial organization of the replication forks. Thus, we reconcile observations from complementary experimental techniques probing bacterial chromosome organization.

Rigidity Dictates Spontaneous Helix Formation of Thermoresponsive Colloidal Chains in Poor Solvent.

The formation of helical motifs typically requires specific directional interactions. Here, we demonstrate that isotropic interparticle attraction can drive self-assembly of colloidal chains into thermo-reversible helices, for chains with a critical level of backbone rigidity. We prepare thermoresponsive colloidal chains by cross-linking PNIPAM microgel-coated polystyrene colloids ("monomers"), aligned in an AC electric field. We control the chain rigidity by varying cross-linking time. Above the LCST of PNIPAM, there is an effective attraction between monomers so that the colloidal chains are in a bad solvent. On heating, the chains decrease in size. For the most rigid chains, the decrease is modest and is not accompanied by a change in shape. Much less rigid chains form relatively compact structures, resulting in a large increase in the local monomer density. Unusually, chains with intermediate rigidity spontaneously assemble into helical structures. The chain helicity increases with temperature and plateaus above the collapse transition temperature of the microgel particles. We simulate a minimal model that captures the spontaneous emergence of the helical conformations of the polymeric chain and provides insight into this shape transition. Our work suggests that a purely mechanical instability for semiflexible filaments can drive helix formation, without the need to invoke directional interactions.

Transient helix formation in charged semiflexible polymers without confinement effects.

Switching on generic interactions e.g. the Coulomb potential or other long ranged spherically symmetric repulsive interactions between monomers of bead-spring model of a semi-flexible polymer induce instabilities in a semiflexible polymer chain to form transient helical structures. Our proposed mechanism could explain the spontaneous emergence of helical order in stiff (bio-) polymers as a chain gets charged from a neutral state. But since the obtained helical structures dissolve away with time, hydrogen bonding (or other additional mechanisms), would be required to form stabilized helical structures as observed in nature (such as in biological macro-molecules). The emergence of the helix is independent of the molecular details of the monomer constituent. The key factors which control the emergence of the helical structure is the persistence length and the charge density. We have avoided using torsional potentials to obtain the transient helical structures. Moreover, we can drive the semiflexible polymer to form helices in a recurring manner by periodically increasing and decreasing the effective charge of the monomers. If the two polymer ends are tethered to two surfaces separated by a distance equal to the contour length of the polymeric chain, which could be in the range 10 nm-µ, the life time of the helical structures formed is increased.

Differential tissue stiffness of body column facilitates locomotion of Hydra on solid substrates.

The bell-shaped members of the Cnidaria typically move around by swimming, whereas the Hydra polyp can perform locomotion on solid substrates in an aquatic environment. To address the biomechanics of locomotion on rigid substrates, we studied the 'somersaulting' locomotion in Hydra We applied atomic force microscopy to measure the local mechanical properties of Hydra's body column and identified the existence of differential Young's modulus between the shoulder region versus rest of the body column at 3:1 ratio. We show that somersaulting primarily depends on differential tissue stiffness of the body column and is explained by computational models that accurately recapitulate the mechanics involved in this process. We demonstrate that perturbation of the observed stiffness variation in the body column by modulating the extracellular matrix polymerization impairs the 'somersault' movement. These results provide a mechanistic basis for the evolutionary significance of differential extracellular matrix properties and tissue stiffness.

Hierarchical and synergistic self-assembly in composites of model wormlike micellar-polymers and nanoparticles results in nanostructures with diverse morphologies.

Using Monte Carlo simulations, we investigate the self-assembly of model nanoparticles inside a matrix of model equilibrium polymers (or matrix of wormlike micelles) as a function of the polymeric matrix density and the excluded volume parameter between polymers and nanoparticles. In this paper, we show morphological transitions in the system architecture via synergistic self-assembly of nanoparticles and the equilibrium polymers. In a synergistic self-assembly, the resulting morphology of the system is a result of the interaction between the nanoparticles and the polymers and corresponding re-organization of both the assemblies. This is different from the polymer templating method. We report the morphological transition of nanoparticle aggregates from percolating network-like structures to non-percolating clusters as a result of the change in the excluded volume parameter between nanoparticles and polymeric chains. Corresponding to the change in the self-assembled structures of nanoparticles, the matrix of equilibrium polymers also simultaneously shows a transition from a dispersed state to a percolating network-like structure formed by the clusters of polymeric chains. We show that the shape anisotropy of the nanoparticle clusters formed is governed by the polymeric density resulting in rod-like, sheet-like or other anisotropic nanoclusters. It is also shown that the pore shape and the pore size of the porous network of nanoparticles can be changed by changing the minimum approaching distance between nanoparticles and polymers. We provide a theoretical understanding of why various nanostructures with very different morphologies are obtained.

Bacterial chromosome organization. I. Crucial role of release of topological constraints and molecular crowders.

We showed in our previous studies that just 3% cross-links (CLs), at special points along the contour of the bacterial DNA, help the DNA-polymer to get organized at micron length scales [T. Agarwal et al., J. Phys.: Condens. Matter 30, 034003 (2018) and T. Agarwal et al., EPL (Europhys. Lett.) 121, 18004 (2018)]. In this work, we investigate how does the release of topological constraints help in the "organization" of the DNA-polymer. Furthermore, we show that the chain compaction induced by the crowded environment in the bacterial cytoplasm contributes to the organization of the DNA-polymer. We model the DNA chain as a flexible bead-spring ring polymer, where each bead represents 1000 base pairs. The specific positions of the CLs have been taken from the experimental contact maps of the bacteria Caulobacter crescentus and Escherichia coli. We introduce different extents of ease of release of topological constraints in our model by systematically changing the diameter of the monomer bead. It varies from the value where chain crossing can occur freely to the value where chain crossing is disallowed. We also study the role of compaction of the chain due to molecular crowders by introducing an "effective" weak Lennard-Jones attraction between the monomers. Using Monte Carlo simulations, we show that the release of topological constraints and the crowding environment play a crucial role to obtain a unique organization of the polymer.

Bacterial chromosome organization. II. Few special cross-links, cell confinement, and molecular crowders play the pivotal roles.

Using a coarse-grained bead-spring model of bacterial chromosomes of Caulobacter crescentus and Escherichia coli, we show that just 33 and 38 effective cross-links in 4017 and 4642 monomer chains at special positions along the chain contour can lead to the large-scale organization of the DNA polymer, where confinement effects of the cell walls play a key role in the organization. The positions of the 33/38 cross-links along the chain contour are chosen from the Hi-C contact map of bacteria C. crescentus and E. coli. We represent 1000 base pairs as a coarse-grained monomer in our bead-spring flexible ring polymer model of the DNA polymer. Thus, 4017/4642 beads on a flexible ring polymer represent the C. crescentus/E. coli DNA polymer with 4017/4642 kilo-base pairs. Choosing suitable parameters from Paper I, we also incorporate the role of compaction of the polymer coil due to the presence of molecular crowders and the ability of the chain to release topological constraints. We validate our prediction of the organization of the bacterial chromosomes with available experimental data and also give a prediction of the approximate positions of different segments within the cell. In the absence of confinement, the minimal number of effective cross-links required to organize the DNA chains of 4017/4642 monomers was 60/82 [Agarwal et al., Europhys. Lett. 121, 18004 (2018) and Agarwal et al., J. Phys.: Condens. Matter 30, 034003 (2018)].

Self assembled linear polymeric chains with tuneable semiflexibility using isotropic interactions.

We propose a two-body spherically symmetric (isotropic) potential such that particles interacting by the potential self-assemble into linear semiflexible polymeric chains without branching. By suitable control of the potential parameters, we can control the persistence length of the polymer and can even introduce a controlled number of branches. Thus we show how to achieve effective directional interactions starting from spherically symmetric potentials. The self-assembled polymers have an exponential distribution of chain lengths akin to what is observed for worm-like micellar systems. On increasing particle density, the polymeric chains self-organize to an ordered line-hexagonal phase where every chain is surrounded by six parallel chains, the transition is first order. On further increase in monomer density, the order is destroyed and we get a branched gel-like phase. This potential can be used to model semi-flexible equilibrium polymers with tunable semiflexibility and excluded volume. The use of the potential is computationally cheap and hence can be used to simulate and probe equilibrium polymer dynamics with long chains. The potential also gives a plausible method of tuning colloidal interactions in experiments such that one can obtain self-assembling polymeric chains made up of colloids and probe polymer dynamics using an optical microscope. Furthermore, we show how a modified potential leads to the observation of an intermediate nematic phase of self-assembled chains in between the low density disordered phase and the line-ordered hexagonal phase.

Role of special cross-links in structure formation of bacterial DNA polymer.

Using data from contact maps of the DNA-polymer of Escherichia coli (E. Coli) (at kilobase pair resolution) as an input to our model, we introduce cross-links between monomers in a bead-spring model of a ring polymer at very specific points along the chain. Via suitable Monte Carlo simulations, we show that the presence of these cross-links leads to a particular organization of the chain at large (micron) length scales of the DNA. We also investigate the structure of a ring polymer with an equal number of cross-links at random positions along the chain. We find that though the polymer does get organized at the large length scales, the nature of the organization is quite different from the organization observed with cross-links at specific biologically determined positions. We used the contact map of E. Coli bacteria which has around 4.6 million base pairs in a single circular chromosome. In our coarse-grained flexible ring polymer model, we used 4642 monomer beads and observed that around 80 cross-links are enough to induce the large-scale organization of the molecule accounting for statistical fluctuations caused by thermal energy. The length of a DNA chain even of a simple bacterial cell such as E. Coli is much longer than typical proteins, hence we avoided methods used to tackle protein folding problems. We define new suitable quantities to identify the large scale structure of a polymer chain with a few cross-links.

Network formation and gelation in telechelic star polymers.

We investigate the efficiency of gelation and network formation in telechelic star polymer melt, where the tips of polymer arms are dipoles while the rest of the monomers are uncharged. Our work is motivated by the experimental observations [A. Kulkarni et al., Macromolecules 48, 6580 (2015)] in which rheological studies of telechelic star polymers of poly-(L-lactide), a bio-degradable polymer, showed a drastic increase in elastic properties (up to 2000 times) compared to corresponding star polymers without the telechelic arm ends. In contrast to previous studies, we avoid using effective attractive Lennard-Jones potentials or dipolar potentials to model telechelic interactions. Instead we use explicit Coulomb positive and negative charges at the tip of polymer-arms of our bead-spring model of star polymers. By our simulations we show that the dipoles at the tip of star arms aggregate together to form clusters of dipoles. Each cluster has contributions from several stars, and in turn each star contributes to several clusters. Thus the entire polymer melt forms a connected network. Network forming tendencies decrease with a decrease of the value of the effective charge constituting the dipole: this can be experimentally realized by choosing a different ionomer for the star tip. We systematically varied the value of dipole charges, the fraction of star-arms with dipoles at the tip, and the length of the arms. The choice of explicit charges in our calculations enables us to make better quantitative predictions about the onset of gelation; moreover we get qualitatively distinct results about structural organization of dipoles within a dipole-cluster.

Hierarchical self-assembly: Self-organized nanostructures in a nematically ordered matrix of self-assembled polymeric chains.

We report many different nanostructures which are formed when model nanoparticles of different sizes (diameter σn) are allowed to aggregate in a background matrix of semiflexible self-assembled polymeric wormlike micellar chains. The different nanostructures are formed by the dynamical arrest of phase-separating mixtures of micellar monomers and nanoparticles. The different morphologies obtained are the result of an interplay of the available free volume, the elastic energy of deformation of polymers, the density (chemical potential) of the nanoparticles in the polymer matrix, and, of course, the ratio of the size of self-assembling nanoparticles and self-avoidance diameter of polymeric chains. We have used a hybrid semi-grand-canonical Monte Carlo simulation scheme to obtain the (nonequilibrium) phase diagram of the self-assembled nanostructures. We observe rodlike structures of nanoparticles which get self-assembled in the gaps between the nematically ordered chains, as well as percolating gel-like network of conjoined nanotubes. We also find a totally unexpected interlocked crystalline phase of nanoparticles and monomers, in which each crystal plane of nanoparticles is separated by planes of perfectly organized polymer chains. We identified the condition which leads to such interlocked crystal structure. We suggest experimental possibilities of how the results presented in this paper could be used to obtain different nanostructures in the laboratory.

Self-assembly of monodisperse CdS nanocylinders with an axial pore.

We present our conclusions of the investigation of the self-assembly and growth of an array of CdS nanotubes: a consequence of a fine balance of directed motion, diffusion, and aggregation of reacting Cd(+2) and S(-2) ions. In a previous communication [Kiruthiga and Chatterji, J. Chem. Phys. 138, 024905 (2013)], we identified the mechanism of an unexpected growth of a CdS nanocylinder of uniform radial cross section from the end of a nanochannel. Furthermore, the cylinder had a pore along the axis but was closed at one end. This unique phenomenon of self-assembly of very monodisperse CdS nanocylinders had been observed in a rather simple experiment where two chambers containing 0.1M CdCl(2) and 0.1M Na(2)S solutions were joined by an array of anodized aluminium oxide (AAO) nanochannels [Varghese and Datta, Phys. Rev. E 85, 056104 (2012)]. Interestingly, the growth of CdS nanotubes was observed only in the Na(2)S chamber. The primary focus of our previous study was on identifying the principles governing the growth of a single nanotube at the exit point of a single AAO nanochannel. In this communication, we identify factors affecting the self-assembly of a nanotube in the presence of other similar neighboring nanotubes growing out from an array of closely spaced AAO nanochannel exit points, a scenario closer to the experimental situation. Our model is not Cd(+2) or S(-2) specific, thus our conclusions suggest that the experimental scheme can be extended to the self-assembly of a general class of reacting-diffusing A and B ions with A (in this case Cd(+2) ions) selectively migrating out from a nanochannel. In particular, we note that after the initial growth of nanotubes for a period of time, there can arise a severe deficiency of B ions (S(-2)) near the AAO-nanochannel exits. The low concentration of B near the nanochannel exits impedes further growth of uniform CdS nanotubes. We further identify the parameters which can be tuned to obtain an improved crop of monodisperse nanotubes. Thereby we predict the necessary characteristics of reacting systems which can be self-assembled using suitable adaptations of experiments used to grow CdS cylinders.

Modeling steady-state dynamics of macromolecules in exponential-stretching flow using multiscale molecular-dynamics-multiparticle-collision simulations.

We introduce a method to obtain steady-state uniaxial exponential-stretching flow of a fluid (akin to extensional flow) in the incompressible limit, which enables us to study the response of suspended macromolecules to the flow by computer simulations. The flow field in this flow is defined by v(x) = εx, where v(x) is the velocity of the fluid and ε is the stretch flow gradient. To eliminate the effect of confining boundaries, we produce the flow in a channel of uniform square cross section with periodic boundary conditions in directions perpendicular to the flow, but simultaneously maintain uniform density of fluid along the length of the tube. In experiments a perfect elongational flow is obtained only along the axis of symmetry in a four-roll geometry or a filament-stretching rheometer. We can reproduce flow conditions very similar to extensional flow near the axis of symmetry by exponential-stretching flow; we do this by adding the right amounts of fluid along the length of the flow in our simulations. The fluid particles added along the length of the tube are the same fluid particles which exit the channel due to the flow; thus mass conservation is maintained in our model by default. We also suggest a scheme for possible realization of exponential-stretching flow in experiments. To establish our method as a useful tool to study various soft matter systems in extensional flow, we embed (i) spherical colloids with excluded volume interactions (modeled by the Weeks-Chandler potential) as well as (ii) a bead-spring model of star polymers in the fluid to study their responses to the exponential-stretched flow and show that the responses of macromolecules in the two flows are very similar. We demonstrate that the variation of number density of the suspended colloids along the direction of flow is in tune with our expectations. We also conclude from our study of the deformation of star polymers with different numbers of arms f that the critical flow gradient ε(c) at which the star undergoes the coil-to-stretch transition is independent of f for f = 2,5,10, and 20.

Nanochannel induced advection with reaction-diffusion dynamics of reacting ions leads to self-assembly of ordered cylindrical structures.

We present a theoretical model, which elucidates the physical principles involved in the formation of very uniform CdS nanocylinders of different radii by combining the physics of flow, diffusion, self assembly, and aggregation of constituent particles. Very recent experiments report that when 0.1M solutions of CdCl(2) and Na(2)S were allowed to mix through some anodised aluminium oxide (AAO) nanochannels, one observes the growth of an array of CdS nano-cylinders on only one end of the AAO template [A. Varghese and S. Datta, Phys. Rev. E 85, 056104 (2012)]. These cylinders have a pore along the center of the cylinder but closed at one end. The reaction happens only in the Na(2)S chamber, and growth of cylinders of uniform size and shape continues as long as the supply of the reactant molecules (CdCl(2) and Na(2)S) is maintained. To try to understand the physics of the observed phenomenon, we propose a model where the Cd(+2) ions exit the AAO-nanochannel to enter Na(2)S chamber with a finite velocity; these ions then react with the diffusing S(-2) ions to form CdS, which then self assemble to form cylinders of uniform width and cross-section. The flow of Cd(+2) out of the AAO nano-channel is the key symmetry breaking feature, which facilitates the formation of uniform cylindrical structures of CdS instead of a CdS precipitate. Since our model does not crucially depend on the chemical details of the reaction, this mechanism can be extended to self-assemble other structures of relevance.

Conformational and dynamical properties of ultra-soft colloids in semi-dilute solutions under shear flow.

We investigate structural and dynamical properties of ultra-soft colloids in dilute and semi-dilute solutions by hybrid mesoscale simulations under linear shear flow. In particular, the influence of functionality on these properties is addressed. Our study combines molecular dynamics simulations for the solute with the multiparticle collision dynamics approach for the coarse-grained solvent. The star polymers exhibit large conformational and orientational changes in shear flow, which we characterize by the radius of gyration tensor and the alignment angle. These quantities show a universal dependence on a concentration- and functionality-dependent Weissenberg number with slight deviations at high shear rates. Moreover, the star polymers display a rotational dynamics with a shear-rate- and concentration-dependent rotation frequency. We attribute the concentration dependence to the screening of hydrodynamic interactions in semi-dilute star-polymer solutions.

The role of effective charges in the electrophoresis of highly charged colloids.

We study the variation of electrophoretic mobility µ of highly charged spherical colloidal macroions for varying surface charge density σ on the colloid using computer simulations of the primitive model for charged colloids. Hydrodynamic interactions between ions are incorporated by coupling the primitive model of charged colloids to the lattice Boltzmann model (LB) of the fluid. In the highly charged regime, the mobility µ of the colloid is known to decrease with the increase of bare charge Q of the colloid; the aim of this paper is to investigate the cause of this. We have identified that the two main factors contributing to the decrease of µ are counterion charge condensation on the highly charged colloid and an increase in effective friction of the macroion-counterion complex due to the condensed counterions. Thus the established O'Brien and White theory, which identified the dipolar force originating from distortion of the electric double layer as the cause of decreasing µ, seems to break down for the case of highly charged colloids with σ in the range of 30-400 µC cm (- 2). To arrive at our conclusions, we counted the number of counterions q0 moving along with the spherical macroion. We observe in our simulations that q0 increases with the increase of bare charge Q, such that the effective charge Qeff = Q - q0 remains approximately constant. Interestingly for our nanometer-sized charged colloid, we observe that, if surface charge density σ of the colloid is increased by decreasing the radius RM of the colloid but fixed bare charge Q, the effective charge Q - q0 decreases with the increase of σ. This behavior is qualitatively different when σ is increased by increasing Q keeping RM fixed. Our observations address a controversy about the effective charge of a strongly charged macroion: some studies claim that effective charge is independent of the bare charge (Alexander et al 1984 J. Chem. Phys. 80 5776; Trizac et al 2003 Langmuir 19 4027) whereas others claim that Qeff decreases with increasing bare charge Q (dos Santos 2009 J. Chem. Phys. 130 124110; Diehl et al 2004 J. Chem. Phys. 121 12100; Groot et al 1991 J. Chem. Phys. 95 9191) at relatively high values of σ.

Electrophoretic properties of highly charged colloids: a hybrid molecular dynamics/lattice Boltzmann simulation study.

Using computer simulations, the electrophoretic motion of a positively charged colloid (macroion) in an electrolyte solution is studied in the framework of the primitive model. In this model, the electrolyte is considered as a system of negatively and positively charged microions (counterions and coions, respectively) that are immersed into a structureless medium. Hydrodynamic interactions are fully taken into account by applying a hybrid simulation scheme, where the charged ions (i.e., macroion and electrolyte), propagated via molecular dynamics, are coupled to a lattice Boltzmann (LB) fluid. In a recent electrophoretic experiment by Martin-Molina et al. [J. Phys. Chem. B 106, 6881 (2002)], it was shown that, for multivalent salt ions, the mobility mu initially increases with charge density sigma, reaches a maximum, and then decreases with further increase of sigma. The aim of the present work is to elucidate the behavior of mu at high values of sigma. Even for the case of monovalent microions, a decrease of mu with sigma is found. A dynamic Stern layer is defined that includes all the counterions that move with the macroion while subjected to an external electrical field. The number of counterions in the Stern layer, q(0), is a crucial parameter for the behavior of mu at high values of sigma. In this case, the mobility mu depends primarily on the ratio q(0)/Q (with Q the valency of the macroion). The previous contention that the increase in the distortion of the electric double layer (EDL) with increasing sigma leads to the lowering of mu does not hold for high sigma. In fact, it is shown that the deformation of the EDL decreases with the increase of sigma. The role of hydrodynamic interactions is inferred from direct comparisons to Langevin simulations where the coupling to the LB fluid is switched off. Moreover, systems with divalent counterions are considered. In this case, at high values of sigma the phenomenon of charge inversion is found.

Combining molecular dynamics with Lattice Boltzmann: a hybrid method for the simulation of (charged) colloidal systems.

We present a hybrid method for the simulation of colloidal systems that combines molecular dynamics (MD) with the Lattice Boltzmann (LB) scheme. The LB method is used as a model for the solvent in order to take into account the hydrodynamic mass and momentum transport through the solvent. The colloidal particles are propagated via MD and they are coupled to the LB fluid by viscous forces. With respect to the LB fluid, the colloids are represented by uniformly distributed points on a sphere. Each such point [with a velocity V(r) at any off-lattice position r] is interacting with the neighboring eight LB nodes by a frictional force F = xi0(V(r)-u(r)), with xi0 being a friction coefficient and u(r) being the velocity of the fluid at the position r. Thermal fluctuations are introduced in the framework of fluctuating hydrodynamics. This coupling scheme has been proposed recently for polymer systems by Ahlrichs and Dunweg [J. Chem. Phys. 111, 8225 (1999)]. We investigate several properties of a single colloidal particle in a LB fluid, namely, the effective Stokes friction and long-time tails in the autocorrelation functions for the translational and rotational velocity. Moreover, a charged colloidal system is considered consisting of a macroion, counterions, and coions that are coupled to a LB fluid. We study the behavior of the ions in a constant electric field. In particular, an estimate of the effective charge of the macroion is yielded from the number of counterions that move with the macroion in the direction of the electric field.