Search

Skin wounds often form scar tissue during healing. Early intervention with tissue-engineered materials and cell therapies may promote scar-free healing. Exosomes and extracellular vesicles (EV) secreted by mesenchymal stromal cells (MSC) are believed to have high regenerative capacity. EV bioactivity is preserved after lyophilization and storage to enable use in remote and typically resource-constrained environments. We developed a bioprinted bandage containing reconstituted EVs that can be fabricated at the point-of-need. An alginate/carboxymethyl cellulose (CMC) biomaterial ink was prepared, and printability and mechanical properties were assessed with rheology and compression testing. Three-dimensional printed constructs were evaluated for Young's modulus relative to infill density and crosslinking to yield material with stiffness suitable for use as a wound dressing. We purified EVs from human MSC-conditioned media and characterized them with nanoparticle tracking analysis and mass spectroscopy, which gave a peak size of 118 nm and identification of known EV proteins. Fluorescently labeled EVs were mixed to form bio-ink and bioprinted to characterize EV release. EV bandages were bioprinted on both a commercial laboratory bioprinter and a custom ruggedized 3D printer with bioprinting capabilities, and lyophilized EVs, biomaterial ink, and thermoplastic filament were deployed to an austere Arctic environment and bioprinted. This work demonstrates that EVs can be bioprinted with an alginate/CMC hydrogel and released over time when in contact with a skin-like substitute. The technology is suitable for operational medical applications, notably in resource-limited locations, including large-scale natural disasters, humanitarian crises, and combat zones.

A key step of hepatitis B virus (HBV) replication is the selective packaging of pregenomic RNA (pgRNA) by core protein (Cp) dimers, forming a nucleocapsid where the reverse transcriptional viral DNA replication takes place. One approach in the development of new anti-HBV drugs is to disrupt the assembly of HBV nucleocapsids by misdirecting Cp dimers to assemble morphologically normal capsids devoid of pgRNA. In this study, we built upon our previous discovery of benzamide-derived HBV capsid assembly modulators by exploring fused bicyclic scaffolds with an exocyclic amide that is β, γ to the fused ring, and identified 1,2,3,4-tetrahydroquinoxaline derived phenyl ureas as a novel scaffold. Structure-activity relationship studies showed that a favorable hydrophobic substitution can be tolerated at the 2-position of the 1,2,3,4-tetrahydroquinoxaline core, and the resulting compound 88 demonstrated comparable or improved antiviral potencies in mouse and human hepatocyte-derived HBV-replicating cell lines compared to our previously reported benzamide compound, 38017 (8). In addition, a novel bis-urea series based on 1,2,3,4-tetrahydroquinoxaline was also found to inhibit HBV DNA replication with sub-micromolar EC 50 values. The mode of action of these compounds is consistent with specific inhibition of pgRNA encapsidation into nucleocapsids in hepatocytes., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Masson SAS. All rights reserved.)

The analysis of thermally labile and high-boiling point compounds by gas chromatography (GC) can be a challenge. One technique to overcome these challenges is low-pressure GC, which uses the vacuum produced from the mass spectrometer and wide-bore columns to elute compounds at significantly lower temperatures. While GC-MS is a powerful technique, comprehensive two-dimensional gas chromatography (GC × GC), allows for resolution of compounds that would typically coelute using GC. In this study, a pesticide standard mixture (8270 MegaMix Standard) was analyzed using a conventional GC × GC-TOFMS configuration (0.25 mm inner diameter (i.d.) to a 0.18 mm i.d. column) and low-pressure GC × GC-TOFMS configuration (0.53 mm i.d. to a 0.53 mm i.d. column). Elution temperatures, sensitivity, and peak capacity were investigated for both configurations. Compounds eluted an average of 30 °C less on the low-pressure GC × GC-TOFMS configuration compared to the conventional GC × GC-TOFMS configuration. Moreover, the compounds were separated in ∼13 min on the low-pressure GC × GC-TOFMS as opposed to 33 min for conventional GC × GC-TOFMS. However, due to the wide-bore columns and faster runtimes the low-pressure GC × GC-TOFMS had a lower, β corrected 2D peak capacity, n c,β,2D , of 1260 while the conventional GC × GC-TOFMS was 3588. Interestingly, both configurations yielded a similar peak capacity production of 93 peaks/min and 107 peaks/min for low-pressure and conventional GC × GC-TOFMS, respectively. A "real world" sample of diesel fuel was tested on the low-pressure and conventional GC × GC-TOFMS configurations and similar results were obtained compared to the pesticide standard mix except the peak capacity production of the low-pressure GC × GC-TOFMS configuration was higher than that of the conventional GC × GC-TOFMS method., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier B.V. All rights reserved.)

