Next Article in Journal
Effect of the Interface/Interphase on the Water Ingress Properties of Joints with PBT-GF30 and Aluminum Substrates Using Silicone Adhesive
Next Article in Special Issue
Analysing Sources of Error in Total Internal Reflection Microscopy (TIRM) Experiments and Data Analysis
Previous Article in Journal
Designing Multistimuli-Responsive Anisotropic Bilayer Hydrogel Actuators by Integrating LCST Phase Transition and Photochromic Isomerization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Brownian Motion in Optical Tweezers, a Comparison between MD Simulations and Experimental Data in the Ballistic Regime

by
Krzysztof Zembrzycki
1,*,
Sylwia Pawłowska
2,
Filippo Pierini
1,* and
Tomasz Aleksander Kowalewski
1
1
Department of Biosystem and Soft Matter, Institute of Fundamental Technological Research, Polish Academy of Sciences, ul. Pawinskiego 5B, 02-106 Warsaw, Poland
2
Faculty of Electronics, Telecommunications and Informatics, Gdańsk University of Technology, ul. G. Narutowicza 11/12, 80-233 Gdańsk, Poland
*
Authors to whom correspondence should be addressed.
Polymers 2023, 15(3), 787; https://doi.org/10.3390/polym15030787
Submission received: 3 January 2023 / Revised: 27 January 2023 / Accepted: 31 January 2023 / Published: 3 February 2023

Abstract

:
The four most popular water models in molecular dynamics were studied in large-scale simulations of Brownian motion of colloidal particles in optical tweezers and then compared with experimental measurements in the same time scale. We present the most direct comparison of colloidal polystyrene particle diffusion in molecular dynamics simulations and experimental data on the same time scales in the ballistic regime. The four most popular water models, all of which take into account electrostatic interactions, are tested and compared based on yielded results and resources required. Three different conditions were simulated: a freely moving particle and one in a potential force field with two different strengths based on 1 pN/nm and 10 pN/nm. In all cases, the diameter of the colloidal particle was 50 nm. The acquired data were compared with experimental measurements performed using optical tweezers with position capture rates as high as 125 MHz. The experiments were performed in pure water on polystyrene particles with a 1 μm diameter in special microchannel cells.

Graphical Abstract

1. Introduction

The knowledge of interactions of colloidal particles suspended in water with their surrounding medium is crucial to understanding their behavior in crowded environments, such as cells and other biological materials [1,2], and diffusion is one of the most important processes taking place everywhere where there are fluids or gasses involved [3,4,5]. Often, molecular dynamics is used for simulating such conditions. Unfortunately, water molecule behavior is complex, and there are many theoretical models that try to reproduce it [6,7,8,9]. Some are not suitable for simulating the diffusion of particles in water, but many are, and choosing the correct one can be an arduous task. Another critical factor is the time scale of the simulation. Since molecular dynamics is computationally intensive, only relatively small systems or short time scales are actively computed. Unfortunately, the process of diffusion during a very small-time frame is very different from what is generally observed on macro scales; the movement of the particle becomes more ballistic than random, which means that below a characteristic relaxation time τ, the diffusing particle retains some of its momentum. This is called the ballistic regime, and the characteristic time can be calculated using the following Equations (1) and (2) [10]:
τ = m + m f 2 γ
γ = 6 π η R
where:
m—particle mass
mf—mass of the fluid displaced by the particle when it moves
γ—friction coefficient
η—dynamic viscosity of the fluid
R—particle radius
A 1 μm diameter polystyrene colloidal particle has a relaxation time τ of about 93 ns, and during this time, it travels around 1 Å. If we reduce the diameter to 50 nm, the relaxation time goes down to 0.21 ns. Knowing the behavior of colloids in such short time scales has become very important in recent studies [11,12]. Another critical factor is that a biological system, similar to most others, is full of charged particles or particles that are polarized, such as water itself; therefore, electrostatic interactions will play a major role in the system.
Four water molecules were selected for comparison, all of which include electrostatic interactions of charged atoms: SPC (simple point charge) [13,14], SPCE (extended simple point charge) [15,16], TIP3P (transferable intermolecular potential 3 point) [17,18], and TIP4P (transferable intermolecular potential 4 point) [19,20]. Those were chosen because they are designed to be used in conditions close to a normal temperature and pressure or those common in biological systems. All of them can be effectively described as rigid pairs with potential composed of Lenanrd–Jones (LJ) and Coulombic terms, with all but TIP4P having three sites of interaction and the latter having a fourth mass-less point charge. On the other hand, they differ in how the water particle is described, in terms of parameters, and can vary quite substantially in computational complexity, which is reflected by the time required for a given simulation. All of them are used extensively to model the behavior of biological systems [21,22].
Since their discovery in 1970, optical tweezers have emerged as a tool to trap and manipulate nano and micrometer-sized material using a highly focused laser beam [23,24,25]. It is an extremely sensitive and precise instrument capable of manipulating objects and detecting their position with sub-nanometer precision and measuring forces with femtonewton (10–15 N) accuracy [26,27]. Now, they are finding applications in many fields of science, such as biology and chemistry, where they are used extensively in studying the unfolding of proteins [28,29,30]. Thanks to the advancements in the field electronics, it is now possible to further increase the resolution of optical tweezers in the time domain as well and observe the phenomena that happen on even smaller time scales in the range of ns [31,32]. Therefore, this technique is the best choice to measure the motion of colloidal particles with both time and spatial scales that match those achieved in numerical simulations.
An experimental study of the diffusion of polystyrene particles trapped in optical tweezers was conducted. This technique allows very precise measurements with high temporal and spatial resolution, even in the ballistic regime [33,34,35]. This allows a comparison of the results of numerical simulations with experimental results more directly and therefore can verify which water models are best suited for this range of applications.

