Monitoring the effects of strong cosolvent hexafluoroisopropanol in investigation of the tetrameric structure and stability of K-channel KcsA

Adsorption of small chain alcohols into lipid membranes significantly changes the conformational states of intrinsic membrane proteins. In this study, the effects of membrane-active strong cosolvent hexafluoroisopropanol (HFIP) on the intrinsic tetrameric stability of potassium channel KcsA were investigated. Presence of acidic phosphatidylglycerol (PG) in non-bilayer phosphatidylethanolamine (PE) or bilayer phosphatidylcholine (PC) significantly increased the tetrameric stability compared to zwitterionic pure PC bilayers. The stabilizing effect of PG in both lipid bilayers was completely abolished upon deletion of the membrane-anchored N-terminus. Tryptophan fluorescence and circular dichroism experiments indicated that HFIP destabilizes the tetramer possibly via drastic changes in the lateral pressure profile close to the membrane-water interface. The data suggest that HFIP disturbs the ionic, H-bonding and hydrophobic interactions among KcsA subunits where N-terminus presumably plays a crucial role in determining the channel proper folding and tetrameric structure via ionic/H-bond interactions between the helix dipole and the membrane lipids.     

Disruption of a tight cluster surrounding tyrosine 131 in the native conformation of antithrombin III activates it for factor Xa inhibition

The native conformation of antithrombin III (ATIII) is a poor inhibitor of its coagulation pathway target enzymes because of the partial insertion of its reactive center loop (RCL) in its central A beta-sheet. This study focused on tyrosine 131, which is located at the helix D-sheet A interface, adjacent to the ATIII pentasaccharide and heparin cofactor-binding sites and some 17A away from the RCL insertion. Crystallographic structures show that the Tyr(131) ring is buried in native ATIII and then becomes exposed when pentasaccharide binds to the inhibitor and activates it. This change suggested that Tyr(131) might serve as a switch for ATIII conformational activation. The hypothesis is supported by results from this study, which progressively removed atoms from the Tyr(131) side chain. Rates of heparin-independent Y131L and Y131A factor Xa inhibition were 25 and 29 times faster than for the control and Y131F, suggesting that Tyr(131) ring interactions with neighboring helix D and strand 2A residues shift the uncatalyzed native-to-activated conformational equilibrium toward the RCL-inserted state. Thermal denaturation experiments showed Y131A and Y131L were less stable than the control and Y131F, implying an increased tendency toward A-sheet mobility in these genetically activated molecules. Thus, the tight Tyr(131)-Asn(127)-Leu(130)-Leu(140)-Ser(142) cluster at the helix D-strand 2A interface of native antithrombin contributes significantly to the stability of the ground state conformation, and tyrosine 131 serves as a heparin-responsive molecular switch during the allosteric activation of ATIII anticoagulant activity.     

Conformational changes in phospholipase A2 upon binding to micellar interfaces in the absence and presence of competitive inhibitors. A 1H and 15N NMR study.

An NMR study has been made of porcine pancreatic phospholipase A2 (PLA) in three environments: free in solution, in a binary complex with dodecylphosphocholine micelles, and in a ternary complex with a micelle and the substrate-like inhibitor (R)-1-octyl-2-(N-dodecanoylamino)-2-deoxyglycero-3-phosph oglycol. 1H and 15N chemical shifts, amide exchange rates, and NOE intensities are compared for the enzyme in different environments. From these data, structural differences are found for the N-terminal part, the end of the surface loop at residues Tyr69 and Thr70, and the active site residue His48, and also for the Ca-binding loop (residues 28-32). Specifically, when binding to a micelle, the side chains of residues Ala1, Trp3, and Tyr69, as well as all protons of Thr70, are found to be closer together. After subsequent introduction of the competitive inhibitor, further changes are found for these residues. The N-terminus is flexible in PLA free in solution, in contrast with the crystal structures where it adopts an alpha-helical conformation. According to the NMR data, this helix is rigidly formed only in the ternary complex. Furthermore, in the ternary complex, the N-terminal amino group and the exchangeable hydrogen at N3 of the ring of His48 are observed. We propose that PLA is activated in two steps. An initial conformational change occurs upon binding to a micellar interface. The catalytically active conformation of the enzyme, which has an extensive network of hydrogen bonds, is formed only when binding a substrate or competitive inhibitor at a lipid-water interface.     

