Free energy determinants of binding the rRNA substrate and small ligands to ricin A-chain

Toxicity and Efficacy Evaluation of Soluble Recombinant Ricin Vaccine

BACKGROUND

Ricin, a toxin extracted from the seeds of Ricinus communis, is classified as a ribosome-inactivating protein. The A-subunit of ricin shows RNA N-glycosidase activity that cleaves ribosomal RNA (rRNA) and exhibits toxicity by inhibiting protein synthesis and inducing vascular leak syndrome.

METHODS

In this study, we created a truncated version of the previously developed R51 ricin vaccine (RTA 1-194 D75C Y80C) through in silico analysis.

RESULTS

The resulting R51-3 vaccine showed a more-than-six-fold increase in soluble protein expression when compared to R51, with over 85% solubility. In a pilot toxicity test, no toxicity was observed in hematological and biochemical parameters in BALB/c mice and New Zealand white rabbits following five repeated administrations of R51-3. Furthermore, R51-3 successfully protected mice and rabbits from a 20 × LD50 ricin challenge after three intramuscular injections spaced 2 weeks apart. Similarly, monkeys that received three injections of R51-3 survived a 60 µg/kg ricin challenge.

CONCLUSIONS

These findings support R51-3 as a promising candidate antigen for ricin vaccine development.

Thermochemical and Quantum Descriptor Calculations for Gaining Insight into Ricin Toxin A (RTA) Inhibitors

In this work, we performed a study to assess the interactions between the ricin toxin A (RTA) subunit of ricin and some of its inhibitors using modern semiempirical quantum chemistry and ONIOM quantum mechanics/molecular mechanics (QM/MM) methods. Two approaches were followed (calculation of binding enthalpies, ΔH _bind, and reactivity quantum chemical descriptors) and compared with the respective half-maximal inhibitory concentration (IC₅₀) experimental data, to gain insight into RTA inhibitors and verify which quantum chemical method would better describe RTA-ligand interactions. The geometries for all RTA-ligand complexes were obtained after running classical molecular dynamics simulations in aqueous media. We found that single-point energy calculations of ΔH _bind with the PM6-DH+, PM6-D3H4, and PM7 semiempirical methods and ONIOM QM/MM presented a good correlation with the IC₅₀ data. We also observed, however, that the correlation decreased significantly when we calculated ΔH _bind after full-atom geometry optimization with all semiempirical methods. Based on the results from reactivity descriptors calculations for the cases studied, we noted that both types of interactions, molecular overlap and electrostatic interactions, play significant roles in the overall affinity of these ligands for the RTA binding pocket.

Unraveling the Roots of Selectivity of Peptide Affinity Reagents for Structurally Similar Ribosomal Inactivating Protein Derivatives

Peptide capture agents have become increasingly useful tools for a variety of sensing applications due to their ease of discovery, stability, and robustness. Despite the ability to rapidly discover candidates through biopanning bacterial display libraries and easily mature them to Protein Catalyzed Capture (PCC) agents with even higher affinity and selectivity, an ongoing challenge and critical selection criteria is that the peptide candidates and final reagent be selective enough to replace antibodies, the gold-standard across immunoassay platforms. Here, we have discovered peptide affinity reagents against abrax, a derivative of abrin with reduced toxicity. Using on-cell Fluorescence Activated Cell Sorting (FACS) assays, we show that the peptides are highly selective for abrax over RiVax, a similar derivative of ricin originally designed as a vaccine, with significant structural homology to abrax. We rank the newly discovered peptides for strongest affinity and analyze three observed consensus sequences with varying affinity and specificity. The strongest (Tier 1) consensus was FWDTWF, which is highly aromatic and hydrophobic. To better understand the observed selectivity, we use the XPairIt peptide-protein docking protocol to analyze binding location predictions of the individual Tier 1 peptides and consensus on abrax and RiVax. The binding location profiles on the two proteins are quite distinct, which we determine is due to differences in pocket size, pocket environment (including hydrophobicity and electronegativity), and steric hindrance. This study provides a model system to show that peptide capture candidates can be quite selective for a structurally similar protein system, even without further maturation, and offers an in silico method of analysis for understanding binding and down-selecting candidates.

