PubMed Health⌕ Search

Biomedical subjects

J García de la Torre

Publications and source records attributed to J García de la Torre.

At least 19 recordsLinked to original sources
PubMed2006

Dynamic electro-optic properties of macromolecules and nanoparticles in solution: a review of computational and simulation methodologies.

This paper reviews some theories, and computational and simulation procedures available for the calculation of the time-course of electro-optic properties of particles in solution. For rigid particles, the time evolution of the properties is directly related to their rotational diffusion; therefore, the computational procedures for the calculation of hydrodynamic properties find a direct application in electro-optics. Several of such computational procedures, based on bead models, are reviewed. For flexible particles, the simultaneous effects the external field and the flexibility can be treated with Brownian dynamics simulation. We illustrate the various procedures, with applications to rigid bent rods and flexible, wormlike or hinged rods, trying to show how the absence or presence of flexibility, and its kind, influences the dynamic electro-optic properties, which are therefore valuable sources of information about the conformation of macromolecules and nanoparticles.

Birefringence↗
PubMed2005

Transient electric birefringence of wormlike macromolecules in electric fields of arbitrary strength: a computer simulation study.

We have studied the birefringence decay of linear models of macromolecules for two different types of flexibility, the broken-rod chain and the wormlike chain, using a computer simulation of a transient electric birefringence experiment. We have paid particular attention to the influence of the intensity of the orienting field, including two orienting mechanisms, the induced dipole, and the permanent dipole. We have compared wormlike and broken-rod models of the same radius of gyration, finding that they present a different decay curve under the influence of the same intensity of the field. We have seen that these differences are due to the faster relaxation times (smaller in the wormlike chain model) and amplitudes, because, regardless of the type of flexibility, the overall size of a molecule (measured by the radius of gyration) essentially determines the longest relaxation time. We have also analyzed how the relaxation process is affected by the degree of flexibility, the orientation mechanisms, and the intensity of the field. Studying a different aspect, we have paid attention to the deformation of a molecule in a transient electric birefringence experiment as a source of information. In this work we have developed equations to characterize this deformation in terms of one of the components of the gyration tensor, if a dynamic light scattering experiment under the influence of an electric field could be performed. To develop this work we have simulated the Brownian dynamics of the different models, relaxing after the removal of an orienting external electric field of arbitrary strength. A comparison with other methods such a the rigid body treatment or the correlation analysis of Brownian trajectories has also been included. We have seen that differences between the two Brownian dynamics methods are small and that the rigid-body treatment is only an acceptable approximation to obtain the longest relaxation time.

Journal Article↗
PubMed2005

Hydrodynamic models and computational methods for NMR relaxation.

Interpretation of NMR relaxation data of macromolecules is based on the analysis of their dynamic behavior in solution. For quasirigid molecules, in addition to a minor, separable contribution from local mobility, the main contribution corresponds to the overall rotational diffusion of the complete molecule. Therefore, theoretical descriptions and computational methodologies for hydrodynamic calculations, which yield the full, anisotropic rotational diffusion tensor of rigid molecules, are extremely helpful in the analysis of NMR relaxation. Recent approaches allow realistic predictions of the rotational diffusion tensor from structures at atomic detail. This enables measured relaxation rates and structural models to be compared. Such a comparison (1) provides an independent test of the structural model, (2) provides a framework for the interpretation of local motion, even for highly anisotropic systems, (3) provides a simple method for the detection of additional sources of relaxation, such as chemical exchange, and (4) provides a sensitive method for the detection of nonspecific aggregation or oligomer formation. Although hydrodynamic calculations usually assume a rigid structure, Brownian dynamics simulations extend their range of applications to flexible multidomain structures. Hydrodynamic applications are not restricted to globular proteins. Small DNA fragments, which could be otherwise considered cylindrical objects, can also be treated with atomic detail using the same methodology used for proteins.

Computational Biology↗
PubMed2003

Calculation of the solution properties of flexible macromolecules: methods and applications.

While the prediction of hydrodynamic properties of rigid particles is nowadays feasible using simple and efficient computer programs, the calculation of such properties and, in general, the dynamic behavior of flexible macromolecules has not reached a similar situation. Although the theories are available, usually the computational work is done using solutions specific for each problem. We intend to develop computer programs that would greatly facilitate the task of predicting solution behavior of flexible macromolecules. In this paper, we first present an overview of the two approaches that are most practical: the Monte Carlo rigid-body treatment, and the Brownian dynamics simulation technique. The Monte Carlo procedure is based on the calculation of properties for instantaneous conformations of the macromolecule that are regarded as if they were instantaneously rigid. We describe how a Monte Carlo program can be interfaced to the programs in the HYDRO suite for rigid particles, and provide an example of such calculation, for a hypothetical particle: a protein with two domains connected by a flexible linker. We also describe briefly the essentials of Brownian dynamics, and propose a general mechanical model that includes several kinds of intramolecular interactions, such as bending, internal rotation, excluded volume effects, etc. We provide an example of the application of this methodology to the dynamics of a semiflexible, wormlike DNA.