Solution structure of porcine pancreatic phospholipase A2.

The lipolytic enzyme phospholipase A2 (PLA2) is involved in the degradation of high-molecular weight phospholipid aggregates in vivo. The enzyme has very high catalytic activities on aggregated substrates compared with monomeric substrates, a phenomenon called interfacial activation. Crystal structures of PLA2s in the absence and presence of inhibitors are identical, from which it has been concluded that enzymatic conformational changes do not play a role in the mechanism of interfacial activation. The high-resolution NMR structure of porcine pancreatic PLA2 free in solution was determined with heteronuclear multidimensional NMR methodology using doubly labeled 13C, 15N-labeled protein. The solution structure of PLA2 shows important deviations from the crystal structure. In the NMR structure the Ala1 alpha-amino group is disordered and the hydrogen bonding network involving the N-terminus and the active site is incomplete. The disorder observed for the N-terminal region of PLA2 in the solution structure could be related to the low activity of the enzyme towards monomeric substrates. The NMR structure of PLA2 suggests, in contrast to the crystallographic work, that conformational changes do play a role in the interfacial activation of this enzyme.     

Alteration of hydrogen bonding in the vicinity of histidine 48 disrupts millisecond motions in RNase A

The motion of amino acid residues on the millisecond (ms) time scale is involved in the tight regulation of catalytic function in numerous enzyme systems. Using a combination of mutational, enzymological, and relaxation-compensated (15)N Carr-Purcell-Meiboom-Gill (CPMG) methods, we have previously established the conformational significance of the distant His48 residue and the neighboring loop 1 in RNase A function. These studies suggested that RNase A relies on an intricate network of hydrogen bonding interactions involved in propagating functionally relevant, long-range ms motions to the catalytic site of the enzyme. To further investigate the dynamic importance of this H-bonding network, this study focuses on the individual replacement of Thr17 and Thr82 with alanine, effectively altering the key H-bonding interactions that connect loop 1 and His48 to the rest of the protein. (15)N CPMG dispersion studies, nuclear magnetic resonance (NMR) chemical shift analysis, and NMR line shape analysis of point mutants T17A and T82A demonstrate that the evolutionarily conserved single H-bond linking His48 to Thr82 is essential for propagating ms motions from His48 to the active site of RNase A on the time scale of catalytic turnover, whereas the T17A mutation increases the off rate and conformational exchange motions in loop 1. Accumulating evidence from our mutational studies indicates that residues experiencing conformational exchange in RNase A can be grouped into two separate clusters displaying distinct dynamical features, which appear to be independently affected by mutation. Overall, this study illuminates how tightly controlled and finely tuned ms motions are in RNase A, suggesting that designed modulation of protein motions may be possible.     

Linkage between the intramembrane H-bond network around aspartic acid 83 and the cytosolic environment of helix 8 in photoactivated rhodopsin

Understanding the coupling between conformational changes in the intramembrane domain and at the membrane-exposed surface of the bovine photoreceptor rhodopsin, a prototypical G protein-coupled receptor (GPCR), is crucial for the elucidation of molecular mechanisms in GPCR activation. Here, we have combined Fourier transform infrared (FTIR) and fluorescence spectroscopy to address the coupling between conformational changes in the intramembrane region around the retinal and the environment of helix 8, a putative cytosolic surface switch region in class-I GPCRs. Using FTIR/fluorescence cross-correlation we show specifically that surface alterations monitored by emission changes of fluorescein bound to Cys316 in helix 8 of rhodopsin are highly correlated with (i) H-bonding to Asp83 proximal of the retinal Schiff base but not to Glu122 close to the beta-ionone and (ii) with a metarhodopsin II (MII)-specific 1643 cm(-1) IR absorption change, indicative of a partial loss of secondary structure in helix 8 upon MII formation. These correlations are disrupted by limited C-terminal proteolysis but are maintained upon binding of a transducin alpha-subunit (G(talpha))-derived peptide, which stabilizes the MII state. Our results suggest that additional C-terminal cytosolic loop contacts monitored by an amide II absorption at 1557 cm(-1) play a functionally crucial role in keeping helix 8 in the position in which its environment is strongly coupled to the retinal-binding site near the Schiff base. In the intramembrane region, this coupling is mediated by the H-bonding network that connects Asp83 to the NPxxY(x)F motif preceding helix 8.     

