Ionic interactions in crystalline bovine pancreatic ribonuclease A

Protein Thermodynamic Properties, Crystallisation, and the Hofmeister Series

The Hofmeister series is a series of ions ordered according to their ability to precipitate proteins. It has also been linked to a host of (bio)chemical phenomena. Several attempts over the years to correlate the series to the varying success of different salts in crystallising proteins have been largely inconclusive. A correlation, based on published data and crystallisation conditions for several proteins, is proposed here between some thermodynamic properties of proteins and the position in the Hofmeister series of the salts from which they preferentially crystallise. Namely, a high ratio between the entropic or enthalpic protein-solvent interactions contribution to thermodynamic stability and the total thermodynamic stability of a given protein, indicate the protein's high propensity to crystallise in solutions of highly kosmotropic salts. Low such ratios on the other hand, indicate that chaotropic salts can be equally successful, i. e. that the protein in question is rather indifferent to the Hofmeister character of the salt. Testing various model proteins for crystallisation against screens containing salts found at different points on the Hofmeister series, as well as further bibliographic analysis, have yielded results that appear to largely corroborate this hypothesis. These conclusions may conceivably be used as a crystallisation predictive tool.

pH-dependent reaction triggering in PmHMGR crystals for time-resolved crystallography

Serial crystallography and time-resolved data collection can readily be employed to investigate the catalytic mechanism of Pseudomonas mevalonii 3-hydroxy-3-methylglutaryl (HMG)-coenzyme-A (CoA) reductase (PmHMGR) by changing the environmental conditions in the crystal and so manipulating the reaction rate. This enzyme uses a complex mechanism to convert mevalonate to HMG-CoA using the co-substrate CoA and cofactor NAD+. The multi-step reaction mechanism involves an exchange of bound NAD+ and large conformational changes by a 50-residue subdomain. The enzymatic reaction can be run in both forward and reverse directions in solution and is catalytically active in the crystal for multiple reaction steps. Initially, the enzyme was found to be inactive in the crystal starting with bound mevalonate, CoA, and NAD+. To observe the reaction from this direction, we examined the effects of crystallization buffer constituents and pH on enzyme turnover, discovering a strong inhibition in the crystallization buffer and a controllable increase in enzyme turnover as a function of pH. The inhibition is dependent on ionic concentration of the crystallization precipitant ammonium sulfate but independent of its ionic composition. Crystallographic studies show that the observed inhibition only affects the oxidation of mevalonate but not the subsequent reactions of the intermediate mevaldehyde. Calculations of the pKa values for the enzyme active site residues suggest that the effect of pH on turnover is due to the changing protonation state of His381. We have now exploited the changes in ionic inhibition in combination with the pH-dependent increase in turnover as a novel approach for triggering the PmHMGR reaction in crystals and capturing information about its intermediate states along the reaction pathway.

Cross-Linked Crystals of Dirhodium Tetraacetate/RNase A Adduct Can Be Used as Heterogeneous Catalysts

Due to their unique coordination structure, dirhodium paddlewheel complexes are of interest in several research fields, like medicinal chemistry, catalysis, etc. Previously, these complexes were conjugated to proteins and peptides for developing artificial metalloenzymes as homogeneous catalysts. Fixation of dirhodium complexes into protein crystals is interesting to develop heterogeneous catalysts. Porous solvent channels present in protein crystals can benefit the activity by increasing the probability of substrate collisions at the catalytic Rh binding sites. Toward this goal, the present work describes the use of bovine pancreatic ribonuclease (RNase A) crystals with a pore size of 4 nm (P3₂21 space group) for fixing [Rh₂(OAc)₄] and developing a heterogeneous catalyst to perform reactions in an aqueous medium. The structure of the [Rh₂(OAc)₄]/RNase A adduct was investigated by X-ray crystallography: the metal complex structure remains unperturbed upon protein binding. Using a number of crystal structures, metal complex accumulation over time, within the RNase A crystals, and structures at variable temperatures were evaluated. We also report the large-scale preparation of microcrystals (∼10-20 μm) of the [Rh₂(OAc)₄]/RNase A adduct and cross-linking reaction with glutaraldehyde. The catalytic olefin cyclopropanation reaction and self-coupling of diazo compounds by these cross-linked [Rh₂(OAc)₄]/RNase A crystals were demonstrated. The results of this work reveal that these systems can be used as heterogeneous catalysts to promote reactions in aqueous solution. Overall, our findings demonstrate that the dirhodium paddlewheel complexes can be fixed in porous biomolecule crystals, like those of RNase A, to prepare biohybrid materials for catalytic applications.

Acetonitrile allows indirect replacement of nondeuterated lipid detergents by deuterated lipid detergents for the nuclear magnetic resonance study of detergent-soluble proteins

