Biomedical subjects

Stephen C Harvey

Publications and source records attributed to Stephen C Harvey.

At least 19 recordsLinked to original sources

PubMed2007

Structural and thermodynamic principles of viral packaging.

Packaging of genetic material inside a capsid is one of the major processes in the lifecycle of bacteriophages. To establish the basic principles of packing double-stranded DNA into a phage, we present a low-resolution model of bacteriophage varphi29 and report simulations of DNA packaging. The simulations show excellent agreement with available experimental data, including the forces of packaging and the average structures seen in cryo-electron microscopy. The conformation of DNA inside the bacteriophage is primarily determined by the shape of the capsid and the elastic properties of DNA, but the energetics of packaging are dominated by electrostatic repulsions and the large entropic penalty associated with DNA confinement. In this slightly elongated capsid, the DNA assumes a folded toroidal conformation, rather than a coaxial spool. The model can be used to study packaging of other bacteriophages with different shapes under a range of environmental conditions.

Bacillus Phages↗

PubMed2006

DNA-directed assembly of polyanilines: modified cytosine nucleotides transfer sequence programmability to a conjoined polymer.

A series of polyaniline (PANI) oligomers was constructed from monomer units covalently linked to duplex DNA through N-(2-aminoethyl) groups bonded through cytosines. DNA oligomers containing the aniline monomers were treated with horseradish peroxidase (HRP) and H2O2 under conditions known to cause polymerization of aniline. No change in the absorption spectrum of the DNA was observed for samples containing fewer than four contiguous aniline groups. However, for oligomers containing four, five, or six aniline units, treatment with HRP and H2O2 led to the appearance of absorption features characteristic of the conducting "proton doped" emeraldine oxidation state of PANI. Molecular modeling shows that the DNA is distorted in the region of the PANI, but flanking regions of the DNA maintain their B-form structure. These findings provide a method to exploit the self-recognition, self-assembly, and sequence programmability of DNA for the formation of conducting polymers.

Aniline Compounds↗

PubMed2006

Novel changes in discoidal high density lipoprotein morphology: a molecular dynamics study.

ApoA-I is a uniquely flexible lipid-scavenging protein capable of incorporating phospholipids into stable particles. Here we report molecular dynamics simulations on a series of progressively smaller discoidal high density lipoprotein particles produced by incremental removal of palmitoyloleoylphosphatidylcholine via four different pathways. The starting model contained 160 palmitoyloleoylphosphatidylcholines and a belt of two antiparallel amphipathic helical lipid-associating domains of apolipoprotein (apo) A-I. The results are particularly compelling. After a few nanoseconds of molecular dynamics simulation, independent of the starting particle and method of size reduction, all simulated double belts of the four lipidated apoA-I particles have helical domains that impressively approximate the x-ray crystal structure of lipid-free apoA-I, particularly between residues 88 and 186. These results provide atomic resolution models for two of the particles produced by in vitro reconstitution of nascent high density lipoprotein particles. These particles, measuring 95 angstroms and 78 angstroms by nondenaturing gradient gel electrophoresis, correspond in composition and in size/shape (by negative stain electron microscopy) to the simulated particles with molar ratios of 100:2 and 50:2, respectively. The lipids of the 100:2 particle family form minimal surfaces at their monolayer-monolayer interface, whereas the 50:2 particle family displays a lipid pocket capable of binding a dynamic range of phospholipid molecules.

Apolipoprotein A-I↗

PubMed2006

The RNA Ontology Consortium: an open invitation to the RNA community.

