PubMed Health⌕ Search

Biomedical subjects

S C Harvey

Publications and source records attributed to S C Harvey.

At least 19 recordsLinked to original sources
PubMed1992

A comprehensive classification of nucleic acid structural families based on strand direction and base pairing.

We propose a classification of DNA structures formed from 1 to 4 strands, based only on relative strand directions, base to strand orientation and base pairing geometries. This classification and its associated notation enable all nucleic acids to be grouped into structural families and bring to light possible structures which have not yet been observed experimentally. It also helps in understanding transitions between families and can assist in the design of multistrand structures.

Base Composition↗
PubMed1992

Analyzing the normal mode dynamics of macromolecules by the component synthesis method.

This paper presents a general method for studying the harmonic dynamics of large biomolecules and molecular complexes. The performance and accuracy of the method applied to a number of molecules are also reported. The basic approach of the method is to divide a macromolecule into a number of smaller components. The local normal modes of the components are first calculated by treating individual components and the interactions between nearest neighboring components. The physical displacements of all atoms are then represented in the local normal mode space, in which a selected range of high-frequency local modes is neglected. The equation of motion of the molecule in the local normal mode space will then have a smaller dimension, and consequently the normal modes of the whole structure, particularly for large molecules, can be solved much more easily. The normal modes of two polypeptides--(Ala)6 and (Ala)12--and a double-helical DNA--d(ATATA).d(TATAT)--are analyzed with this method. Reductions on the dimensions of harmonic dynamic equations for these molecules have been made, with the fraction of the deleted high-frequency modes ranging from 1/2 to 5/6. The calculated low-frequency normal modes are found to be very accurate as compared to the exact solutions by standard procedure. The major advantage of the present approach on macromolecule harmonic dynamics is that the reduction on the dimensionality of the eigenvalue problems can be varied according to the size of molecules, so the method can be easily applied to large macromolecules with controlled accuracy.

DNA↗
PubMed1992

Teaching macromolecular modeling.

Training newcomers to the field of macromolecular modeling is as difficult as is training beginners in x-ray crystallography, nuclear magnetic resonance, or other methods in structural biology. In one or two lectures, the most that can be conveyed is a general sense of the relationship between modeling and other structural methods. If a full semester is available, then students can be taught how molecular structures are built, manipulated, refined, and analyzed on a computer. Here we describe a one-semester modeling course that combines lectures, discussions, and a laboratory using a commercial modeling package. In the laboratory, students carry out prescribed exercises that are coordinated to the lectures, and they complete a term project on a modeling problem of their choice. The goal is to give students an understanding of what kinds of problems can be attacked by molecular modeling methods and which problems are beyond the current capabilities of those methods.

Biochemistry↗
PubMed1992

Molecular modeling to predict the structural and biological effects of mutations in a highly conserved histone mRNA loop sequence.

The 3'-end of histone mRNAs contains a highly conserved sequence motif which is believed to form a 6 base pair stem and a 4 base loop. These sequences are involved in both the efficiency of 3'-end formation and stability of the mature histone mRNA. We have modeled four stem basepairs and the loop portion of this structure using the wildtype sequences and several mutant sequences. A structure for the wildtype stem-loop is proposed that is based on energy minimization using a representative wildtype sequence and comparison with structures obtained using naturally occurring mutations which do not alter loop function. A wildtype structure is proposed in which the top basepair of the stem is broken, forming a six base loop. Mutant sequences with altered bases in the loop and in the stem were also modeled. The effect of these mutations on the proposed wildtype structure is discussed and possible biological consequences considered.

Base Sequence↗
PubMed1991

Mean field stochastic boundary molecular dynamics simulation of a phospholipid in a membrane.

Computer simulations of phospholipid membranes have been carried out by using a combined approach of molecular and stochastic dynamics and a mean field based on the Marcelja model. First, the single-chain mean field simulations of Pastor et al. [(1988) J. Chem. Phys. 89, 1112-1127] were extended to a complete dipalmitoylphosphatidylcholine molecule; a 102-ns Langevin dynamics simulation is presented and compared with experiment. Subsequently, a hexagonally packed seven-lipid array was simulated with Langevin dynamics and a mean field at the boundary and with molecular dynamics (and no mean field) in the center. This hybrid method, mean field stochastic boundary molecular dynamics, reduces bias introduced by the mean field and eliminates the need for periodic boundary conditions. As a result, simulations extending to tens of nanoseconds may be carried out by using a relatively small number of molecules to model the membrane environment. Preliminary results of a 20-ns simulation are reported here. A wide range of motions, including overall reorientation with a nanosecond decay time, is observed in both simulations, and good agreement with NMR, IR, and neutron diffraction data is found.

