Search

Many scientists and philosophers of science are troubled by the relative isolation of developmental from evolutionary biology. Reconciling the science of development with the science of heredity preoccupied a minority of biologists for much of the twentieth century, but these efforts were not corporately successful. Mainly in the past fifteen years, however, these previously dispersed integrating programmes have been themselves synthesized and so reinvigorated. Two of these more recent synthesizing endeavours are evolutionary developmental biology (EDB, or "evo-devo") and developmental systems theory (DST). While the former is a bourgeoning and scientifically well-respected biological discipline, the same cannot be said of DST, which is virtually unknown among biologists. In this review, we provide overviews of DST and EDB, summarize their key tenets, examine how they relate to one another and to the study of epigenetics, and survey the impact that DST and EDB have had (and in future should have) on biological theory and practice., (Copyright 2001 John Wiley & Sons, Inc.)

Little is known about the developmental biology of caecilians-tropical, elongate, limbless, mostly fossorial amphibians that are members of the Order Gymnophiona. Ichthyophis kohtaoensis (Family Ichthyophiidae; southeast Asia) is an oviparous species in which maternal care of the clutch is provided. The clutch is laid in a burrow on land, and the embryos develop in their egg membranes, curved around a large yolk mass. Larvae are aquatic and exhibit characteristic features that are not present in the terrestrial adults. Because accurate descriptions of ontogenies and the establishment of standardized stages of embryonic and larval development are useful for both experimental and comparative embryology, a staging table for I.kohtaoensis was developed based on external morphological features. Development from the end of neurulation to metamorphosis was divided into 20 stages. Principal diagnostic features include development of the lateral line organs, formation of three pairs of external gills, development of the eyes, changes in yolk structure, changes in the structure of the cloacal aperture and growth of the tail, including the formation and regression of the tail fin. This study provides a comparison with descriptions of embryonic stages of I.glutinosus and Hypogeophis rostratus and with a recent staging table for the aquatic, viviparous caecilian Typhlonectes compressicauda, the only other caecilians for which reasonably complete ontogenetic information exists in the literature. Comparisons with established staging tables for selected frogs and salamanders are also presented., (Copyright 2000 Wiley-Liss, Inc.)

Using high precision densimetric and ultrasonic measurements, we have determined, at 25 degrees C, the apparent molar volumes PhiV and the apparent molar compressibilities PhiK(S) of four nucleic acid duplexes-namely, the DNA duplex, poly(dIdC)poly(dIdC); the RNA duplex, poly(rA)poly(rU); and the two DNA/RNA hybrid duplexes, poly(rA)poly(dT) and poly(dA)poly(rU). Using available fiber diffraction data on these duplexes, we have calculated the molecular volumes as well as the solvent-accessible surface areas of the constituent charged, polar, and nonpolar atomic groups. We found that the hydration properties of these nucleic acid duplexes do not correlate with the extent and the chemical nature of the solvent-exposed surfaces, thereby suggesting a more specific set of duplex-water interactions beyond general solvation effects. A comparative analysis of our volumetric data on the four duplexes, in conjunction with available structural information, suggests the following features of duplex hydration: (a) The four duplexes exhibit different degrees of hydration, in the order poly(dIdC)poly(dIdC) > poly(dGdC)poly(dGdC) > poly(dAdT)poly(dAdT) approximately poly(dA)poly(dT). (b) Repetitive AT and IC sequences within a duplex are solvated beyond general effects by a spine of hydration in the minor groove, with this sequence-specific water network involving about 8 additional water molecules from the second and, perhaps, even the third hydration layers. (c) Repetitive GC and IC sequences within a duplex are solvated beyond general effects by a "patch of hydration" in the major groove, with this water network involving about 13 additional water molecules from the second and, perhaps, even the third hydration layers. (d) Random sequence, polymeric DNA duplexes, which statistically lack extended regions of repetitive AT, GC, or IC sequences, do not experience such specific enhancements of hydration. Consequently, consistent with our previous observations (T. V. Chalikian, A. P. Sarvazyan, G. E. Plum, and K. J. Breslauer, Biochemistry, 1994, Vol. 33, pp. 2394-2401), duplexes with approximately 50% AT content exhibit the weakest hydration, while an increase or decrease from this AT content causes enhancement of hydration, either due to stronger hydration of the minor groove (an increase in AT content) or due to stronger hydration of the major groove (an increase in GC content). (e) In dilute aqueous solutions, a B-DNA duplex is more hydrated than an A-DNA duplex, a volumetric-based conclusion that is in agreement with previous results obtained on crystals, fibers, and DNA solutions in organic solvent-water mixtures. (f) the A-like, RNA duplex poly(rA)poly(rU) and the structurally similar A-like, hybrid duplex poly(rA)poly(dT), exhibit similar hydration properties, while the structurally distinct A-like, hybrid duplex poly(rA)poly(dT) and non-A-like, hybrid duplex poly(dA)poly(rU) exhibit differential hydration properties, consistent with structural features dictating hydration characteristics. We discuss how volumetric characterizations, in conjunction with structural studies, can be used to describe, define, and resolve the general and sequence/conformation-specific hydration properties of nucleic acid duplexes., (Copyright 1999 John Wiley & Sons, Inc.)