Proposal of a new hydrogen-bonding form to maintain curdlan triple helix

Curdlan and other beta-1,3-D-glucans form right-handed triple helices, and it has been believed that the intermolecular H-bond is present at the center of the helix to maintain the structure. In this H-bond model, three secondary OH groups form an inequilateral hexagonal shape perpendicular to the helix axis. This hexagonal form seems to be characteristic for beta-1,3-D-glucans and is widely accepted. We carried out MOPAC and ab initio calculations for the curdlan helix, and we propose a new intermolecular H-bonding model. In our model, the H-bonds are formed between the O2-atoms on different x-y planes along the curdlan helix, hence the H-bonds are not perpendicular to the helix axis. The new H-bonds are connected along the helix, traversing three curdlan chains to make a left-handed helix. Therefore, the H-bonding array leads to a reverse helix of the main chain. According to our MOPAC calculation, this model is more stable than the previous one. We believe that the continuous H-bonding array is stabilized by cooperative phenomena in the polymeric system.     

Site-directed mutagenesis and X-ray crystallography of two phospholipase A2 mutants: Y52F and Y73F.

Tyr52 and Tyr73 are conserved amino acid residues throughout all vertebrate phospholipases A2. They are part of an extended hydrogen bonding system that links the N-terminal alpha-NH3(+)-group to the catalytic residues His48 and Asp99. These tyrosines were replaced by phenylalanines in a porcine pancreatic phospholipase A2 mutant, in which residues 62-66 had been deleted (delta 62-66PLA2). The mutations did not affect the catalytic properties of the enzyme, nor the folding kinetics. The stability against denaturation by guanidine hydrochloride was decreased, however. To analyse how the enzyme compensates for the loss of the tyrosine hydroxyl group, the X-ray structures of the delta Y52F and delta Y73F mutants were determined. After crystallographic refinement the final crystallographic R-factors were 18.1% for the delta Y52F mutant (data between 7 and 2.3 A resolution) and 19.1% for the delta Y73F mutant (data between 7 and 2.4 A resolution). No conformational changes occurred in the mutants compared with the delta 62-66PLA2, but an empty cavity formed at the site of the hydroxyl group of the former tyrosine. In both mutants the Asp99 side chain loses one of its hydrogen bonds and this might explain the observed destabilization.     

NMR determination of the global structure of the 113Cd derivative of desulforedoxin: investigation of the hydrogen bonding pattern at the metal center.

Desulforedoxin (Dx) is a simple homodimeric protein isolated from Desulfovibrio gigas (Dg) containing a distorted rubredoxin-like center with one iron coordinated by four cysteinyl residues (7.9 kDa with 36 amino acids per monomer). In order to probe the geometry and the H-bonding at the active site of Dx, the protein was reconstituted with 113Cd and the solution structure determined using 2D NMR methods. The structure of this derivative was initially compared with the NMR solution structure of the Zn form (Goodfellow BJ et al., 1996, J Biol Inorg Chem 1:341-353). Backbone amide protons for G4, D5, G13, L11 NH, and the Q14 NH side-chain protons, H-bonded in the X-ray structure, were readily exchanged with solvent. Chemical shift differences observed for amide protons near the metal center confirm the H-bonding pattern seen in the X-ray model (Archer M et al., 1995, J Mol Biol 251:690-702) and also suggest that H-bond lengths may vary between the Fe, Zn, and 113Cd forms. The H-bonding pattern was further probed using a heteronuclear spin echo difference (HSED) experiment; the results confirm the presence of NH-S H-bonds inferred from D2O exchange data and observed in the NMR family of structures. The presence of "H-bond mediated" coupling in Dx indicates that the NH-S H-bonds at the metal center have significant covalent character. The HSED experiment also identified an intermonomer "through space" coupling for one of the L26 methyl groups, indicating its proximity to the 113Cd center in the opposing monomer. This is the first example of an intermonomer "through space" coupling. Initial structure calculations produced subsets of NMR families with the S of C28 pointing away from or toward the L26 methyl: only the subset with the C28 sulfur pointing toward the L26 methyl could result in a "through space" coupling. The HSED result was therefore included in the structure calculations. Comparison of the Fe, Zn, and 113Cd forms of Dx suggests that the geometry of the metal center and the global fold of the protein does not vary to any great extent, although the H-bond network varies slightly when Cd is introduced. The similarity between the H-bonding pattern seen at the metal center in Dx, Rd (including H-bonded and through space-mediated coupling), and many zinc-finger proteins suggests that these H-bonds are structurally vital for stabilization of the metal centers in these proteins.     