Do the A subunits contribute to the differences in the toxicity of Shiga toxin 1 and Shiga toxin 2?

Shiga toxin producing Escherichia coli O157:H7 (STEC) is one of the leading causes of food-poisoning around the world. Some STEC strains produce Shiga toxin 1 (Stx1) and/or Shiga toxin 2 (Stx2) or variants of either toxin, which are critical for the development of hemorrhagic colitis (HC) or hemolytic uremic syndrome (HUS). Currently, there are no therapeutic treatments for HC or HUS. E. coli O157:H7 strains carrying Stx2 are more virulent and are more frequently associated with HUS, which is the most common cause of renal failure in children in the US. The basis for the increased potency of Stx2 is not fully understood. Shiga toxins belong to the AB5 family of protein toxins with an A subunit, which depurinates a universally conserved adenine residue in the α-sarcin/ricin loop (SRL) of the 28S rRNA and five copies of the B subunit responsible for binding to cellular receptors. Recent studies showed differences in the structure, receptor binding, dependence on ribosomal proteins and pathogenicity of Stx1 and Stx2 and supported a role for the B subunit in differential toxicity. However, the current data do not rule out a potential role for the A1 subunits in the differential toxicity of Stx1 and Stx2. This review highlights the recent progress in understanding the differences in the A1 subunits of Stx1 and Stx2 and their role in defining toxicity.

Arginine residues on the opposite side of the active site stimulate the catalysis of ribosome depurination by ricin A chain by interacting with the P-protein stalk

Ricin inhibits protein synthesis by depurinating the α-sarcin/ricin loop (SRL). Ricin holotoxin does not inhibit translation unless the disulfide bond between the A (RTA) and B (RTB) subunits is reduced. Ricin holotoxin did not bind ribosomes or depurinate them but could depurinate free RNA. When RTA is separated from RTB, arginine residues located at the interface are exposed to the solvent. Because this positively charged region, but not the active site, is blocked by RTB, we mutated arginine residues at or near the interface of RTB to determine if they are critical for ribosome binding. These variants were structurally similar to wild type RTA but could not bind ribosomes. Their K(m) values and catalytic rates (k(cat)) for an SRL mimic RNA were similar to those of wild type, indicating that their activity was not altered. However, they showed an up to 5-fold increase in K(m) and up to 38-fold decrease in kcat toward ribosomes. These results suggest that the stalk binding stimulates the catalysis of ribosome depurination by RTA. The mutated arginines have side chains behind the active site cleft, indicating that the ribosome binding surface of RTA is on the opposite side of the surface that interacts with the SRL. We propose that stalk binding stimulates the catalysis of ribosome depurination by orienting the active site of RTA toward the SRL and thereby allows docking of the target adenine into the active site. This model may apply to the translation factors that interact with the stalk.

Disruption of the putative vascular leak peptide sequence in the stabilized ricin vaccine candidate RTA1-33/44-198

