Stephen C. Harvey

A DETAILED MOLECULAR BELT MODEL FOR APOLIPOPROTEIN A-I IN DISCOIDAL HIGH DENSITY LIPOPROTEIN

Apolipoprotein A-I (apoA-I) is the principal protein of high density lipoprotein particles (HDL). ApoA-I contains a globular N-terminal domain (residues 1–43) and a lipid-binding C-terminal domain (residues 44–243). Here we propose a detailed model for the smallest discoidal HDL, consisting of two apoA-I molecules wrapped beltwise around a small patch of bilayer containing 160 lipid molecules. The C-terminal domain of each monomer is ringlike, a curved, planar amphipathic α helix with an average of 3.67 residues per turn, and with the hydrophobic surface curved toward the lipids. We have explored all possible geometries for forming the dimer of stacked rings, subject to the hypothesis that the optimal geometry will maximize intermolecular salt bridge interactions. The resulting model is an antiparallel arrangement with an alignment matching that of the (nonplanar) crystal structure of lipid-free apoA-I.

The Amphipathic α Helix: A Multifunctional Structural Motif in Plasma Apolipoproteins

Publisher Summary The dominant structural motif of the peripheral apolipoproteins is the amphipathic helix, which is responsible for the reversible association of these proteins with lipids, as well as for many biological functions mediated by these apolipoproteins. This chapter reviews the different classes of amphipathic helices, using a combination of powerful computer programs to develop a comparison database and to analyze these structures. It also discusses their evolutionary origins, physical-chemical properties, X-ray structure determination, and conformational analysis. Although the structures of these lipoprotein classes are similar, they differ in relative proportion of lipids, in the apolipoprotein: lipid ratio and in the apolipoprotein species. The amphipathic α helix plays a pivotal role in the structure and functions of the exchangeable apolipoproteins. Site-directed mutagenesis and other molecular biology-based techniques are available for probing the structural motif. The location and properties of the amphipathic helices in apolipoproteins and the results are compared with recently developed and ever-expanding computer methods for the location and characterization. A variety of structure-function studies, including the activation of lipoprotein lipase, receptor recognition, lecithin-cholesterol acyltransferase (LCAT) activation, and antiviral and anti-inflammatory activities are also discussed.

The flying ice cube: Velocity rescaling in molecular dynamics leads to violation of energy equipartition

This article describes an unexpected phenomenon encountered during MD simulations: velocity rescaling using standard protocols can systematically change the proportion of total kinetic energy (KE) found in motions associated with the various degrees of freedom. Under these conditions, the simulation violates the principle of equipartition of energy, which requires a mean kinetic energy of RT/2 in each degree of freedom. A particularly pathological form of this problem occurs if one does not periodically remove the net translation of (and rotation about) the center of mass. In this case, almost all of the kinetic energy is converted into these two kinds of motion, producing a system with almost no kinetic energy associated with the internal degrees of freedom. We call this phenomenon “the flying ice cube.” We present a mathematical analysis of a simple diatomic system with two degrees of freedom, to document the origin of the problem. We then present examples from three kinds of MD simulations, one being an in vacuo simulation on a diatomic system, one involving a low resolution model of DNA in vacuo, and the third using a traditional all‐atom DNA model with full solvation, periodic boundary conditions, and the particle mesh Ewald method for treating long‐range electrostatics. Finally, we discuss methods for avoiding the problem. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 726–740, 1998

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.

Use of photoaffinity crosslinking and molecular modeling to analyze the global architecture of ribonuclease P RNA.

Bacterial ribonuclease P (RNase P), an endonuclease involved in tRNA maturation, is a ribonucleoprotein containing a catalytic RNA. The secondary structure of this ribozyme is well established, but comparatively little is understood about its 3‐D structure. In this analysis, orientation and distance constraints between elements within the Escherichia coli RNase P RNA‐pre‐tRNA complex were determined by intra‐ and intermolecular crosslinking experiments. A molecular mechanics‐based RNA structure refinement protocol was used to incorporate the distance constraints indicated by crosslinking, along with the known secondary structure of RNase P RNA and the tertiary structure of tRNA, into molecular models. Seven different structures that satisfy the constraints equally well were generated and compared by superposition to estimate helix positions and orientations. Manual refinement within the range of conformations indicated by the molecular mechanics analysis was used to derive a model of RNase P RNA with bound substrate pre‐tRNA that is consistent with the crosslinking results and the available phylogenetic comparisons.

Structure and function of apolipoprotein A-I and high-density lipoprotein.

Structural biology and molecular modeling have provided intriguing insights into the atomic details of the lipid-associated structure of the major protein component of HDL, apo A-I. For the first time, an atomic resolution map is available for future studies of the molecular interactions of HDL in such biological processes as ABC1-regulated HDL assembly, LCAT activation, receptor binding, reverse lipid transport and HDL heterogeneity. Within the context of this paradigm, the current review summarizes the state of HDL research.

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.

Static contributions to the persistence length of DNA and dynamic contributions to DNA curvature

Long molecules of DNA have the statistical properties of a worm-like coil. Deviations from linearity occur both because of small dynamic bends induced by thermal motion and from a random distribution of static bends. The latter originate in the different conformations of each of the possible base pair sequences. In this paper a statistical theory of the persistence length of DNA is developed which includes both static and dynamic effects for each base pair sequence, as well as the sequence-dependent correlations of bending angles. The result applies to a generic DNA, i.e., the average over an ensemble of all possible sequences. The theory is also applied to the generation of the average properties of curved DNAs by an analytic method that includes dynamic averaging as well as correlated bends. These results provide information which supplements that obtained by others using Monte Carlo methods. The additivity relation 1/P = 1/P(S) + 1/P(d) proposed by Trifonov et al., where P is the persistence length and P(S) and P(d) are the persistence lengths arising from purely static and dynamic effects, respectively, has been verified to be accurate to better than 0.5%. This is true for both a simplified model and one that includes a complete set of static bends at all base pair sequences.

FLEXIBLE DNA : GENETICALLY UNSTABLE CTG.CAG AND CGG.CCG FROM HUMAN HEREDITARY NEUROMUSCULAR DISEASE GENES

The properties of duplex CTG·CAG and CGG·CCG, which are involved in the etiology of several hereditary neurodegenerative diseases, were investigated by a variety of methods, including circularization kinetics, apparent helical repeat determination, and polyacrylamide gel electrophoresis. The bending moduli were 1.13 × 10−19 erg·cm for CTG and 1.27 × 10−19 erg·cm for CGG, ∼40% less than for random B-DNA. Also, the persistence lengths of the triplet repeat sequences were ∼60% the value for random B-DNA. However, the torsional moduli and the helical repeats were 2.3 × 10−19 erg·cm and 10.4 base pairs (bp)/turn for CTG and 2.4 × 10−19 erg·cm and 10.3 bp/turn for CGG, respectively, all within the range for random B-DNA. Determination of the apparent helical repeat by the band shift assay indicated that the writhe of the repeats was different from that of random B-DNA. In addition, molecules of 224–245 bp in length (64–71 triplet repeats) were able to form topological isomers upon cyclization. The low bending moduli are consistent with predictions from crystallographic variations in slide, roll, and tilt. No unpaired bases or non-B-DNA structures could be detected by chemical and enzymatic probe analyses, two-dimensional agarose gel electrophoresis, and immunological studies. Hence, CTG and CGG are more flexible and highly writhed than random B-DNA and thus would be expected to act as sinks for the accumulation of superhelical density.

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.