Dynamics of Cleft Closure of the GluA2 Ligand-binding Domain in the Presence of Full and Partial Agonists Revealed by Hydrogen-Deuterium Exchange

The majority of excitatory neurotransmission in the CNS is mediated by tetrameric AMPA receptors. Channel activation begins with a series of interactions with an agonist that binds to the cleft between the two lobes of the ligand-binding domain of each subunit. Binding leads to a series of conformational transitions, including the closure of the two lobes of the binding domain around the ligand, culminating in ion channel opening. Although a great deal has been learned from crystal structures, determining the molecular details of channel activation, deactivation, and desensitization requires measures of dynamics and stabilities of hydrogen bonds that stabilize cleft closure. The use of hydrogen-deuterium exchange at low pH provides a measure of the variation of stability of specific hydrogen bonds among agonists of different efficacy. Here, we used NMR measurements of hydrogen-deuterium exchange to determine the stability of hydrogen bonds in the GluA2 (AMPA receptor) ligand-binding domain in the presence of several full and partial agonists. The results suggest that the stabilization of hydrogen bonds between the two lobes of the binding domain is weaker for partial than for full agonists, and efficacy is correlated with the stability of these hydrogen bonds. The closure of the lobes around the agonists leads to a destabilization of the hydrogen bonding in another portion of the lobe interface, and removing an electrostatic interaction in Lobe 2 can relieve the strain. These results provide new details of transitions in the binding domain that are associated with channel activation and desensitization.     

The influence of oligodeoxyribonucleotide phosphorothioate pyrimidine strands on triplex stability.

The analogues of the homopyrimidine oligonucleotide dT15 has been synthesized. The analogues, contains phosphorothioate bonds of a mixture of diastereoisomers or one of the two stereoisomer (either Rp or Sp). The analogues were mixed under conditions conductive to the formation of triple stranded assemblies. The mixtures were characterized by their thermal stabilities (Tm values) and CD spectra.     

Positive contribution of hydration structure on the surface of human lysozyme to the conformational stability

Water molecules make a hydration structure with the network of hydrogen bonds, covering on the surface of proteins. To quantitatively estimate the contribution of the hydration structure to protein stability, a series of hydrophilic mutant human lysozymes (Val to Ser, Tyr, Asp, Asn, and Arg) modified at three different positions on the surface, which are located in the alpha-helix (Val-110), the beta-sheet (Val-2), and the loop (Val-74), were constructed. Their thermodynamic parameters of denaturation and crystal structures were examined by calorimetry and by x-ray crystallography at 100 K, respectively. The introduced polar residues made hydrogen bonds with protein atoms and/or water molecules, sometimes changing the hydration structure around the mutation site. Changes in the stability of the mutant proteins can be evaluated by a unique equation that considers the conformational changes resulting from the substitutions. Using this analysis, the relationship between the changes in the stabilities and the hydration structures for mutant human lysozymes substituted on the surface could be quantitatively estimated. The analysis indicated that the hydration structure on protein surface plays an important role in determining the conformational stability of the protein.     

Hydrogen bonding in helical polypeptides from molecular dynamics simulations and amide hydrogen exchange analysis: alamethicin and melittin in methanol.

