The solution properties, including hydrodynamic quantities and the radius of gyration, of globular proteins are calculated from their detailed, atomic-level structure, using bead-modeling methodologies described in our previous article (, Biophys. J. 76:3044-3057). We review how this goal has been pursued by other authors in the past. Our procedure starts from a list of atomic coordinates, from which we build a primary hydrodynamic model by replacing nonhydrogen atoms with spherical elements of some fixed radius. The resulting particle, consisting of overlapping spheres, is in turn represented by a shell model treated as described in our previous work. We have applied this procedure to a set of 13 proteins. For each protein, the atomic element radius is adjusted, to fit all of the hydrodynamic properties, taking values close to 3 A, with deviations that fall within the error of experimental data. Some differences are found in the atomic element radius found for each protein, which can be explained in terms of protein hydration. A computational shortcut makes the procedure feasible, even in personal computers. All of the model-building and calculations are carried out with a HYDROPRO public-domain computer program.

No abstract

No abstract

We propose a partially flexible, double-helical model for describing the conformational and dynamic properties of DNA. In this model, each nucleotide is represented by one element (bead), and the known geometrical features of the double helix are incorporated in the equilibrium conformation. Each bead is connected to a few neighbor beads in both strands by means of stiff springs that maintain the connectivity but still allow for some extent of flexibility and internal motion. We have used Brownian dynamics simulation to sample the conformational space and monitor the overall and internal dynamics of short DNA pieces, with up to 20 basepairs. From Brownian trajectories, we calculate the dimensions of the helix and estimate its persistence length. We obtain translational diffusion coefficient and various rotational relaxation times, including both overall rotation and internal motion. Although we have not carried out a detailed parameterization of the model, the calculated properties agree rather well with experimental data available for those oligomers.

A Brownian dynamics simulation of a saturated hydrocarbon chain with simple mean-field potentials, namely anchorage, orientation and enclosing, reproducing a biological membrane environment is presented. The simulation was performed for a time equivalent to 1.4 micros thanks to the simplicity of our model. The results are compared with those obtained for a hydrocarbon chain simulated in the absence of the membrane potentials but with confinement. With the appropriate choice of parameters, equilibrium properties, such as deuterium order parameter, chain length, tilt angle and geometry, and dynamic properties, such as dihedral angle transition rate, rotational and translational diffusion, recovered from our simulations, correctly reproduced, are consistent with hydrocarbon-derived molecule experimental results and simulation results obtained from other more complex studies.