Normal mode theory for the Brownian dynamics of a weakly bending rod: Comparison with Brownian dynamics simulations

A normal mode theory is developed for the Brownian dynamics of weakly bending rods with preset hydrodynamic interactions. The rod is replaced by a chain of contiguous spheres whose radius is chosen to yield the appropriate uniform translational and rotational diffusion coefficients. Despite the inclusion of preset hydrodynamic interactions in the dynamical operator, its normal modes are not coupled by the potential energy, so their amplitudes remain pairwise “orthogonal” under equilibrium averaging. The uniform translational and rotational diffusion coefficients obtained from Langevin theory are shown to be identical to those obtained from the Kirkwood algorithm, despite their rather different appearance. An expression is given for the mean squared angular displacement 〈Δ_xm(t)²〉 of the mth bond vector around the instantaneous x axis (perpendicular to the end-to-end vector z). Necessary algorithms are presented for the numerical evaluation of all quantities. The normal mode theory is compared with Brownian dynamics simulations for the same model by examining 3〈Δ_xm(t)²〉 for the central bond vector of rods comprising 10 and 30 subunits with various persistence lengths. The normal mode theory works very well for all times for L/P ? 0.6, where P = κ/k_BT is the persistence length and κ is the bending rigidity. With increasing flexibility, the domain of validity of the normal mode theory is restricted to shorter times, where violations of the weak bending approximation are less severe. However, increasing the length of the rod from 10 to 30 subunits yields improved agreement with the simulations for the same and even longer times. This latter effect is tentatively attributed to the greater fluctuating tension in the longer chains, which acts to retard the rotational relaxation in the simulations, but is not taken into account in the present normal mode theory.     

Dynamic light scattering from thin rigid rods: Anisotropy of translational diffusion of tobacco mosaic virus

An Exact theoretical expression for the apparent diffusion coefficient D_app(K) of a thin rigid rod with arbitrary anisotropy of its translational diffusion diffusion coefficient is derived from the first cumulant of its dynamic structure factor. D_app(K) is predicted to reach a limiting plateau value at extermely large values of KL, where K is the scattering vector and L the rod length. Howerver, that limiting plateau value is approached only very slowly along a quasi-plateau with a very gradual slope. Dynamic light-scattering studies have been performed on tobacco mosaic virus from K² = (0.4–20) × 10¹⁰ cm^?2 using 632-8-nm laser radiation. The present data yield D₀ = (4.19 ± 0.10) × 10^?8 cm²/s (corrected to 20,w conditions) and, with literature data to establish L = 2980 Å and the rotational diffusion coefficient D_R = 318s^?1, yield also Δ ≡ D_∥ ? D_⊥ = (1.79 ± 0.38) × 10^?8 cm²/s. The experimental data closely follow the curve of D_app(K) vs K² calcuated for these parameters. The present value of D₀ substantially exceeds all previous dynamic light-scattering values, but is in good aggreement with previous sedimentation data, which were confirmed for the presemt sample. The anisotropy ratio Δ/D₀ = 0.43 ± 0.09 is in accord with theoretical predictions based on the modified Kirkwood algorithm, despite the fact the D₀ lies significantly below its corresponding theoretical value. The present data largely predlude the possibility that both D₀ and Δ/D₀ could simultaneously match their theoretical predictions. We present a detailed comparison of the experimental data with the calculations of Tirado and Garcia de la Torre based on the modified Kirkwood algorithm and with the Broersma formulas.     

Monte Carlo study of cycloamylose: chain conformation, radius of gyration, and diffusion coefficient

Cyclic (1 --> 4)-alpha-D-glucan chains with or without excluded volume have been collected from a huge number (about 10(7)) of linear amylosic chains generated by the Monte Carlo method with a conformational energy map for maltose, and their mean-square radii of gyration and translational diffusion coefficients D (based on the Kirkwood formula) have been computed as functions of x (the number of glucose residues in a range from 7 to 300) and the excluded-volume strength represented by the effective hard-core radius. Both /x and D in the unperturbed state weakly oscillate for x < 30 and the helical nature of amylose appears more pronouncedly in cyclic chains than in linear chains. As x increases, these properties approach the values expected for Gaussian rings. Though excluded-volume effects on them are always larger in cycloamylose than in the corresponding linear amylose, the ratios of and the hydrodynamic radius of the former to the respective properties of the latter in good solvents can be slightly lower than or comparable to the (asymptotic) Gaussian-chain values when x is not sufficiently large. An interpolation expression is constructed for the relation between the gyration-radius expansion factors for linear and cyclic chains from the present Monte Carlo data and the early proposed asymptotic relation with the aid of the first-order perturbation theories.     