2. Materials and Methods

In this study, use was made of both numerical simulations and experimental measurements of a spherical polystyrene bead moving in water. For time steps Δt>>τ, this is called Brownian motion. These fluctuations are stochastic in nature, and the general equation for the position as a function time is shown in the following differential Equation (3) [36]:
x ˙ t = 2 k B T γ W t
where:
x—position
t—time
kB—Boltzmann constant
T—temperature
W(t)—white noise
The white noise term is not a standard function, but rather a process characterized by the following properties:
  • its mean <W(t)> is equal to 0 for all values of t
  • Mean Squared Displacement (MSD) <W2(t)> is equal to 1 for all values of t
  • W(ti) and W(tj) are independent for all i ≠ j
If an external harmonic potential is applied to the particle, as generated by an optical trap, with the force if defined by Equation (4),
F x = k x
where:
F—force
k—force constant
then Equation (3) for the particle position takes the following form (5):
x ˙ t = k γ x t + 2 k B T γ W t
The following Equation (6) ties the MSD to the diffusion coefficient, where N is a constant that depends on the degrees of freedom and takes vales of 2, 4, and 6 for 1, 2, and 3 degrees of freedom respectively.
D = M S D N d t
The Stokes–Einstein Equation (7) can be used to calculate the theoretical diffusion of a colloidal particle in a stationary fluid without external force fields.
D = k B T γ
As stated earlier, the above equations are true tor time steps Δt >> τ, where the MSD is proportional to t. In the ballistic regime, the MSD is proportional to t2 [34]. It is important to note however that in the Euclidean space, the positions of the particle in each exes are statistically independent; therefore, it is possible to compare diffusion coefficients obtained by from sources with different dimensions or degrees of freedom. For the purpose of this article, that is to compare MD water models, we will be calculating the diffusion coefficient, keeping in mind that in short time scales and in the presence of the trapping potential, it is not equal to the Stokes–Einstein relation. This is performed because the diffusion coefficient, as described by Equation (6), normalizes the MSD to the time step, making it easier to compare results obtained for different measurements.

2.1. Simulations

The LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator, version 22 August 2018) [37] was used to perform simulations using all of the selected models [38,39]. It is one of the most commonly used simulators in many fields of science. The simulation consisted of a 2-dimensional rectangular box, 350 × 350 nm in size, filled with water, with a colloidal bead in the center. The box boundary in the simulation plane was periodic, meaning that particles interact across the boundary and can exit one end of the box and re-enter the other. The third unused dimension was fixed to 0. All molecules were able to freely rotate within the plane of simulation. Initially, the water molecules were placed on a rectangular grid with 3.1 Å spacing, which is equal to their normal mean distance [40]. The potential describing the molecule interactions was a standard Lennard–Jones given by the following Equation (8):
E r = 4 ϵ σ r 12 σ r 6
where:
E—potential energy
r—distance between atoms
ε—depth of the potential well
σ—distance at which the intermolecular potential is zero
To calculate the LJ interaction parameters between different types of atoms, the Lorentz–Berthelot rules where applied, given by the following Equations (9) and (10)
q i j = q i + q j 2
ϵ i j = ϵ i ϵ j
The Coulombic pairwise interactions were given by the following Equation (11):
E i j r = C q i q j ϵ r
where
Eij—potential energy the between i-type and j-type atom
r—distance between atoms
ε—dielectric constant
q—charge of atom
C—energy conversion constant
The exact parameters of the molecule depend on the model and are shown in Table 1.
In all cases, the cutoff distance, above which the interaction forces were not calculated, was set to 10 Å for LJ and 15 Å for Coulombic interactions. For calculating the electrostatic interactions, the Particle–Particle Particle–Mesh (PPPM) method was set and the default velocity–verlet method was employed for LJ interactions. A round disc was placed in the middle of the simulation box, making the colloid particle. Since polystyrene is often used to make colloidal particles, the disc was made to simulate this material. Its composition was carbon and hydrogen in the exact same proportions as in a polystyrene chain. The atoms were placed in ring patterns, filling the entire disc, and a spacing corresponding with polystyrene, 1.4 Å between carbon atoms and 1.1 Å between hydrogen, and 1.8 Å spacing between carbon rings so that the overall density of the colloid would be the same as water in the simulations, which is 1 g/cm3. A real polystyrene bead used in the experiments has a slightly higher density of 1.005 g/cm3, but this difference is insignificant in the time scales discussed in this paper. The overall diameter of the disc was 50 nm, and it was configured to be one rigid body and interacted with water with the same LJ and Coulombic potential with parameters listed in Table 1. The size of the void space between colloid and water atoms was set so that the overall mass density of the system was also equal to 1 g/cm3, and the overall atom count was almost 4 million, with more than 1.25 million water molecules. To simulate an optical trap, a harmonic-type force was applied to the center carbon atom of the particle. At the beginning of each simulation, the water molecules were assigned velocities corresponding to a temperature of 298.0 K, generated with a random number generator. During the whole run, an NVE thermostat was applied to the whole system