Centrifugation, Density Gradient↗
PubMed2002

Joint determination by Brownian dynamics and fluorescence quenching of the in-depth location profile of biomolecules in membranes.

The in-depth molar distribution function of fluorophores is revealed by a new methodology for fluorescence quenching data analysis in membranes. Brownian dynamics simulation was used to study the in-depth location profile of quenchers. A Lorentzian profile was reached. Since the Stern-Volmer equation is valid at every depth in the membrane for low quencher concentrations, the molar distribution of the fluorophore (also regarded as a Lorentzian) can be achieved. The average location and the broadness of the fluorophore distribution can be calculated. The importance of the knowledge of the location width is demonstrated and discussed, since this parameter reveals important conclusions on structural features of the interaction of membranes with probes and biomolecules (e.g., conformational freedom in proteins), as well as photophysical properties (e.g., differential fluorophore quantum yields). Subsequent use of this methodology by the reader does not, necessarily, involve the performance of simulations and is not limited to the use of Lorentzian function distributions.

Cell Membrane↗
PubMed2002

Calculation of hydrodynamic properties of small nucleic acids from their atomic structure.

Hydrodynamic properties (translational diffusion, sedimentation coefficients and correlation times) of short B-DNA oligonucleotides are calculated from the atomic-level structure using a bead modeling procedure in which each non-hydrogen atom is represented by a bead. Using available experimental data of hydrodynamic properties for several oligonucleotides, the best fit for the hydrodynamic radius of the atoms is found to be approximately 2.8 A. Using this value, the predictions for the properties corresponding to translational motion and end-over-end rotation are accurate to within a few percent error. Analysis of NMR correlation times requires accounting for the internal flexibility of the double helix, and allows an estimation of approximately 0.85 for the Lipari-Szabo generalized order parameter. Also, the degree of hydration can be determined from hydrodynamics, with a result of approximately 0.3 g (water)/g (DNA). These numerical results are quite similar to those found for globular proteins. If the hydrodynamic model for the short DNA is simply a cylindrical rod, the predictions for overall translation and rotation are slightly worse, but the NMR correlation times and the degree of hydration, which depend more on the cross-sectional structure, are more severely affected.

Base Sequence↗
PubMed2002

Hydrodynamic properties of rigid macromolecules composed of ellipsoidal and cylindrical subunits.

A procedure is devised for the calculation of hydrodynamic properties of rigid macromolecules composed subunits that are modeled as ellipsoids of revolution and cylinders. Owing to the axial symmetry of these shapes, smooth shell models can be constructured for the subunit structure. The bead shell model so constructed is employed for the calculation of the properties. A computer program, HYDROSUB, has been written implementing both the model building and the hydrodynamic calculation. A detailed example of the use of this methodology is presented for the case of the solution properties of the human antibody molecule immunoglobulin G3 (IgG3). Finally, hints are given on other uses and applications of the procedure.

Biophysical Phenomena↗
PubMed2001

An analytical solution to the problem of the orientation of rigid particles by planar obstacles. Application to membrane systems and to the calculation of dipolar couplings in protein NMR spectroscopy.

Nonspherical particles or molecules experience an ordering effect in the presence of obstacles due to the restrictions they place on the orientation of those molecules that are in their proximity. Obstacles may be the limits of a membrane in which the molecule is embedded, oriented mesoscopic systems such as bicelles, or membrane fragments used to induce weak protein alignment in a magnetic field. The overall shape of most proteins can be described to a good approximation by an ellipsoidal particle. Here we describe and solve analytically the problem of the orientation of ellipsoidal particles by planar obstacles. Simple expressions are derived for the orientational distribution function and the order parameter. These expressions allow the analytical calculation of the residual dipolar couplings for a protein of known three-dimensional structure oriented by steric effects. The results are in good agreement with experiment and with the results of previously described simulations. However, they are obtained analytically in a fraction of the time and therefore open the possibility to include the optimization of the overall shape in the determination of three-dimensional structures using residual dipolar coupling constraints. The equations derived are general and can also be applied to problems of a completely different nature. In particular, previous equations describing the orientation of particles embedded in membranes are verified and generalized here.

Membrane Proteins↗
PubMed2001

Hydration from hydrodynamics. General considerations and applications of bead modelling to globular proteins.