Vitetta and colleagues identified and characterized a putative vascular leak peptide (VLP) consensus sequence in recombinant ricin toxin A-chain (RTA) that contributed to dose-limiting human toxicity when RTA was administered intravenously in large quantities during chemotherapy. We disrupted this potentially toxic site within the more stable RTA1-33/44-198 vaccine immunogen and determined the impact of these mutations on protein stability, structure and protective immunogenicity using an experimental intranasal ricin challenge model in BALB/c mice to determine if the mutations were compatible. Single amino acid substitutions at the positions corresponding with RTA D75 (to A, or N) and V76 (to I, or M) had minor effects on the apparent protein melting temperature of RTA1-33/44-198 but all four variants retained greater apparent stability than the parent RTA. Moreover, each VLP(−) variant tested provided protection comparable with that of RTA1-33/44-198 against supralethal intranasal ricin challenge as judged by animal survival and several biomarkers. To understand better how VLP substitutions and mutations near the VLP site impact epitope structure, we introduced a previously described thermal stabilizing disulfide bond (R48C/T77C) along with the D75N or V76I substitutions in RTA1-33/44-198. The D75N mutation was compatible with the adjacent stabilizing R48C/T77C disulfide bond and the T_m was unaffected, whereas the V76I mutation was less compatible with the adjacent disulfide bond involving C77. A crystal structure of the RTA1-33/44-198 R48C/T77C/D75N variant showed that the structural integrity of the immunogen was largely conserved and that a stable immunogen could be produced from E. coli. We conclude that it is feasible to disrupt the VLP site in RTA1-33/44-198 with little or no impact on apparent protein stability or protective efficacy in mice and such variants can be stabilized further by introduction of a disulfide bond.

Structure/function studies on two type 1 ribosome inactivating proteins: Bouganin and lychnin

The three-dimensional structures of two type 1 RIPs, bouganin and lychnin, has been solved. Their adenine polynucleotide glycosylase activity was also determined together with other known RIPs: dianthin 30, PAP-R, momordin I, ricin A chain and saporin-S6. Saporin-S6 releases the highest number of adenine molecules from rat ribosomes, and poly(A), while its efficiency is similar to dianthin 30, bouganin and PAP-R on herring sperm DNA. Measures of the protein synthesis inhibitory activity confirmed that saporin-S6 is the most active. The overall structure of bouganin and lychnin is similar to the other considered RIPs and the typical RIP fold is conserved. The superimpositioning of their C(alpha) atoms highlights some differences in the N-terminal and C-terminal domains. A detailed structural analysis indicates that the efficiency of saporin-S6 on various polynucleotides can be ascribed to a negative electrostatic surface potential at the active site and several exposed positively charged residues in the region around that site. These two conditions, not present at the same time in other examined RIPs, could guarantee an efficient interaction with the substrate and an efficient catalysis.

Defense Against Biological Weapons (Biodefense)

Biological warfare (germ warfare) is defined as the use of any disease-causing organism or toxin(s) found in nature as weapons of war with the intent to destroy an adversary. Though rare, the use of biological weapons has occurred throughout the centuries.

Calculation of absolute ligand binding free energy to a ribosome-targeting protein as a function of solvent model

A comparative analysis is provided of the effect of different solvent models on the calculation of a potential of mean force (PMF) for determining the absolute binding affinity of the small molecule inhibitor pteroic acid bound to ricin toxin A-chain (RTA). Solvent models include the distance-dependent dielectric constant, several different generalized Born (GB) approximations, and a hybrid explicit/GB-based implicit solvent model. We found that the simpler approximation of dielectric screening and a GB model, with Born radii fitted to a switching-window dielectric-boundary surface Poisson solvent model, severely overpredicted the binding affinity as compared to the experimental value, estimated to range from -4.4 to -6.0 kcal/mol. In contrast, GB models that are parametrized to fit the Lee-Richards molecular surface performed much better, predicting binding free energy within 1-3 kcal/mol of experimental estimates. However, the predicted free-energy profiles of these GB models displayed alternative binding modes not observed in the crystal structure. Finally, the most rigorous and computationally costly approach in this work, which used a hybrid explicit/implicit solvent model, correctly determined a binding funnel in the PMF near the crystallographic bound state and predicted an absolute binding affinity that was 2 kcal/mol more favorable than the estimated experimental binding affinity.

Structural Re-Alignment in an Immunogenic Surface Region of Ricin a Chain