2.2. Experiment

The optical trap [41,42] used in the experiments was developed earlier [43], but for the purpose of these measurements, the electronic front end had to be modified to increase its bandwidth. The detector was a quadrant avalanche photodiode (QA4000, First Sensor AG, Berlin, Germany), and the signals were digitized using a PC acquisition card (OctaveCSE8349, GaGe Vitrek LLC, Lockport, IL, USA), with a 125 MHz maximum sample rate, which equals a 8 ns time between individual samples, which is in the same order of magnitude as the overall simulated time. Experiments were conducted in microchannels made from PDMS, using a photolithography-based technique, with a geometry that consisted of a main channel and a series of round wells in which the measurements made were connected to the main channel by a single smaller one. This solved the problem of particle movement being influenced by stray fluid flow. All measurements were performed above 5 μm from the bottom channel wall in order to ensure that the effect of slip was negligible. Before the experiment, each PDMS microchannel was placed in a vacuum for 15 min and immediately filled with water. That caused any air bubbles in the wells to be removed. Polystyrene particles used in the measurements were 1 μm in diameter (R1000 Thermo Fisher Scientific, Waltham, MA, USA) and were dispersed in distilled water with a 1:10,000 concentration. The optical tweezers were calibrated to the same particle type. A fresh solution was made before every experiment to prevent aggregation and sedimentation.

3. Results and Discussion

