Global and local motions in ribonuclease A: a molecular dynamics study

The understanding of protein dynamics is one of the major goals of structural biology. A direct link between protein dynamics and function has been provided by x-ray studies performed on ribonuclease A (RNase A) (B. F. Rasmussen et al., Nature, 1992, Vol. 357, pp. 423-424; L. Vitagliano et al., Proteins: Structure, Function, and Genetics, 2002, Vol. 46, pp. 97-104). Here we report a 3 ns molecular dynamics simulation of RNase A in water aimed at characterizing the dynamical behavior of the enzyme. The analysis of local and global motions provides interesting insight on the dynamics/function relationship of RNase A. In agreement with previous crystallographic reports, the present study confirms that the RNase A active site is constituted by rigid (His12, Asn44, Thr45) and flexible (Lys41, Asp83, His119, Asp121) residues. The analysis of the global motions, performed using essential dynamics, shows that the two beta-sheet regions of RNase A move coherently in opposite directions, thus modifying solvent accessibility of the active site, and that the mixed alpha/3(10)-helix (residues 50-60) behaves as a mechanical hinge during the breathing motion of the protein. These data demonstrate that this motion, essential for RNase A substrate binding and release, is an intrinsic dynamical property of the ligand-free enzyme.     

Subtle functional collective motions in pancreatic-like ribonucleases: from ribonuclease A to angiogenin

The analysis of the dynamic behavior of enzymes is fundamental to structural biology. A direct relationship between protein flexibility and biological function has been shown for bovine pancreatic ribonuclease (RNase A) (Rasmussen et al., Nature 1992;357:423-424). More recently, crystallographic studies have shown that functional motions in RNase A involve the enzyme beta-sheet regions that move concertedly on substrate binding and release (Vitagliano et al., Proteins 2002;46:97-104). These motions have been shown to correspond to intrinsic dynamic properties of the native enzyme by molecular dynamics (MD) simulations. To unveil the occurrence of these collective motions in other members of pancreatic-like superfamily, we carried out MD simulations on human angiogenin (Ang). Essential dynamics (ED) analyses performed on the trajectories reveal that Ang exhibits collective motions similar to RNase A, despite the limited sequence identity (33%) of the two proteins. Furthermore, we show that these collective motions are also present in ensembles of experimentally determined structures of both Ang and RNase A. Finally, these subtle concerted beta-sheet motions were also observed for other two members of the pancreatic-like superfamily by comparing the ligand-bound and ligand-free structures of these enzymes. Taken together, these findings suggest that pancreatic-like ribonucleases share an evolutionary conserved dynamic behavior consisting of subtle beta-sheet motions, which are essential for substrate binding and release.     

Eukaryotic ribonucleases P/MRP: the crystal structure of the P3 domain

Ribonuclease (RNase) P is a site‐specific endoribonuclease found in all kingdoms of life. Typical RNase P consists of a catalytic RNA component and a protein moiety. In the eukaryotes, the RNase P lineage has split into two, giving rise to a closely related enzyme, RNase MRP, which has similar components but has evolved to have different specificities. The eukaryotic RNases P/MRP have acquired an essential helix‐loop‐helix protein‐binding RNA domain P3 that has an important function in eukaryotic enzymes and distinguishes them from bacterial and archaeal RNases P. Here, we present a crystal structure of the P3 RNA domain from Saccharomyces cerevisiae RNase MRP in a complex with RNase P/MRP proteins Pop6 and Pop7 solved to 2.7 Å. The structure suggests similar structural organization of the P3 RNA domains in RNases P/MRP and possible functions of the P3 domains and proteins bound to them in the stabilization of the holoenzymes' structures as well as in interactions with substrates. It provides the first insight into the structural organization of the eukaryotic enzymes of the RNase P/MRP family.     

Classification and annotation of the relationship between protein structural change and ligand binding