Detergent-soluble proteins (DSPs) are commonly dissolved in lipid buffers for NMR experiments, but the huge lipid proton signal prevents recording of high-quality spectra. The use of costly deuterated lipids is thus required to replace nondeuterated ones. With conventional methods, detergents like dodecylphosphocholine (DPC) cannot be fully exchanged due to their high binding affinity to hydrophobic proteins. We propose an original and simple protocol which combines the use of acetonitrile, dialysis and lyophilization to disrupt the binding of lipids to the protein and allow their indirect replacement by their deuterated equivalents, while maintaining the native structure of the protein. Moreover, by this protocol, the detergent-to-protein molar ratio can be controlled as it challenges the protein structure. This protocol was applied to solubilize the Vpx protein that was followed upon addition of DPC-d₃₈ by ¹ H-¹⁵ N SOFAST-HMQC spectra and the best detergent-to-DSPs molar ratio was obtained for structural studies.

Do soft anions promote protein denaturation through binding interactions? A case study using ribonuclease A

It has long been known that large soft anions like bromide, iodide and thiocyanate are protein denaturing agents, but their mechanism of action is still unclear. In this work we have investigated the protein denaturing properties of these anions using Ribonuclease A (RNase A) as a model protein system. Salt-induced perturbations to the protein folding free energy were determined using differential scanning calorimetry and the results demonstrate that the addition of sodium iodide and sodium thiocyanate significantly decreases the melting temperature of the protein. In order to account for this reduction in protein stability, we show that the introduction of salts that contain soft anions to the aqueous solvent perturbs the protein unfolding free energy through three mechanisms: (a) screening Coulomb interactions that exist between charged protein residues, (b) Hofmeister effects, and (c) specific anion binding to CH and CH2 moieties in the protein polypeptide backbone. Using the micellization of 1,2-hexanediol as a ruler for hydrophobicity, we have devised a practical methodology that separates the Coulomb and Hofmeister contributions of salts to the protein unfolding free energy. This allowing us to isolate the contribution of soft anion binding interactions to the unfolding process. The analysis shows that binding contributions have the largest magnitude, confirming that it is the binding of soft anions to the polypeptide backbone that is the main promoter of protein unfolding.

The ammonium sulfate inhibition of human angiogenin

In this study, we investigate the inhibition of human angiogenin by ammonium sulfate. The inhibitory potency of ammonium sulfate for human angiogenin (IC50 = 123.5 ± 14.9 mm) is comparable to that previously reported for RNase A (119.0 ± 6.5 mm) and RNase 2 (95.7 ± 9.3 mm). However, analysis of two X-ray crystal structures of human angiogenin in complex with sulfate anions (in acidic and basic pH environments, respectively) indicates an entirely distinct mechanism of inhibition. While ammonium sulfate inhibits the ribonucleolytic activity of RNase A and RNase 2 by binding to the active site of these enzymes, sulfate anions bind only to peripheral substrate anion-binding subsites of human angiogenin, and not to the active site.

The first crystal structure of human RNase 6 reveals a novel substrate-binding and cleavage site arrangement

We describe the first human RNase 6 crystal structure in complex with sulfate anions. Kinetic analysis, site-directed mutagenesis and molecular dynamics simulations identified novel substrate recognition and cleavage sites.

Human RNase 6 is a cationic secreted protein that belongs to the RNase A superfamily. Its expression is induced in neutrophils and monocytes upon bacterial infection, suggesting a role in host defence. We present here the crystal structure of RNase 6 obtained at 1.72 Å (1 Å=0.1 nm) resolution, which is the first report for the protein 3D structure and thereby setting the basis for functional studies. The structure shows an overall kidney-shaped globular fold shared with the other known family members. Three sulfate anions bound to RNase 6 were found, interacting with residues at the main active site (His¹⁵, His¹²² and Gln¹⁴) and cationic surface-exposed residues (His³⁶, His³⁹, Arg⁶⁶ and His⁶⁷). Kinetic characterization, together with prediction of protein–nucleotide complexes by molecular dynamics, was applied to analyse the RNase 6 substrate nitrogenous base and phosphate selectivity. Our results reveal that, although RNase 6 is a moderate catalyst in comparison with the pancreatic RNase type, its structure includes lineage-specific features that facilitate its activity towards polymeric nucleotide substrates. In particular, enzyme interactions at the substrate 5′ end can provide an endonuclease-type cleavage pattern. Interestingly, the RNase 6 crystal structure revealed a novel secondary active site conformed by the His³⁶–His³⁹ dyad that facilitates the polynucleotide substrate catalysis.

Platinated oligomers of bovine pancreatic ribonuclease: Structure and stability

The reaction between cis-diamminedichloroplatinum(II) (CDDP), cisplatin, a common anticancer drug, and bovine pancreatic ribonuclease (RNase A), induces extensive protein aggregation, leading to the formation of one dimer, one trimer and higher oligomers whose yields depend on cisplatin/protein ratio. Structural and functional properties of the purified platinated species, together with their spontaneous dissociation and thermally induced denaturation, have been characterized. Platinated species preserve a significant, although reduced, ribonuclease activity. The high resistance of the dimers against dissociation and the different thermal unfolding profiles suggest a quaternary structure different from those of the well-known swapped dimers of RNase A.

