HYDRONMR: prediction of NMR relaxation of globular proteins from atomic-level structures and hydrodynamic calculations.

The heteronuclear NMR relaxation of globular proteins depends on the anisotropic rotational diffusion tensor. Using our previous developments for prediction of hydrodynamic properties of arbitrarily shaped particles, by means of bead models, we have constructed a computational procedure to calculate the rotational diffusion tensor and other properties of proteins from their detailed, atomic-level structure. From the atomic coordinates file used to build the bead model, the orientation of the pertinent dipoles can be extracted and combined with the hydrodynamic information to predict, for each residue in the protein, the relaxation times. All of these developments have been implemented in a computer program, HYDRONMR, which will be of public domain.

Calculation of hydrodynamic properties of globular proteins from their atomic-level structure.

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.

Conformation and dynamic properties of a saturated hydrocarbon chain confined in a model membrane: a Brownian dynamics simulation.

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.

Simulation of the distribution and diffusion of a rigid amphipathic particle embedded in a model membrane.

We simulate, by Brownian dynamics, the distribution, orientation and diffusion of a rigid molecule, represented as a dumbbell, with amphipathic nature, embedded in a model membrane. The significant features of a biological membrane are reproduced by means of a Maier-Saupe orienting potential, an enclosing potential and a lipophobic potential. We also evaluate the equilibrium quantities, such as order parameter, and dynamic features, such as rotational or translational diffusivity, of the embedded molecule in terms of the system parameters and compare the obtained results with those obtained from model independent theory.

Simulation of the conformation and dynamics of a double-helical model for DNA.

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.

Molecular dynamics simulation of a dye molecule in the interior of a bilayer: 1,6-diphenyl-1,3,5-hexatriene in dipalmitoylphosphatidylcholine.

A molecular dynamics simulation was carried out for a dipalmitoylphosphatidylcholine (DPPC) membrane in its liquid crystalline state containing different concentrations of the dye molecule 1,6-diphenyl-1,3,5-hexatriene (DPH). From a numerical analysis of the trajectories, we obtained information concerning structural changes of the membrane due to the presence of the probe and some hydrodynamic information concerning the probe itself. The hydrodynamic properties regarding dye molecules that have been reported in this article are: rotational and translational diffusion coefficient and relaxation times. From this analysis, we estimated a range of values of 0.6-0.9 cP for the micro-viscosity in the mid-membrane. These simulations also afforded us some information regarding structural changes in the membrane as a consequence of the presence of the fluorescent dyes at different concentrations. Thus, the disorder inside the membrane, the surface area per lipid and thickness of the membrane were also investigated.

Distribution and diffusivity of a hydrophobic probe molecule in the interior of a membrane: theory and simulation.

We propose a simple model for the distribution of position and orientation and the diffusion of a hydrophobic probe molecule embedded in a membrane. The molecule experiences both a Maier-Saupe orienting potential as well as an enclosing potential of repulsion from the membrane walls. A statistical thermodynamics treatment of the model provides predictions of the location and orientation of the molecule within the membrane. In particular, we evaluate the order parameter of the molecule in terms of the model constants. The diffusivity of the probe is studied by Brownian dynamics simulation. For rotational diffusion, we check an available analytical approximate treatment that allows for the prediction of the dynamics in terms of equilibrium quantities. We also pay attention to quantities related to the initial and mean reorientational rate of the probe. For translational diffusion, we use the simulation results to analyze some general aspects of lateral and transversal diffusion.

Medical social work in the 'Centro Nacional de Paraplejicos', Toledo, Spain.

The Medical Social Service in the 'Centro Nacional de Rehabilitación de Parapléjicos' of Toledo is briefly explained. Until 30 April 1975, 141 patients were treated in this hospital and the pertinent statistics are shown. The main social problems were housing, pension and job. The Medical Social Service makes use of all its resources to solve these problems.