The causal relationship between protein structural change and ligand binding was classified and annotated for 839 nonredundant pairs of crystal structures in the Protein Data Bank—one with and the other without a bound low-molecular-weight ligand molecule. Protein structural changes were first classified into either domain or local motions depending on the size of the moving protein segments. Whether the protein motion was coupled with ligand binding was then evaluated based on the location of the ligand binding site and by application of the linear response theory of protein structural change. Protein motions coupled with ligand binding were further classified into either closure or opening motions. This classification revealed the following: (i) domain motions coupled with ligand binding are dominated by closure motions, which can be described by the linear response theory; (ii) local motions frequently accompany order-disorder or α-helix-coil conformational transitions; and (iii) transferase activity (Enzyme Commission   number 2) is the predominant function among coupled domain closure motions. This could be explained by the closure motion acting to insulate the reaction site of these enzymes from environmental water.     

Interactions of a Pop5/Rpp1 heterodimer with the catalytic domain of RNase MRP

Ribonuclease (RNase) MRP is a multicomponent ribonucleoprotein complex closely related to RNase P. RNase MRP and eukaryotic RNase P share most of their protein components, as well as multiple features of their catalytic RNA moieties, but have distinct substrate specificities. While RNase P is practically universally found in all three domains of life, RNase MRP is essential in eukaryotes. The structural organizations of eukaryotic RNase P and RNase MRP are poorly understood. Here, we show that Pop5 and Rpp1, protein components found in both RNase P and RNase MRP, form a heterodimer that binds directly to the conserved area of the putative catalytic domain of RNase MRP RNA. The Pop5/Rpp1 binding site corresponds to the protein binding site in bacterial RNase P RNA. Structural and evolutionary roles of the Pop5/Rpp1 heterodimer in RNases P and MRP are discussed.     

Association and dissociation kinetics of colicin E3 and immunity protein 3: convergence of theory and experiment

The rapid binding of cytotoxic colicin E3 by its cognate immunity protein Im3 is essential in safeguarding the producing cell. The X-ray structure of the E3/Im3 complex shows that the Im3 molecule interfaces with both the C-terminal ribonuclease (RNase) domain and the N-terminal translocation domain of E3. The association and dissociation rates of the RNase domain and Im3 show drastically different sensitivities to ionic strength, as previously rationalized for electrostatically enhanced diffusion-limited protein-protein associations. Relative to binding to the RNase domain, binding to full-length E3 shows a comparable association rate but a significantly lower dissociation rate. This outcome is just what was anticipated by a theory for the binding of two linked domains to a protein. The E3/Im3 system thus provides a powerful paradigm for the interplay of theory and experiment.     

Finding biologically relevant protein domain interactions: conserved binding mode analysis

Proteins evolved through the shuffling of functional domains, and therefore, the same domain can be found in different proteins and species. Interactions between such conserved domains often involve specific, well-determined binding surfaces reflecting their important biological role in a cell. To find biologically relevant interactions we developed a method of systematically comparing and classifying protein domain interactions from the structural data. As a result, a set of conserved binding modes (CBMs) was created using the atomic detail of structure alignment data and the protein domain classification of the Conserved Domain Database. A conserved binding mode is inferred when different members of interacting domain families dock in the same way, such that their structural complexes superimpose well. Such domain interactions with recurring structural themes have greater significance to be biologically relevant, unlike spurious crystal packing interactions. Consequently, this study gives lower and upper bounds on the number of different types of interacting domain pairs in the structure database on the order of 1000-2000. We use CBMs to create domain interaction networks, which highlight functionally significant connections by avoiding many infrequent links between highly connected nodes. The CBMs also constitute a library of docking templates that may be used in molecular modeling to infer the characteristics of an unknown binding surface, just as conserved domains may be used to infer the structure of an unknown protein. The method's ability to sort through and classify large numbers of putative interacting domain pairs is demonstrated on the oligomeric interactions of globins.     

Structural plasticity and Mg2+ binding properties of RNase P P4 from combined analysis of NMR residual dipolar couplings and motionally decoupled spin relaxation