The simulations were split into three groups based on the applied harmonic potential simulating the optical tweezers. First, it was set to zero, corresponding to a freely moving particle; next, was a force constant of 1 pN/nm, and lastly, 10 pN/nm, but only for the two most basic models, SPC and SPCE. Figure 1 shows a scaled-down version of the simulation box.
Due to computational complexity and real-time constraints, the TIP4P model was simulated only with a freely moving particle, without simulations that included a harmonic potential, similar to that in the case of other models discussed here. Each group consisted of four simulation runs, identical in parameters, the only difference being the random number generator seed that generated the initial water velocities. The value of the time step was different for each model and condition and was fine-tuned to the highest value achievable for the simulation to be stable. High values resulted in lost atoms, that is, atoms that escaped the simulation box or the temperature of the system exceeding the set value. The velocity of the colloid particle was set to 0 at the start of the run, and its position was printed every 1000 steps, and the time step value varied depending on the water model used. In all cases, at least 10 M steps were computed, which corresponds to about 30 ns of real time in most cases, depending on the actual time step setting. The equilibration phase was very short, and each system reached a steady state after about 10 k time steps. Each simulation was run on an identical computer equipped with two Intel Xeon X5650 processors with six (physical) cores running at 2.67 GHz and 24 GB of RAM. Table 2 contains the summarized model performance and resource usage.
When an external force field was applied, in some cases, the time step value had to be modified for the simulation to be stable. However, in all cases, it did not impact the time step calculation time or memory usage. Figure 2 shows a typical particle movement, one from each model. Due to variations in the time step value, the overall simulated real time varied depending on the model and whether external force field was applied.
From the computational perspective, both SPC and SPCE models behaved very similarly, having the same memory usage, and the same time step value and application of additional force field did not impact their performance. The only major difference is that the SPC model showed results slightly faster. On the other hand, the TIP3P and TIP4P behaved very differently depending on whether the additional force field was applied. In the latter case, its performance was not much worse than that of the SPC(E) models, but when an external force field was applied, the maximum time step value decreased dramatically by an order of magnitude. This significantly increased the need for resources required to compute a given amount of real time and was the primary reason why in this experiment, the TIP4P model was used only to simulate a freely moving colloidal particle.
Figure 3 shows calculated diffusion coefficients using Equation (12), which is a special case of Equation (6):
D = d s 2 4 d t
where:
ds—particle displacement
dt—time of the displacement, 1000 time steps in our case
Most considered models give an order of magnitude of higher diffusion than the theoretical value of 1.45 × 10−11 m2/s; only TIP4P is very close but is the most complicated model. As expected, with force applied, the diffusion drops, and the particle is confined close to the center of the simulated optical tweezers, with all models giving similar results. In the performed simulations, although very similar to those of SPCE, the SPC model had a faster computation time and yielded similar results with exactly the same memory usage. Employing the TIP3P model was more time-consuming and slightly more memory-consuming while giving similar results. It also had a lower time step value; thus, it required more computation time to simulate the same amount of real time. The model that produced the most accurate result was TIP4P, but on the other hand, it was significantly more resource-consuming, taking about 3.3 times the time to compute a time step and requiring twice the memory.
Figure 3. Calculated average diffusion coefficients obtained from simulation data. For the given colloidal particle, the theoretical diffusion coefficient without external force, calculated using the Stokes–Einstein equation, is 9.81 × 10−12 m2/s.
Figure 3. Calculated average diffusion coefficients obtained from simulation data. For the given colloidal particle, the theoretical diffusion coefficient without external force, calculated using the Stokes–Einstein equation, is 9.81 × 10−12 m2/s.
Polymers 15 00787 g003
A series of experimental measurements was made on 1 μm polystyrene beads suspended in water trapped in optical tweezers with a stiffness of 1 pN/nm, the same as in the simulations. Multiple different particles were measured in order to eliminate errors caused by particle size differences. Every measurement was made in a special experimental chamber with only one inlet, the particle being dragged in by the optical tweezers. This technique ensured that there was no fluid movement caused by residue pressure, fluctuations in temperature, and mechanical vibrations that could cause additional forces influencing the particle and thus cause measurement errors. Figure 4 shows a picture of one such chamber.
Each particle was measured twice, first with a low sample rate of 1k samples/s and the second time with a maximum acquisition rate of 125 M samples/s. Care was taken to ensure that only one particle was within the tweezers. While hard to see under the microscope, the diffraction pattern of the laser light back-scattered from the objects within the optical trap is very different that when more than one particle is trapped. Results, calculated using the same formula as with simulations (6), are shown in Figure 5. The theoretical diffusion was calculated using the Stokes–Einstein Equation (13):
D = k B T 6 π η R
where:
T—absolute temperature, 298 °C in our case
η—dynamic viscosity, 0.89 mPa s−1 in our case
R—radius of the particle, 0.5 μm in our case
Figure 5. Measured diffusion of 1 μm polystyrene particle in an optical tweezers system with a 1 pN/nm trap stiffness, along with theoretical calculations for a freely moving particle.
Figure 5. Measured diffusion of 1 μm polystyrene particle in an optical tweezers system with a 1 pN/nm trap stiffness, along with theoretical calculations for a freely moving particle.
Polymers 15 00787 g005
As observed, the measured diffusion with a slow sample rate is smaller than that of a theoretical calculation for a freely moving particle. For comparison, the diffusion measured with a fast sample rate is higher. This is because the time scale is in the ballistic regime, where the mean square displacement has a nonlinear relationship with time. This increase in the measured diffusion coefficient in the ballistic regime is in agreement with the numerical simulations discussed in the previous paragraph. Although the particle used in the experiments is much bigger than that used in the simulations, as can be seen in Equation (7), the diffusion coefficient is independent of mass and is inversely proportional to the radius. Therefore, it is possible to qualitatively compare the obtained results.
The obtained results of simulations can be compared to the work of Pekka and Lennart [44] where the authors compared the MD simulations of water molecules simulated with SPC, SPCE, and TIP3P models. Some of the simulation parameters are also very similar to those used in this work, such as the time step (2 fs), temperature (298 K), and cutoff distance (12 Å). Results obtained by Pekka and Lennart differ based on the model used, almost by an order of magnitude, and all give higher values than those obtained by experiments, as was the case in this study. Yet, there are also major differences between the conducted simulations, mainly the scale; there, the box contained around 1k atoms compared to almost 8M in this work. Moreover, here, the diffusion of a polystyrene particle, with an external force applied, was simulated compared to pure water molecules.
In the work by Pancorbo, Rubio, and Domínguez-García [45] the authors studied Brownian motion in optical tweezers using 2D simulations, very similar to this work. The disk had a diameter of 1.9 μm, at a temperature of 295.9 K, at trap. Their results agree to within 1% of theoretical values, which is a very good result. Yet, there are many differences between their work and this. Firstly, the simulation was not molecular dynamics-based and the time is in the rage of milliseconds, so it is not possible to verify different water models and simulate Brownian dynamics in the ballistic regime. In addition, there is no mention of electrostatic interactions between water and the molecule simulated.
To the best of our knowledge, this is the first attempt at a large-scale MD simulation of Brownian motion of a charged colloidal particle in different water models. Yet, as is always the case, there are some limitations of the employed methodology. Firstly, everything was performed in 2D. It would be even better to perform 3D simulations, which would take into account the non-symmetric potential of the optical trap, which is different along the optical axis of the lens to the potential in the plane perpendicular to the optical axis. Another drawback is that the studied water models are one of the simplest ones. There are other more sophisticated ones, such as TIP5P [46], which may yield more accurate results. One of the most evident differences between the experiment and the simulation is the size difference between the colloids. In the MD, it was 50 nm, whereas in the experiment, it was 1 μm. While this issue is addressed in this work, it would be good to perform simulations and experimental measurements with exactly the same sized colloids. While it is possible to detect nanometer-sized objects with optical tweezers, it has not been performed with high acquisition rates due to the a very small signal available form such a small object.