What's in your buffer? Solute altered millisecond motions detected by solution NMR

To date, little work has been conducted on the relationship between solute and buffer molecules and conformational exchange motion in enzymes. This study uses solution NMR to examine the effects of phosphate, sulfate, and acetate in comparison to MES- and HEPES-buffered references on the chemical shift perturbation and millisecond, chemical, or conformational exchange motions in the enzyme ribonuclease A (RNase A), triosephosphate isomerase (TIM) and HisF. The results indicate that addition of these solutes has a small effect on (1)H and (15)N chemical shifts for RNase A and TIM but a significant effect for HisF. For RNase A and TIM, Carr-Purcell-Meiboom-Gill relaxation dispersion experiments, however, show significant solute-dependent changes in conformational exchange motions. Some residues show loss of millisecond motions relative to the reference sample upon addition of solute, whereas others experience an enhancement. Comparison of exchange parameters obtained from fits of dispersion data indicates changes in either or both equilibrium populations and chemical shifts between conformations. Furthermore, the exchange kinetics are altered in many cases. The results demonstrate that common solute molecules can alter observed enzyme millisecond motions and play a more active role than what is routinely believed.

Protein flexibility in docking and surface mapping

Structure-based drug design has become an essential tool for rapid lead discovery and optimization. As available structural information has increased, researchers have become increasingly aware of the importance of protein flexibility for accurate description of the native state. Typical protein-ligand docking efforts still rely on a single rigid receptor, which is an incomplete representation of potential binding conformations of the protein. These rigid docking efforts typically show the best performance rates between 50 and 75%, while fully flexible docking methods can enhance pose prediction up to 80-95%. This review examines the current toolbox for flexible protein-ligand docking and receptor surface mapping. Present limitations and possibilities for future development are discussed.

Histidine in continuum electrostatics protonation state calculations

A modification to the standard continuum electrostatics approach to calculate protein pK(a)s, which allows for the decoupling of histidine tautomers within a two-state model, is presented. Histidine with four intrinsically coupled protonation states cannot be easily incorporated into a two-state formalism, because the interaction between the two protonatable sites of the imidazole ring is not purely electrostatic. The presented treatment, based on a single approximation of the interrelation between histidine's charge states, allows for a natural separation of the two protonatable sites associated with the imidazole ring as well as the inclusion of all protonation states within the calculation.

The crystal structure of ribonuclease A in complex with thymidine-3'-monophosphate provides further insight into ligand binding

Thymidine-3'-monophosphate (3'-TMP) is a competitive inhibitor analogue of the 3'-CMP and 3'-UMP natural product inhibitors of bovine pancreatic ribonuclease A (RNase A). Isothermal titration calorimetry experiments show that 3'-TMP binds the enzyme with a dissociation constant (K(d)) of 15 microM making it one of the strongest binding members of the five natural bases found in nucleic acids (A, C, G, T, and U). To further investigate the molecular properties of this potent natural affinity, we have determined the crystal structure of bovine pancreatic RNase A in complex with 3'-TMP at 1.55 A resolution and we have performed NMR binding experiments with 3'-CMP and 3'-TMP. Our results show that binding of 3'-TMP is very similar to other natural and non-natural pyrimidine ligands, demonstrating that single nucleotide affinity is independent of the presence or absence of a 2'-hydroxyl on the ribose moiety of pyrimidines and suggesting that the pyrimidine binding subsite of RNase A is not a significant contributor of inhibitor discrimination. Accumulating evidence suggests that very subtle structural, chemical, and potentially motional variations contribute to ligand discrimination in this enzyme.

Structure of bovine pancreatic ribonuclease complexed with uridine 5'-monophosphate at 1.60 A resolution

Bovine pancreatic ribonuclease A (RNase A) was crystallized from a mixture of small molecules containing basic fuchsin, tobramycin and uridine 5'-monophosphate (U5P). Solution of the crystal structure revealed that the enzyme was selectively bound to U5P, with the pyrimidine ring of U5P residing in the pyrimidine-binding site at Thr45. The structure was refined to an R factor of 0.197 and an R(free) of 0.253.

Selectivity in the post-translational, transglutaminase-dependent acylation of lysine residues

Transglutaminases (TGs) are known to exhibit remarkable specificities not only for the Q (or Gln) sites but also for the K (or Lys) sites of proteins with which they react. To gain further insight into K-site specificity, we examined the reactions of dansyl-epsilon-aminocaproyl-GlnGlnIleVal with three chemically and structurally well-characterized proteins (bovine pancreatic ribonuclease A, bovine pancreatic trypsin inhibitor, and chicken egg white lysozyme), as catalyzed by TG2, a biologically important post-translational enzyme. The substrates represent a total of 20 potential surface sites for acylation by the fluorescent Gln probe, yet only two of the lysine side chains reacted with TG2. While the K1 site of ribonuclease and the K15 site of the trypsin inhibitor could be readily acylated by the enzyme, none of the lysines in lysozyme were modified. The findings lead us to suggest that the selection of lysine residues by TG2 is not encoded in the primary amino acid sequence surrounding the target side chain but depends primarily on its being positioned in an accessible segment of the protein structure.