We compared structure alignments generated by several protein structure comparison programs to determine whether existing methods would satisfactorily align residues at a highly conserved position within an immunogenic loop in ribosome inactivating proteins (RIPs). Using default settings, structure alignments generated by several programs (CE, DaliLite, FATCAT, LGA, MAMMOTH, MATRAS, SHEBA, SSM) failed to align the respective conserved residues, although LGA reported correct residue-residue (R-R) correspondences when the beta-carbon (Cb) position was used as the point of reference in the alignment calculations. Further tests using variable points of reference indicated that points distal from the beta carbon along a vector connecting the alpha and beta carbons yielded rigid structural alignments in which residues known to be highly conserved in RIPs were reported as corresponding residues in structural comparisons between ricin A chain, abrin-A, and other RIPs. Results suggest that approaches to structure alignment employing alternate point representations corresponding to side chain position may yield structure alignments that are more consistent with observed conservation of functional surface residues than do standard alignment programs, which apply uniform criteria for alignment (i.e. alpha carbon (Ca) as point of reference) along the entirety of the peptide chain. We present the results of tests that suggest the utility of allowing user-specified points of reference in generating alternate structural alignments, and we present a web server for automatically generating such alignments: http://as2ts.llnl.gov/AS2TS/LGA/lga_pdblist_plots.html.

GBr6NL: a generalized Born method for accurately reproducing solvation energy of the nonlinear Poisson-Boltzmann equation

The nonlinear Poisson-Boltzmann (NLPB) equation can provide accurate modeling of electrostatic effects for nucleic acids and highly charged proteins. Generalized Born methods have been developed to mimic the linearized Poisson-Boltzmann (LPB) equation at substantially reduced cost. The computer time for solving the NLPB equation is approximately fivefold longer than for the LPB equation, thus presenting an even greater obstacle. Here we present the first generalized Born method, GBr(6)NL, for mimicking the NLPB equation. GBr(6)NL is adapted from GBr(6), a generalized Born method recently developed to reproduce the solvation energy of the LPB equation [Tjong and Zhou, J. Phys. Chem. B 111, 3055 (2007)]. Salt effects predicted by GBr(6)NL on 55 proteins overall deviate from NLPB counterparts by 0.5 kcal/mol from ionic strengths from 10 to 1000 mM, which is approximately 10% of the average magnitudes of the salt effects. GBr(6)NL predictions for the salts effects on the electrostatic interaction energies of two protein:RNA complexes are very promising.

Do electrostatic interactions destabilize protein-nucleic acid binding?

The negatively charged phosphates of nucleic acids are often paired with positively charged residues upon binding proteins. It was thus counter-intuitive when previous Poisson-Boltzmann (PB) calculations gave positive energies from electrostatic interactions, meaning that they destabilize protein-nucleic acid binding. Our own PB calculations on protein-protein binding have shown that the sign and the magnitude of the electrostatic component are sensitive to the specification of the dielectric boundary in PB calculations. A popular choice for the boundary between the solute low dielectric and the solvent high dielectric is the molecular surface; an alternative is the van der Waals (vdW) surface. In line with results for protein-protein binding, in this article, we found that PB calculations with the molecular surface gave positive electrostatic interaction energies for two protein-RNA complexes, but the signs are reversed when the vdW surface was used. Therefore, whether destabilizing or stabilizing effects are predicted depends on the choice of the dielectric boundary. The two calculation protocols, however, yielded similar salt effects on the binding affinity. Effects of charge mutations differentiated the two calculation protocols; PB calculations with the vdW surface had smaller deviations overall from experimental data.

Circular DNA and DNA/RNA hybrid molecules as scaffolds for ricin inhibitor design