Sedimentation coefficient of polyoma virus DNA

The sedimentation coefficient of the twisted circular form of polyoma virus DNA is calculated from the Kirkwood sedimentation–diffusion equation, the structure being assumed to be a rigid double superhelix. Agreement with the experimental sedimentation coefficient can be obtained, with the use of an experimental value for the number of superhelical turns, when the pitch of the superhelix is intermediate between its minimal and maximal possible values. Another model, which has been proposed for polyoma DNA at low ionic strengths, may be visualized as a superhelical structure wound about a torus. Calculations of sedimentation coefficients for this model agree qualitatively with experimental data at ionic strengths Below 10^?2M.     

Optical activity of aromatic chromophores. I. O, m, and p-tyrosine

The optical activity and the conformational energy of the amino acids tyrosine, o-tyrosine, and m-tyrosine have been calculated as a function of molecular conformation. A new graphic technique, rotatory strength–conformational energy maps, was developed for the presentation of these calculations. Experimental circular dichroism spectra were determined by utilization of a new experimental technique that involves complete computer control of a Cary spectropolarimeter. This permitted repetitive scanning with signal averaging over extended time periods and resulted in a greatly enhanced signal to noise ratio. On the basis of these calculations, it, was concluded that the Kirkwood “coupled oscillator” mechanism is capable of explaining the observed optical activity of the aromatic chromophores of the molecules that were investigated.     

Transport properties of oligomeric subunit structures

The translational friction coefficients, rotational friction coefficient, and intrinsic viscosity of rigid regular structures composed of up to eight identical spherical subunits have been accurately calculated. The aim of this calculation is to interpret the hydrodynamic properties of oligomeric subunit proteins. To avoid the well-known failure of the theory in the evaluation of rotational coefficients and intrinsic viscosities, each subunit is hydrodynamically modeled as a polyhedral array of smaller spheres. The analysis of several alternatives suggests that a cubic array is the best choice. The reliability of this strategy is checked by comparison of the calculated values for all the transport properties of a sphere and the translational friction coefficients of a dimer with their exact values. Finally, the hydrodynamic properties of a number of subunit structures with varying number of subunits and different geometries are tabulated.     

On the Approximation of Solvent Effects on the Conformation and Dynamics of Cyclosporin A by Stochastic Dynamics Simulation Techniques

Abstract

The molecular simulation technique of stochastic dynamics (SD) is tested by application to the immunosuppressive drug cyclosporin A (CPA). Two stochastic dynamics simulations are performed, one (SD_CCl4 ) with atomic friction coefficients proportional to the viscosity of the nonpolar solvent CCl₄, and one (SD_H2O) with atomic friction coefficients corresponding to an aqueous solution. The atomic friction coefficients are also taken proportional to an approximate expression for the atomic accessible surface area. The properties of both stochastic dynamics simulations are compared to those of two full molecular dynamics (MD) simulations of cyclosporin A, one in a box with 591 CCl₄ molecules, and one in a box with 632 H₂O molecules.

The properties of cyclosporin A as found in the molecular dynamics simulation in CCl₄ are well reproduced by the SD_CCl4 simulation. This indicates that the neglect of a mean force reresenting the average solvent effects on the solute is justified in the case of nonpolar solvents. For polar solvents, like water, this mean force may not be neglected. The SD_H2O simulation of cyclosporin A clearly fails to reproduce the amount of hydrogen bonding found in the molecular dynamics stimulation of cyclosporin A in water.

A comparison with a molecular dynamics simulation of cyclosporin A in vacuo shows that both the SD_CCl4 and the SD_H2O simulation come closer to the properties of the molecular dynamics simulations in CCl₄ and in H₂O than a molecular dynamics simulation in vacuo.     

Intensity-fluctuation spectroscopy and tRNA conformation. II. Changes of size and shape of tRNA in the melting process