Ribonuclease A homologues of the zebrafish: polymorphism, crystal structures of two representatives and their evolutionary implications

The widespread and functionally varied members of the ribonuclease A (RNase A) superfamily provide an excellent opportunity to study evolutionary forces at work on a conserved protein scaffold. Representatives from the zebrafish are of particular interest as the evolutionary distance from non-ichthyic homologues is large. We conducted an exhaustive survey of available zebrafish DNA sequences and found significant polymorphism among its four known homologues. In an extension of previous nomenclature, the variants have been named RNases ZF-1a–c,-2a–d,-3a–e and-4. We present the first X-ray crystal structures of zebrafish ribonucleases, RNases ZF-1a and-3e at 1.35-and 1.85 Å resolution, respectively. Structure-based clustering with ten other ribonuclease structures indicates greatest similarity to mammalian angiogenins and amphibian ribonucleases, and supports the view that all present-day ribonucleases evolved from a progenitor with three disulphide bonds. In their details, the two structures are intriguing melting-pots of features present in ribonucleases from other vertebrate classes. Whereas in RNase ZF-1a the active site is obstructed by the C-terminal segment (as observed in angiogenin), in RNase ZF-3e the same region is open (as observed in more catalytically efficient homologues). The progenitor of present-day ribonucleases is more likely to have had an obstructive C terminus, and the relatively high similarity (late divergence) of RNases ZF-1 and-3 infers that the active site unblocking event has happened independently in different vertebrate lineages.

A new crystal form of bovine pancreatic RNase A in complex with 2'-deoxyguanosine-5'-monophosphate

The structure of bovine pancreatic RNase A has been determined in complex with 2'-deoxyguanosine-5'-monophosphate (dGMP) at 1.33 A resolution at room temperature in a previously unreported unit cell belonging to space group P3(1). There are two molecules of nucleotide per enzyme molecule, one of which lies in the active-site cleft in the productive binding mode, whilst the guanine base of the other dGMP occupies the pyrimidine-specific binding site in a nonproductive mode such that it forms hydrogen bonds to the phosphate group of the first dGMP. This is the first RNase A structure containing a guanine base in the B2 binding site. Each dGMP molecule is involved in intermolecular interactions with adjacent RNase A molecules in the lattice and the two nucleotides appear to direct the formation of the crystal lattice. Because GMP may be produced during degradation of RNA, this association could represent an inhibitor complex and thereby affect the observed enzyme kinetics.

A novel cross-linked RNase A dimer with enhanced enzymatic properties

A new cross-linked ribonuclease A (RNase A) dimer composed of monomeric units covalently linked by a single amide bond between the side-chains of Lys(66) and Glu(9) is described. The dimer was prepared in the absence of water by incubating a lyophilized preparation of RNase, sealed under vacuum, in an oven at 85 degrees C. It was determined that the in vacuo procedure does not induce any significant conformational changes to the overall structure of RNase A, yet the amide cross-link has an increased acid lability, indicating that it is exposed and conformationally strained. Examination of X-ray crystallographic structures indicates that Lys(66) and Glu(9) are not close enough for the in vacuo dimer to adopt any of the known domain-swapped conformations. Therefore, the in vacuo RNase A dimer appears to be a novel dimeric structure. The in vacuo RNase A dimer also exhibits a twofold increase in activity over monomeric RNase A on a per monomer basis. This doubling of enzymatic activity was shown using dsRNA and ssRNA as substrates. In addition to this enhanced ability to degrade RNA, the dimer is not inhibited by the cellular ribonuclease inhibitor protein (cRI).

The binding of 3′-N-piperidine-4-carboxyl-3′-deoxy-ara-uridine to ribonuclease A in the crystal

The binding of a moderate inhibitor, 3'-N-piperidine-4-carboxyl-3'-deoxy-ara-uridine, to ribonuclease A has been studied by X-ray crystallography at 1.7A resolution. Two inhibitor molecules are bound in the central RNA binding cavity of RNase A exploiting interactions with residues from peripheral binding sites rather than from the active site of the enzyme. The uracyl moiety of the first inhibitor molecule occupies the purine-preferring site of RNase A, while the rest of the molecule projects to the solvent. The second inhibitor molecule binds with the carboxyl group at the pyrimidine recognition site and the uridine moiety exploits interactions with RNase A residues Lys66, His119 and Asp121. Comparative structural analysis of the 3'-N-piperidine-4-carboxyl-3'-deoxy-ara-uridine complex with other RNase A-ligand complexes provides a structural explanation of its potency. The crystal structure of the RNase A-3'-N-piperidine-4-carboxyl-3'-deoxy-ara-uridine complex provides evidence of a novel ligand-binding pattern in RNase A for 3'-N-aminonucleosides that was not anticipated by modelling studies, while it also suggests ways to improve the efficiency and selectivity of such compounds to develop pharmaceuticals against pathologies associated with RNase A homologues.