Synthesizing solids in molten fluxes enables the rapid diffusion of soluble species at temperatures lower than in solid-state reactions, leading to crystal formation of kinetically stable compounds. In this study, we demonstrate the effectiveness of mixed hydroxide and halide fluxes in synthesizing complex Sr/Ag/Se in mixed LiOH/LiCl. We have accessed a series of two-dimensional Sr(Ag 1-x Li x ) 2 Se 2 layered phases. With increased LiOH/LiCl ratio or reaction temperature, Li partially substituted Ag to form solid solutions of Sr(Ag 1-x Li x ) 2 Se 2 with x up to 0.45. In addition, a new type of intergrowth compound [Sr 3 Se 2 ][(Ag 1-x Li x ) 2 Se 2 ] was synthesized upon further reaction of Sr(Ag 1-x Li x ) 2 Se 2 with SrSe. Both Sr(Ag 1-x Li x ) 2 Se 2 and [Sr 3 Se 2 ][(Ag 1-x Li x ) 2 Se 2 ] exhibit a direct band gap, which increases with increasing Li substitution (x). Therefore, the band gap of Sr(Ag 1-x Li x ) 2 Se 2 can be precisely tuned via fine-tuning x that is controlled by only the flux ratio and temperature., (© 2023 Wiley-VCH GmbH.)

This paper investigates the electrochemical behavior of p -aminophenol (PAP) on commercially available carbon screen-printed electrodes (CSPEs) and gold screen-printed electrodes (GSPEs) at neutral and basic pHs for the development of inexpensive immunoassays. The electrochemical oxidative signal from PAP results from its adsorption to the electrode. The formation of self-assembled monolayers on gold electrodes prevented PAP adsorption but also reduced its oxidative current, confirming that adsorption increases signal production. On bare electrodes, PAP adsorption results in oxidative current variability depending on the electroactive surface area of the screen-printed electrode. This variability could not be remedied by cleaning and reusing the same GSPE. Decreasing the PAP concentration to 3.8 μM greatly improved the consistency of the measurements, suggesting that the adsorption of PAP is concentration-dependent. Multiple PAP oxidations on the same electrode caused polymerization, limiting PAP in continuous monitoring applications. Infrared and Raman spectroscopy allow the distinction between adsorbed PAP and electropolymerized PAP on the surface of a gold wafer. The results from this study suggest that the use of PAP production in immunoassays with SPEs must be fine-tuned, and electrodes must be cleaned or disposed of between measurements.

Introduction: There are no reports in the literature describing risk factors for failure of nonoperative treatment of patients with posterior labral tears on magnetic resonance imaging (MRI). The purpose of this study is to identify risk factors for failure of nonoperative treatment in patients with an isolated posterior glenoid labral tear identified on MRI only. Patients with posterior labral tears on MRI who fail to improve with nonoperative treatment likely share a constellation of clinical history, physical exam, and radiographic findings., Methods: One hundred and fifty-nine active duty military service members under the age of 40 with a posterior labral tear seen on MRI and who were clinically evaluated by a musculoskeletal trained physician were identified. We retrospectively evaluated their records ensuring a minimum of 2 years follow-up after MRI to identify surgical intervention for the posterior labral tear during this time period. Patients were stratified into two groups, those treated with any combination of nonoperative modalities and those treated with posterior labral repair surgery during the 2 years after the MRI. The electronic medical records were reviewed for clinical presentation and physical exam results. We measured multiple radiographic parameters, including glenoid version, size of the tear, and bone loss on MRI. Qualitative and quantitative data were compared between groups using Fisher's exact test and Student's t-test, respectively. This study was conducted under institutional review board approval., Results: Of the 157 patients' shoulders in our study, 52% (n = 82) of patients with posterior labral tears underwent nonoperative treatment while 48% (n = 75) underwent surgery. The significant risk factors associated with surgery were a history of a specific injury, primary presenting complaint of instability, patient reported history of subluxation, inability to trust their shoulder with overhead activity, decreased strength with weight lifting, positive posterior load/shift exam, positive anterior apprehension, increased osseous glenoid retroversion, increased humeral head subluxation ratio, and anterior labral height (P < 0.05). Patients with a chief complaint of pain were much more likely to succeed with nonoperative treatment while those with instability underwent surgery more often. Ten (12.5%) of the surgical procedures included an anterior and posterior labral repair/stabilization procedure., Conclusion: Patients with an MRI confirmed posterior labral tear, which present with subjective complaints and physical exam maneuvers consistent with instability, appear less likely to be treated nonoperatively. Increased glenoid retroversion and posterior humeral head subluxation may also predispose patients toward surgical treatment. Additionally, posterior labral tears may extend into the anterior labrum more frequently than is recognized on MRI., (Published by Oxford University Press on behalf of the Association of Military Surgeons of the United States 2020. This work is written by US Government employees and is in the public domain in the US.)