Computer Simulation↗
PubMed1990

Nuclear magnetic resonance and molecular modeling studies on O-beta-D-galactopyranosyl-(1----4)-O-beta-D-xylopyranosyl-(1----0)-L-se rine, a carbohydrate-protein linkage region fragment from connective tissue proteoglycans.

The solution conformation of O-beta-D-galactopyranosyl-(1----4)-O-beta-D-xylopyranosyl-(1----0)-L-ser ine (GXS), a carbohydrate-protein linkage region fragment from connective tissue proteoglycans, was investigated by two-dimensional NMR spectroscopy and molecular modeling calculations. Specifically, the 1H and 13C resonances were assigned by 2D-COSY and by 1H-13C heteronuclear correlation spectroscopy methods. 2D-NOESY was used to generate distance constraints between the galactose and xylose and between the xylose and serine residues. The 1H vicinal coupling constants for the sugars and the serine were also determined. A general molecular modeling methodology suitable for complex carbohydrates was developed. This methodology employed molecular dynamics and energy minimization procedures together with the application of inter-residue spatial constraints across the linkages derived from 2D-NOESY. The first step in this methodology is the generation of a wide variety of starting conformations that span the (phi, psi) space for each linkage. In the present study, nine such conformations were constructed for each linkage using the torsion angles phi and psi corresponding to the gauche+, gauche-, and trans configurations across each of the two bonds constituting the linkage. These conformations were subjected to a combined molecular dynamics/energy minimization refinement using the NOESY derived constraints as pseudoenergy functions. Families of conformations for the whole molecule were then constructed from the structures derived for each linkage. Characterization of GXS using this methodology identified a single family of conformations that are consistent with the solution phase NMR data on this molecule.

Computer Graphics↗
PubMed1990

Electrical potential of transfer RNAs: codon-anticodon recognition.

Calculations of the electrostatic potentials were made around yeast elongator phenylalanine, aspartate tRNAs, and yeast initiator methionine tRNA in aqueous solution at physiological ionic strength. The calculations were carried out with a finite difference algorithm for solving the nonlinear Poisson-Boltzmann equation that incorporates the screening effects of the electrolyte, the exclusion of ions by the molecule, the molecular shape, and the different polarizabilities of the solvent and the tRNA. The initiator tRNA is surrounded by uniformly spaced contours of negative potential. The elongator tRNAs are also surrounded by a similar contour pattern except in the anticodon region where there is a pronounced "hole" in the potential surface. This hole is caused by an invagination of the potential contours, which also results in an increase in the local field strength. The effect of this hole is that the anticodon region in the elongator tRNAs is the least negative, or conversely the most positive, region of the molecule. This hole, which is not found when simple Coulombic potentials are used, is due both to the structure of the elongator tRNA anticodon loops and to the different polarizabilities of the solvent and tRNA. The existence of the potential hole in elongator tRNAs may account in part for their ability to associate with other negatively charged macromolecules, in particular mRNA. Moreover, it suggests that the anticodon loop of elongator tRNAs is the energetically most favorable point of approach of mRNA to tRNA.

Anticodon↗
PubMed1990

Prediction of the three-dimensional structure of Escherichia coli 30S ribosomal subunit: a molecular mechanics approach.

We introduce a computer-assisted procedure for folding large RNA chains into three-dimensional conformations consistent with their secondary structure and other known experimental constraints. The RNA chain is modeled using pseudoatoms at different levels of detail--from a single pseudoatom per helix to a single pseudoatom for each nucleotide. A stepwise procedure is used, starting with a simple representation of the macromolecule that is refined and then extrapolated into higher resolution for further refinement. The procedure is capable of folding different random-walk chains by using energy minimization, allowing generation of a range of conformations consistent with given experimental data. We use this procedure to generate several possible conformations of the 16S RNA in the 30S ribosomal subunit of Escherichia coli by using secondary structure and the neutron-scattering map of the 21 proteins in the small subunit. The RNA chain is modeled using a single pseudoatom per helix. RNA-RNA and RNA-protein crosslinks, reported in current literature, are included in our model. Footprinting data for different ribosomal proteins in the 16S RNA are also used. Several conformations of the 16S RNA are generated and compared to predict gross structural features of the 30S subunit as well as to identify regions of the 16S RNA that cannot be well-defined with current experimental data.