Chemically accurate protein structures: validation of protein NMR structures by comparison of measured and predicted pKa values

A new method is presented for evaluating the quality of protein structures obtained by NMR. This method exploits the dependence between measurable chemical properties of a protein, namely pKa values of acidic residues, and protein structure. The accurate and fast empirical computational method employed by the PROPKA program ( http://www.propka.chem.uiowa.edu) allows the user to test the ability of a given structure to reproduce known pKa values, which in turn can be used as a criterion for the selection of more accurate structures. We demonstrate the feasibility of this novel idea for a series of proteins for which both NMR and X-ray structures, as well as pKa values of all ionizable residues, have been determined. For the 17 NMR ensembles used in this study, this criterion is shown effective in the elimination of a large number of NMR structure ensemble members.

Very fast empirical prediction and rationalization of protein pKa values

A very fast empirical method is presented for structure-based protein pKa prediction and rationalization. The desolvation effects and intra-protein interactions, which cause variations in pKa values of protein ionizable groups, are empirically related to the positions and chemical nature of the groups proximate to the pKa sites. A computer program is written to automatically predict pKa values based on these empirical relationships within a couple of seconds. Unusual pKa values at buried active sites, which are among the most interesting protein pKa values, are predicted very well with the empirical method. A test on 233 carboxyl, 12 cysteine, 45 histidine, and 24 lysine pKa values in various proteins shows a root-mean-square deviation (RMSD) of 0.89 from experimental values. Removal of the 29 pKa values that are upper or lower limits results in an RMSD = 0.79 for the remaining 285 pKa values.

The binding of IMP to ribonuclease A

The binding of inosine 5' phosphate (IMP) to ribonuclease A has been studied by kinetic and X-ray crystallographic experiments at high (1.5 A) resolution. IMP is a competitive inhibitor of the enzyme with respect to C>p and binds to the catalytic cleft by anchoring three IMP molecules in a novel binding mode. The three IMP molecules are connected to each other by hydrogen bond and van der Waals interactions and collectively occupy the B1R1P1B2P0P(-1) region of the ribonucleolytic active site. One of the IMP molecules binds with its nucleobase in the outskirts of the B2 subsite and interacts with Glu111 while its phosphoryl group binds in P1. Another IMP molecule binds by following the retro-binding mode previously observed only for guanosines with its nucleobase at B1 and the phosphoryl group in P(-1). The third IMP molecule binds in a novel mode towards the C-terminus. The RNase A-IMP complex provides structural evidence for the functional components of subsite P(-1) while it further supports the role inferred by other studies to Asn71 as the primary structural determinant for the adenine specificity of the B2 subsite. Comparative structural analysis of the IMP and AMP complexes highlights key aspects of the specificity of the base binding subsites of RNase A and provides a structural explanation for their potencies. The binding of IMP suggests ways to develop more potent inhibitors of the pancreatic RNase superfamily using this nucleotide as the starting point.

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.

The structural and functional studies of His119 and His12 in RNase A via chemical modification

The histidyl residues of bovine pancreatic ribonuclease A (RNase A) play a crucial role in enzymatic activity. Diethylpyrocarbonate (DEPC) is a potent inhibitor of RNase A, and its precise sites of action on the imidazole rings of the four histidyl residues of RNase A are not clearly defined. We have used a multidisciplinary approach including enzyme assay, calculation of accessible surface area (ASA), isoelectric pH gradient technique, fluorescence investigations, circular dichroism spectroscopy, differential scanning calorimetry, and 1H NMR analysis to study the sites of DEPC interaction with the imidazole rings of the four histidyl residues. Our results demonstrate that among the histidyl residues of RNase A, His48 is not accessible to react with DEPC. However, the sequential carbethoxylation of the imidazole rings of His119, His105, and His12 occurs on the nitrogen atoms of Ndelta, Nepsilon, and Nepsilon, respectively. Carbethoxylation of His119 was followed by conversion of the A conformation to the B conformation in the active site. However, the carbethoxylation of His12 was accompanied by a second spatial rotation of the corresponding imidazole ring in the active site to adopt a new conformation. These conformation changes are accompanied by subsequent decrements in the thermal stability of the protein. Therefore, these findings reinforce the important structural roles of the spatial positions for His119 and His12 in the active site of RNase A.

High-resolution crystal structures of ribonuclease A complexed with adenylic and uridylic nucleotide inhibitors. Implications for structure-based design of ribonucleolytic inhibitors

The crystal structures of bovine pancreatic ribonuclease A (RNase A) in complex with 3',5'-ADP, 2',5'-ADP, 5'-ADP, U-2'-p and U-3'-p have been determined at high resolution. The structures reveal that each inhibitor binds differently in the RNase A active site by anchoring a phosphate group in subsite P1. The most potent inhibitor of all five, 5'-ADP (Ki = 1.2 microM), adopts a syn conformation (in contrast to 3',5'-ADP and 2',5'-ADP, which adopt an anti), and it is the beta- rather than the alpha-phosphate group that binds to P1. 3',5'-ADP binds with the 5'-phosphate group in P1 and the adenosine in the B2 pocket. Two different binding modes are observed in the two RNase A molecules of the asymmetric unit for 2',5'-ADP. This inhibitor binds with either the 3' or the 5' phosphate groups in subsite P1, and in each case, the adenosine binds in two different positions within the B2 subsite. The two uridilyl inhibitors bind similarly with the uridine moiety in the B1 subsite but the placement of a different phosphate group in P1 (2' versus 3') has significant implications on their potency against RNase A. Comparative structural analysis of the RNase A, eosinophil-derived neurotoxin (EDN), eosinophil cationic protein (ECP), and human angiogenin (Ang) complexes with these and other phosphonucleotide inhibitors provides a wealth of information for structure-based design of inhibitors specific for each RNase. These inhibitors could be developed to therapeutic agents that could control the biological activities of EDN, ECP, and ANG, which play key roles in human pathologies.

Predictive crystallization of ribonuclease A via rapid screening of osmotic second virial coefficients

Important progress has been made in recent years toward developing a molecular-level understanding of protein phase behavior in terms of the osmotic second virial coefficient, a thermodynamic parameter that characterizes pairwise protein interactions. Yet there has been little practical application of this knowledge to the field of protein crystallization, largely because of the difficult and time-consuming nature of traditional techniques for characterizing protein interactions. Self-interaction chromatography has recently been proposed as a highly efficient method for measuring the osmotic second virial coefficient. The utility of the technique is examined in this work by characterizing virial coefficients for ribonuclease A under 59 solution conditions using several crystallization additives, including PEG, sodium chloride, ammonium sulfate, and propanol. The virial coefficient measurements show some counterintuitive trends and shed light on the previous difficulties in crystallizing ribonuclease A. Crystallization experiments at the corresponding solution conditions were conducted by using ultracentrifugal crystallization. Using this methodology, ribonuclease A crystals were obtained under conditions for which the virial coefficients fell within the "crystallization slot." Crystallographic characterization showed that the crystals diffract to high resolution. Metastable crystals were also obtained for conditions outside, but near, the "crystallization slot," and they could also be frozen and used to collect structural information.

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.

Crystallographic and functional studies of a modified form of eosinophil-derived neurotoxin (EDN) with novel biological activities

The crystal structure of a post-translationally modified form of eosinophil-derived neurotoxin (EDN) with four extra residues on its N terminus ((-4)EDN) has been solved and refined at atomic resolution (1 A). Two of the extra residues can be placed unambiguously, while the density corresponding to two others is poor. The modified N terminus appears to influence the position of the catalytically important His129, possibly explaining the diminished catalytic activity of this variant. However, (-4)EDN has been shown to be cytotoxic to a Kaposi's sarcoma tumor cell line and other endothelial cell lines. Analysis of the structure and function suggests that the reason for cytotoxicity is most likely due to cellular recognition by the N-terminal extension, since the intrinsic activity of the enzyme is not sufficient for cytotoxicity and the N-terminal extension does not affect the conformation of EDN.

Structures of the two 3D domain-swapped RNase A trimers

When concentrated in mildly acidic solutions, bovine pancreatic ribonuclease (RNase A) forms long-lived oligomers including two types of dimer, two types of trimer, and higher oligomers. In previous crystallographic work, we found that the major dimeric component forms by a swapping of the C-terminal beta-strands between the monomers, and that the minor dimeric component forms by swapping the N-terminal alpha-helices of the monomers. On the basis of these structures, we proposed that a linear RNase A trimer can form from a central molecule that simultaneously swaps its N-terminal helix with a second RNase A molecule and its C-terminal strand with a third molecule. Studies by dissociation are consistent with this model for the major trimeric component: the major trimer dissociates into both the major and the minor dimers, as well as monomers. In contrast, the minor trimer component dissociates into the monomer and the major dimer. This suggests that the minor trimer is cyclic, formed from three monomers that swap their C-terminal beta-strands into identical molecules. These conclusions are supported by cross-linking of lysyl residues, showing that the major trimer swaps its N-terminal helix, and the minor trimer does not. We verified by X-ray crystallography the proposed cyclic structure for the minor trimer, with swapping of the C-terminal beta-strands. This study thus expands the variety of domain-swapped oligomers by revealing the first example of a protein that can form both a linear and a cyclic domain-swapped oligomer. These structures permit interpretation of the enzymatic activities of the RNase A oligomers on double-stranded RNA.

Conformational strictness required for maximum activity and stability of bovine pancreatic ribonuclease A as revealed by crystallographic study of three Phe120 mutants at 1.4 A resolution

The replacement of Phe120 with other hydrophobic residues causes a decrease in the activity and thermal stability in ribonuclease A (RNase A). To explain this, the crystal structures of wild-type RNase A and three mutants--F120A, F120G, and F120W--were analyzed up to a 1.4 A resolution. Although the overall backbone structures of all mutant samples were nearly the same as that of wild-type RNase A, except for the C-terminal region of F120G with a high B-factor, two local conformational changes were observed at His119 in the mutants. First, His119 of the wild-type and F120W RNase A adopted an A position, whereas those of F120A and F120G adopted a B position, but the static crystallographic position did not reflect either the efficiency of transphosphorylation or the hydrolysis reaction. Second, His119 imidazole rings of all mutant enzymes were deviated from that of wild-type RNase A, and those of F120W and F120G appeared to be "inside out" compared with that of wild-type RNase A. Only approximately 1 A change in the distance between N(epsilon2) of His12 and N(delta1) of His119 causes a drastic decrease in k(cat), indicating that the active site requires the strict positioning of the catalytic residues. A good correlation between the change in total accessible surface area of the pockets on the surface of the mutant enzymes and enthalpy change in their thermal denaturation also indicates that the effects caused by the replacements are not localized but extend to remote regions of the protein molecule.

Mapping the ribonucleolytic active site of eosinophil-derived neurotoxin (EDN). High resolution crystal structures of EDN complexes with adenylic nucleotide inhibitors

Eosinophil-derived neurotoxin (EDN), a basic ribonuclease found in the large specific granules of eosinophils, belongs to the pancreatic RNase A family. Although its physiological function is still unclear, it has been shown that EDN is a neurotoxin capable of inducing the Gordon phenomenon in rabbits. EDN is also a potent helminthotoxin and can mediate antiviral activity of eosinophils against isolated virions of the respiratory syncytial virus. EDN is a catalytically efficient RNase sharing similar substrate specificity with pancreatic RNase A with its ribonucleolytic activity being absolutely essential for its neurotoxic, helminthotoxic, and antiviral activities. The crystal structure of recombinant human EDN in the unliganded form has been determined previously (Mosimann, S. C., Newton, D. L., Youle, R. J., and James, M. N. G. (1996) J. Mol. Biol. 260, 540-552). We have now determined high resolution (1.8 A) crystal structures for EDN in complex with adenosine-3',5'-diphosphate (3',5'-ADP), adenosine-2',5'-di-phosphate (2',5'-ADP), adenosine-5'-diphosphate (5'-ADP) as well as for a native structure in the presence of sulfate refined at 1.6 A. The inhibition constant of these mononucleotides for EDN has been determined. The structures present the first detailed picture of differences between EDN and RNase A in substrate recognition at the ribonucleolytic active site. They also provide a starting point for the design of tight-binding inhibitors, which may be used to restrain the RNase activity of EDN.