We have measured the translational (D^T) and rotational (D^R) diffusion coefficients of bulk tRNA from baker's yeast during the thermal unfolding process by means of photon-correlation spectroscopy. It should be noted that our estimate of the rotational diffusion coefficient represented, for the first time, measurements on a small macromolecule in solution by the photoelectron time-of-arrival technique with a delay-time resolution of 1 nsec. The melting curves expressed in terms of δD^T vs temperature were consistent with the literature data in revealing the melting steps and their dependence on NaCl concentration. Additionally, it was possible to prove the existence of an intermediate, more compact structure during the initial steps of the thermal unfolding process. We found that the temperature ranges over which this intermediate structure appears depend strongly on salt concentration. By utilizing both translational and rotational diffusion coefficients and Perrin's equations for ellipsoids of revolution, we have computed the values of the equivalent length and width of tRNA molecules in solution at four different temperatures for NaCl concentrations of 0.2, 0.5, and 1M. The approximate model of ellipsoids of revolution also permits us to obtain an estimate of the radius of gyration, which is in very good agreement with literature data measured by means of small-angle x-ray scattering. Furthermore, we have measured the shape and size changes of tRNA with varying NaCl concentrations at room temperatures (25°C). The molecule becomes smaller and more spherical when NaCl concentration increases. As a result of partial melting at 70°C, the macromolecule is surprisingly elongated with an approximate axial ratio of 8:1 and has dimensions of about 180/22Å. Such information on conformational changes by a simultaneous determination of rotational and translational diffusion coefficients illustrates the potential of this approach, not available by other methods.     

Hydrogen Bonding and Stacking of DNA Bases: A Review of Quantum-chemical ab initio Studies

Abstract

Ab initio quantum-chemical calculations with inclusion of electron correlation made since 1994 (such reliable calculations were not feasible before) significantly modified our view on interactions of nucleic acid bases. These calculations allowed to perform the first reliable comparison of the strength of stacked and hydrogen bonded pairs of nucleic acid bases, and to characterize the nature of the base-base interactions. Although hydrogen-bonded complexes of nucleobases are primarily stabilized by the electrostatic interaction, the dispersion attraction is also important. The stacked pairs are stabilized by dispersion attraction, however, the mutual orientation of stacked bases is determined rather by the electrostatic energy. Some popular theories of stacking were ruled out: The theory based on attractive interactions of polar exocyclic groups of bases with delocalized electrons of the aromatic rings (Bugg et al., Biopolymers 10, 175 (1971).), and the II-II interactions model (C.A. Hunter, J. Mol. Biol. 230, 1025 (1993)). The calculations demonstrated that amino groups of nucleobases are very flexible and intrinsically nonplanar, allowing hydrogen-bond-like interactions which are oriented out of the plane of the nucleobase. Many H-bonded DNA base pairs are intrinsically nonplanar. Higher-level ab initio calculations provide a unique set of reliable and consistent data for parametrization and verification of empirical potentials. In this article, we present a short survey of the recent calculations, and discuss their significance and limitations. This summary is written for readers which are not experts in computational quantum chemistry.     

Electrostatic effects on polynucleotide transitions. I. Behavior at neutral pH

An approximate analytical expression for the electrostatic free energy of a polynucleo-tide in any of its possible ordered or random conformations is derived by integration of the screened-Coulomb potential energy function over all charge pairs in the structure. The electrostatic free energy of any form is found to be a linear function of the logarithm of the monovalent counterion concentration, in the range of low salt concentrations. Hence the electrostatic free energy difference between ordered and disordered forms in a polynucleotide structural transition is a linear function of the logarithm of the monovalent counterion concentration. A free energy balance applied to a two-state model for the transition then yields a linear dependence of the transition temperature T_m upon the logarithm of the counterion concentration. Calculation of the quantity dT_m/d log M, where M is the monovalent counterion concentration, shows it to be a characteristic constant for a given transition, with a magnitude and sign proportional to the charge density difference between the ordered and disordered forms. Use of any one of several alternate, simple assumptions yields predicted dT_m/d log M values in good agreement with experimental data for various polynucleotide transitions.     

Hydrodynamic properties of macromolecular complexes. IV. Intrinsic viscosity theory,with applications to once-broken rods and multisubunit proteins