We demonstrate that differential scanning calorimetry (DSC) can be used to yield high-resolution melting profiles for DNA plasmids that agree in all major features with the corresponding plasmid melting profiles derived using more traditional optical techniques. We further demonstrate that by combining information derived from both calorimetric and optical melting profiles one can glean insights that are unavailable from either melting curve alone. By using both optical and calorimetric observables, we show how one can resolve, identify, and measure the thermodynamic properties of particular sequences/domains of interest within a plasmid. We also show that complementary DSC and optical melting studies on plasmids with and without specifically designed inserts can provide fundamental advantages over the corresponding melting studies on other model system constructs for thermodynamically characterizing nucleic acid sequences/structures., (Copyright 1999 John Wiley & Sons, Inc.)

Motivation: MELTSIM is a windows-based statistical mechanical program for simulating melting curves of DNAs of known sequence and genomic dimensions under different conditions of ionic strength with great accuracy. The program is useful for mapping variations of base compositions of sequences, conducting studies of denaturation, establishing appropriate conditions for hybridization and renaturation, determinations of sequence complexity, and sequence divergence., Results: Good agreement is achieved between experimental and calculated melting curves of plasmid, bacterial, yeast and human DNAs. Denaturation maps that accompany the calculated curves indicate non-coding regions have a significantly lower (G+C) composition than coding regions in all species examined. Curves of partially sequenced human DNA suggest the current database may be heavily biased with coding regions, and excluding large (A+T)-rich elements., Availability: MELTSIM 1.0 is available at: //www.uml.edu/Dept/Chem/UMLBIC/Apps/MEL TSIM/MELTSIM-1.0-Win/meltsim. zip. Melting curve plots in this paper were made with GNUPLOT 3.5, available at: http://www.cs.dartmouth.edu/gnuplot&#95;inf o.html Contact : blake@maine.maine.edu;

The slime mold, Dictyostelium discoideum, possesses an (A+T) rich eukaryotic genome that is being sequenced in the Human Genome Project. High resolution melting curves of isolated total and fractionated nuclear D. discoideum DNA(AX3 strain) were determined experimentally and are compared to melting curves calculated from GENBANK sequences (1.59% of genome) by the statistical thermodynamics program MELTSIM (1), parameterized for long DNA sequences (2,3). The lower and upper temperature limits of calculated melting agree well with the observed melting of total DNA. The experimental curve is unusual in that it contains a number of sharp peaks. MELTSIM allowed us to calculate positional denaturation maps of D. discoideum GENBANK sequence documents containing the 26S, 5.8S and 17S rDNA gene sequences, a major satellite DNA and repetitive sequence family present in 100-200 copies/nucleus. These denaturation maps contain subtransitions that correspond with a number of the experimentally observed peaks, some of which we show to correspond with rDNA gene enriched CsCl gradient fractions of D. discoideum DNA. MELTSIM calculated curves of coding, intron and flanking sequences indicate that both intron and flanking sequences are extremely (A+T) rich and account for most of the low temperature melting. There is no temperature overlap between thermal stabilities of these sequence domains and those of coding DNA. The latter must satisfy triplet codon constraints of higher (G+C) content. These large stability property differences enable a denaturation mapping feature of MELTSIM to clearly distinguish exon positions from those of introns and flanking DNA in long D. discoideum gene containing sequences.

Tij and Delta Hij for stacking of pair i upon j in DNA have been obtained over the range 0.034-0.114 M Na+from high-resolution melting curves of well-behaved synthetic tandemly repeating inserts in recombinant pN/MCS plasmids. Results are consistent with neighbor-pair thermodynamic additivity, where the stability constant, sij , for different domains of length N depend quantitatively on the product of stability constants for each individual pair in domains, sijN . Unit transition enthalpies with average errors less than +/-5%, were determined by analysis of two-state equilibria associated with the melting of internal domains and verified from variations of Tij with [Na+]. Enthalpies increase with Tij , in close agreement with the empirical function: Delta Hij = 52.78@ Tij - 9489, and in parallel with a smaller increase in Delta Sij . Delta Hij and Delta Sij are in good agreement with the results of an extensive compilation of published Delta Hcal and Delta Scal for synthetic and natural DNAs. Neighbor-pair additivity was also observed for (dA@dT)-tracts at melting temperatures; no evidence could be detected of the familiar and unusual structural features that characterize tracts at lower temperatures. The energetic effects of loops were determined from the melting behavior of repeating inserts installed between (G+C)-rich barrier domains in the pN/MCS plasmids. A unique set of values for the cooperativity, loop exponent and stiffness parameters were found applicable to internal domains of all sizes and sequences. Statistical mechanical curves calculated with values of Tij([Na+]) , Delta Hij and these loop parameters are in good agreement with observation.

