Similarity of protein-RNA interfaces based on motif analysis.

We have developed a method for determination of the similarity of pairs of protein-RNA complexes, which we refer to as SIMA (Similarity by Identity and Motif Alignment). The key element in the SIMA method is the description of the protein-RNA interface in terms of motifs (salt bridges, aromatic stacking interactions, nonaromatic stacks, hydrophobic interactions, and hydrogen-bonded motifs), in addition to single hydrogen bonds and van der Waals contacts. Based on a pairwise scoring function combining motif alignment with identity of the protein and RNA sequences, we define a SIMA score for any pair of protein-RNA complexes. A positive score indicates similarity between the complexes. We used the SIMA method to identify 284 nonredundant binary protein-RNA complexes out of 776 such complexes in 382 nonribosomal protein-RNA structure files obtained from the RCSB database. SIMA allows rapid and quantitative comparison of protein-RNA interfaces and may be useful for interface classification with potential functional and evolutionary implications.

Site-directed spin labeling measurements of nanometer distances in nucleic acids using a sequence-independent nitroxide probe.

In site-directed spin labeling (SDSL), local structural and dynamic information is obtained via electron paramagnetic resonance (EPR) spectroscopy of a stable nitroxide radical attached site-specifically to a macromolecule. Analysis of electron spin dipolar interactions between pairs of nitroxides yields the inter-nitroxide distance, which provides quantitative structural information. The development of pulse EPR methods has enabled such distance measurements up to 70 A in bio-molecules, thus opening up the possibility of SDSL global structural mapping. This study evaluates SDSL distance measurement using a nitroxide (designated as R5) that can be attached, in an efficient and cost-effective manner, to a phosphorothioate backbone position at arbitrary DNA or RNA sequences. R5 pairs were attached to selected positions of a dodecamer DNA duplex with a known NMR structure, and eight distances, ranging from 20 to 40 A, were measured using double electron-electron resonance (DEER). The measured distances correlated strongly (R2 = 0.98) with the predicted values calculated based on a search of sterically allowable R5 conformations in the NMR structure, thus demonstrating accurate distance measurements using R5. Furthermore, distance measurement in a 42 kD DNA was demonstrated. The results establish R5 as a sequence-independent probe for global structural mapping of DNA and DNA-protein complexes.

The role of positively charged amino acids and electrostatic interactions in the complex of U1A protein and U1 hairpin II RNA.

Previous kinetic investigations of the N-terminal RNA recognition motif (RRM) domain of spliceosomal protein U1A, interacting with its RNA target U1 hairpin II, provided experimental evidence for a 'lure and lock' model of binding in which electrostatic interactions first guide the RNA to the protein, and close range interactions then lock the two molecules together. To further investigate the 'lure' step, here we examined the electrostatic roles of two sets of positively charged amino acids in U1A that do not make hydrogen bonds to the RNA: Lys20, Lys22 and Lys23 close to the RNA-binding site, and Arg7, Lys60 and Arg70, located on 'top' of the RRM domain, away from the RNA. Surface plasmon resonance-based kinetic studies, supplemented with salt dependence experiments and molecular dynamics simulation, indicate that Lys20 predominantly plays a role in association, while nearby residues Lys22 and Lys23 appear to be at least as important for complex stability. In contrast, kinetic analyses of residues away from the RNA indicate that they have a minimal effect on association and stability. Thus, well-positioned positively charged residues can be important for both initial complex formation and complex maintenance, illustrating the multiple roles of electrostatic interactions in protein-RNA complexes.

Kinetic analysis of the role of the tyrosine 13, phenylalanine 56 and glutamine 54 network in the U1A/U1 hairpin II interaction.