Molecular dynamics simulations of ion channel peptides alamethicin and melittin, solvated in methanol at 27 degrees C, were run with either regular alpha-helical starting structures (alamethicin, 1 ns; melittin 500 ps either with or without chloride counterions), or with the x-ray crystal coordinates of alamethicin as a starting structure (1 ns). The hydrogen bond patterns and stabilities were characterized by analysis of the dynamics trajectories with specified hydrogen bond angle and distance criteria, and were compared with hydrogen bond patterns and stabilities previously determined from high-resolution NMR structural analysis and amide hydrogen exchange measurements in methanol. The two alamethicin simulations rapidly converged to a persistent hydrogen bond pattern with a high level of 3(10) hydrogen bonding involving the amide NH's of residues 3, 4, 9, 15, and 18. The 3(10) hydrogen bonds stabilizing amide NH's of residues C-terminal to P2 and P14 were previously proposed to explain their high amide exchange stabilities. The absence, or low levels of 3(10) hydrogen bonds at the N-terminus or for A15 NH, respectively, in the melittin simulations, is also consistent with interpretations from amide exchange analysis. Perturbation of helical hydrogen bonding in the residues before P14 (Aib10-P14, alamethicin; T11-P14, melittin) was characterized in both peptides by variable hydrogen bond patterns that included pi and gamma hydrogen bonds. The general agreement in hydrogen bond patterns determined in the simulations and from spectroscopic analysis indicates that with suitable conditions (including solvent composition and counterions where required), local hydrogen-bonded secondary structure in helical peptides may be predicted from dynamics simulations from alpha-helical starting structures. Each peptide, particularly alamethicin, underwent some large amplitude structural fluctuations in which several hydrogen bonds were cooperatively broken. The recovery of the persistent hydrogen bonding patterns after these fluctuations demonstrates the stability of intramolecular hydrogen-bonded secondary structure in methanol (consistent with spectroscopic observations), and is promising for simulations on extended timescales to characterize the nature of the backbone fluctuations that underlie amide exchange from isolated helical polypeptides.     

Simulation Studies on the Stabilities of Aggregates Formed by Fibril-Forming Segments of α-Synuclein

Abstract

We performed molecular dynamics simulations for various oligomers with different β-sheet conformations consisting of α-Synuclein 71–82 residues using an all atom force field and explicit water model. Tetramers of antiparallel β-sheet are shown to be stable, whereas parallel sheets are highly unstable due to the repulsive interactions between bulky and polar side chains as well as the weaker backbone hydrogen bonds. We also investigated the stabilities of double antiparallel β-sheets stacked with asymmetric and symmetric geometries. Our results show that this 12 amino acid residue peptide can form stable β-sheet conformers at 320K and higher temperatures. The backbone hydrogen bonds in β-sheet and the steric packing between hydrophobic side chains between β-sheets are shown to give conformational stabilities.     

Stabilities of disulfide bond intermediates in the folding of apamin.

Apamin is an 18-residue bee venom peptide with the sequence CNCKAPETALCARRCQQH-amide and contains 2 disulfide bonds connecting C-1 to C-11 and C-3 to C-15. In the folding of reduced, unfolded apamin to native apamin with two disulfide bonds, the one-disulfide folding intermediate states are not populated to significant levels. To study the properties of the one-disulfide intermediates, we have synthesized two peptide models to mimic the one-disulfide intermediates, Apa-1 and Apa-2, in which two cysteines in the sequence have been replaced by alanines. These peptides can form only one of the native disulfide bonds, C-1 to C-11 in the case of Apa-1 and C-3 to C-15 in the case of Apa-2. The stabilities of these disulfide bonds have been measured as a function of pH, concentration of urea, and temperature, in order to understand which contributions stabilize the disulfide-bonded structures. Using oxidized and reduced glutathione, the equilibrium constants for forming the disulfide bonds at 25 degrees C and pH 7.0 are 0.018 M for Apa-1 and 0.033 M for Apa-2 and show little dependence on pH or temperature. Both disulfide bonds are destabilized slightly (by approximately a factor of 2) between 0 and 8 M urea. Circular dichroism spectra indicate that although both Apa-1 and Apa-2 exhibit some structure, Apa-2 exhibits more than Apa-1. The results suggest that in the folding of apamin, the one-disulfide intermediate containing the C-3 to C-15 disulfide bond, as in Apa-2, is favored slightly. Secondary structure provides modest stabilization to this intermediate.     

Role of disulfide bonds in goose-type lysozyme