Ricin Toxin A-chain (RTA) catalyzes the hydrolytic depurination of A4324, the first adenosine of the GAGA tetra-loop portion of 28S eukaryotic ribosomal RNA. Truncated stem-loop versions of the 28S rRNA are RTA substrates. Here, we investigate circular DNA and DNA/RNA hybrid GAGA sequence oligonucleotides as minimal substrates and inhibitor scaffolds for RTA catalysis. Closing the 5'- and 3'-ends of a d(GAGA) tetraloop creates a substrate with 92-fold more activity with RTA (kcat/Km) than that for the d(GAGA) linear form. Circular substrates have catalytic rates (kcat) comparable to and exceeding those of RNA and DNA stem-loop substrates, respectively. RTA inhibition into the nanomolar range has been achieved by introducing an N-benzyl-hydroxypyrrolidine (N-Bn) transition state analogue at the RTA depurination site in a circular GAGA motif. The RNA/DNA hybrid oligonucleotide cyclic GdAGA provides a new scaffold for RTA inhibitor design, and cyclic G(N-Bn)GA is the smallest tight-binding RTA inhibitor (Ki = 70 nM). The design of such molecules that lack the base-paired stem-loop architecture opens new chemical synthetic approaches to RTA inhibition.

Improved formulation of a recombinant ricin A-chain vaccine increases its stability and effective antigenicity

Ricin is a potent toxin associated with bioterrorism for which no vaccine or specific countermeasures are currently available. A stable, non-toxic and immunogenic recombinant ricin A-chain vaccine (RTA 1-33/44-198) has been developed by protein engineering. We identified optimal formulation conditions for this vaccine under which it remained stable and potent in storage for up to 18 months, and resisted multiple rounds of freeze-thawing without stabilizing co-solvents. Reformulation from phosphate buffer to succinate buffer increased adherence of the protein to aluminum hydroxide adjuvant from 15 to 91%, with a concomitant increase of nearly threefold in effective antigenicity in a mouse model. Using Fourier-transform infrared spectroscopy, we examined the secondary structure of the protein while it was adhered to aluminum hydroxide. Adjuvant adsorption produced only a small apparent change in secondary structure, while significantly stabilizing the protein to thermal denaturation. The vaccine therefore may be safely stored in the presence of adjuvant. Our results suggest that optimization of adherence of a protein antigen to aluminum adjuvant can be a useful route to increasing both stability and effectiveness, and support a role for a "depot effect" of adjuvant.

Improved stability of a protein vaccine through elimination of a partially unfolded state

Ricin is a potent toxin presenting a threat as a biological weapon. The holotoxin consists of two disulfide-linked polypeptides: an enzymatically active A chain (RTA) and a galactose/N-acetylgalactosamine-binding B chain. Efforts to develop an inactivated version of the A chain as a vaccine have been hampered by limitations of stability and solubility. Previously, recombinant truncated versions of the 267-amino-acid A chain consisting of residues 1-33/44-198 or 1-198 were designed by protein engineering to overcome these limits and were shown to be effective and nontoxic as vaccines in mice. Herein we used CD, dynamic light scattering, fluorescence, and Fourier-transform infrared spectroscopy to examine the biophysical properties of these proteins. Although others have found that recombinant RTA (rRTA) adopts a partially unfolded, molten globule-like state at 45 degrees C, rRTA 1-33/44-198 and 1-198 are significantly more thermostable, remaining completely folded at temperatures up to 53 degrees C and 51 degrees C, respectively. Deleting both an exposed loop region (amino acids 34-43) and the C-terminal domain (199-267) contributed to increased thermostability. We found that chemically induced denaturation of rRTA, but not the truncated variants, proceeds through at least a three-state mechanism. The intermediate state in rRTA unfolding has a hydrophobic core accessible to ANS and an unfolded C-terminal domain. Removing the C-terminal domain changed the mechanism of rRTA unfolding, eliminating a tendency to adopt a partially unfolded state. Our results support the conclusion that these derivatives are superior candidates for development as vaccines against ricin and suggest an approach of reduction to minimum essential domains for design of more thermostable recombinant antigens.

The effect of mutations surrounding and within the active site on the catalytic activity of ricin A chain