High-resolution derivative melting was used to obtain detailed distributions of local (G + C) contents in a number of ruminant DNAs. Profiles over low (G + C) regions [20-36% (G + C)] are congruent for all ruminants. This region represents 45-50% of the nuclear DNA content and primarily contains intergenic and intron sequences. The high (G + C) region, where most coding sequences are found [38-68% (G + C)], is marked by satellite bands denoting the presence of transcriptionally inert, tandemly repetitive sequence families. These bands can be analyzed for the abundance, base composition, and sequence divergence of satellite families with relatively high precision. Band patterns are unique to each species; even closely related species can be readily distinguished by their base distribution profiles. Variations in nuclear DNA contents in ruminants, determined by flow cytometry, are primarily due to variations in abundances of these repetitive sequence families. Thus, A. alces (moose) is found to have 8.85 +/- 0.2 pg DNA/cell, 25% more than the average in ruminants, while the base distribution curve indicates the presence of an unusually abundant satellite of 52.6% (G + C). The size (1 kb) and sequence of this satellite corresponds to satellite-I of other cervids, and in consequence it is designated Alces-I. The sequence of a cloned repeat of Alces-I has a length of 968 bp, a (G + C) content of 52.6%, and contributes 35%, or almost 3 million copies to the nuclear DNA, exceeding by approximately 300% the average array size of this repeat family in related cervids. In situ hybridization indicates the repeat is distributed throughout centromeric regions of all 62 acrocentric autosomes. Alces-I has much greater-than-expected numbers of GG, GA, and AG and far fewer numbers of TA and CG duplets, characteristics of all tandem repeats. The sequence is judged to be orthologous with satellite-I sequences from Rangifer tarandus (caribou), Capreolus capreolus (roe deer), Muntiacus muntjac (Chinese muntjac) and Muntiacus reevesi (Indian muntjac), as well as Antilocapra americana (pronghorn), and the bovids Bos taurus and Ovis aries. A tentative tree for the five cervids is in excellent agreement with one proposed on the basis of morphological characteristics. Differences from a consensus sequence indicate transversions exceed transitions by almost twofold, suggesting that substitutions occur randomly, or nearly so.

Formamide lowers melting temperatures (Tm) of DNAs linearly by 2.4-2.9 degrees C/mole of formamide (C(F)) depending on the (G+C) composition, helix conformation and state of hydration. The inherent cooperativity of melting is unaffected by the denaturant. dTm/dC(F)for 11 plasmid domains of 0.23 < (G+C)<0.71 generally fit to a linear dependence on (G+C)-content, which, however, is consistent with a (G+C)-independent alteration in the apparent equilibrium constant for thermally induced helix <--> coil transitions. Results indicate that formamide has a destabilizing effect on the helical state, and that sequence-dependent variations in hydration patterns are primarily responsible for small variations in sensitivity to the denaturant. The average unit transition enthalpy delta H(m)[see text for complete expression], exhibits a biphasic dependence on formamide concentration. The initial drop of -0.8 kcal/mol bp at low formamide concentrations is attributable to a delta delta H(m)[see text for complete expression], for exchange of solvent in the vicinity of the helix: displacement by formamide of weakly bound hydrate or counterion. The phenomenological effects are equivalent to lowering the bulk counterion concentration. Poly(dA.dT) exhibits a much lower sensitivity to formamide, due to the specific pattern of tightly bound, immobilized water bridges that buttress the helix from within the narrow minor groove. Tracts of three (A.T)-pairs behave normally, but tracts of six exhibit the same level of reduced sensitivity as the polymer, suggesting a conformational shift as tracts are elongated beyond some critical length [McCarthy J.G. and Rich,A. (1991) Nucleic Acids Res. 19, 3421-3429].