The ultrahigh resolution crystal structure of ribonuclease A containing an isoaspartyl residue: hydration and sterochemical analysis

Crystals of the deamidated form of bovine pancreatic ribonuclease which contains an isoaspartyl residue in position 67 diffract to 0. 87 A at 100 K. We have refined the crystallographic model using anisotropic displacement parameters for all atoms to a conventional crystallographic residual R=0.101 for all observed reflections in the resolution range 61.0-0.87 A. The ratio observations/parameters is 7.2 for the final model. This structure represents one of the highest resolution protein structures to date and interestingly, it is the only example containing more than one molecule in the asymmetric unit with a resolution better than 1.0 A. The non-crystallographic symmetry has been used as a validation check of the geometrical parameters and it has allowed an estimate for an upper limit of errors associated with this high resolution model. In the present structure it was possible to obtain a more accurate picture of the active site whose electron density was not clearly interpretable in the previous 1.9 A resolution structure. In particular, the P1 site is alternatively occupied either by a sulphate anion or by a water molecule network. Most of hydrogen atoms were visible in the electron density maps, including those involved in C(alpha)-H(alpha).O interactions. Analysis of protein-solvent interactions has revealed the occurrence of an extensive cluster of water molecules, predominantly arranged in pentagonal fused rings and surrounding hydrophobic moiety of side-chains. Finally, in spite of the limited sample of residues, we have detected a clear dependence of backbone N-C(alpha)-C angle on residue conformation. This correlation can be fruitfully used as a valuable tool in protein structure validation.