Contribution of the hydrogen-bond network involving a tyrosine triad in the active site to the structure and function of a highly proficient ketosteroid isomerase from Pseudomonas putida biotype B

Delta(5)-3-Ketosteroid isomerase from Pseudomonas putida biotype B is one of the most proficient enzymes catalyzing an allylic isomerization reaction at rates comparable to the diffusion limit. The hydrogen-bond network (Asp99... Wat504...Tyr14...Tyr55...Tyr30) which links the two catalytic residues, Tyr14 and Asp99, to Tyr30, Tyr55, and a water molecule in the highly apolar active site has been characterized in an effort to identify its roles in function and stability. The DeltaG(U)(H2O) determined from equilibrium unfolding experiments reveals that the elimination of the hydroxyl group of Tyr14 or Tyr55 or the replacement of Asp99 with leucine results in a loss of conformational stability of 3.5-4.4 kcal/mol, suggesting that the hydrogen bonds of Tyr14, Tyr55, and Asp99 contribute significantly to stability. While decreasing the stability by about 6.5-7.9 kcal/mol, the Y55F/D99L or Y30F/D99L double mutation also reduced activity significantly, exhibiting a synergistic effect on k(cat) relative to the respective single mutations. These results indicate that the hydrogen-bond network is important for both stability and function. Additionally, they suggest that Tyr14 cannot function efficiently alone without additional support from the hydrogen bonds of Tyr55 and Asp99. The crystal structure of Y55F as determined at 1.9 A resolution shows that Tyr14 OH undergoes an alteration in orientation to form a new hydrogen bond with Tyr30. This observation supports the role of Tyr55 OH in positioning Tyr14 properly to optimize the hydrogen bond between Tyr14 and C3-O of the steroid substrate. No significant structural changes were observed in the crystal structures of Y30F and Y30F/Y55F, which allowed us to estimate approximately the interaction energies mediated by the hydrogen bonds Tyr30...Tyr55 and Tyr14...Tyr55. Taken together, our results demonstrate that the hydrogen-bond network provides the structural support that is needed for the enzyme to maintain the active-site geometry optimized for both function and stability.     

A conformational equilibrium in a protein fragment caused by two consecutive capping boxes: 1H-, 13C-NMR, and mutational analysis.

The conformational properties of an 18 residues peptide spanning the entire sequence, L1KTPA5QFDAD10ELRAA15MKG, of the first helix (A-helix) of domain 2 of annexin I, were thoroughly investigated. This fragment exhibits several singular features, and in particular, two successive potential capping boxes, T3xxQ6 and D8xxE11. The former corresponds to the native hydrogen bond network stabilizing the alpha helix N-terminus in the protein; the latter is a non-native capping box able to break the helix at residue D8, and is observed in the domain 2 partially folded state. Using 2D-NMR techniques, we showed that two main populations of conformers coexist in aqueous solution. The first corresponds to a single helix extending from T3 to K17. The second corresponds to a broken helix at residue Ds. Four mutants, T3A, F7A, D8A, and E11A, were designed to further analyze the role of key amino acids in the equilibrium between the two ensembles of conformers. The sensitivity of NMR parameters to account for the variations in the populations of conformers was evaluated for each peptide. Our data show the delta13Calpha chemical shift to be the most relevant parameter. We used it to estimate the population ratio in the various peptides between the two main ensembles of conformers, the full helix and the broken helix. For the WT, E11A, and F7A peptides, these ratios are respectively 35/65, 60/40, 60/40. Our results were compared to the data obtained from helix/coil transition algorithms.     

Conserved residues in the extracellular loops of short-wavelength cone visual pigments