We report the results of a theoretical study, combining the results of sequence analysis and integral equation structural methods for nucleic acids in aqueous solutions, on the effects of nearest neighbors on the (T.G) mispair in solution, for 12 nearest neighbor contexts. Attempts have been made to classify the structural and energetic effects of the 5' and 3' neighbors with respect to the observed spontaneous mutation rates in vertebrates. It is found that 5' nearest neighbor is probably the most critical structural factor in facilitating or discouraging mutations. Local conformational states correlate with discrimination of bases to be excised in mispairs. Our study confirms the role of the flexibility of the DNA molecule in governing the rates of spontaneous mutations.

The (G + C) distribution and the presence and amounts of repetitive sequence families in the white-tailed deer (Odocoileus virginianus) have been examined. The distribution ranges from 20 to 70% (G+C) and shows four distinct repeat families. A 0.7-kb family, DII, corresponds to satellite II in domestic bovids--ox, sheep, and goat--and was singled out for detailed characterization. DII has a prototypic repeat of 67% (G + C), consists of 25,000 tandem copies, and contributes 1.7% to the genomic DNA. Sequencing and electrophoretic analysis indicate a repeat length of 691 bp. These characteristics are similar to those of the bovid satellite II families as well as to those of other cervids that we have examined. The intraspecific sequence divergence within this family has a variance of only 2.5 +/- 0.3%.

It has been shown that the frequency versus size distribution of A and T overlapping and non-overlapping homopolymer tracts of N > 5 in D. discoideum gene flanking and intron regions are significantly greater than in coding regions(1). In the present report, we demonstrate, that a spatial periodicity exists in long A and T tracts (N > 10) in long flanking sequences by scored alignments of those tracts (N > 10) with the nucleosomal repeat. A tract spacing was found at 185-190 bp that corresponds to a maximum alignment score. This is exactly the average spacing of D. discoideum nucleosomes determined experimentally. A majority of A and T tracts in flanking sequences are often spaced by short DNA stretches and the total length of adjacent A and T tracts plus the interrupting short DNA stretch corresponds closely to the average experimentally measured nucleosomal linker DNA size in D. discoideum-42 bp. These data suggest a model which has A and T runs of N > 10 bp in flanking DNA of D. discoideum organized in a regular phase with nonhomopolymer sequences along the DNA. This model has functional implications for A and T tracts, suggesting that they are found in nucleosomal linker DNA regions of chromatin during some necessary portion(s) of the life of the cell.

The pattern of 20,200 point substitutions in the 16 unique neighbor-pair environments has been determined from aligned gene/pseudogene sequences in the current database of human DNA sequences. Substitution rates, representing averages over those for different regions of the genome, are distributed over a 60-fold range with strong biases in particular neighbor-pair environments. The rates for substitutions involving the CG doublet are the most rapid overall, where changes of the C.G pair vary over a tenfold range depending on the type of substitution and the 5' neighbor-pair. In general, the rates are fastest in alternating purine-pyrimidine sequences and slowest in purine.pyrimidine tracts, suggesting that the frequencies of one or both key molecular misadventures that can occur during replication, dNTP misinsertion and transient misalignment, may be associated with structural alternations and flexibility of the backbone. By contrast, purine.pyrimidine tracts are less flexible, less prone to substitution, and therefore their proportions accumulate in sequences over time. Characteristic biases of the content and arrangement of oligonucleotide strings or tuples in all sequence elements, but particularly in non-coding regions, appear to be due to the pattern of different neighbor-dependent substitution rates. Computer simulations of numerous replicative cycles have been carried out with substitutions occurring on the same schedule found in this study for pseudogenes. Statistical analyses of tuple frequencies at periodic intervals during the simulation experiment indicate that sequences slowly change in lexical complexity toward a quasi-equilibrium state that corresponds to that for introns.

The results of this theoretical study combining sequence analysis and minimization with integral equation liquid structural methods indicate that the local sequence context of a T-G wobble mismatch influences the local conformation of the helix, and that conformational alterations are correlated with mutational activity. Studies on the mismatch in four different 5' and 3' neighbor contexts indicate that the nature of the 5' base to the thymine of the mispair is probably the single most critical factor in determining the structural features that facilitate or discourage mutations. When cytosine is the 5' neighbor, the helix adopts a mostly BII conformation, whereas a 5' guanine preserves the canonical BI. Structures that vary little from the BI structure on the incorporation of the mismatch have sequences that correspond to lower rates of transition, whereas those with mostly BII conformations, have sequences with high mutation rates. Subtle variations in stacking patterns around the mismatch precipitate a structural Domino-effect, with a variety of changes in conformation. The helix opens at the mismatch with increased roll angle and propeller twist, causing the thymine to migrate into the major groove and the guanine into the minor groove, exposing the heteroatomic groups to the solvent in the major and minor grooves, respectively, and allowing for some unusual hydrogen bonds. These alterations show a tentative correlation with mutation rates, implying that stacking and structure around the mismatch are important features in the discrimination by proofreading activities of canonical W-C and wobble mismatch base pairs during replication-repair. Variations in the C1'-C1' distances, high propeller twists, changes in the electrostatic complementarity leading to unusual hydrogen bonding patterns probably all correlate with detectability.