The P4 helix is an essential element of ribonuclease P (RNase P) that is believed to bind catalytically important metals. Here, we applied a combination of NMR residual dipolar couplings (RDCs) and a recently introduced domain-elongation strategy for measuring "motionally decoupled" relaxation data to characterize the structural dynamics of the P4 helix from Bacillus subtilis RNase P. In the absence of divalent ions, the two P4 helical domains undergo small amplitude (approximately 13 degrees) collective motions about an average interhelical angle of 10 degrees. The highly conserved U7 bulge and helical residue C8, which are proposed to be important for substrate recognition and metal binding, are locally mobile at pico- to nanosecond timescales and together form the pivot point for the collective domain motions. Chemical shift mapping reveals significant association of Mg2+ ions at the P4 major groove near the flexible pivot point at residues (A5, G22, G23) previously identified to bind catalytically important metals. The Mg2+ ions do not, however, significantly alter the structure or dynamics of P4. Analysis of results in the context of available X-ray structures of the RNA component of RNase P and structural models that include the pre-tRNA substrate suggest that the internal motions observed in P4 likely facilitate adaptive changes in conformation that take place during folding and substrate recognition, possibly aided by interactions with Mg2+ ions. Our results add to a growing view supporting the existence of functionally important internal motions in RNA occurring at nanosecond timescales.     

Modular construction for function of a ribonucleoprotein enzyme: the catalytic domain of Bacillus subtilis RNase P complexed with B. subtilis RNase P protein

The bacterial RNase P holoenzyme catalyzes the formation of the mature 5′-end of tRNAs and is composed of an RNA and a protein subunit. Among the two folding domains of the RNase P RNA, the catalytic domain (C-domain) contains the active site of this ribozyme. We investigated specific binding of the Bacillus subtilis C-domain with the B.subtilis RNase P protein and examined the catalytic activity of this C-domain–P protein complex. The C-domain forms a specific complex with the P protein with a binding constant of ~0.1 µM. The C-domain–P protein complex and the holoenzyme are equally efficient in cleaving single-stranded RNA (~0.9 min^–1 at pH 7.8) and substrates with a hairpin–loop 3′ to the cleavage site (~40 min^–1). The holoenzyme reaction is much more efficient with a pre-tRNA substrate, binding at least 100-fold better and cleaving 10–500 times more efficiently. These results demonstrate that the RNase P holoenzyme is functionally constructed in three parts. The catalytic domain alone contains the active site, but has little specificity and affinity for most substrates. The specificity and affinity for the substrate is generated by either the specificity domain of RNase P RNA binding to a T stem–loop-like hairpin or RNase P protein binding to a single-stranded RNA. This modular construction may be exploited to obtain RNase P-based ribonucleoprotein complexes with altered substrate specificity.     

Crystallographic and modeling studies of RNase III suggest a mechanism for double-stranded RNA cleavage.

BACKGROUND: Aquifex aeolicus Ribonuclease III (Aa-RNase III) belongs to the family of Mg(2+)-dependent endonucleases that show specificity for double-stranded RNA (dsRNA). RNase III is conserved in all known bacteria and eukaryotes and has 1-2 copies of a 9-residue consensus sequence, known as the RNase III signature motif. The bacterial RNase III proteins are the simplest, consisting of two domains: an N-terminal endonuclease domain, followed by a double-stranded RNA binding domain (dsRBD). The three-dimensional structure of the dsRBD in Escherichia coli RNase III has been elucidated; no structural information is available for the endonuclease domain of any RNase III. RESULTS: We present the crystal structures of the Aa-RNase III endonuclease domain in its ligand-free form and in complex with Mn(2+). The structures reveal a novel protein fold and suggest a mechanism for dsRNA cleavage. On the basis of structural, genetic, and biological data, we have constructed a hypothetical model of Aa-RNase III in complex with dsRNA and Mg(2+) ion, which provides the first glimpse of RNase III in action. CONCLUSIONS: The functional Aa-RNase III dimer is formed via mainly hydrophobic interactions, including a "ball-and-socket" junction that ensures accurate alignment of the two monomers. The fold of the polypeptide chain and its dimerization create a valley with two compound active centers at each end of the valley. The valley can accommodate a dsRNA substrate. Mn(2+) binding has significant impact on crystal packing, intermolecular interactions, thermal stability, and the formation of two RNA-cutting sites within each compound active center.     

Swapping the domains of exoribonucleases RNase II and RNase R: conferring upon RNase II the ability to degrade ds RNA