Models for the binding of the sarcin-ricin loop (SRL) of 28S ribosomal RNA to ricin A chain (RTA) suggest that several surface exposed arginine residues surrounding the active site cleft make important interactions with the RNA substrate. The data presented in this study suggest differing roles for these arginyl residues. Substitution of Arg48 or Arg213 with Ala lowered the activity of RTA 10-fold. Furthermore, substitution of Arg213 with Asp lowered the activity of RTA 100-fold. The crystal structure of this RTA variant showed it to have an unaltered tertiary structure, suggesting that the positively charged state of Arg213 is crucial for activity. Substitution of Arg258 with Ala had no effect on activity, although substitution with Asp lowered activity 10-fold. Substitution of Arg134 prevented expression of folded protein, suggesting a structural role for this residue. Several models have been proposed for the binding of the SRL to the active site of RTA in which the principal difference lies in the conformation of the second 'G' in the target GAGA motif in the 28S rRNA substrate. In one model, the sidechain of Asn122 is proposed to make interactions with this G, whereas another model proposes interactions with Asp75 and Asn78. Site-directed mutagenesis of these residues of RTA favours the first of these models, as substitution of Asn78 with Ser yielded an RTA variant whose activity was essentially wild-type, whereas substitution of Asn122 reduced activity 37.5-fold. Substitution of Asp75 failed to yield significant folded protein, suggesting a structural role for this residue.

Therapy and prophylaxis of inhaled biological toxins

This review highlights the current lack of therapeutic and prophylactic treatments for use against inhaled biological toxins, especially those considered as potential biological warfare (BW) or terrorist threats. Although vaccine development remains a priority, the use of rapidly deployable adjunctive therapeutic or prophylactic drugs could be life-saving in severe cases of intoxication or where vaccination has not been possible or immunity not established. The current lack of such drugs is due to many factors. Thus, methods involving molecular modelling are limited by the extent to which the cellular receptor sites and mode of action and structure of a toxin need to be known. There is also our general lack of knowledge of what effect individual toxins will have when inhaled into the lungs - whether and to what extent the action will be cell specific and cytotoxic or rather an acute inflammatory response requiring the use of immunomodulators. Possible sources of specific high-affinity toxin antagonists being investigated include monoclonal antibodies, selected oligonucleotides (aptamers) and derivatized dendritic polymers (dendrimers). The initial selection of suitable agents of these kinds can be made using cytotoxicity assays involving cultured normal human lung cells and a range of suitable indicators. The possibility that a mixture of selected antibody, aptamer or dendrimer-based materials for one or more toxins could be delivered simultaneously as injections or as inhaled aerosol sprays should be investigated.

Can the calculation of ligand binding free energies be improved with continuum solvent electrostatics and an ideal-gas entropy correction?

The prediction of a ligand binding constant requires generating three-dimensional structures of the complex concerned and reliably scoring these structures. Here, the scoring problem is investigated by examining benzamidine-like inhibitors of trypsin, a system for which errors in the structures are small. Precise and consistent binding free energies for the inhibitors are determined experimentally for this test system. To examine possible improvement of scoring methods, we test the suitability of continuum electrostatics to account for solvation effects and use an ideal-gas entropy correction to account for the changes in the degrees of freedom of the ligand. The small observed root-mean-square deviation of 0.55 kcal/mol of the calculated relative to the experimental values indicates that the essentials of the binding process have been captured. Even though all six ligands make the same salt bridge and H-bonds to the protein, the electrostatic contribution varies among the ligands by as much as 2 kcal/mol. Moreover, although the ligands are rigid and similar in size, the entropic terms also significantly affect the relative binding affinities (by up to 2.7 kcal/mol). The present approach to solvation and entropy may allow the ranking of the ligands to be considerably improved at a cost that makes the method applicable to the optimization of lead compounds or to the screening of small collections of ligands.