The role of the extracellular loop region of a short-wavelength sensitive pigment, Xenopus violet cone opsin, is investigated via computational modeling, mutagenesis, and spectroscopy. The computational models predict a complex H-bonding network that stabilizes and connects the EC2-EC3 loop and the N-terminus. Mutations that are predicted to disrupt the H-bonding network are shown to produce visual pigments that do not stably bind chromophore and exhibit properties of a misfolded protein. The potential role of a disulfide bond between two conserved Cys residues, Cys(105) in TM3 and Cys(182) in EC2, is necessary for proper folding and trafficking in VCOP. Lastly, certain residues in the EC2 loop are predicted to stabilize the formation of two antiparallel β-strands joined by a hairpin turn, which interact with the chromophore via H-bonding or van der Waals interactions. Mutations of conserved residues result in a decrease in the level of chromophore binding. These results demonstrate that the extracellular loops are crucial for the formation of this cone visual pigment. Moreover, there are significant differences in the structure and function of this region in VCOP compared to that in rhodopsin.     

The Mechanism of Allosteric Inhibition of Protein Tyrosine Phosphatase 1B

As the prototypical member of the PTP family, protein tyrosine phosphatase 1B (PTP1B) is an attractive target for therapeutic interventions in type 2 diabetes. The extremely conserved catalytic site of PTP1B renders the design of selective PTP1B inhibitors intractable. Although discovered allosteric inhibitors containing a benzofuran sulfonamide scaffold offer fascinating opportunities to overcome selectivity issues, the allosteric inhibitory mechanism of PTP1B has remained elusive. Here, molecular dynamics (MD) simulations, coupled with a dynamic weighted community analysis, were performed to unveil the potential allosteric signal propagation pathway from the allosteric site to the catalytic site in PTP1B. This result revealed that the allosteric inhibitor compound-3 induces a conformational rearrangement in helix α7, disrupting the triangular interaction among helix α7, helix α3, and loop11. Helix α7 then produces a force, pulling helix α3 outward, and promotes Ser190 to interact with Tyr176. As a result, the deviation of Tyr176 abrogates the hydrophobic interactions with Trp179 and leads to the downward movement of the WPD loop, which forms an H-bond between Asp181 and Glu115. The formation of this H-bond constrains the WPD loop to its open conformation and thus inactivates PTP1B. The discovery of this allosteric mechanism provides an overall view of the regulation of PTP1B, which is an important insight for the design of potent allosteric PTP1B inhibitors.     

Influence of the 2'-hydroxyl group conformation on the stability of A-form helices in RNA

The 2'-hydroxyl group plays fundamental roles in both the structure and the function of RNA, and is the major determinant of the conformational and thermodynamic differences between RNA and DNA. Here, we report a conformational analysis of 2'-OH groups of the HIV-2 TAR RNA by means of NMR scalar coupling measurements in solution. Our analysis supports the existence of a network of water molecules spanning the minor groove of an RNA A-form helix, as has been suggested on the basis of a high-resolution X-ray study of an RNA duplex. The 2'-OH protons of the lower stem nucleotides of the TAR RNA project either towards the O3' or towards the base, where the 2'-OH group can favorably participate in H-bonding interactions with a water molecule situated in the nucleotide base plane. We observe that the k(ex) rate of the 2'-OH proton with the bulk solvent anti-correlates with the base-pair stability, confirming the involvement of the 2'-OH group in a collective network of H-bonds, which requires the presence of canonical helical secondary structure. The methodology and conformational analysis presented here are broadly applicable and facilitate future studies aimed to correlate the conformation of the 2'-OH group with both the structure and the function of RNA and RNA-ligand complexes.     

Conformationally linked interaction in the cytochrome bc(1) complex between inhibitors of the Q(o) site and the Rieske iron-sulfur protein