D. discoideum, the slime mold, is one of the most AT rich eukaryotic genomes known. In this paper we examine this organism's database for overlapping N-tuples of high frequency and find A and T tracts possess among the highest frequencies in flanking sequences but not in coding sequences. We examined both overlapping and non-overlapping frequencies of the A, T, G and C homopolymer tracts of 2 < N < 6. Overlapping (dG).(dC) and (dA).(dT) tracts occur at greater frequencies than expected, based on random occurrence. Long (dA).(dT) tracts of N > 10 occur at well above expected frequencies in flanking and intron regions, while (dG).(dC) tracts above N = 5 are rarely found. Some of the implications of these findings for tract origins in slip-strand replication and for chromatin structure are discussed.

The numbers and local sequence environments of the two types of substitution mutation plus additions and deletions have been obtained directly in this study from differences between a large number of extant primate gene and pseudogene sequences. A total of 3786 mutations were scored in regions where similarities between pseudogene and corresponding gene sequences is greater than or equal to 85%, comprising approximately 30% of the pseudogene database of 80,584 bp. The pattern of mutations obtained in this fashion is almost identical to that obtained by Li et al. (1984) using a slightly different, more direct approach and with a smaller database. When mutations were scored, the neighbor pairs on the 5' and 3' sides were also noted, leading to a large 16 x 12 matrix of transitions and transversions. Biases of varying magnitude are found in the rates of substitution of the same base pair in different local sequence environments. The overall order for the effect of the 5' neighbor on the rates of substitution mutation of a pyrimidine is A greater than C much greater than T greater than G, and G greater than A greater than T greater than C for the 3' neighbor; where these results represent the average of substitution rates for the complement purine with complement neighbors of bases ordered above. The order for the 3' neighbor is essentially the same for the two transitions and most of the four transversions as well; however, the order for the 5' neighbor is more variable. The overall rate for the C.G----T.A transition is not unusual, however the presence of a 3' neighboring G.C pair boosts the rate substantially, presumably due to specific cytosine methylation of the CG doublet in primate DNAs. The rate of the T.A----C.G transition is also well above average when the 3' neighbor is an A.T, and to a lesser extent a G.C, pair. The latter bias is typical in that it reflects the association of alternating pyrimidine-purine sequences with increasing mutation rates. The substitution of the pyrimidine in a 5'purine-pyrimidine-purine3' sequence generally occurs much faster than in a pyrimidine tract and points to the local conformation as a major determining factor of the substitution rate. An apparent inverse relationship is found between starting and product doublet frequencies of base pairs undergoing mutations with specific 3' neighbors, indicating that differences in intrinsic substitution rates of base pairs with specific neighbors are a key factor in producing the familiar biases of nearest-neighbor frequencies.

Variations in base mono- and dipoles result in variations in stacking energies for the 10 unique neighbor pairs in DNA. Stacking energies for pair M on N, expressed as TMN, were derived by matrix decomposition of a large set of linear algebraic expressions relating the measured Tm for subtransitions emanating from large polymeric DNAs, and the fractional neighbor frequencies, fMN, for the domains responsible for the transitions, Tm = sigma fMNTMN. Tm were determined for subtransitions that dissociate in approximately all-or-none fashion in high resolution melting profiles of partially deleted and recombinant forms of pBR322 DNA. Three different analytical maneuvers were undertaken to resolve subtransitions: site-specific cleavage of domains; deletion of domains; and addition of domains. Three dozen domains of widely divergent, quasi-random neighbor frequencies were identified and assigned, resulting in a unique set of values for TMN with standard deviation, sigma = +/- 0.23 degree C. The average difference between calculated and experimental Tm for domains is only +/- 0.17 degree C, indicating that the thermodynamic properties of these domains are not in any way unusual. Assuming delta S to be constant for all pairs, the corresponding delta HMN are found to have a precision of +/- 10 calories.mol-1 and an accuracy of +/- 606 calories.mol-1. TMN used to calculate melting curves by statistical mechanical analysis of sequences of the different plasmid specimens in this study were in quantitative agreement with observed curves for most sequences. These TMN differ significantly from those determined previously and also correlate poorly with values determined by quantum chemical analysis. Stabilities of neighbor pairs, expressed as the difference in free energy between that for a given pair (MN) and that for the average of like pairs (M, N), depend on the relationship of stacked purines and pyrimidines as follows. delta delta Gpu-py(-466 cal) greater than delta delta Gpu-pu(+52 cal) greater than delta delta Gpy-pu(+335 cal) Differences between experimental Tm and Tm calculated with TMN for the isolated neighbor pairs in the B-conformation are useful in the identification of altered structures and unusual modes of dissociation of helixes. A significantly higher Tm is observed for the highly biased repeated sequence synthetic helixes dA.dT, d(AGC).d(GCT), and d(GAT).d(ATC), reflecting auxiliary sources of stability such as bifurcated hydrogen bonds and/or altered structures for these helixes.