Free-energy contributions to complex formation between botulinum neurotoxin type B and synaptobrevin fragment

Free-energy terms that contribute to complex formation between the catalytic domain of botulinum neurotoxin type B (BoNT/B-L(C)) and a 36-residue synaptobrevin fragment were estimated by using a combination of microscopic simulations and continuum methods. The complex for a non-hydrolyzed substrate was calculated by optimizing an energy function applied to the X-ray co-crystal structure of BoNT/B-L(C) bound with reaction products from a cleaved synaptobrevin peptide, refined to high crystallographic thermal factors. The estimated absolute binding affinity of the simulation structure is in good qualitative agreement with the experimental free energy of Michaelis complex formation, given the approximations of the model calculations. The simulation structure revealed significant complex stabilization from the hydrophobic effect, while the electrostatic cost of releasing water molecules from the interface determined to be highly unfavorable. By partitioning the total electrostatic and hydrophobic terms into residue free-energy contributions, a binding-affinity 'signature' for synaptobrevin was developed from the optimized conformation. The results demonstrate the effect of substrate length on complex formation and identify a peripheral high-affinity binding site near the N-terminal region that might initiate cooperative activation responsible for the large minimal substrate length requirement. The so-called SNARE motif is observed to contribute negligible free energy of binding.

Calculations of free-energy contributions to protein-RNA complex stabilization

The problem of calculating binding affinities of protein-RNA complexes is addressed by analyzing a computational strategy of modeling electrostatic free energies based on a nonlinear Poisson-Boltzmann (NLPB) model and linear response approximation (LRA). The underlying idea is to treat binding as a two-step process. Solutions to the NLPB equation calculate free energies arising from electronic polarizability and the LRA is constructed from molecular dynamics simulations to model reorganization free energies due to conformational transitions. By implementing a consistency condition of requiring the NLPB model to reproduce the solute-solvent free-energy transitions determined by the LRA, a "macromolecule dielectric constant" (epsilon(m)) for treating reorganization is obtained. The applicability of this hybrid approach was evaluated by calculating the absolute free energy of binding and free-energy changes for amino acid substitutions in the complex between the U1A spliceosomal protein and its cognate RNA hairpin. Depending on the residue substitution, epsilon(m) varied from 3 to 18, and reflected dipolar reorientation not included in the polarization modeled by epsilon(m) = 2. Although the changes in binding affinities from substitutions modeled strictly at the implicit level by the NLPB equation with epsilon(m) = 4 reproduced the experimental values with good overall agreement, substitutions problematic to this simple treatment showed significant improvement when solved by the NLPB-LRA approach.

Electrostatic effects on the free-energy balance in folding a ribosome-inactivating protein

Electrostatics of globular proteins provides structural integrity as well as specificity of biological function. This dual role is particularly striking for ricin A-chain (RTA), an N-glycosidase which hydrolyzes a single adenine base from a conserved region of rRNA. The reported X-ray crystallographic structure of the RTA mutant E177A demonstrated a remarkable rescue of charge balance in the active site, achieved by the rotation of a second glutamic acid (Glu-208) into the vacated space. To understand this conformational reorganization, molecular-dynamics simulations were applied to estimate relative free energies that govern the thermodynamic stability of E177A together with mutants E177Q and E177D. The simulations anticipate that while E177A is a non-conservative substitution, the protein is more stable than the other two mutants. However, the structural plasticity of the RTA active site is not obtained penalty-free, rather E177A among the mutants shows the largest unfavorable net change in the electrostatic contribution to folding. Of the E177A folded state, reorganization of Glu-208 lowers the electrostatic cost of the free-energy change, yet interestingly, protein interactions oppose the rotational shift, while solvent effects favor the transition.

Hypoxic vasodilation in porcine coronary artery is preferentially inhibited by organ culture