The effect of hydration on hydrodynamic properties of globular proteins can be expressed in terms of two quantities: the delta (g/g) parameter and the thickness of the hydration layer. The two paradigms on hydration that originate these alternative measures are described and compared. For the numerical calculation of hydrodynamic properties, from which estimates of hydration can be made, we employ the bead modelling with atomic resolution implemented in programs HYDROPRO and HYDRONMR. As typical, average values, we find 0.3 g/g and a thickness of only approximately 1.2 A. However, noticeable differences in this parameter are found from one protein to another. We have made a numerical analysis, which leaves apart marginal influences of modelling imperfections by simulating properties of a spherical protein. This analysis confirms that the errors that one can attribute to the experimental quantities suffice to explain the observed fluctuations in the hydration parameters. However, for the main purpose of predicting protein solution properties, the above mentioned typical values may be safely used. Particularly for atomic bead modelling, a hydrodynamic radius of approximately 3.2 A yields predictions in very good agreement with experiments.

Models, Chemical↗
PubMed2001

HYDROMIC: prediction of hydrodynamic properties of rigid macromolecular structures obtained from electron microscopy images.

We have developed a procedure for the prediction of hydrodynamic coefficients and other solution properties of macromolecules and macromolecular complexes whose volumes have been generated from electron microscopy images. Starting from the structural files generated in the three-dimensional reconstructions of such molecules, it is possible to construct a hydrodynamic model for which the solution properties can be calculated. We have written a computer program, HYDROMIC, that implements all the stages of the calculation. The use of this procedure is illustrated with a calculation of the solution properties of the volume of the cytosolic chaperonin CCT, obtained from cryoelectron microscopy images.

Biophysical Phenomena↗
PubMed2000

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.

Hydrocarbons↗
PubMed2000

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.

Algorithms↗
PubMed1999

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.

Diffusion↗
PubMed1999

Calculation of NMR relaxation, covolume, and scattering-related properties of bead models using the SOLPRO computer program.

The hydrodynamic properties of macromolecules and bioparticles, represented by bead models, can be calculated using methods implemented in the computer routine HYDRO. Recently, a new computer routine, SOLPRO, has been presented for the calculation of various SOLution PROperties. These include (1) time-dependent electro-optic and spectroscopic properties related to rotational diffusion, (2) non-dynamic properties like scattering curves, and (3) dimensionless quantities that combine two or more solution properties in a form which depends on the shape of the macromolecule but not on its size. In the present work we describe the inclusion of more of those types of properties in a new version of SOLPRO. Particularly, we describe the calculation of relaxation rates in nuclear magnetic resonance (NMR). For dipolar coupling, given the direction of the dipole the program calculates values of the spectral density, from which the NMR relaxation times can be obtained. We also consider scattering-related properties, namely the distribution of distances for the bead model, which is directly related to the angular dependence of scattered intensity, and the particle's longest distance. We have devised and programmed a procedure to calculate the covolume of the bead model, related to the second virial coefficient and, in general, to the concentration dependence of solution properties. Various shape-dependent dimensionless quantities involving the covolume are calculated. In this paper we also discuss some aspects, namely bead overlapping and hydration, that are not explicitly included in SOLPRO, but should be considered by the user.

Mathematical Computing↗
PubMed1999

Calculation of hydrodynamic properties of macromolecular bead models with overlapping spheres.

For the calculation of hydrodynamic properties of rigid macromolecules using bead modelling, models with overlapping beads of different sizes are used in some applications. The hydrodynamic interaction tensor between unequal overlapping beads is unknown, and an oversimplified treatment with the Oseen tensor may introduce important errors. Here we discuss some aspects of the overlapping problem, and explore an ad hoc form of the interaction tensor, proposed by Zipper and Durchschlag. We carry out a systematic numerical study of the hydrodynamic properties of a two-spheres model, showing how the Zipper-Durchschlag correction removes efficiently the numerical instabilities, and predicts the correct limits.

Biophysical Phenomena↗
PubMed1999

Hydrodynamic properties of rigid particles: comparison of different modeling and computational procedures.

The hydrodynamic properties of rigid particles are calculated from models composed of spherical elements (beads) using theories developed by Kirkwood, Bloomfield, and their coworkers. Bead models have usually been built in such a way that the beads fill the volume occupied by the particles. Sometimes the beads are few and of varying sizes (bead models in the strict sense), and other times there are many small beads (filling models). Because hydrodynamic friction takes place at the molecular surface, another possibility is to use shell models, as originally proposed by Bloomfield. In this work, we have developed procedures to build models of the various kinds, and we describe the theory and methods for calculating their hydrodynamic properties, including approximate methods that may be needed to treat models with a very large number of elements. By combining the various possibilities of model building and hydrodynamic calculation, several strategies can be designed. We have made a quantitative comparison of the performance of the various strategies by applying them to some test cases, for which the properties are known a priori. We provide guidelines and computational tools for bead modeling.

Biophysical Phenomena↗