The Tm of internal loop-forming (dA.dT)N domains in pBR322 DNA has been measured over a tenfold range of [Na+]. The slopes SN = dTm/d log [Na+] are linear and decrease in magnitude with decreasing loop size N, signaling a reduction in Na+ released during the transition of these domains to the coil state. Values of SN decrease linearly with increasing N-1 in accordance with the expectation of a simple model for the occurrence of a gradient of long-range electrostatic forces at helix-coil boundaries, and extrapolate almost precisely to the value of S infinity observed for (dA.dT) infinity. These results indicate (1) less counterion is released per phosphate residue from the finite loop than from the infinite-sized loop, and (2) the difference in binding is constant for each boundary formed and independent of the size of the loop within the range examined: approximately 350 base pair (bp) greater than N greater than 71 bp. The slope of the dependence of SN on N-1 indicates the region of higher charge density at the boundary extends at least 18 A into the coil and probably 40-50 A before dropping to a value characteristic of the unperturbed coil. The free energy for excess counterion binding at boundaries can be expressed by -delta G/RT = 10.47 log[Na+] + 5.234 When the loop entropy function in a statistical mechanical algorithm for the dissociation of DNA is weighted by this quantity, calculated Tm are seen to vary by only +/- 0.09 degrees C from observed.

Oligonucleotide and codon frequencies have been determined in published sequences of E. coli DNA totaling 103,100bp with 18,459 reading frame trinucleotides; corresponding to 2.5% of the total genome. Dinucleotide frequencies are in excellent agreement with those determined by nearest neighbor chemical analysis, indicating the computer count of a limited sampling to be a good representation of the overall frequencies in total genomic DNA. The distinctive nonrandom codon pattern is found to be uniformly distributed and contributes to a distinctive nonrandom oligonucleotide pattern; enabling correlations between frequency levels to be extended beyond reading frame sequences. Correlation analysis indicates a surprisingly high degree of correlation everywhere in the genome. Coefficients of correlation between oligonucleotide frequencies overall and those in specific segments vary as follows: primary strands of individual coding sequences greater than 0.9 greater than lambda DNA greater than noncoding, non-RNA greater than phi X174 DNA greater than complementary strands greater than RNA genes congruent to 0.6 greater than transposon-insertion elements greater than T7DNA much greater than eukaryotic sequences congruent to 0. It is concluded that this high degree of oligonucleotide and codon correspondence in E. coli reflects the widespread distribution of remnants of an early and slowly changing codon pattern that has been continually dispersed by duplication-divergence processes, leading to the present genome.

Codon usage tables have been produced for E. coli, yeast, human, and mouse. The nonrandom employment of codons allows assignment of probability values to trinucleotides in any DNA sequence. These values represent the probability that a given trinucleotide is used as a codon in the organism from which the table is derived. For the graphical delineation of coding areas in DNA sequences, a probability is assigned to each trinucleotide equal to its frequency in the codon table. Averaging and smoothing procedures then greatly enhance the detectability of areas of high average codon probability and better represent the mean codon probability. These manipulations increase graphical clarity without altering the overall magnitude of probabilities. Averaging introduces an error of less than 0.5% between "raw" and smoothed data. This graphical delineation of coding sequences does not depend on the presence of punctuation, ribosomal binding sites, etc: moreover the delineation of introns and exons is also possible.

Fifty-three gene sequences from E. coli containing 18,288 reading frame triplets have been characterized according to the nature and level of average codon preference. The distribution of average preferences is bimodal, with approximately half the genes using an average of only 36 codons, and the remainder just 42 codons. There is a high correlation between the level of codon bias, the tRNA population and the abundance of protein product, indicating biased patterns are exploited by the cell for the production of widely different levels of gene product. This relationship is especially striking in genes involved in the production of components for transcription and translation. Overall, the genes for these processes generate some five-fold more protein than the average in the genome, and use about five fewer codons. The very high codon bias found in the RNA polymerase gene thus provides a simple, autogenous mechanism for the coordinate synthesis of these components and RNA polymerase. A surprisingly high level of codon probability is also found in triplets of the complement of coding sequences. This is apparently due to the evolutionary dispersion of coding sequences and/or the requirement for increased levels of secondary structure in messenger RNAs.