The role of the two disulfide bonds (Cys4-Cys60 and Cys18-Cys29) in the activity and stability of goose-type (G-type) lysozyme was investigated using ostrich egg-white lysozyme as a model. Each of the two disulfide bonds was deleted separately or simultaneously by substituting both Cys residues with either Ser or Ala. No remarkable differences in secondary structure or catalytic activity were observed between the wild-type and mutant proteins. However, thermal and guanidine hydrochloride unfolding experiments revealed that the stabilities of mutants lacking one or both of the disulfide bonds were significantly decreased relative to those of the wild-type. The destabilization energies of mutant proteins agreed well with those predicted from entropic effects in the denatured state. The effects of deleting each disulfide bond on protein stability were found to be approximately additive, indicating that the individual disulfide bonds contribute to the stability of G-type lysozyme in an independent manner. Under reducing conditions, the thermal stability of the wild-type was decreased to a level nearly equivalent to that of a Cys-free mutant (C4S/C18S/C29S/C60S) in which all Cys residues were replaced by Ser. Moreover, the optimum temperature of the catalytic activity for the Cys-free mutant was downshifted by about 20 degrees C as compared with that of the wild-type. These results indicate that the formation of the two disulfide bonds is not essential for the correct folding into the catalytically active conformation, but is crucial for the structural stability of G-type lysozyme.     

Cellular quality control screening to identify amino acid pairs for substituting the disulfide bonds in immunoglobulin fold domains

We are interested in determining which amino acid pairs can be substituted for the disulfide (S-S) bonds in proteins without disrupting their native structures under physiological conditions. In this study, we focused on the intradomain S-S bonds in Ig fold domains and aimed to determine a simple rule for replacement of their S-S bonds. The cysteines of four different Ig fold domains were mutated randomly, and the amino acid pairs substituted for the S-S bonds were screened by the method utilizing a cellular quality control system. Among the 36 selected mutants, 31 were natively folded without S-S bonds, as judged from the cooperativity of thermal unfolding. In addition, the selected mutant llama heavy chain antibodies retained antigen-binding affinity. At least two of the pairs Ala:Ala, Ala:Val, Val: Ala, and Val:Val were found in the selected mutants for all four different Ig fold domains, and they were stably folded at 30 degrees C. This suggests that examination of these four pairs could be enough to obtain natively folded Ig fold domains without S-S bonds.     

Effects of Turn Stability and Side-Chain Hydrophobicity on the Folding of β-Structures

Key elements of β-structure folding include hydrophobic core collapse, turn formation, and assembly of backbone hydrogen bonds. In the present folding simulations of several β-hairpins and β-sheets (peptide 1, protein G B1 domain peptide, TRPZIP2, TRPZIP4, 20mer, and 20mer^DP6D), the folding free-energy landscape as a function of several reaction coordinates corresponding to the three key elements indicates apparent dependence on turn stability and side-chain hydrophobicity, which demonstrates different folding mechanisms of similar β-structures of varied sequences. Turn stability is found to be the key factor in determining the formation order of the three structural elements in the folding of β-structures. Moreover, turn stability and side-chain hydrophobicity both affect the stability of backbone hydrogen bonds. The three-stranded β-sheets fold through a three-state transition in which the formation of one hairpin always takes precedence over the other. The different stabilities of two anti-parallel hairpins in each three-stranded β-sheet are shown to correlate well with the different levels of their hydrophobic interactions.     

Folding mechanism of beta-hairpins studied by replica exchange molecular simulations

The folding process of trpzip2 beta-hairpin is studied by the replica exchange molecular dynamics (REMD) and normal MD simulations, aiming to understand the folding mechanism of this unique small, stable, and fast folder, as well as to reveal the general principles in the folding of beta-hairpins. According to our simulations, the TS ensemble is mainly characterized by a largely formed turn and the interaction between the inner pair of hydrophobic core residues. The folding is a zipping up of hydrogen bonds. However, the nascent turn has to be stabilized by the partially formed hydrophobic core to cross the TS. Thus our folding picture is in essence a blend of hydrogen bond-centric and hydrophobic core-centric mechanism. Our simulations provide a direct evidence for a very recent experiment (Du et al., Proc Natl Acad Sci USA 2004;101:15915-15920), which suggests that the turn formation is the rate-limiting step for beta-hairpin folding and the unfolding is mainly determined by the hydrophobic interactions. Besides, the relationship between hydrogen bond stabilities and their relative importance in folding are investigated. It is found that the hydrogen bonds with higher stabilities need not play more important roles in the folding process, and vice versa.     

Role of intermolecular disulfide bonds of the organic solvent-stable PST-01 protease in its organic solvent stability