RNase II and RNase R are the two E. coli exoribonucleases that belong to the RNase II super family of enzymes. They degrade RNA hydrolytically in the 3' to 5' direction in a processive and sequence independent manner. However, while RNase R is capable of degrading structured RNAs, the RNase II activity is impaired by dsRNAs. The final end-product of these two enzymes is also different, being 4 nt for RNase II and 2 nt for RNase R. RNase II and RNase R share structural properties, including 60% of amino acid sequence similarity and have a similar modular domain organization: two N-terminal cold shock domains (CSD1 and CSD2), one central RNB catalytic domain, and one C-terminal S1 domain. We have constructed hybrid proteins by swapping the domains between RNase II and RNase R to determine which are the responsible for the differences observed between RNase R and RNase II. The results obtained show that the S1 and RNB domains from RNase R in an RNase II context allow the degradation of double-stranded substrates and the appearance of the 2 nt long end-product. Moreover, the degradation of structured RNAs becomes tail-independent when the RNB domain from RNase R is no longer associated with the RNA binding domains (CSD and S1) of the genuine protein. Finally, we show that the RNase R C-terminal Lysine-rich region is involved in the degradation of double-stranded substrates in an RNase II context, probably by unwinding the substrate before it enters into the catalytic cavity.     

Dynamic properties of the N-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two 3-dimensional domain-swapped dimers with different quaternary structures. One dimer is characterized by the swapping of the C-terminal region (C-Dimer) and presents a rather loose structure. The other dimer (N-Dimer) exhibits a very compact structure with exchange of the N-terminal helix. Here we report the results of a molecular dynamics/essential dynamics (MD/ED) study carried out on the N-Dimer. This investigation, which represents the first MD/ED analysis on a three-dimensional domain-swapped enzyme, provides information on the dynamic properties of the active site residues as well as on the global motions of the dimer subunits. In particular, the analysis of the flexibility of the active site residues agrees well with recent crystallographic and site-directed mutagenesis studies on monomeric RNase A, thus indicating that domain swapping does not affect the dynamics of the active sites. A slight but significant rearrangement of N-Dimer quaternary structure, favored by the formation of additional hydrogen bonds at subunit interface, has been observed during the MD simulation. The analysis of collective movements reveals that each subunit of the dimer retains the functional breathing motion observed for RNase A. Interestingly, the breathing motion of the two subunits is dynamically coupled, as they open and close in phase. These correlated motions indicate the presence of active site intercommunications in this dimer. On these bases, we propose a speculative mechanism that may explain negative cooperativity in systems preserving structural symmetry during the allosteric transitions.     

Specific binding of a Pop6/Pop7 heterodimer to the P3 stem of the yeast RNase MRP and RNase P RNAs

Pop6 and Pop7 are protein subunits of Saccharomyces cerevisiae RNase MRP and RNase P. Here we show that bacterially expressed Pop6 and Pop7 form a soluble heterodimer that binds the RNA components of both RNase MRP and RNase P. Footprint analysis of the interaction between the Pop6/7 heterodimer and the RNase MRP RNA, combined with gel mobility assays, demonstrates that the Pop6/7 complex binds to a conserved region of the P3 domain. Binding of these proteins to the MRP RNA leads to local rearrangement in the structure of the P3 loop and suggests that direct interaction of the Pop6/7 complex with the P3 domain of the RNA components of RNases MRP and P may mediate binding of other protein components. These results suggest a role for a key element in the RNase MRP and RNase P RNAs in protein binding, and demonstrate the feasibility of directly studying RNA-protein interactions in the eukaryotic RNases MRP and P complexes.     

A comparison of BRCT domains involved in nonhomologous end-joining: introducing the solution structure of the BRCT domain of polymerase lambda