New approaches to rational drug design

Rational drug design has emerged as a powerful technique. In this review, three new developments in rational drug design are explored. These developments include new methods to find binding sites for small molecules on the surface of a protein, the suggestion that the protein environment may change the shape of a protein sufficiently to alter drug design, and the use of data emerging from structural genomics in drug design. Although these are three new and distinct areas, the insights derived from these studies suggest a reason for the observation that similar drugs do not always bind to a protein in the same manner.

What is the true structure of liganded haemoglobin?

Does the crystal structure of a protein accurately represent its structure in solution? Or does the crystallization process perturb the structure significantly? Although aware of the problem, most crystallographers would argue that the highly solvated and weakly held lattice in protein crystals is, in general, unlikely to shift ordered parts of the molecule. In the case of conformationally flexible proteins, however, there is the possibility that one form might be favoured over another. Several lines of evidence suggest that this might be the case for the crystal structure of liganded Hb, although conflicting data exist.

Monoacylation of ribonuclease A enables its transport across an in vitro model of the blood-brain barrier

A major challenge in correcting disorders affecting the central nervous system is to induce blood-brain barrier (BBB) crossing of exogenous biological compounds such as proteins or specific nucleic acid sequences. Fatty acids, due to their high membrane affinity and low toxicity, are good potential candidates to promote this barrier crossing when covalently bound to proteins. In this paper, we report that regiospecific monoacylation of ribonuclease A (RNase A) enables its transport across an in vitro model of the BBB. Myristoylated, palmitoylated and stearoylated RNases A were prepared using reversed micelles as microreactors. All the purified acylated RNases A kept their original enzymatic activity. A single fatty acid moiety was linked to RNase A through the alpha-amino group of its N-terminal lysine as shown by powerful analytical techniques. The ability of monoacylated RNases A to cross an in vitro model of the BBB is strictly dependent on the acyl chain length, which must be at least 16 carbon atoms long.