The PST-01 protease is secreted by the organic solvent-tolerant microorganism Pseudomonas aeruginosa PST-01 and is stable in the presence of various organic solvents. Therefore, the PST-01 strain and the PST-01 protease are very useful for fermentation and reactions in the presence of organic solvents, respectively. The organic solvent-stable PST-01 protease has two disulfide bonds (between Cys-30 and Cys-58 and between Cys-270 and Cys-297) in its molecule. Mutant PST-01 proteases in which one or both of the disulfide bonds were deleted were constructed by site-directed mutagenesis, and the effect of the disulfide bonds on the activity and the various stabilities was investigated. The disulfide bond between Cys-270 and Cys-297 in the PST-01 protease was found to be essential for its activity. The disulfide bond between Cys-30 and Cys-58 played an important role in the organic solvent stability of the PST-01 protease.     

Role of Intermolecular Disulfide Bonds of the Organic Solvent-Stable PST-01 Protease in Its Organic Solvent Stability

The PST-01 protease is secreted by the organic solvent-tolerant microorganism Pseudomonas aeruginosa PST-01 and is stable in the presence of various organic solvents. Therefore, the PST-01 strain and the PST-01 protease are very useful for fermentation and reactions in the presence of organic solvents, respectively. The organic solvent-stable PST-01 protease has two disulfide bonds (between Cys-30 and Cys-58 and between Cys-270 and Cys-297) in its molecule. Mutant PST-01 proteases in which one or both of the disulfide bonds were deleted were constructed by site-directed mutagenesis, and the effect of the disulfide bonds on the activity and the various stabilities was investigated. The disulfide bond between Cys-270 and Cys-297 in the PST-01 protease was found to be essential for its activity. The disulfide bond between Cys-30 and Cys-58 played an important role in the organic solvent stability of the PST-01 protease.     

Stability engineering of the human antibody repertoire

Human monoclonal antibodies often display limited thermodynamic and colloidal stabilities. This behavior hinders their production, and places limitations on the development of novel formulation conditions and therapeutic applications. Antibodies are highly diverse molecules, with much of the sequence variation observed within variable domain families and, in particular, their complementarity determining regions. This has complicated the development of comprehensive strategies for the stability engineering of the human antibody repertoire. Here we provide an overview of the field, and discuss recent advances in the development of robust and aggregation resistant antibody therapeutics.     

Mass spectrometric localization of carbon-carbon double bonds: A critical review of recent methods

Current methods are reviewed for determining the position of double bonds in fatty acids, and other unsaturated organic compounds, using mass spectrometry. ‘On-site’ and ‘remote-site’ derivatization methods are described, and their advantages and disadvantages for mass spectrometric analysis discussed. Chemical transformation of double bonds by methoxylation, silyloxylation or deuteration, together with electron impact (EI), chemical ionization (CI) or collisionally induced decomposition (CID) techniques in combination with fast atom bombardment (FAB) or CI, are found to be most suitable for polyunsaturated fatty acids (PUFA): for the analysis of less unsaturated compounds on a submicrogram scale those methods are most promising which either do not involve derivatization of the double bonds, or give derivatives in quantitative yields. Effects of mass spectrometer geometry and operating conditions are also considered.     

Hydrogen bond stabilities in membrane-reconstituted alamethicin from amide-resolved hydrogen-exchange measurements.