The aim of the RNA Ontology Consortium (ROC) is to create an integrated conceptual framework-an RNA Ontology (RO)-with a common, dynamic, controlled, and structured vocabulary to describe and characterize RNA sequences, secondary structures, three-dimensional structures, and dynamics pertaining to RNA function. The RO should produce tools for clear communication about RNA structure and function for multiple uses, including the integration of RNA electronic resources into the Semantic Web. These tools should allow the accurate description in computer-interpretable form of the coupling between RNA architecture, function, and evolution. The purposes for creating the RO are, therefore, (1) to integrate sequence and structural databases; (2) to allow different computational tools to interoperate; (3) to create powerful software tools that bring advanced computational methods to the bench scientist; and (4) to facilitate precise searches for all relevant information pertaining to RNA. For example, one initial objective of the ROC is to define, identify, and classify RNA structural motifs described in the literature or appearing in databases and to agree on a computer-interpretable definition for each of these motifs. To achieve these aims, the ROC will foster communication and promote collaboration among RNA scientists by coordinating frequent face-to-face workshops to discuss, debate, and resolve difficult conceptual issues. These meeting opportunities will create new directions at various levels of RNA research. The ROC will work closely with the PDB/NDB structural databases and the Gene, Sequence, and Open Biomedical Ontology Consortia to integrate the RO with existing biological ontologies to extend existing content while maintaining interoperability.

Databases, Genetic↗

PubMed2006

A structural model for the large subunit of the mammalian mitochondrial ribosome.

Protein translation is essential for all forms of life and is conducted by a macromolecular complex, the ribosome. Evolutionary changes in protein and RNA sequences can affect the 3D organization of structural features in ribosomes in different species. The most dramatic changes occur in animal mitochondria, whose genomes have been reduced and altered significantly. The RNA component of the mitochondrial ribosome (mitoribosome) is reduced in size, with a compensatory increase in protein content. Until recently, it was unclear how these changes affect the 3D structure of the mitoribosome. Here, we present a structural model of the large subunit of the mammalian mitoribosome developed by combining molecular modeling techniques with cryo-electron microscopic data at 12.1A resolution. The model contains 93% of the mitochondrial rRNA sequence and 16 mitochondrial ribosomal proteins in the large subunit of the mitoribosome. Despite the smaller mitochondrial rRNA, the spatial positions of RNA domains known to be involved directly in protein synthesis are essentially the same as in bacterial and archaeal ribosomes. However, the dramatic reduction in rRNA content necessitates evolution of unique structural features to maintain connectivity between RNA domains. The smaller rRNA sequence also limits the likelihood of tRNA binding at the E-site of the mitoribosome, and correlates with the reduced size of D-loops and T-loops in some animal mitochondrial tRNAs, suggesting co-evolution of mitochondrial rRNA and tRNA structures.

Animals↗

PubMed2005

Exploring the flexibility of ribosome recycling factor using molecular dynamics.

Ribosome recycling factor is proposed to be flexible, and that flexibility is believed to be important to its function. Here we use molecular dynamics to test the flexibility of Escherichia coli RRF (ecRRF) with and without decanoic acid bound to a hydrophobic pocket between domains 1 and 2, and Thermus thermophilus RRF (ttRRF) with and without a mutation in the hinge between domains 1 and 2. Our simulations show that the structure of ecRRF rapidly goes from having an interdomain angle of 124 degrees to an angle of 98 degrees independently of the presence of decanoic acid. The simulations also show that the presence or absence of decanoic acid leads to changes in ecRRF flexibility. Simulations of wild-type and mutant ttRRF (R32G) show that mutating Arg-32 to glycine decreases RRF flexibility. This was unexpected because the range of dihedral angles for arginine is limited relative to glycine. Furthermore, the interdomain angle of wild-type T. thermophilus goes from 81 degrees to 118 degrees whereas the R32G mutant remains very close to the crystallographic angle of 78 degrees . We propose that this difference accounts for the fact that mutant ttRRF complements an RRF deficient strain of E. coli whereas wild-type ttRRF does not. When the ensemble of RRF structures is modeled into the ribosomal crystal structure, a series of overlaps is found that corresponds with regions where conformational changes have been found in the cryoelectron microscopic structure of the RRF/ribosome complex, and in the crystal structure of a cocomplex of RRF with the 50S subunit. There are also overlaps with the P-site, suggesting that RRF flexibility plays a role in removing the deacylated P-site tRNA during termination of translation.

Binding Sites↗

PubMed2005

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↗

PubMed2004

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↗

PubMed2004

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↗

PubMed2004

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↗

PubMed2003

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↗

PubMed2003

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↗

PubMed2003

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↗

PubMed2003

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↗

PubMed2003

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↗

PubMed2003

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↗

PubMed2002

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↗

PubMed2002

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↗