The A protein of the U1 small nuclear ribonucleoprotein particle, interacting with its stem-loop RNA target (U1hpII), is frequently used as a paradigm for RNA binding by recognition motif domains (RRMs). U1A/U1hpII complex formation has been proposed to consist of at least two steps: electrostatically mediated alignment of both molecules followed by locking into place, based on the establishment of close-range interactions. The sequence of events between alignment and locking remains obscure. Here we examine the roles of three critical residues, Tyr13, Phe56 and Gln54, in complex formation and stability using Biacore. Our mutational and kinetic data suggest that Tyr13 plays a more important role than Phe56 in complex formation. Mutational analysis of Gln54, combined with molecular dynamics studies, points to Arg52 as another key residue in association. Based on our data and previous structural and modeling studies, we propose that electrostatic alignment of the molecules is followed by hydrogen bond formation between the RNA and Arg52, and the sequential establishment of interactions with loop bases (including Tyr13). A quadruple stack, sandwiching two bases between Phe56 and Asp92, would occur last and coincide with the rearrangement of a C-terminal helix that partially occludes the RRM surface in the free protein.

The effect of cholesterol and monosialoganglioside (GM1) on the release and aggregation of amyloid beta-peptide from liposomes prepared from brain membrane-like lipids.

In order to investigate the influence of cholesterol (Ch) and monosialoganglioside (GM1) on the release and subsequent deposition/aggregation of amyloid beta peptide (Abeta)-(1-40) and Abeta-(1-42), we have examined Abeta peptide model membrane interactions by circular dichroism, turbidity measurements, and transmission electron microscopy (TEM). Model liposomes containing Abeta peptide and a lipid mixture composition similar to that found in the cerebral cortex membranes (CCM-lipid) have been prepared. In all, four Abeta-containing liposomes were investigated: CCM-lipid; liposomes with no GM1 (GM1-free lipid); those with no cholesterol (Ch-free lipid); liposomes with neither cholesterol nor GM1 (Ch-GM1-free lipid). In CCM liposomes, Abeta was rapidly released from membranes to form a well defined fibril structure. However, for the GM1-free lipid, Abeta was first released to yield a fibril structure about the membrane surface, then the membrane became disrupted resulting in the formation of small vesicles. In Ch-free lipid, a fibril structure with a phospholipid membrane-like shadow formed, but this differed from the well defined fibril structure seen for CCM-lipid. In Ch-GM1-free lipid, no fibril structure formed, possibly because of membrane solubilization by Abeta. The absence of fibril structure was noted at physiological extracellular pH (7.4) and also at liposomal/endosomal pH (5.5). Our results suggest a possible role for both Ch and GM1 in the membrane release of Abeta from brain lipid bilayers.

Computation of DNA backbone conformations.

We present an algorithm for the computation of 2'-deoxyribose-phosphodiester backbone conformations that are stereochemically compatible with a given arrangement of nucleic acid bases in a DNA structure. The algorithm involves the sequential computation of 2'-deoxyribose and phosphodiester conformers (collectively referred to as a backbone 'segment'), beginning at the 5'-end of a DNA strand. Computation of the possible segment conformations is achieved by the initial creation of a fragment library, with each fragment representing a set of bond lengths, bond angles and torsion angles. Following exhaustive searching of sugar conformations, each segment conformation is reduced to a single vector, defined by a specific distance, angle and torsion angle, that allows calculation of the O(1)' position. A given 'allowed' conformation of a backbone segment is determined based on its compatibility with the base positions and with the position of the preceding backbone segment. Initial computation of allowable segment conformations of a strand is followed by the determination of continuous backbone solutions for the strand, beginning at the 3'-end. The algorithm is also able to detect repeating segment conformations that arise in structures containing geometrically repeating dinucleotide steps. To illustrate the utility and properties of the algorithm, we have applied it to a series of experimental DNA structures. Regardless of the conformational complexity of these structures, we are able to compute backbone conformations for each structure. Hence, the algorithm, which is currently implemented within a new computer program NASDAC (Nucleic Acids: Structure, Dynamics and Conformation), should have generally applicability to the computation of DNA structures.