Three of the four family X polymerases, DNA polymerase lambda, DNA polymerase mu, and TdT, have been associated with repair of double-strand DNA breaks by nonhomologous end-joining. Their involvement in this DNA repair process requires an N-terminal BRCT domain that mediates interaction with other protein factors required for recognition and binding of broken DNA ends. Here we present the NMR solution structure of the BRCT domain of DNA polymerase lambda, completing the structural portrait for this family of enzymes. Analysis of the overall fold of the polymerase lambda BRCT domain reveals structural similarity to the BRCT domains of polymerase mu and TdT, yet highlights some key sequence and structural differences that may account for important differences in the biological activities of these enzymes and their roles in nonhomologous end-joining. Mutagenesis studies indicate that the conserved Arg57 residue of Pol lambda plays a more critical role for binding to the XRCC4-Ligase IV complex than its structural homolog in Pol mu, Arg43. In contrast, the hydrophobic Leu60 residue of Pol lambda contributes less significantly to binding than the structurally homologous Phe46 residue of Pol mu. A third leucine residue involved in the binding and activity of Pol mu, is nonconservatively replaced by a glutamine in Pol lambda (Gln64) and, based on binding and activity data, is apparently unimportant for Pol lambda interactions with the NHEJ complex. In conclusion, both the structure of the Pol lambda BRCT domain and its mode of interaction with the other components of the NHEJ complex significantly differ from the two previously studied homologs, Pol mu and TdT.     

Ternary borate–nucleoside complex stabilization by ribonuclease A demonstrates phosphate mimicry

Phosphate esters play a central role in cellular energetics, biochemical activation, signal transduction and conformational switching. The structural homology of the borate anion with phosphate, combined with its ability to spontaneously esterify hydroxyl groups, suggested that phosphate ester recognition sites on proteins might exhibit significant affinity for nonenzymatically formed borate esters. ¹¹B NMR studies and activity measurements on ribonuclease A (RNase A) in the presence of borate and several cytidine analogs demonstrate the formation of a stable ternary RNase A·3′-deoxycytidine–2′-borate ternary complex that mimics the complex formed between RNase A and a 2′-cytidine monophosphate (2′-CMP) inhibitor. Alternatively, no slowly exchanging borate resonance is observed for a ternary RNase A, borate, 2′-deoxycytidine mixture, demonstrating the critical importance of the 2′-hydroxyl group for complex formation. Titration of the ternary complex with 2′-CMP shows that it can displace the bound borate ester with a binding constant that is close to the reported inhibition constant of RNase A by 2′-CMP. RNase A binding of a cyclic cytidine-2′,3′-borate ester, which is a structural homolog of the cytidine-2′,3′-cyclic phosphate substrate, could also be demonstrated. The apparent dissociation constant for the cytidine-2′,3′-borate·RNase A complex is 0.8 mM, which compares with a Michaelis constant of 11 mM for cytidine-2′,3′-cyclic phosphate at pH 7, indicating considerably stronger binding. However, the value is 1,000-fold larger than the reported dissociation constant of the RNase A complex with uridine–vanadate. These results are consistent with recent reports suggesting that in situ formation of borate esters that mimic the corresponding phosphate esters supports enzyme catalysis.     

Structural plasticity of peptidyl-prolyl isomerase sFkpA is a key to its chaperone function as revealed by solution NMR

Intramolecular dynamics of periplasmic chaperone FkpA-deltaCT (sFkpA) and its complexes with partially structured substrates are studied by NMR in solution. The backbone amide 15N relaxation of sFkpA reveals flexibility in the relative orientation between the dimerization domain and two juxtaposed catalytic domains identified in the X-ray structure of sFkpA. This flexibility is attributed to the structural plasticity within the long alpha-helical arm (helix III) consisting of residues 84 and 91. Residual dipolar couplings (RDCs) indicate an absence of fixed orientation between the sFkpA domains. The substrate binding surface of sFkpA is defined on the X-ray structure by mapping of chemical shift perturbations introduced by complexation of sFkpA with its corresponding protein substrates: partially folded RNase A S-protein and reduced carboxymethylated bovine alpha-lactalbumin (RCM-la). A comparison of 15N relaxation of apo-sFkpA and its complex with RNase A S-protein indicates an increased rigidity within the long alpha-helix III and decreased interdomain mobility of the complex. We speculate that these dynamic properties may play a key role in the chaperone activity of sFkpA, since ability to bind different substrates potentially requires structural adaptations of the chaperone protein. We show that binding of sFkpA to RNase A S-protein greatly reduces the population of aggregated oligomeric species of RNase A S-protein. Finally, a molecular model, the so-called "mother's arms" model, is proposed to illustrate the mechanism of chaperone activity by FkpA.     