4. Conclusions

This work compares the behavior of different water models in molecular dynamics simulations and verifies those data against experimental measurements. Four models, SPC, SPCE, TPI3P, and TIP4P, were evaluated in regards to their resource usage and from a qualitative perspective, all yielded satisfying results, with the SPC being the most economical, but from a quantitative point of view, the TIP4P proved to be the most accurate. Both the SPC and SPCE models yielded very similar results and required the same amount of memory, with SPCE being 14% slower on the same machine. Moreover, the addition of an external force field did not impact their overall performance. As expected, the slowest model out of the tested four was the TIP4P, being about 70% slower than SPC and requiring 100% more memory. Further, the TIP3P and TIP4P models required a reduction of the time step value by an order of magnitude when an external force was added, dramatically increasing the amount of time required to simulate a given amount of real time. This was not the case with the SPC and SPCE models, meaning that the use of TIP3P and TIP4P models can cause significant drawbacks in many applications requiring long-term simulations. All tested models had the same behavior with the particle diffusion increased in the ballistic regime.
Experimental measurements, made using optical tweezers with a 1 pN/nm trap stiffness, show an increased diffusion coefficient in the ballistic regime compared to that of the standard one measured with a relatively slow sample rate. This expected behavior corresponds well with numerical data, in which the diffusion in the ballistic regime is orders of magnitude higher in larger timescales. The data obtained in our study can prove crucial, as The Brownian motion of confined particles is often used in calibration procedures of optical traps. Therefore if the measurements were taken with a high enough acquisition speed, the calibration would be incorrect. In recent years a new technique of measuring fluid viscosity through observations of diffusion in the ballistic regime has been reported [11], and the results of this study should help in better predicting the results of simulations of such systems.
Possible future improvements on the subject of this work would be to change the number of dimensions in MD simulations to three, increase the size of the colloid to better match experimental capabilities, and finally employ a more sophisticated water model, such as the mentioned TIP5P. While all this requires significantly more computational power, due to advances in technology, the required computers will become more accessible.

Author Contributions