The modified Q cycle mechanism accounts for the proton and charge translocation stoichiometry of the bc(1) complex, and is now widely accepted. However the mechanism by which the requisite bifurcation of electron flow at the Q(o) site reaction is enforced is not clear. One of several proposals involves conformational gating of the docking of the Rieske ISP at the Q(o) site, controlled by the stage of the reaction cycle. Effects of different Q(o)-site inhibitors on the position of the ISP seen in crystals may reflect the same conformational mechanism, in which case understanding how different inhibitors control the position of the ISP may be a key to understanding the enforcement of bifurcation at the Q(o) site (Table?1). Here we examine the available structures of cytochrome bc(1) with different Q(o)-site inhibitors and different ISP positions to look for clues to this mechanism. The effect of ISP removal on binding affinity of the inhibitors stigmatellin and famoxadone suggest a "mutual stabilization" of inhibitor binding and ISP docking, however this thermodynamic observation sheds little light on the mechanism. The cd(1) helix of cytochrome b moves in such a way as to accommodate docking when inhibitors favoring docking are bound, but it is impossible with the current structures to say whether this movement of α-cd(1) is a cause or result of ISP docking. One component of the movement of the linker between E and F helices also correlates with the type of inhibitor and ISP position, and seems to be related to the H-bonding pattern of Y279 of cytochrome b. An H-bond from Y279 to the ISP, and its possible modulation by movement of F275 in the presence of famoxadone and related inhibitors, or its competition with an alternate H-bond to I269 of cytochrome b that may be destabilized by bound famoxadone, suggest other possible mechanisms. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.     

Analysis of protein loop closure. Two types of hinges produce one motion in lactate dehydrogenase

As shown in previous crystallographic investigations, upon binding lactate and NAD, lactate dehydrogenase undergoes a large conformational change that results in a surface loop moving roughly 10 A to cover the active site. In addition, there are appreciable movements (approximately 2 A) of five helices and three other loops. We demonstrate by a new fitting procedure that the loop moves on two hinges separated by a relatively rigid type II turn. The first hinge has few steric constraints on it, and its motion can be well accounted for by large changes in two torsion angles, i.e. as in a classic hinge motion. In contrast, the second hinge, which is part of a helix connected to the end of the loop, has many more constraints on it and distributes its deformation over more torsion angles. This novel motion involves the helix stretching and splitting into alpha-helical and 3(10)-helical components and substantial side-chain repacking in the sense of "cogs hopping between grooves" at its interface with the end of a neighboring helix. The loop is stabilized by five transverse (across loop) hydrogen bonds. These are preserved, through the conformational change and through 17 lactate dehydrogenase sequences, more than the longitudinal hydrogen bonds down the sides of the loop. Through a network of contacts, many of them conserved hydrophobic residues, the motion of the loop is propagated outward to structures that have no direct contact with the ligands. These moving structures are on the surface of the protein, and the whole protein can be subdivided into concentric shells of increasing mobility.     

Phospholipase A2 engineering. Structural and functional roles of highly conserved active site residues tyrosine-52 and tyrosine-73.

Site-directed mutagenesis was used to probe the structural and functional roles of two highly conserved residues, Tyr-52 and Tyr-73, in interfacial catalysis by bovine pancreatic phospholipase A2 (PLA2, overproduced in Escherichia coli). According to crystal structures, the side chains of these two active site residues form H-bonds with the carboxylate of the catalytic residue Asp-99. Replacement of either or both Tyr residues by Phe resulted in only very small changes in catalytic rates, which suggests that the hydrogen bonds are not essential for catalysis by PLA2. Substitution of either Tyr residue by nonaromatic amino acids resulted in substantial decreases in the apparent kcat toward 1,2-dioctanoyl-sn-glycero-3-phosphocholine (DC8PC) micelles and the v(o) (turnover number at maximal substrate concentration, i.e., mole fraction = 1) toward 1,2-dimyristoyl-sn-glycero-3-phosphomethanol (DC14PM) vesicles in scooting mode kinetics [Berg, O. G., Yu, B.-Z., Rogers, J., & Jain, M. K. (1991) Biochemistry 30, 7283-7297]. The Y52V mutant was further analyzed in detail by scooting mode kinetics: the E to E* equilibrium was examined by fluorescence; the dissociation constants of E*S, E*P, and E*I (KS*, KP*, and KI*, respectively) in the presence of Ca2+ were measured by protection of histidine-48 modification and by difference UV spectroscopy; the Michaelis constant KM* was calculated from initial rates of hydrolysis in the absence and presence of competitive inhibitors; and the turnover number under saturating conditions (kcat, which is a theoretical value since the enzyme may not be saturated at the interface) was calculated from the vo and KM* values. The results indicated little perturbation in the interfacial binding step (E to E*) but ca. 10-fold increases in KS*, KP*, KI*, and KM* and a less than 10-fold decrease in kcat. Such changes in the function of Y52V are not due to global conformational changes since the proton NMR properties of Y52V closely resemble those of wild-type PLA2; instead, it is likely to be caused by perturbed enzyme-substrate interactions at the active site. Tyr-73 appears to play an important structural role. The conformational stability of all Tyr-73 mutants decreased by 4-5 kcal/mol relative to that of the wild-type PLA2. The proton NMR properties of Y73A suggested significant conformational changes and substantially increased conformational flexibility. These detailed structural and functional analyses represent a major advancement in the structure-function study of an enzyme involved in interfacial catalysis.     