We have extended our previous theories of the translational and rotational frictional properties of multisubunit complexes to calculate the intrinsic viscosity of such structures. Our theory is similar to those recently construced by McCammon and Deutch, and by Nakajima and Wada, in that it uses a modified hydrodynamic interaction tensor and solves the system of simultaneous interaction equations by digital computation rather than by successive approximations. However, there are some differences in the formulation and averaging of these equations. Extensive numerical comparison is made between this theory and others that are available—associated with the names of Hearst and Tagami, Abdel-Khalik and Bird, and Tsuda—using as a basis exact results for prolate ellipsoids of revolution. For large axial ratios, only our theory asymptotically approaches the correct limit; but for small axial ratios, only the Tsuda “shell-model” theory is adequate, because the other theories neglect the preponderant influence of the sphere located at the center of rotation. Intrinsic viscosities, translational frictional coefficients, and Scheraga-Mandelkern β parameters, are tabulated for a large number of polygonal and polyhedral subunit structures, with up to eight elements, using both our theory and Tsuda's. Particular application is made to hemerythrin and aspartate transcarbamylase. Finally, the viscosities and friction coefficients o once-broken rods are calculated and compared with an approximate theory by Wilenski.     

Calculation of translational friction and intrinsic viscosity. I. General formulation for arbitrarily shaped particles.

A general method for calculating translational friction and intrinsic viscosity is developed through exploiting relations between hydrodynamics and electrostatics. An approximate relation xi = 6 pi eta 0C between the translational friction coefficient xi of a particle (eta 0: solvent viscosity) and its capacitance C was derived previously. This involved orientationally preaveraging the Oseen tensor, but the result was found to be very accurate. Based on preaveraging, we find that the intrinsic viscosity [eta] of a particle can be estimated from its polarizability alpha through [eta] = 3/4 alpha + 1/4 Vp, where Vp is the volume of the particle. Both the capacitance and the polarizability can be obtained in a single calculation using the boundary-element technique. An efficient approach is thus found for estimating [eta], a quantity that is very useful in practice because of its sensitivity to particle shape but is notoriously difficult to calculate. Illustrative calculations on ellipsoids, cylinders, and dumbbells demonstrate both the accuracy of the approximate relations and the efficiency of the present method.     

Diffusion coefficients of segmentally flexible macromolecules with two spherical subunits

The translational and rotational diffusion coefficients have been calculated for a simple, segmentally flexible model: the hinged dumbbell (HD). In the HD, two spherical subunits are attached to an universal joint by means of frictionless connectors. In addition to the case in which hydrodynamic interactions are neglected (NI), we have also considered two more cases, including hydrodynamic interaction by means of the Kirkwood-Riseman approximate treatment (KR) and using accurate procedure based in the series expansions for the two-sphere diffusion tensor (SE). Expressions for the friction coefficients of the HD are given for the three cases, and the diffusion coefficients are evaluted inverting the 9 × 9 resistance matrix, for two HDs with different dimensions. The KR treatment, which includes a contribution from the finite volume of the subunits, is shown to be an excellent approximation to the more rigorous procedure. In the NI case for rotation, the various coefficients present different deviations with respect to the SE results. A rough estimate of the errors of the NI relaxation times indicates that they may be smaller than 15% for a HD with identical beads. However, the influence of hydrodynamic interaction should be more important for the rotational diffusivity of a small sphere attached to a larger one. The error of the NI result for the translational diffusion coefficient is of about 25% for the two HDs.     

Effect of drop interactions on the performance of hydrocarbon fermentors

A drop-interaction model devised by the authors, which successfully predicts experimentally observed drop-size distributions in stirred batch vessels, has been used to assess the importance of droplet coalescence and redispersion in hydrocarbon fermentation yields. Calculations, which have been made with both n-alkane and gas–oil fermentations, have been compared with previous calculations based on simpler and, hence, less realistic droplet interaction mechanisms. The comparison shows that the differences in yield are not spectacularly different.     

Translational diffusion coefficient of cycloamylose in aqueous sodium hydroxide