Cancer-related mutations in BRCA1-BRCT cause long-range structural changes in protein-protein binding sites: a molecular dynamics study

Cancer-associated mutations in the BRCT domain of BRCA1 (BRCA1-BRCT) abolish its tumor suppressor function by disrupting interactions with other proteins such as BACH1. Many cancer-related mutations do not cause sufficient destabilization to lead to global unfolding under physiological conditions, and thus abrogation of function probably is due to localized structural changes. To explore the reasons for mutation-induced loss of function, the authors performed molecular dynamics simulations on three cancer-associated mutants, A1708E, M1775R, and Y1853ter, and on the wild type and benign M1652I mutant, and compared the structures and fluctuations. Only the cancer-associated mutants exhibited significant backbone structure differences from the wild-type crystal structure in BACH1-binding regions, some of which are far from the mutation sites. Backbone differences of the A1708E mutant from the liganded wild type structure in these regions are much larger than those of the unliganded wild type X-ray or molecular dynamics structures. These BACH1-binding regions of the cancer-associated mutants also exhibited increases in their fluctuation magnitudes compared with the same regions in the wild type and M1562I mutant, as quantified by quasiharmonic analysis. Several of the regions of increased fluctuation magnitude correspond to correlated motions of residues in contact that provide a continuous path of fluctuating amino acids in contact from the A1708E and Y1853ter mutation sites to the BACH1-binding sites with altered structure and dynamics. The increased fluctuations in the disease-related mutants suggest an increase in vibrational entropy in the unliganded state that could result in a larger entropy loss in the disease-related mutants upon binding BACH1 than in the wild type. To investigate this possibility, vibrational entropies of the A1708E and wild type in the free state and bound to a BACH1-derived phosphopeptide were calculated using quasiharmonic analysis, to determine the binding entropy difference DeltaDeltaS between the A1708E mutant and the wild type. DeltaDeltaS was determined to be -4.0 cal mol(-1) K(-1), with an uncertainty of 2 cal mol(-1) K(-1); that is, the entropy loss upon binding the peptide is 4.0 cal mol(-1) K(-1) greater for the A1708E mutant, corresponding to an entropic contribution to the DeltaDeltaG of binding (-TDeltaDeltaS) 1.1 kcal mol(-1) more positive for the mutant. The observed differences in structure, flexibility, and entropy of binding likely are responsible for abolition of BACH1 binding, and illustrate that many disease- related mutations could have very long-range effects. The methods described here have potential for identifying correlated motions responsible for other long-range effects of deleterious mutations.     

2-5A induces a conformational change in the ankyrin-repeat domain of RNase L

RNase L is responsible for the 2-5A host defense system, an RNA degradation pathway present in cells of higher vertebrates that functions in both the antiviral and anticellular activities of interferon. The activity of RNase L is tightly regulated and is exerted only in the presence of 2-5A. The postulated mechanism of its regulation is as follows: the N-terminal half ankyrin-repeat domain masks the C-terminal half nuclease domain in the absence of 2-5A. On binding 2-5A at the ankyrin-repeat domain, RNase L forms a homodimer and removes the ankyrin-repeat domain from the nuclease domain to become the active form. A conformational change in the ankyrin-repeat domain is a key step in this hypothetical mechanism, but there is as yet no evidence for such a change. To clarify the events induced by 2-5A binding, we established procedures for expression and purification of the ankyrin-repeat domain of human RNase L. Fluorescence spectra of the protein showed clear difference in the presence and absence of 2-5A. The alterations in the spectra supported conformational changes of the protein. Time-resolved anisotropy measurements indicated that 2-5A binding led to a significant decrease in the rotational radius of the protein. In addition, 2-5A provided the domain with resistance to protease digestion as a result of a conformational change. These results indicated that the ankyrin-repeat domain of RNase L constricts its structure by binding of 2-5A. This observation suggests a revised model of the 2-5A-induced activation of RNase L.