Conceptualization, K.Z. and T.A.K.; methodology, K.Z. and S.P.; validation, K.Z.; investigation, K.Z., S.P. and F.P.; resources, K.Z. and F.P.; writing—original draft preparation, K.Z.; writing—review and editing, K.Z., S.P. and F.P.; supervision, F.P.; funding acquisition, K.Z. and F.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Polish NCN grant No. 2015/19/N/ST8/03931. This study was partially supported by the LIDER program funded by the National Center for Research and Development (NCBiR) under NCBiR Grant No. LIDER/28/0067/L-7/15/NCBR/2016.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ellis, R.J. Macromolecular crowding: An important but neglected aspect of the intracellular environment. Curr. Opin. Struct. Biol. 2001, 11, 114–119. [Google Scholar] [CrossRef] [PubMed]
  2. Andrews, S.S. Effects of surfaces and macromolecular crowding on bimolecular reaction rates. Phys. Biol. 2020, 17, 045001. [Google Scholar] [CrossRef] [PubMed]
  3. Hall, D.; Hoshino, M. Effects of macromolecular crowding on intracellular diffusion from a single particle perspective. Biophys. Rev. 2010, 2, 39–53. [Google Scholar] [CrossRef] [PubMed]
  4. Minton, A.P. Holobiochemistry: The effect of local environment upon the equilibria and rates of biochemical reactions. Int. J. Biochem. 1990, 10, 1063–1067. [Google Scholar] [CrossRef]
  5. Novak, I.L.; Kraikivski, P.; Slepchenko, B.M. Diffusion in cytoplasm: Effects of excluded volume due to internal membranes and cytoskeletal structures. Biophys. J. 2009, 97, 758–767. [Google Scholar] [CrossRef]
  6. Berezhkovskii, A.M.; Szabo, A. Theory of crowding effects on bimolecular reaction rates. J. Phys. Chem. B 2016, 120, 5998–6002. [Google Scholar] [CrossRef]
  7. Kim, J.S.; Yethiraj. A. Effect of macromolecular crowding on reaction rates: A computational and theoretical study. Biophys. J. 2009, 96, 1330–1334. [Google Scholar] [CrossRef]
  8. Andrews, S.S.; Addy, N.J.; Brent, R.; Arkin, A.P. Detailed simulations of cell biology with Smoldyn 2.1. PLoS Comput. Biol. 2010, 16, e1000705. [Google Scholar] [CrossRef]
  9. Wieczorek, G.; Zielenkiewicz, P. Influence of Macromolecular Crowding on Protein-Protein Association Rates—a Brownian Dynamics Study. Biophys. J. 2008, 95, 5030–5036. [Google Scholar] [CrossRef]
  10. Kheifets, S.; Simha, A.; Melin, K.; Li, T.; Raizen, M.G. Observation of Brownian Motion in Liquids at Short Times: Instantaneous Velocity and Memory Loss. Science 2014, 343, 1493–1496. [Google Scholar] [CrossRef]
  11. Madsen, L.S.; Waleed, M.; Casacio, C.A.; Terrasson, A.; Stilgoe, A.B.; Taylor, M.A.; Bowen, W.P. Ultrafast viscosity measurement with ballistic optical tweezers. Nat. Photon 2021, 15, 386–392. [Google Scholar] [CrossRef]
  12. Rubtsov, I.V.; Burin, A.L. Ballistic and diffusive vibrational energy transport in molecules. J. Chem. Phys. 2019, 150, 020901. [Google Scholar] [CrossRef]
  13. Berendsen, H.J.C.; Postma, J.P.M.; Van Gunsteren, W.F.; Hermans, J. Interaction models for water in relation to protein hydration. In Intermolecular Forces; Pullman, B., Ed.; Springer: Dordrecht, The Netherlands, 1981; Volume 14, pp. 331–342. [Google Scholar] [CrossRef]
  14. Fitzgerald, C.; Hosny, N.A.; Tong, H.; Seville, P.C.; Gallimore, P.J.; Davidson, N.M.; Athanasiadis, A.; Botchway, S.W.; Ward, A.D.; Kalberer, M.; et al. Fluorescence lifetime imaging of optically levitated aerosol: A technique to quantitatively map the viscosity of suspended aerosol particles. Phys. Chem. Chem. Phys. 2016, 18, 21710–21719. [Google Scholar] [CrossRef]
  15. Berendsen, H.J.C.; Grigera, J.R.; Straatsma, T.P. The missing term in effective pair potentials. J. Phys. Chem. 1987, 91, 6269–6271. [Google Scholar] [CrossRef]
  16. Ding, H.; Kollipara, P.S.; Lin, L.; Zheng, Y. Atomistic modeling and rational design of optothermal tweezers for targeted applications. Nano Res. 2021, 14, 295–303. [Google Scholar] [CrossRef]
  17. Price, D.J.; Brooks, C.L. A modified TIP3P water potential for simulation with Ewald summation. J. Chem. Phys. 2004, 121, 10096–10103. [Google Scholar] [CrossRef]
  18. Carney, S.P.; Ma, W.; Whitley, K.D.; Jia, H.; Lohman, T.M.; Luthey-Schulten, Z.; Chemla, Y.R. Kinetic and structural mechanism for DNA unwinding by a non-hexameric helicase. Nat. Commun. 2021, 12, 1–14. [Google Scholar] [CrossRef]
  19. Horn, H.W.; Swope, W.C.; Pitera, J.W.; Madura, J.D.; Dick, T.J.; Hura, G.L.; Head-Gordon, T. Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. J. Chem. Phys. 2004, 120, 9665–9678. [Google Scholar] [CrossRef]
  20. Jakobi, A.J.; Mashaghi, A.; Tans, S.J.; Huizinga, E.G. Calcium modulates force sensing by the von Willebrand factor A2 domain. Nat. Commun. 2011, 2, 385. [Google Scholar] [CrossRef] [Green Version]
  21. Ferrario, V.; Pleiss, J. Simulation of protein diffusion: A sensitive probe of protein–solvent interactions. J. Biomol. Struct. Dyn. 2019, 37, 1534–1544. [Google Scholar] [CrossRef]
  22. Zhang, H.; Yin, C.; Jiang, Y.; van der Spoel, D. Force Field Benchmark of Amino Acids: I. Hydration and Diffusion in Different Water Models. J. Chem. Inf. Model. 2018, 58, 1037–1052. [Google Scholar] [CrossRef] [PubMed]
  23. Ashkin, A. Acceleration and Trapping of Particles by Radiation Pressure. Phys. Rev. Lett. 1970, 24, 156–159. [Google Scholar] [CrossRef]
  24. Ashkin, A.; Dziedzic, J.M.; Bjorkholm, J.E.; Chu, S. Observation of a single-beam gradient force optical trap for dielectric particles. Opt. Lett. 1986, 11, 288–290. [Google Scholar] [CrossRef] [PubMed]
  25. Bolognesi, G.; Friddin, M.S.; Salehi-Reyhani, A.; Barlow, N.E.; Brooks, N.J.; Ces, O.; Elani, Y. Sculpting and fusing biomimetic vesicle networks using optical tweezers. Nat. Commun. 2018, 9, 1882. [Google Scholar] [CrossRef]
  26. Kolbow, J.D.; Lindquist, N.C.; Ertsgaard, C.T.; Yoo, D.; Oh, S. Nano-Optical Tweezers: Methods and Applications for Trapping Single Molecules and Nanoparticles. ChemPhysChem 2021, 22, 1409–1420. [Google Scholar] [CrossRef]
  27. Killian, J.L.; Inman, J.T.; Wang, M.D. High-Performance Image-Based Measurements of Biological Forces and Interactions in a Dual Optical Trap. ACS Nano 2018, 12, 11963–11974. [Google Scholar] [CrossRef]
  28. Murugesapillai, D.; McCauley, M.J.; Maher, L.J.; Williams, M.C. Single-molecule studies of high-mobility group B architectural DNA bending proteins. Biophys. Rev. 2016, 9, 17–40. [Google Scholar] [CrossRef]
  29. Jagannathan, B.; Marqusee, S. Protein folding and unfolding under force. Biopolymers 2013, 99, 860–869. [Google Scholar] [CrossRef]
  30. Kreysing, M.; Ott, D.; Schmidberger, M.J.; Otto, O.; Schürmann, M.; Martín-Badosa, E.; Whyte, G.; Guck, J. Dynamic operation of optical fibres beyond the single-mode regime facilitates the orientation of biological cells. Nat. Commun. 2014, 5, 5481. [Google Scholar] [CrossRef] [Green Version]
  31. Soni, G.V.; Jonsson, M.P.; Dekker, C. Periodic Modulations of Optical Tweezers Near Solid-State Membranes. Small 2012, 9, 679–684. [Google Scholar] [CrossRef]
  32. Rahman, S.; Torun, R.; Zhao, Q.; Boyraz, O. Electronic control of optical tweezers using space-time-wavelength mapping. J. Opt. Soc. Am. B 2016, 33, 313–319. [Google Scholar] [CrossRef]
  33. Huang, R.; Chavez, I.; Taute, K.M.; Lukić, B.; Jeney, S.; Raizen, M.G.; Florin, E.-L. Direct observation of the full transition from ballistic to diffusive Brownian motion in a liquid. Nat. Phys. 2011, 7, 576–580. [Google Scholar] [CrossRef]
  34. Li, T.; Kheifets, S.; Medellin, D.; Raizen, M.G. Measurement of the Instantaneous Velocity of a Brownian Particle. Science 2010, 328, 1673–1675. [Google Scholar] [CrossRef]
  35. Lu, D.; Labrador-Páez, L.; Ortiz-Rivero, E.; Frades, P.; Antoniak, M.A.; Wawrzyńczyk, D.; Nyk, M.; Brites, C.D.S.; Carlos, L.D.; Solé, J.A.G.; et al. Exploring Single-Nanoparticle Dynamics at High Temperature by Optical Tweezers. Nano Lett. 2020, 20, 8024–8031. [Google Scholar] [CrossRef]
  36. Volpe, G.; Volpe, G. Simulation of a Brownian particle in an optical trap. Am. J. Phys. 2013, 81, 224–230. [Google Scholar] [CrossRef]
  37. LAMMPS Molecular Dynamics Simulator. Available online: https://www.lammps.org/ (accessed on 22 September 2022).
  38. Sun, H.; Ferasat, K.; Nowak, P.; Gravelle, L.; Gaffran, N.; Anderson, C.; Sirola, T.; Pintar, O.; Lievers, W.B.; Kim, I.Y.; et al. Implementing a non-local lattice particle method in the open-source large-scale atomic/molecular massively parallel simulator. Model. Simul. Mater. Sci. Eng. 2022, 30, 054001. [Google Scholar] [CrossRef]
  39. Murashima, T.; Urata, S.; Li, S. Coupling finite element method with large scale atomic/molecular massively parallel simulator (LAMMPS) for hierarchical multiscale simulations. Eur. Phys. J. B 2019, 92, 211. [Google Scholar] [CrossRef]
  40. Perkins, S.J. Protein volumes and hydration effects. The calculations of partial specific volumes, neutron scattering matchpoints and 280-nm absorption coefficients for proteins and glycoproteins from amino acid sequences. JBIC J. Biol. Inorg. Chem. 1986, 157, 169–180. [Google Scholar] [CrossRef]
  41. Sato, S.; Inaba, H. Optical trapping and manipulation of microscopic particles and biological cells by laser beams. Opt. Quantum Electron. 1996, 28, 1–16. [Google Scholar] [CrossRef]
  42. Hansen, P.M.; Bhatia, V.K.; Harrit, N.; Oddershede, L. Expanding the Optical Trapping Range of Gold Nanoparticles. Nano Lett. 2005, 5, 1937–1942. [Google Scholar] [CrossRef]
  43. Pierini, F.; Zembrzycki, K.; Nakielski, P.; Pawłowska, S.; A Kowalewski, T. Atomic force microscopy combined with optical tweezers (AFM/OT). Meas. Sci. Technol. 2016, 27, 025904. [Google Scholar] [CrossRef]
  44. Pekka, M.; Lennart, N. Structure and Dynamics of the TIP3P, SPC, and SPC/E Water Models at 298 K. J. Phys. Chem. A 2001, 105, 9954–9960. [Google Scholar]
  45. Pancorbo, M.; Rubio, M.A.; Domínguez-García, P. Brownian dynamics simulations to explore experimental microsphere diffusion with optical tweezers. Procedia Comput. Sci. 2017, 108, 166–174. [Google Scholar] [CrossRef]
  46. Mahoney, M.W.; Jorgensen, W.L. A five-site model for liquid water and the reproduction of the density anomaly by rigid, nonpolarizable potential functions. J. Chem. Phys. 2000, 112, 8910. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Scaled-down (computed simulation consisted of almost 4 million atoms) image of a simulation box filled with atoms. Water molecules were placed on a rectangular grid with 3.1 Å spacing, and in the center, was a solid particle composed of carbon and hydrogen atoms placed on a circular pattern.