A microcomputer-controlled spectrophotometer is described for obtaining the base composition of melting domains in DNA from derivative melting curves. Values have been determined for the differential molar extinction coefficients for the A-T and G-C base pair at the three wavelengths most useful for spectral analysis of base composition, 260, 270 and 282 nm. The average RMS error for these values was 29 l(mol X cm)-1 for the melting of 14 DNA specimens ranging in base composition from 0-0.72 F(G + C). A precision of approximately 1% in base composition of domains is possible. Such analysis is useful for confirming or establishing assignments of domains to particular subtransitional features in high resolution melting curves.

Methods of computer analysis for the recurrence of symmetrical and repetitive elements in large numbers of DNA sequences are described, together with derivations of appropriate quantitative criteria for the evaluation of the statistical significance of these elements in DNAs of different base composition. Examples of some extraordinary variations in the occurrence of symmetrical and repetitive elements are provided, many of which are new. Special consideration is devoted to a determination of the statistical significance of a two-fold palindrome at the origin of replication. A computer search of 14 independently determined DNA sequences containing an origin of replication locus indicates each contains a large two-fold palindrome. The average length of this palindrome is 28 +/- 6 base pairs, of which 22 contribute to the palindromic symmetry. The probability of occurrence of such a palindrome is only 1/26000, while the probability of occurrence in all 14 different species is (1/26000).

The mean (G + C) composition (51.0%) and standard deviation (+/- 3.8%) of published DNA sequences accounting for 10% of the E. coli genome is in excellent agreement with the principal overall distribution determined by high resolution melting. While differences in base and neighbor characteristics are small and uniform throughout all regions of the genome, it is found that the (G + C) content of sequences varies in segmented fashion within boundaries corresponding to coding (53% G + C) and noncoding (46% G + C) regions; with variances in the latter being six-fold greater than in coding regions. The variance in different regions shows a strong negative dependence on (G + C) content of the region, reflecting the condition that A-T and G-C base pairs are preferred neighbors of A-T and C-G pairs, respectively; with the bias increasing with decreasing (G + C) content. Neighbor analysis indicates the most extreme positive biases occur in AA, TT, GC and CG throughout all regions, but particularly in noncoding regions. Extraordinary numbers of oligomeric strings of (A)n, etc., are the further consequence of this bias. These and other characteristics point to the existence of inherent biases in neighbor frequencies levied during replication or repair, and which reflect, in turn, neighbor influences during mutation. The bias in codon usage noted by Grantham and others is seen here as due, in part, to the adaptation of coding sequences to this microenvironment through selection among synonymous codons so as to preserve inherent neighbor biases.

Derivative melting profiles of DNA have been obtained directly by recording the difference in absorbance between two identical solutions maintained at a small constant temperature differential. This deltaA is monitored continuously with increasing temperature in a ratio recording spectrophotometer. Resolution of complex hyperfine structure in the profiles of small homogeneous viral DNAs appears to be significantly better than has been produced by various numerical methods of differentiation. In addition, a spectral method has been modified that permits easy analysis for DNA base composition from the ratio of derivative melting curves obtained at 282 and 260 nm. Eight bacterial and three vertebrate DNAs have been analyzed for total base composition from the product of the instantaneous base composition at small temperature intervals (0.05 degrees C) throughout the entire melting region and the integrated area of the 282 nm profile. The results are in excellent agreement with values determined by traditional methods.

In the present investigations the molecular unfolding of yeast tRNA(Gly) has been studied by a combination of nuclear magnetic resonance spectroscopy, melting techniques, and relaxation kinetics. From these studies the following pathway of unfolding was found. In a coupled melting transition the tertiary, the DHU, and the anticodon structure are disrupted. This is followed by the melting of the acceptor arm, while the T psi C arm, which only contains G-C pairs, melts out last. Interestingly, during the first melting transition a new structure not belonging to the original cloverleaf structure is formed. The thermodynamic and kinetic parameters of the melting transitions were determined and are discussed in relation to earlier work. The present nuclear magnetic resonance (NMR) experiments as well as earlier studies show that the ring current calculations based on the cloverleaf structure provide a good first-order interpretation of the NMR spectra of tRNA.