A method is offered for obtaining minimum energy configurations of DNA minicircles constrained by one or more DNA-binding proteins. The minicircles are modeled as elastic rods, while the presence of bound protein is implied by rigidly fixing portions of these chains. The configurations of the geometrically constrained circular rods are sampled stochastically and optimized according to a simple elastic energy model of nicked DNA. The shapes of the minimum energy structures identified after a simulated annealing process are analyzed in terms of relative protein orientation and writhing number. The procedure is applied to minicircles 500 base pairs in length, bound to two evenly spaced DNA-wrapping proteins. The presence of histone octamers is suggested by rigidly fixing the two protein-bound portions of each minicircle as small superhelices similar in dimension to nucleosomal DNA. The folded minimum energy forms of sample chains with different degrees of protein wrapping are noteworthy in themselves in that they offer a new resolution to the well-known minichromosome linking number paradox and point to future minicircle simulations of possible import.

A combination of detailed energy minimization and molecular dynamics studies of closed circular DNA offers here new information that may be relevant to the dynamics of short DNA chains and/or low superhelical densities. We find a complex dependence of supercoiled DNA energies and geometries on the linking number difference delta Lk as physiological superhelical densities (magnitude of sigma approximately 0.06) are approached. The energy minimization results confirm and extend predictions of classical elasticity theory for the equilibria of elastic rods. The molecular dynamics results suggest how these findings may affect the dynamics of supercoiled DNA. The minimization reveals sudden higher order configurational transitions in addition to the well-known catastrophic buckling from the circle to the figure-8. The competition among the bending, twisting, and self-contact forces leads to different families of supercoiled forms. Some of those families begin with configurations of near-zero twist. This offers the intriguing possibility that nicked DNA may relax to low-twist forms other than the circle, as generally assumed. Furthermore, for certain values of delta Lk, more than one interwound DNA minimum exists. The writhing number as a function of delta Lk is discontinuous in some ranges; it exhibits pronounced jumps as delta Lk is increased from zero, and it appears to level off to a characteristic slope only at higher values of delta Lk. These findings suggest that supercoiled DNA may undergo systematic rapid interconversions between different minima that are both close in energy and geometry. Our molecular dynamics simulations reveal such transitional behavior. We observe the macroscopic bending and twisting fluctuations of interwound forms about the global helix axis as well as the end-over-end tumbling of the DNA as a rigid body. The overall mobility can be related to magnitude of sigma and to the bending, twisting, and van der Waals energy fluctuations. The general character of molecular motions is thus determined by the types of energy minima found at a given delta Lk. Different time scales may be attributed to each type of motion: The overall chain folding occurs on a time scale almost an order of magnitude faster than the end-over-end tumbling. The local bending and twisting of individual chain residues occur at an even faster rate, which in turn correspond to several cycles of local variations for each large-scale bending and straightening motion of the DNA.

The base-to-base virtual bond treatment of nucleic acids used in statistical mechanical calculations of polynucleotide chain properties has been refined by incorporating the six parameters that relate the positions and orientations of sequential rigid bodies. The scheme allows for the sequence-dependent bending, twisting, and displacement of base pairs as well as for asymmetry in the angular and translational fluctuations of individual residues. Expressions are developed for the generator matrices required for the computation, as a function of chain length, of various parameters measuring the overall mean extension and shape of the DNA. Quantities of interest include the end-to-end vector r, the square of the end-to-end distance r2, the square radius of gyration s2, the center-of-gravity vector g, the second moments of inertia Sx2, and the higher moments of r and g. The matrix expressions introduced in the 1960s by Flory and co-workers for the determination of configuration-dependent polymer chain averages are decomposed into their translational and orientational contributions so that the methods can be extended to the rigid body analysis of chemical moieties. The new expressions permit, for the first time, examination of the effects of sequence-dependent translations, such as the lateral sliding of residues in A- and B-helices and the vertical opening of base pairs in drug-DNA complexes, on the average extension and shape of the long flexible double helix. The approach is in the following paper using conformational energy estimates of the base sequence-dependent flexibility of successive B-DNA base pairs.