Figure 1. Scaled-down (computed simulation consisted of almost 4 million atoms) image of a simulation box filled with atoms. Water molecules were placed on a rectangular grid with 3.1 Å spacing, and in the center, was a solid particle composed of carbon and hydrogen atoms placed on a circular pattern.
Polymers 15 00787 g001
Figure 2. Sample traces of particles for models: SPC (a); SPCE (b); TIP3P, in this case, there was no simulation with a 10 pN/nm trap stiffness (c); and TIP4P, with only a freely moving particle (d).
Figure 2. Sample traces of particles for models: SPC (a); SPCE (b); TIP3P, in this case, there was no simulation with a 10 pN/nm trap stiffness (c); and TIP4P, with only a freely moving particle (d).
Polymers 15 00787 g002
Figure 4. Picture of the experimental chamber, with one inlet, filled with pure water, in which the measurements on single particles were made. The inlet is at the top. It is part of a bigger microchannel chip made from PDMS, each chip containing around a dozen of such wells, each connected to one big channel (above the inlet, not shown here). Dark spots on the image are from dust on the optics and camera sensor.
Figure 4. Picture of the experimental chamber, with one inlet, filled with pure water, in which the measurements on single particles were made. The inlet is at the top. It is part of a bigger microchannel chip made from PDMS, each chip containing around a dozen of such wells, each connected to one big channel (above the inlet, not shown here). Dark spots on the image are from dust on the optics and camera sensor.
Polymers 15 00787 g004
Table 1. Parameters of atoms used in each simulation.
Table 1. Parameters of atoms used in each simulation.
Water ModelSPCSPCETIP3PTIP4P
O-H distance1 Å1 Å0.9572 Å0.9572 Å
H-H angle109.47 O109.47 O104.52 O104.52 O
O charge−0.82 e−0.8476 e−0.83 e−1.04844 e
H charge0.41 e0.4238 e0.415 e0.52422 e
ϵ of O0.1553 kcal mol−10.1553 kcal mol−10.102 kcal mol−10.16275 kcal mol−1
σ of O3.166 Å3.166 Å3.188 Å3.16435 Å
ϵ of C0.644 kcal mol−10.644 kcal mol−10.644 kcal mol−10.644 kcal mol−1
σ of C3.554 Å3.554 Å3.554 Å3.554 Å
ϵ of H0 kcal mol−10 kcal mol−10 kcal mol−10 kcal mol−1
σ of H0 Å0 Å0 Å0 Å
O—angle degrees, superscript. e—electron charge, normal, no superscript etc.
Table 2. Model performance and resource usage.
Table 2. Model performance and resource usage.
ModelComputation Time
[Time Steps/s]
Memory Usage [GB]Time Step Value with No Force [fs]Time Step Value with Force Applied [fs]
SPC0.1799.83.253.25
SPCE0.1549.83.253.25
TIP3P0.144103 0.75
TIP4P0.05319.53-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zembrzycki, K.; Pawłowska, S.; Pierini, F.; Kowalewski, T.A. Brownian Motion in Optical Tweezers, a Comparison between MD Simulations and Experimental Data in the Ballistic Regime. Polymers 2023, 15, 787. https://doi.org/10.3390/polym15030787

AMA Style

Zembrzycki K, Pawłowska S, Pierini F, Kowalewski TA. Brownian Motion in Optical Tweezers, a Comparison between MD Simulations and Experimental Data in the Ballistic Regime. Polymers. 2023; 15(3):787. https://doi.org/10.3390/polym15030787

Chicago/Turabian Style

Zembrzycki, Krzysztof, Sylwia Pawłowska, Filippo Pierini, and Tomasz Aleksander Kowalewski. 2023. "Brownian Motion in Optical Tweezers, a Comparison between MD Simulations and Experimental Data in the Ballistic Regime" Polymers 15, no. 3: 787. https://doi.org/10.3390/polym15030787

APA Style

Zembrzycki, K., Pawłowska, S., Pierini, F., & Kowalewski, T. A. (2023). Brownian Motion in Optical Tweezers, a Comparison between MD Simulations and Experimental Data in the Ballistic Regime. Polymers, 15(3), 787. https://doi.org/10.3390/polym15030787

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop