PubMed Health⌕ Search

Biomedical subjects

Stephen C Harvey

Publications and source records attributed to Stephen C Harvey.

17 recordsLinked to original sources

Assembly of lipoprotein particles containing apolipoprotein-B: structural model for the nascent lipoprotein particle.

Apolipoprotein B (apoB) is the major protein component of large lipoprotein particles that transport lipids and cholesterol. We have developed a detailed model of the first 1000 residues of apoB using standard sequence alignment programs (ClustalW and MACAW) and the MODELLER6 package for three-dimensional homology modeling. The validity of the apoB model was supported by conservation of disulfide bonds, location of all proline residues in turns and loops, and conservation of the hydrophobic faces of the two C-terminal amphipathic beta-sheets, betaA (residues 600-763) and betaB (residues 780-1000). This model suggests a lipid-pocket mechanism for initiation of lipoprotein particle assembly. In a previous model we suggested that microsomal triglyceride transfer protein might play a structural role in completion of the lipid pocket. We no longer think this likely, but instead propose a hairpin-bridge mechanism for lipid pocket completion. Salt-bridges between four tandem charged residues (717-720) in the turn of the hairpin-bridge and four tandem complementary residues (997-1000) at the C-terminus of the model lock the bridge in the closed position, enabling the deposition of an asymmetric bilayer within the lipid pocket.

Amino Acid Sequence↗

Double belt structure of discoidal high density lipoproteins: molecular basis for size heterogeneity.

We recently proposed an all-atom model for apolipoprotein (apo) A-I in discoidal high-density lipoprotein in which two monomers form stacked antiparallel helical rings rotationally aligned by interhelical salt-bridges. The model can be derived a priori from the geometry of a planar bilayer disc that constrains the hydrophobic face of a continuous amphipathic alpha helix in lipid-associated apoA-I to a plane inside of an alpha-helical torus. This constrains each apoA-I monomer to a novel conformation, that of a slightly unwound, curved, planar amphipathic alpha 11/3 helix (three turns per 11 residues). Using non-denaturing gradient gel electrophoresis, we show that dimyristoylphosphocholine discs containing two apoA-I form five distinct particles with maximal Stokes diameters of 98 A (R2-1), 106 A (R2-2), 110 A (R2-3), 114 A (R2-4) and 120 A (R2-5). Further, we show that the Stokes diameters of R2-1 and R2-2 are independent of the N-terminal 43 residues (the flexible domain) of apoA-I, while the flexible domain is necessary and sufficient for the formation of the three larger complexes. On the basis of these results, the conformation of apoA-I on the R2-2 disc can be modeled accurately as an amphipathic helical double belt extending the full length of the lipid-associating domain with N and C-terminal ends in direct contact. The smallest of the discs, R2-1, models as the R2-2 conformation with an antiparallel 15-18 residue pairwise segment of helixes hinged off the disc edge. The conformations of full-length apoA-I on the flexible domain-dependent discs (R2-3, R2-4 and R2-5) model as the R2-2 conformation extended on the disc edge by one, two or three of the 11-residue tandem amphipathic helical repeats (termed G1, G2 and G3), respectively, contained within the flexible domain. Although we consider these results to favor the double belt model, the topographically very similar hairpin-belt model cannot be ruled out entirely.

Apolipoprotein A-I↗

Packaging double-helical DNA into viral capsids.

DNA packaging in bacteriophage P4 has been examined using a molecular mechanics model with a reduced representation containing one pseudoatom per turn of the double helix. The model is a discretized version of an elastic continuum model. The DNA is inserted piecewise into the model capsid, with the structure being reoptimized after each piece is inserted. Various optimization protocols were investigated, and it was found that molecular dynamics at a very low temperature (0.3 K) produces the optimal packaged structure. This structure is a concentric spool, rather than the coaxial spool that has been commonly accepted for so many years. This geometry, which was originally suggested by Hall and Schellman in 1982 (Biopolymers Vol. 21, pp. 2011-2031), produces a lower overall elastic energy than coaxial spooling.

Bacteriophages↗

Nodavirus coat protein imposes dodecahedral RNA structure independent of nucleotide sequence and length.

The nodavirus Flock house virus (FHV) has a bipartite, positive-sense RNA genome that is packaged into an icosahedral particle displaying T=3 symmetry. The high-resolution X-ray structure of FHV has shown that 10 bp of well-ordered, double-stranded RNA are located at each of the 30 twofold axes of the virion, but it is not known which portions of the genome form these duplex regions. The regular distribution of double-stranded RNA in the interior of the virus particle indicates that large regions of the encapsidated genome are engaged in secondary structure interactions. Moreover, the RNA is restricted to a topology that is unlikely to exist during translation or replication. We used electron cryomicroscopy and image reconstruction to determine the structure of four types of FHV particles that differed in RNA and protein content. RNA-capsid interactions were primarily mediated via the N and C termini, which are essential for RNA recognition and particle assembly. A substantial fraction of the packaged nucleic acid, either viral or heterologous, was organized as a dodecahedral cage of duplex RNA. The similarity in tertiary structure suggests that RNA folding is independent of sequence and length. Computational modeling indicated that RNA duplex formation involves both short-range and long-range interactions. We propose that the capsid protein is able to exploit the plasticity of the RNA secondary structures, capturing those that are compatible with the geometry of the dodecahedral cage.

Animals↗

Incorporation of aminoacyl-tRNA into the ribosome as seen by cryo-electron microscopy.

Aminoacyl-tRNAs (aa-tRNAs) are delivered to the ribosome as part of the ternary complex of aa-tRNA, elongation factor Tu (EF-Tu) and GTP. Here, we present a cryo-electron microscopy (cryo-EM) study, at a resolution of approximately 9 A, showing that during the incorporation of the aa-tRNA into the 70S ribosome of Escherichia coli, the flexibility of aa-tRNA allows the initial codon recognition and its accommodation into the ribosomal A site. In addition, a conformational change observed in the GTPase-associated center (GAC) of the ribosomal 50S subunit may provide the mechanism by which the ribosome promotes a relative movement of the aa-tRNA with respect to EF-Tu. This relative rearrangement seems to facilitate codon recognition by the incoming aa-tRNA, and to provide the codon-anticodon recognition-dependent signal for the GTPase activity of EF-Tu. From these new findings we propose a mechanism that can explain the sequence of events during the decoding of mRNA on the ribosome.

Cryoelectron Microscopy↗

Study of the structural dynamics of the E coli 70S ribosome using real-space refinement.

Cryo-EM density maps showing the 70S ribosome of E. coli in two different functional states related by a ratchet-like motion were analyzed using real-space refinement. Comparison of the two resulting atomic models shows that the ribosome changes from a compact structure to a looser one, coupled with the rearrangement of many of the proteins. Furthermore, in contrast to the unchanged inter-subunit bridges formed wholly by RNA, the bridges involving proteins undergo large conformational changes following the ratchet-like motion, suggesting an important role of ribosomal proteins in facilitating the dynamics of translation.

Bacterial Proteins↗

A structural model for the assembly of the 30S subunit of the ribosome.

The order in which proteins bind to 16S rRNA, the assembly map, was determined by Nomura and co-workers in the early 1970s. The assembly map shows the dependencies of binding of successive proteins but fails to address the relationship of these dependencies to the three-dimensional folding of the ribosome. Here, using molecular mechanics techniques, we rationalize the order of protein binding in terms of ribosomal folding. We determined the specific contacts between the ribosomal proteins and 16S rRNA from a crystal structure of the 30S subunit (1FJG). We then used these contacts as restraints in a rigid body Monte-Carlo simulation with reduced-representation models of the RNA and proteins. Proteins were added sequentially to the RNA in the order that they appear in the assembly map. Our results show that proteins nucleate the folding of the head, platform, and body domains, but they do not strongly restrict the orientations of the domains relative to one another. We also examined the contributions of individual proteins to the formation of binding sites for sequential proteins in the assembly process. Binding sites for the primary binding proteins are generally more ordered in the naked RNA than those for other proteins. Furthermore, we examined one pathway in the assembly map and found that the addition of early binding proteins helps to organize the RNA around the binding sites of proteins that bind later. It appears that the order of assembly depends on the degree of pre-organization of each protein's binding site at a given stage of assembly, and the impact that the binding of each protein has on the organization of the remaining unoccupied binding sites.

Binding Sites↗

Motifs in nucleic acids: molecular mechanics restraints for base pairing and base stacking.

In building and refining nucleic acid structures, it is often desirable to enforce particular base pairing and/or base stacking interactions. Energy-based modeling programs with classical molecular mechanics force fields do not lend themselves to the easy imposition of penalty terms corresponding to such restraints, because the requirement that two bases lie in or near the same plane (pairing) or that they lie in parallel planes (stacking) cannot be easily expressed in terms of traditional interactions involving two atoms (bonds), three atoms (angles), or four atoms (torsions). Here we derive expressions that define a collection of pseudobonds and pseudoangles through which molecular mechanics restraints for base pairing and stacking can be imposed. We have implemented these restraints into the JUMNA package for modeling DNA and RNA structures. JUMNA scripts can specify base pairing with a variety of standard geometries (Watson-Crick, Hoogsteen, wobble, etc.), or with user-defined geometries; they can also specify stacking arrangements. We have also implemented "soft-core" functions to modify van der Waals and electrostatic interactions to avoid steric conflicts in particularly difficult refinements where two backbones need to pass through one another. Test cases are presented to show the utility of the method. The restraints could be adapted for implementation in other molecular mechanics packages.

Base Pairing↗

The calcium-saturated cTnI/cTnC complex: structure of the inhibitory region of cTnI.

The contiguous inhibitory and regulatory regions of troponin I in the heterotrimeric troponin complex play a critical role in Ca(2+) activation of striated muscle. Knowledge of the structure of this critical region within the complex will enhance efforts toward understanding regulatory mechanisms. Toward this goal, we have used simulated annealing to study the structure of the inhibitory and regulatory regions of cardiac muscle troponin I in the calcium-saturated complex formed between cardiac troponin C and cardiac troponin I. We have incorporated distances determined experimentally by Förster resonance energy transfer in the full-length complex, rather than using peptides derived from cTnI. For these models, we assume a helix-loop-helix conformation for the inhibitory region. We have found several structures that satisfy the experimental constraints fairly well. Although it is not possible to eliminate any of these models at this time, future studies with additional experimental restraints will yield insights on the mechanisms of calcium regulation in cardiac muscle.

Amino Acid Sequence↗

Probing the importance of tRNA anticodon: human immunodeficiency virus type 1 (HIV-1) RNA genome complementarity with an HIV-1 that selects tRNA(Glu) for replication.

The initiation of human immunodeficiency virus type 1 (HIV-1) reverse transcription occurs at the primer binding site (PBS) that is complementary to the 3'-terminal nucleotides of tRNA(3)(Lys). Why all known strains of HIV-1 select tRNA(3)(Lys) for replication is unknown. Previous studies on the effect of altering the PBS of HIV-1 on replication identified an HIV-1 with a PBS complementary to tRNA(Glu). Since the virus was not initially designed to use tRNA(Glu), the virus had selected tRNA(Glu) from the intracellular pool of tRNA for use in replication. Further characterization of HIV-1 that uses tRNA(Glu) may provide new insights into the preference for tRNA(3)(Lys). HIV-1 constructed with the PBS complementary to tRNA(Glu) was more stable than HIV-1 with the PBS complementary to tRNA(Met) or tRNA(His); however, all of these viruses eventually reverted back to using tRNA(3)(Lys) following growth in SupT1 cells or peripheral blood mononuclear cells (PBMCs). New HIV-1 mutants with nucleotides in U5 complementary to the anticodon of tRNA(Glu) remained stable when grown in SupT1 cells or PBMCs, although the mutants grew more slowly than the wild-type virus. Sequence analysis of the U5 region and the PBS revealed additional mutations predicted to further promote tRNA-viral genome interaction. The results support the importance of the tRNA anticodon-genome interaction in the selection of the tRNA primer and highlight the fact that unique features of tRNA(3)(Lys) are exploited by HIV-1 for selection as the reverse transcription primer.

Anticodon↗

Investigation of viral DNA packaging using molecular mechanics models.

A simple molecular mechanics model has been used to investigate optimal spool-like packing conformations of double-stranded DNA molecules in viral capsids with icosahedral symmetry. The model represents an elastic segmented chain by using one pseudoatom for each ten basepairs (roughly one turn of the DNA double helix). Force constants for the various terms in the energy function were chosen to approximate known physical properties, and a radial restraint was used to confine the DNA into a sphere with a volume corresponding to that of a typical bacteriophage capsid. When the DNA fills 90% of the spherical volume, optimal packaging is obtained for coaxially spooled models, but this result does not hold when the void volume is larger. When only 60% of the spherical volume is filled with DNA, the lowest energy structure has two layers, with a coiled core packed at an angle to an outer coaxially spooled shell. This relieves bending strain associated with tight curvature near the poles in a model with 100% coaxial spooling. Interestingly, the supercoiling density of these models is very similar to typical values observed in plasmids in bacterial cells. Potential applications of the methodology are also discussed.

DNA, Viral↗

Molecular dynamics simulations on discoidal HDL particles suggest a mechanism for rotation in the apo A-I belt model.

Apolipoprotein A-I (apo A-I) is the major protein component of high-density lipoprotein (HDL) particles. Elevated levels of HDL in the bloodstream have been shown to correlate strongly with a reduced risk factor for atherosclerosis. Molecular dynamics simulations have been carried out on three separate model discoidal high-density lipoprotein particles (HDL) containing two monomers of apo A-I and 160 molecules of palmitoyloleoylphosphatidylcholine (POPC), to a time-scale of 1ns. The starting structures were on the basis of previously published molecular belt models of HDL consisting of the lipid-binding C-terminal domain (residues 44-243) wrapped around the circumference of a discoidal HDL particle. Subtle changes between two of the starting structures resulted in significantly different behavior during the course of the simulation. The results provide support for the hypothesis of Segrest et al. that helical registration in the molecular belt model of apo A-I is modulated by intermolecular salt bridges. In addition, we propose an explanation for the presence of proline punctuation in the molecular belt model, and for the presence of two 11-mer helical repeats interrupting the otherwise regular pattern of 22-mer helical repeats in the lipid-binding domain of apo A-I.

1,2-Dipalmitoylphosphatidylcholine↗

Comparative models for human apolipoprotein A-I bound to lipid in discoidal high-density lipoprotein particles.

We have constructed a series of models for apolipoprotein A-I (apo A-I) bound to discoidal high-density lipoprotein (HDL) particles, based upon the molecular belt model [Segrest, J. P., et al. (1999) J. Biol. Chem. 274, 31755-31758] and helical hairpin models [Rogers, D. P., et al. (1998) Biochemistry 37, 11714-11725], and compared these with picket fence models [Phillips, J. C., et al. (1997) Biophys. J. 73, 2337-2346]. Molecular belt models for discoidal HDL particles with differing diameters are presented, illustrating that the belt model can explain the discrete changes in HDL particle size observed experimentally. Hairpin models are discussed for the binding of apo A-I to discoidal HDL particles with diameters identical to those for the molecular belt model. Two models are presented for the binding of three monomers of apo A-I to a 150 A diameter discoidal HDL particle. In one model, two monomers of apo A-I bind to the exterior of the HDL particle in an antiparallel belt, with a third monomer of apo A-I bound to the disk in a hairpin conformation. In the second model, all three monomers of apo A-I are bound to the discoidal HDL particle in a hairpin conformation. Previously published experimental data for each model are reviewed, with FRET favoring either the belt or hairpin models over the picket fence models for HDL particles with diameters of 105 A. Naturally occurring mutations appear to favor the belt model for the 105 A particles, while the 150 A HDL particles favor the presence of at least one hairpin.

Apolipoprotein A-I↗

Modeling a minimal ribosome based on comparative sequence analysis.

We have determined the three-dimensional organization of ribosomal RNAs and proteins essential for minimal ribosome function. Comparative sequence analysis identifies regions of the ribosome that have been evolutionarily conserved, and the spatial organization of conserved domains is determined by mapping these onto structures of the 30S and 50S subunits determined by X-ray crystallography. Several functional domains of the ribosome are conserved in their three-dimensional organization in the Archaea, Bacteria, Eucaryotic nuclear, mitochondria and chloroplast ribosomes. In contrast, other regions from both subunits have shifted their position in three-dimensional space during evolution, including the L11 binding domain and the alpha-sarcin-ricin loop (SRL). We examined conserved bridge interactions between the two ribosomal subunits, giving an indication of which contacts are more significant. The tRNA contacts that are conserved were also determined, highlighting functional interactions as the tRNA moves through the ribosome during protein synthesis. To augment these studies of a large collection of comparative structural models sampled from all major branches on the phylogenetic tree, Caenorhabditis elegans mitochondrial rRNA is considered individually because it is among the smallest rRNA sequences known. The C.elegans model supports the large collection of comparative structure models while providing insight into the evolution of mitochondrial ribosomes.

Animals↗

Identification of the proteoglycan binding site in apolipoprotein B48.

An initial event in atherosclerosis is the retention of lipoproteins within the intima of the vessel wall. Previously we identified Site B (residues 3359-3369) in apolipoprotein (apo) B100 as the proteoglycan binding sequence in low density lipoproteins (LDLs) and showed that the atherogenicity of apoB-containing lipoproteins is linked to their affinity for artery wall proteoglycans. However, both apoB100- and apoB48-containing lipoproteins are equally atherogenic even though Site B lies in the carboxyl-terminal half of apoB100 and is absent in apoB48. If binding to proteoglycans is a key step in atherogenesis, apoB48-containing lipoproteins must bind to proteoglycans via other proteoglycan binding sites in the amino-terminal 48% of apoB. In vitro studies have identified five clusters of basic amino acids in delipidated apoB48 that bind negatively charged glycosaminoglycans. To determine which of these sites is functional on LDL particles, we analyzed the proteoglycan binding activity of recombinant human LDLs from transgenic mice or rat hepatoma cells. Substitution of neutral amino acids for the basic amino acids in Site B-Ib (residues 84-94) abolished the proteoglycan binding activity of recombinant apoB53. Carboxyl-truncated apoB80 bound biglycan with higher affinity than apoB100 and apoB48. ApoB80 in which Site B was mutated had the same affinity for proteoglycans as apoB48. These data support the hypothesis that the carboxyl terminus of apoB100 "masks" Site B-Ib, the amino-terminal proteoglycan binding site, and that this site is exposed in carboxyl-truncated forms of apoB. The presence of a proteoglycan binding site in the amino-terminal region of apoB may explain why apoB48- and apoB100-containing lipoproteins are equally atherogenic.

Amino Acid Sequence↗

RNAML: a standard syntax for exchanging RNA information.

Analyzing a single data set using multiple RNA informatics programs often requires a file format conversion between each pair of programs, significantly hampering productivity. To facilitate the interoperation of these programs, we propose a syntax to exchange basic RNA molecular information. This RNAML syntax allows for the storage and the exchange of information about RNA sequence and secondary and tertiary structures. The syntax permits the description of higher level information about the data including, but not restricted to, base pairs, base triples, and pseudoknots. A class-oriented approach allows us to represent data common to a given set of RNA molecules, such as a sequence alignment and a consensus secondary structure. Documentation about experiments and computations, as well as references to journals and external databases, are included in the syntax. The chief challenge in creating such a syntax was to determine the appropriate scope of usage and to ensure extensibility as new needs will arise. The syntax complies with the eXtensible Markup Language (XML) recommendations, a widely accepted standard for syntax specifications. In addition to the various generic packages that exist to read and interpret XML formats, an XML processor was developed and put in the open-source MC-Core library for nucleic acid and protein structure computer manipulation.

Databases, Nucleic Acid↗