Structure and stability of the P93G variant of ribonuclease A

The peptide bonds preceding Pro 93 and Pro 114 of bovine pancreatic ribonuclease A (RNase A) are in the cis conformation. The trans-to-cis isomerization of these bonds had been indicted as the slow step during protein folding. Here, site-directed mutagenesis was used to replace Pro 93 or Pro 114 with a glycine residue, and the crystalline structure of the P93G variant was determined by X-ray diffraction analysis to a resolution of 1.7 A. This structure is essentially identical to that of the wild-type protein, except for the 91-94 beta-turn containing the substitution. In the wild-type protein, the beta-turn is of type VIa. In the P93G variant, this turn is of type II with the peptide bond preceding Gly 93 being trans. The thermal stabilities of the P93G and P114G variants were assessed by differential scanning calorimetry and thermal denaturation experiments monitored by ultraviolet spectroscopy. The value of delta deltaGm which reports on the stability lost in the variants, is 1.5-fold greater for the P114G variant than for the P93G variant. The greater stability of the P93G variant is likely due to the relatively facile accommodation of residues 91-94 in a type II turn, which has a preference for a glycine residue in its i + 2 position.

His...Asp catalytic dyad of ribonuclease A: structure and function of the wild-type, D121N, and D121A enzymes

The side chains of histidine and aspartate residues form a hydrogen bond in the active sites of many enzymes. In serine proteases, the His...Asp hydrogen bond of the catalytic triad is known to contribute greatly to catalysis, perhaps via the formation of a low-barrier hydrogen bond. In bovine pancreatic ribonuclease A (RNase A), the His...Asp dyad is composed of His119 and Asp121. Previously, site-directed mutagenesis was used to show that His119 has a fundamental role, to act as an acid during catalysis of RNA cleavage [Thompson, J. E., and Raines, R. T. (1994) J. Am. Chem. Soc. 116, 5467-5468]. Here, Asp121 was replaced with an asparagine or alanine residue. The crystalline structures of the two variants were determined by X-ray diffraction analysis to a resolution of 1.6 A with an R-factor of 0.18. Replacing Asp121 with an asparagine or alanine residue does not perturb the overall conformation of the enzyme. In the structure of D121N RNase A, Ndelta rather than Odelta of Asn121 faces His119. This alignment in the crystalline state is unlikely to exist in solution because catalysis by the D121N variant is not compromised severely. The steady-state kinetic parameters for catalysis by the wild-type and variant enzymes were determined for the cleavage of uridylyl(3'-->5')adenosine and poly(cytidylic acid), and for the hydrolysis of uridine 2',3'-cyclic phosphate. Replacing Asp121 decreases the values of kcat/Km and kcat for cleavage by 10-fold (D121N) and 10(2)-fold (D121A). Replacing Asp121 also decreases the values of kcat/Km and kcat for hydrolysis by 10(0. 5)-fold (D121N) and 10-fold (D121A) but has no other effect on the pH-rate profiles for hydrolysis. There is no evidence for the formation of a low-barrier hydrogen bond between His119 and either an aspartate or an asparagine residue at position 121. Apparently, the major role of Asp121 is to orient the proper tautomer of His119 for catalysis. Thus, the mere presence of a His...Asp dyad in an enzymic active site is not a mandate for its being crucial in effecting catalysis.