Seven cyclic (1 --> 4)-alpha-D-glucan (cycloamylose) samples ranging in weight-average molecular weight from 5 x 10(3) to 1.8 x 10(4) and gamma-cyclodextrin have been studied by sedimentation equilibrium in dimethylsulfoxide (at 25 degrees C) and by dynamic light scattering in 0.5 N aqueous sodium hydroxide (at 25 degrees C), a good solvent for linear amylose. The measured translational diffusion coefficients D in the aqueous NaOH agree fairly closely with previous Monte Carlo results for cyclic (1 --> 4)-alpha-D-glucan chains with excluded volume, when correction is made for the effects of bead diameter and fluctuating hydrodynamic interaction (HI) on the Kirkwood theory on which the computation of D was based. These D data are also explained almost quantitatively by Yamakawa and Fujii's expression for the associated KP ring (based on the Kratky-Porod wormlike chain) with the molecular parameters for linear amylose if the fluctuating HI and excluded-volume effects are taken into account. It is concluded that the translational diffusion behavior of cycloamylose in the aqueous NaOH is consistent with the conformational characteristics derived from the conformational energy of maltose and dilute-solution data for linear amylose.     

Ice friction during speed skating.

During speed skating, the external power output delivered by the athlete is predominantly used to overcome the air and ice frictional forces. Special skates were developed and used to measure the ice frictional forces during actual speed skating. The mean coefficients of friction for the straights and curves were, respectively, 0.0046 and 0.0059. The minimum value of the coefficient of ice friction was measured at an ice surface temperature of about -7 degrees C. It was found that the coefficient of friction increases with increasing speed. In the literature, it is suggested that the relatively low friction in skating results from a thin film of liquid water on the ice surface. Theories about the presence of water between the rubbing surfaces are focused on the formation of water by pressure-melting, melting due to frictional heating and on the 'liquid-like' properties of the ice surface. From our measurements and calculations, it is concluded that the liquid-like surface properties of ice seem to be a reasonable explanation for the low friction during speed skating.     

Transport properties of rigid and flexible macromolecules by brownian dynamics simulation

The method of Ermak and McCammon [Ermak, D. L. & McCammon, J. A. (1978) J. Chem. Phys. 69 , 1352–1360] is used to simulate the Brownian dynamics of a system of identical, interacting beads. In the present study, we use the method to obtain transport coefficients for a variety of rigid and flexible structures modeled as arrays of spherical subunits. Constraints are enforced using the SHAKE algorithm or a modification, SHAKE-HI, that is described for the first time. In SHAKE-HI, hydrodynamic interactions between subunits are taken into account when the constraints are enforced. Use of SHAKE-HI yields transport coefficients that are in perfect agreement with those obtained by other methods. The primary advantage of the present method is its generality. We also propose that multistep Brownian dynamics may be important in simulating actual experiments (such as fluorescence depolarization) on well-defined model systems that possess an arbitrary degree of internal flexibility.     

Finite-size effects of Kirkwood–Buff integrals from molecular simulations

The modelling of thermodynamic properties of liquids from local density fluctuations is relevant to many chemical and biological processes. The Kirkwood–Buff (KB) theory connects the microscopic structure of isotropic liquids with macroscopic properties such as partial derivatives of activity coefficients, partial molar volumes and compressibilities. Originally, KB integrals were formulated for open and infinite systems which are difficult to access with standard Molecular Dynamics (MD) simulations. Recently, KB integrals for finite and open systems were formulated (J Phys Chem Lett. 2013;4:235). From the scaling of KB integrals for finite subvolumes, embedded in larger reservoirs, with the inverse of the size of these subvolumes, estimates for KB integrals in the thermodynamic limit are obtained. Two system size effects are observed in MD simulations: (1) effects due to the size of the simulation box and the size of the finite subvolume embedded in the simulation box, and (2) effects due to computing radial distribution functions (RDF) from a closed and finite system. In this study, we investigate the two effects in detail by computing KB integrals using the following methods: (1) Monte Carlo simulations of finite subvolumes of a liquid with an analytic RDF and (2) MD simulations of a WCA mixture for various simulation box sizes, but at the same thermodynamic state. We investigate the effect of the size of the simulation box and quantify the differences compared to KB integrals computed in the thermodynamic limit. We demonstrate that calculations of KB integrals should not be extended beyond half the size of the simulation box. For finite-size effects related to the RDF, we find that the Van der Vegt correction (J Chem Theory Comput. 2013;9:1347) yields the most accurate results.