Hypoxia (95% N2-5% CO2) elicits an endothelium-independent relaxation (45-80%) in freshly dissected porcine coronary arteries. Paired artery rings cultured at 37 degrees C in sterile DMEM (pH approximately 7.4) for 24 h contracted normally to KCl or 1 microM U-46619. However, relaxation in response to hypoxia was sharply attenuated compared with control (fresh arteries or those stored at 4 degrees C for 24 h). Hypoxic vasorelaxation in organ cultured vessels was reduced at both high and low stimulation, indicating that both Ca2+-independent and Ca2+-dependent components are altered. In contrast, relaxation to G-kinase (sodium nitroprusside) or A-kinase (forskolin and isoproterenol) activation was not significantly affected by organ culture. Additionally, there was no difference in relaxation after washout of the stimulus, indicating that the inhibition is specific to acute hypoxia-induced relaxation. Simultaneous force and intracellular calcium concentration ([Ca2+]i) measurements indicate the reduction in [Ca2+]i concomitant with hypoxia at low stimulus levels in these tissue is abolished by culture. Our results indicate that organ culture at 37 degrees C specifically attenuates hypoxic relaxation in vascular smooth muscle by altering dynamics of [Ca2+]i handling and decreasing a Ca2+-independent component of relaxation. Thus organ culture can be a novel tool for investigating the mechanisms of hypoxia-induced vasodilation.

Simulation of the structure and dynamics of nonhelical RNA motifs

Computer simulation methods are increasingly being used to study possible conformations and dynamics of structural motifs in RNA. Recent results of molecular dynamics simulations and continum solvent studies of RNA structures and RNA-ligand complexes show promising agreement with experimental data. Combined with the ongoing progress in the experimental characterization of RNA structure and thermodynamics, these computational approaches can help to better understand the mechanism of RNA structure formation and the binding of ligands.

Modeling implicit reorganization in continuum descriptions of protein-protein interactions

The determination of free energies that govern protein-protein recognition is essential for a detailed molecular understanding of biological specificity. Continuum models of macromolecular interactions, in which the solvent is treated by an implicit representation and the proteins are treated semi-microscopically, are computationally tractable for estimating free energies, yet many questions remain concerning their accuracy. This article reports a continuum analysis of the free-energy changes underlying the binding of 31 interfacial alanine substitutions of two complexes of the antihen egg white lysozyme (HEL) antibody D1.3 bound with HEL or the antibody E5.2. Two implicit schemes for modeling the effects of protein and solvent relaxation were examined, in which the protein environment was treated as either homogeneous with a "protein dielectric constant" of epsilon(p) = 4 or inhomogeneous, with epsilon(p) = 4 for neutral residues and epsilon(p) = 25 for ionized residues. The results showed that the nonuniform dielectric model reproduced the experimental differences better, with an average absolute error of +/-1.1 kcal/mol, compared with +/-1.4 kcal/mol for the uniform model. More importantly, the error for charged residues in the nonuniform model is +/-0.8 kcal/mol and is nearly half of that corresponding to the uniform model. Several substitutions were clearly problematic in determining qualitative trends and probably required explicit structural reorganization at the protein-protein interface.

Role of Arg163 in the N-glycosidase activity of neo-trichosanthin

Three mutant crystals of neo-trichosanthin (n-TCS), R163K, R163H and R163Q, were obtained by the hanging drop vapor diffusion method. Structure determination indicated that there are no significant differences between the mutants and n-TCS except in the active pocket. All of them were also soaked in sodium citrate buffer (pH 4. 5) containing 20% KCl and 10 mg/ml AMP. Structure determination suggests that in the active pocket of the crystals of R163K and R163H, parallel to the aromatic ring of Tyr70, each mutant possesses an adenine. The relationship between structure and function is discussed. Biochemical analysis reveals that the mutants R163K and R163H have N-glycosidase activity, while R163Q does not. This suggests that R163 is a crucial residue for the enzyme activity of n-TCS, and its role is providing proton.