Escherichia coli↗
PubMed1989

Molecular mechanics model of supercoiled DNA.

We describe a pseudo-atomic model of supercoiled DNA. Each base-pair of the DNA is represented in the model by three particles placed in a plane. The particle triplets are stacked to model stacked base-pairs in double-helical DNA, and closed circular conformations are generated to investigate supercoiling. This model is less detailed than all-atom models, which are too computationally demanding to be used to study supercoiling. On the other hand, this model contains details at the base-pair level and is therefore more elaborate than elastomechanical models. A potential energy function is written in terms of a set of internal co-ordinates defined to resemble a limited number of helical parameters. The modeled helical parameters, helical twist, base-roll, tilt and rise, are the most important parameters of the global shape of DNA. Experimentally measured mechanical properties of DNA are used to define the forces holding the particles together. We then use a procedure incorporating energy minimization and molecular dynamics to locate low energy conformations of the model DNA. The model was found to behave very much like rubber-tubing and elastomechanical models. The conformations and the effects of supercoiling pressure (a number proportional to the degree to which the total twist of the DNA has been altered from its natural value) on these conformations are all very similar to those observed in the latter two models. We also used this model to examine the effects of supercoiling pressure, base-sequence and mechanical properties on the conformations and energies of five sequences. The sequences studied include models of naturally straight DNA and DNA with static or natural bends.

Base Sequence↗
PubMed1989

Differential scanning calorimetric study of the thermal unfolding of myosin rod, light meromyosin, and subfragment 2.

The thermal unfolding of myosin rod, light meromyosin (LMM), and myosin subfragment 2 (S-2) was studied by differential scanning calorimetry (DSC) over the pH range of 6.5-9.0 in 0.5M KCl and either 0.20 M sodium phosphate or 0.15M sodium pyrophosphate. Two rod samples were examined: one was purified by Sephadex G-200 without prior denaturation (native rod), and the other was purified by a cycle of denaturation-renaturation followed by Sephacryl S-200 chromatography (renatured rod). There were clearly distinguishable differences in the calorimetric behavior of these two samples. At pH 7.0 in phosphate the DSC curves of native rod were deconvoluted into six endothermic two-state transitions with melting temperatures in the range of 46-67 degrees C and a total enthalpy of 4346 kJ/mol. Under identical conditions the melting profile of LMM was resolved into five endothermic peaks with transition temperatures in the range of 45-66 degrees C, and the thermal profile of long S-2 was resolved into two endotherms, 46 and 57 degrees C. Transition 4 observed with native rod was present in the deconvoluted DSC curve for long S-2, but absent in the DSC curve for LMM. This transition was identified with the high-temperature transition detected with long S-2 and attributed to the melting of the coiled-coil alpha-helical segment of subfragment 2 (short S-2). The low-temperature transition of long S-2 was attributed to the unfolding of the hinge region. The smallest transition temperatures observed for all three fragments were 45-46 degrees C. It is suggested that the most unstable domain in rod (domain 1) responsible for the 46 degrees C transition includes both the hinge region, which is the C-terminal segment of long S-2, and a short N-terminal segment of LMM. This domain, accounting for 21% of the rod structure, contains the S-2/LMM junction, and upon proteolytic cleavage yields the C-terminal and N-terminal ends of long S-2 and LMM, respectively. Over the pH range of 6.5-7.5, the observed specific heat of denaturation of rod was approximately equal to the sum of the specific heats of LMM and S-2. This finding provides an additional argument for the existence of independent domains in myosin rod.

Animals↗
PubMed1989

Models for two tRNAs bound to successive codons on mRNA on the ribosome.