We have used the elementary generator matrices outlined in the preceding paper to examine the conformational plasticity of the nucleic acid double helix. Here we investigate kinked DNA structures made up of alternating B- and A-type helices and intrinsically curved duplexes perturbed by the intercalation of ligands. We model the B-to-A transition by the lateral translation of adjacent base pairs, and the intercalation of ligands by the vertical displacement of neighboring residues. We report a complete set of average configuration-dependent parameters, ranging from scalars (i.e., persistence lengths) to first- and second-order tensor parameters (i.e., average second moments of inertia), as well as approximations of the associated spatial distributions of the DNA and their angular correlations. The average structures of short chains (of lengths less than 100 base pairs) with local kinks or intrinsically curved sequences are essentially rigid rods. At the smallest chain lengths (10 base pairs), the kinked and curved chains exhibit similar average properties, although they are structurally perturbed compared to the standard B-DNA duplex. In contrast, at lengths of 200 base pairs, the curved and kinked chains are more compact on average and are located in a different space from the standard B- or A-DNA helix. While A-DNA is shorter and thicker than B-DNA in x-ray models, the long flexible A-DNA helix is thinner and more extended on average than its B-DNA counterpart because of more limited fluctuations in local structure. Curved polymers of 50 base pairs or longer also show significantly greater asymmetry than other DNAs (in terms of the distribution of base pairs with respect to the center of gravity of the chain). The intercalation of drugs in the curved DNA straightens and extends the smoothly deformed template. The dimensions of the average ellipsoidal boundaries defining the configurations of the intercalated polymers are roughly double those of the intrinsically curved chain. The altered proportions and orientations of these density functions reflect the changing shape and flexibility of the double helix. The calculations shed new light on the possible structural role of short A-DNA fragments in long B-type duplexes and also offer a model for understanding how GC-specific intercalative ligands can straighten naturally curved DNA. The mechanism is not immediately obvious from current models of DNA curvature, which attribute the bending of the chain to a perturbed structure in repeating tracts of A.T base pairs.

Elasticity theory of naturally curved rods is employed to study the effects of intrinsic curvature on the properties of the equilibrium conformations of supercoiled DNA. The results stand in sharp contrast to those obtained when the molecule is viewed as being straight in its relaxed form. Starting from very fundamental principles of the theory, we show that the torsion of an open segment with a curved duplex axis can vary when the temperature, and along with it, the intrinsic twist is changed. Conversely, an imposed helicity, such as might be associated with binding to a histone, can change the intrinsic twist. It is also shown that another consequence of the presence of naturally curved sequences is that the twist density will, in general, vary with position along the chain in all equilibrium states. Then portions of the molecule will be more or less susceptible to interaction with other agents sensitive to such a variation. Finally, some closed equilibrium global structures uniquely associated with intrinsic curvature are discussed.

A Monte Carlo method has been developed for generating the conformations of short single-stranded DNAs from arbitrary starting states. The chain conformers are constructed from energetically favorable arrangements of the constituent mononucleotides. Minimum energy states of individual dinucleotide monophosphate molecules are identified using a torsion angle minimizer. The glycosyl and acyclic backbone torsions of the dimers are allowed to vary, while the sugar rings are held fixed in one of the two preferred puckered forms. A total of 108 conformationally distinct states per dimer are considered in this first stage of minimization. The torsion angles within 5 kcal/mole of the global minimum in the resulting optimized states are then allowed to vary by +/- 10 degrees in an effort to estimate the breadth of the different local minima. The energies of a total of 2187 (3(7)) angle combinations are examined per local conformational minimum. Finally, the energies of all dinucleotide conformers are scaled so that the populations of differently puckered sugar rings in the theoretical sample match those found in nmr solution studies. This last step is necessitated by limitations in the theoretical methods to predict DNA sugar puckering accurately. The conformer populations of the individual acyclic torsion angles in the composite dimer ensembles are found to be in good agreement with the distributions of backbone conformations deduced from nmr coupling constants and the frequencies of glycosyl conformations in x-ray crystal structures, suggesting that the low energy states are reasonable. The low energy dimer forms (consisting of 150-325 conformational states per dimer step) are next used as variables in a Monte Carlo algorithm, which generates the conformations of single-stranded d(CXnG) chains, where X = A, T and n = 3, 4, 5. The oligonucleotides are built sequentially from the 5' end of the chain using random numbers to select the conformations of overlapping dimer units. The simulations are very fast, involving a total of 10(6) conformations per chain sequence. The potential errors in the buildup procedure are minimized by taking advantage of known rotational interdependences in the sugar-phosphate backbone. The distributions of oligonucleotide conformations are examined in terms of the magnitudes, positions, and orientations of the end-to-end vectors of the chains. The differences in overall flexibility and extension of the oligomers are discussed in terms of the conformations of the constituent dinucleotide steps, while the general methodology is discussed and compared with other nucleic acid model building techniques.