Structural dynamics of precursor and product of the RNA enzyme from the hepatitis delta virus as revealed by molecular dynamics simulations

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis. Site-specific self-cleavage requires cytidine in position 75 (C75). A precursor simulation with unprotonated C75 reveals a rather weak dynamic binding of C75 in the catalytic pocket with spontaneous, transient formation of a H-bond between U-1(O2') and C75(N3). This H-bond would be required for C75 to act as the general base. Upon protonation in the precursor, C75H+ has a tendency to move towards its product location and establish a firm H-bonding network within the catalytic pocket. However, a C75H+(N3)-G1(O5') H-bond, which would be expected if C75 acted as a general acid catalyst, is not observed on the present simulation timescale. The adjacent loop L3 is relatively dynamic and may serve as a flexible structural element, possibly gated by the closing U20.G25 base-pair, to facilitate a conformational switch induced by a protonated C75H+. L3 also controls the electrostatic environment of the catalytic core, which in turn may modulate C75 base strength and metal ion binding. We find that a distant RNA tertiary interaction involving a protonated cytidine (C41) becomes unstable when left unprotonated, leading to disruptive conformational rearrangements adjacent to the catalytic core. A Na ion temporarily compensates for the loss of the protonated hydrogen bond, which is strikingly consistent with the experimentally observed synergy between low pH and high Na+ concentrations in mediating residual self-cleavage of the HDV ribozyme in the absence of divalents.     

Fine-tuning of the stability of β-strands by Y181 in perfringolysin O directs the prepore to pore transition

Perfringolysin O (PFO) is a toxic protein that forms β-barrel transmembrane pores upon binding to cholesterol-containing membranes. The formation of lytic pores requires conformational changes in PFO that lead to the conversion of water-soluble monomers into membrane-bound oligomers. Although the general outline of stepwise pore formation has been established, the underlying mechanistic details await clarification. To extend our understanding of the molecular mechanisms that control the pore formation, we compared the hydrogen-deuterium exchange patterns of PFO with its derivatives bearing mutations in the D3 domain. In the case of two of these mutations F318A, Y181A, known from previous work to lead to a decreased lytic activity, global destabilization of all protein domains was observed in their water-soluble forms. This was accompanied by local changes in D3 β-sheet, including unexpected stabilization of functionally important β1 strand in Y181A. In case of the double mutation (F318A/Y181A) that completely abolished the lytic activity, several local changes were retained, but the global destabilization effects of single mutations were reverted and hydrogen-deuterium exchange (HDX) pattern returned to PFO level. Strong structural perturbations were not observed in case of remaining variants in which other residues of the hydrophobic core of D3 domain were substituted by alanine. Our results indicate the existence in PFO of a well-tuned H-bonding network that maintains the stability of the D3 β-strands at appropriate level at each transformation step. F318 and Y181 moieties participate in this network and their role extends beyond their direct intermolecular interaction during oligomerization that was identified previously.     

Molecular modeling of transmembrane helices 6 and 7 of the heptahelical lutropin receptor

In response to ligand binding and activating mutations, the lutropin receptor undergoes a conformational change to trigger a cellular response. D556 is the most common locus for naturally occurring activating mutations of the lutropin receptor, and a D556A mutant is shown to be constitutively active. A water-mediated proton transfer is postulated as part of the transmembrane signaling mechanism. Using energy minimization and ab initio calculations, a hydrogen bonding network involving a highly constrained water molecule(s) and D556 (helix 6) and N593/N597/Y601 (helix 7) is presented.