We have investigated the structural changes necessary to build a model complex of two molecules of phenylalanine transfer RNA (tRNA(Phe) bound to successive codons in a short segment of a model messenger RNA (mRNA), consisting of U6. We keep the mRNA in an ideal helical conformation, deforming the tRNAs as necessary to eliminate steric overlaps while bringing the two 3' termini together. The resulting model has the two tRNAs oriented relative to one another in a manner that is very similar to a model developed by McDonald and Rein (1) in which the tRNAs maintain their ideal crystallographic conformations and all of the deformations are introduced into the mRNA. Consequently, regardless of how one divides the deformations between the tRNAs and the mRNA it is clear that, on the ribosome, the tRNA in the P site has its "front" side (that side with the variable loop) close to the "back" side of the tRNA in the A site (that side with the D loop). The space between the two molecules must be left free on the ribosome, in order to facilitate the transition from the A site to the P site. A detailed pathway is also proposed for changing the anticodon loop structure from that of the A site to that of the P site. The anticodon loop is always kept in a 3'-stacked conformation, since we find that the shift between the 3'-stacked and 5'-stacked structures proposed by Woese (2) is not feasible.

Codon↗
PubMed1988

DNA stem-loop structures in oligopurine-oligopyrimidine triplexes.

Closed circular DNA containing polypurine-polypyrimidine sequences can adopt a triple helical stem-loop structure under supercoiling pressure. We describe an automated procedure for building model loops and its application to the investigation of the polypyrimidine loop at the end of such a triple helical stem. All possible combinations of 3'-stacked and 5'-stacked structures have been examined for loops containing three, four, five, and six nucleotides. The lowest energy conformation is a four-membered loop with all bases stacked on the strand at the 3' end of the loop. The model predicts that sequences (GA)n, (GGGA)n and (GAAA)n should form the stem-loop structure more easily than (GGA)n and (GAA)n. It is also predicted that when a polypurine-polypyrimidine sequence converts from a double stranded structure to a triple stranded stem-loop, the most favorable conditions are those where an even number of basepairs makes the transition. Experimental tests of these predictions are also described.

DNA↗
PubMed1988

Specificity of mutagenesis by 4-aminobiphenyl. A possible role for N-(deoxyadenosin-8-yl)-4-aminobiphenyl as a premutational lesion.

Mutagenesis by N-acetoxy-N-trifluoroacetyl-4-aminobiphenyl, a reactive form of the human bladder carcinogen 4-aminobiphenyl (ABP), was studied in Escherichia coli virus M13mp10. N-acetoxy-N-trifluoroacetyl-4-ABP-treated DNA containing 140 lesions/duplex genome, when introduced into excision repair-competent cells induced for SOS mutagenic processing, resulted in a 40-fold increase in mutation frequency over background in the lacZ alpha gene fragment. DNA sequence changes were determined for 20 independent mutants. G-C base pairs were the major targets for base pair substitution mutations, although significant mutagenic activity was also observed at certain A-T base pairs. Deletion and frameshift mutations also were found in this sample. The salient feature of this partial "mutational spectrum" was a hotspot that occurred at position 6357 (amino acid 30 of the beta-galactosidase fragment encoded by M13mp10); this A-T to T-A transversion appeared in 6 of the 20 mutants. The property of ABP to mutate A-T base pairs was consistent with the result that N-hydroxy-ABP reverted Salmonella typhimurium strain TA104, which is presumed to revert primarily due to mutations at these sites. The ability of the major carcinogen-DNA adduct formed by ABP in vivo and in vitro, N-(deoxyguanosin-8-yl)-4-aminobiphenyl, to cause base pair substitution mutations was also investigated. This adduct was positioned specifically in the minus strand at position 6270 in duplex M13mp10 DNA. In the presence of the mutagenesis-enhancing plasmid pGW16 and UV induction of SOS mutagenic processing, it was shown that fewer than 0.02% of the adducts resulted in transition or transversion mutations following transfection of DNA into excision-repair competent cells. Similar results were obtained in uvrA and uvrC backgrounds. Although the major adduct did not cause base substitution mutations under these experimental conditions, the contribution of this lesion to the entire spectrum of mutations in the lacZ alpha fragment seems likely.

Amino Acid Sequence↗
PubMed1988

AUGUR: a program to predict, display and analyze the tertiary structure of B-DNA.

AUGUR is a program to predict, display and analyze the three-dimensional structure of B-DNA. The user can choose one of six models to predict the helical parameters of a given sequence. These parameters are then used to generate the coordinates of the DNA model in three-dimensional space (trajectory). The trajectory can be displayed and rotated on a graphics terminal. The trajectory and helical parameters can also be searched for bends and structural homologues.

Computer Graphics↗