We examine the salt dependence of the B-Z transition in DNA by means of the counterion condensation theory adapted to structurally realistic coordinates of the phosphate groups. The ionic contribution to the free energy difference delta G is computed for both the ZI and ZII conformations over broad ranges of NaCl and MgCl2 concentrations and polymer lengths. For the solvent we employ both a constant-dielectric model (dielectric constant set to 78.3) and a dielectric saturation model (distance-dependent dielectric constant). Where comparison can be made, the results for the constant-dielectric model are similar to those obtained by other workers for the same model but with different computational methods. The existence of a low-salt transition, and its location when it does occur, depends strongly on the DNA length and on the dielectric model. The behavior of ZI and ZII are qualitatively similar throughout the entire salt range for the constant-dielectric model, but qualitatively different if dielectric saturation is simulated, as we think is necessary for a realistic description. The ionic delta G, in the presence of dielectric saturation, bears comparison with the high-salt trend of the measured total delta G if "Z-DNA" is predominantly ZI, but not if it is predominantly ZII.

According to counterion condensation theory, one of the contributions to the polyelectrolyte free energy is a pairwise sum of Debye-Hückel potentials between polymer charges that are reduced by condensed counterions. When the polyion model is taken as an infinitely long and uniformly spaced line of charges, a simple closed expression for the summation, combined with entropy-derived mixing contributions, leads to the central result of the theory, a condensed fraction of counterions dependent only on the linear charge density of the polyion and the valence of the counterion, stable against increases of salt up to concentrations in excess of 0.1 M. Here we evaluate the sum numerically for B-DNA models other than the infinite line of B-DNA charges. For a finite-length line there are end effects at low salt. The condensation limit is reached as a flat plateau by increasing the salt concentration. At a fixed salt concentration the condensation limit is reached by increasing the length of the line. At moderate salt even very short B-DNA line-model oligomers have condensed fractions not far from the infinite polymer limit. For a long double-helical array with charge coordinates at the phosphates of B-DNA, the limiting condensed fraction appears to be approached at low salt. In contrast to the results for the line of charges, however, the computed condensed fraction varies strongly with salt in the range of experimentally typical concentrations. Salt invariance is restored, in agreement with both the line model and experimental data, when dielectric saturation is considered by means of a distance-dependent dielectric function. For sufficiently long B-DNA line and helical models, as typical salt concentrations, the counterion binding fraction approaches the polymer limit as a linear function of 1/P, where P is the number of phosphate groups of B-DNA.

The pharmacokinetics of sulfamethazine (SMZ) was studied in 44 healthy volunteers. After an oral dose of 10 mg/kg we found the usual distribution of rapid and slow acetylators. The differences in the t1/2s between the phenotypes were not due to variations in renal clearance. Repeat studies at this dose showed very little intraindividual variation in t1/2s. After a 40 mg/kg dose there was an increase in the t1/2 of SMZ in every subject. The implications of these results on methods of phenotyping are discussed.

A method is offered to model the complex trajectories of closed circular DNA supercoils and knots. The trajectories are approximated by polygons and analytical expressions of the curves are generated from the polygons with B-spline functions. The resulting curves are used to evaluate the writhe and elastic energy of a series of interrelated supercoils, and to generate detailed atomic models of the deformed double helix.

If bone biopsies are performed after the sites are located under radionuclide guidance, the chances of sampling the pathologic tissue are greatly improved. After bone scanning, two to three areas for biopsy are chosen and locating is done with 0.05 ml of Tc-99m in a tuberculin syringe under an Anger camera. Methylene blue-xylocaine is injected into the skin up to the periosteum for marking the site of the biopsy. Using this technique, 27 biopsies in 20 patients have been performed (four were open and 23 were closed needle). Pathology was found in 25 biopsies (92.6%), cancer in 22 (81.5%), benign lesions in three (11.1%), and normal bone in two (7.4%).