Amide-resolved hydrogen-deuterium exchange-rate constants were measured for backbone amides of alamethicin reconstituted in dioleoylphosphatidylcholine vesicles by an exchange-trapping method combined with high-resolution nuclear magnetic resonance spectroscopy. In vesicles containing alamethicin at molar ratios between 1:20 and 1:100 relative to lipid, the exchange-rate constants increased with increasing volume of the D20 buffer in which the vesicles were suspended, indicating that exchange under these conditions is dominated by partitioning of the peptide into the aqueous phase. This was supported by observation of a linear relationship between the exchange-rate constants for amides in membrane-reconstituted alamethicin and those for amides in alamethicin dissolved directly into D2O buffer. Significant protection of amides from exchange with D2O buffer in membrane-reconstituted alamethicin is interpreted in terms of stabilization by helical hydrogen bonding. Under conditions in which amide exchange occurred by partitioning of the peptide into solution, only lower limits for hydrogen-bond stabilities in the membrane were determined; all the potentially hydrogen-bonded amides of alamethicin are at least 1000-fold exchange protected in the membrane-bound state. When partitioning of alamethicin into the aqueous phase was suppressed by hydration of reconstituted vesicles in a limiting volume of water [D2O:dioleoylphosphatidylcholine:alamethicin; 220:1:0.05; (M:M:M)], the exchange-protection factors exhibited helical periodicity with highly exchange-protected, and less well-protected, amides on the nonpolar and polar helix faces, respectively. The exchange data indicate that, under the conditions studied, alamethicin adopts a stable helical structure in DOPC bilayers in which all the potentially hydrogen-bonded amides are stabilized by helical hydrogen bonds. The protection factors define the orientation of the peptide helix with respect to an aqueous phase, which is either the bulk solution or water within parallel or antiparallel transmembrane arrays of reconstituted alamethicin.     

Design and selection of nitrogen-rich bridged di-1,3,5-triazine derivatives with high energy and reduced sensitivity

The heats of formation (HOFs), electronic structures, energetic properties, and thermal stabilities of a series of energetic bridged di-1,3,5-triazine derivatives with different substituents and linkages were studied using density functional theory. It was found that the groups -N(3) and -N=N- are effective structural units for improving the HOF values of the di-1,3,5-triazine derivatives. The effects of the substituents on the HOMO-LUMO gap combine with those of the bridge groups. The calculated detonation velocities and detonation pressures indicate that substituting the -ONO(2), -NF(2), or -N=N- group is very useful for enhancing the detonation performance of these derivatives. Analysis of the bond dissociation energies for several relatively weak bonds suggests that most of the derivatives have good thermal stability. On the whole, the -NH(2), -N(3), -NH-, and -CH=CH- groups are effective structural units for increasing the thermal stabilities of the derivatives. Based on detonation performance and thermal stability, nine of the compounds can be considered potential candidates for high energy density materials with reduced sensitivity.     

Phenylalanyl-tRNA synthetase from E. coli: synergistic coupling between the sites for binding of L-phenylalanine and ATP

Binding of small ligands to the sites binding L-phenylalanine and ATP was measured by fluorescence titration technique. It is found that complex formation is not independent under conditions when both types of ligands are present in solution. The coupling is expressed as a synergistic binding rendering higher stabilities for complexes with ligand couples than expected on basis of separate binding of each ligand. In contrast, the substrate couple L-phenylalanine — ATP does not exhibit synergistic binding.     

Melting and chemical modification of a cyclized self-splicing group I intron: similarity of structures in 1 M Na+, in 10 mM Mg2+, and in the presence of substrate

C IVS is the cyclized form of the intron from the RNA precursor of the Tetrahymena thermophila large subunit (LSU) ribosomal RNA. C IVS was mapped by chemical modification in 1 M Na+, 0.05 M Na+ and 10 mM Mg2+ (Na+/Mg2+), and Na+/Mg2+ with CUCU substrate. The results suggest the secondary structure is similar for all three conditions. Optical melting curves were also measured for C IVS in 1 M Na+ and Na+/Mg2+ and indicate the secondary structures have similar stabilities under both conditions. Computer predictions of secondary structure and stability are in good agreement with observations. The results suggest that many of the approximations used for computer prediction of secondary structure by free energy minimization are reasonable.     

Stable tRNA precursors in HeLa cells.

Two tRNA precursors were isolated from 32P-labeled or unlabeled HeLa cells by two dimensional polyacrylamide gel electrophoresis, and were sequenced. These were the precursors of tRNAMet and tRNALeu, and both contained four extra nucleotides including 5'-triphosphates at their 5'-end and nine extra nucleotides including oligo U at their 3'-end. These RNAs are the first naturally occurring tRNA precursors from higher eukaryotes whose sequences have been determined. In these molecules, several modified nucleosides such as m2G, t6A and ac4C in mature tRNAs were undermodified. Two additional hydrogen bonds were formed in the clover leaf structures of these tRNA precursors. These extra hydrogen bonds may be responsible for the stabilities of these tRNA precursors.