Crystal structure of the E. coli tRNA aminoacyl stem isoacceptor RR-1660 at 2.0 Å resolution

Due to the redundancy of the genetic code there exist six mRNA codons for arginine and several tRNA^Arg isoacceptors which translate these triplets to protein within the context of the mRNA. The tRNA identity elements assure the correct aminoacylation of the tRNA with the cognate amino acid by the aminoacyl-tRNA-synthetases. In tRNA^Arg, the identity elements consist of the anticodon, parts of the D-loop and the discriminator base. The minor groove of the acceptor stem interacts with the arginyl-tRNA-synthetase. We crystallized different Escherichia coli tRNA^Arg acceptor stem helices and solved the structure of the tRNA^Arg isoacceptor RR-1660 microhelix by X-ray structure analysis. The acceptor stem helix crystallizes in the space group P1 with the cell constants a = 26.28, b = 28.92, c = 29.00 Å, α = 105.74, β = 99.01, γ = 97.44° and two molecules per asymmetric unit. The RNA hydration pattern consists of 88 bound water molecules. Additionally, one glycerol molecule is bound within the interface of the two RNA molecules.     

Structure and dynamics of α‐MSH using DRISM integral equation theory and stochastic dynamics

The structural and dynamical features of the hormone α‐MSH in solution have been examined over a 100 ns time scale by using free energy molecular mechanics models at room temperature. The free energy surface has been modeled using methods from integral equation theory and the dynamics by the Langevin equation. A modification of the accessible surface area friction drag model was used to calculate the atomic friction coefficients. The molecule shows a stable β‐turn conformation in the message region and a close interaction between the side chains of His⁶, Phe⁷, and Trp⁹. A salt bridge between Glu⁵ and Arg⁸ was found not to be a preferred interaction, whereas a Glu⁵ and Lys¹¹ salt bridge was not sampled, presumably due to relatively high free energy barriers. The message region was more conformationally rigid than the N‐terminal region. Several structural features observed here agree well with experimental results. The conformational features suggest a receptor–hormone interaction model where the hydrophobic side chains of Phe⁷ and Trp⁹ interact with the transmembrane portion of the MC1 receptor. Also, the positively charged side chain of Arg⁸ and the imidazole side chain of His⁶ may interact with the negatively charged portions of the receptor which may even be on the receptor's extracellular loops. © 1999 John Wiley & Sons, Inc. Biopoly 50: 255–272, 1999     

Dynamics and cooperativity of Trp-cage folding

In this study, multiple independent molecular dynamics (MD) simulations on Trp-cage folding were performed at 300, 325 and 375 K using generalized Born (GB) implicit solvent model. The orientational movement of the side-chain of Trp6 to form a hydrophobic core with 3₁₀-helix was observed. The breaking/formation of a salt bridge between Asp9 and Arg16 was proposed to be the prerequisite for Trp-cage folding/refolding. Our results demonstrate that the cooperation between the salt bridge and the Trp6 orientation leads to a stable tertiary structure of Trp-cage. Analyses on backbone concerted motions at different temperatures indicate that interactions between Trp6 and 3₁₀-helix & Pro18 and between Pro12 and Pro17 & Pro18 are weakened at 375 K but strengthened at lower temperatures, suggesting that they could be the potential driving force of hydrophobic collapse.     

Two Salt Bridges Differentially Contribute to the Maintenance of Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) Channel Function

Previous studies have identified two salt bridges in human CFTR chloride ion channels, Arg³⁵²-Asp⁹⁹³ and Arg³⁴⁷-Asp⁹²⁴, that are required for normal channel function. In the present study, we determined how the two salt bridges cooperate to maintain the open pore architecture of CFTR. Our data suggest that Arg³⁴⁷ not only interacts with Asp⁹²⁴ but also interacts with Asp⁹⁹³. The tripartite interaction Arg³⁴⁷-Asp⁹²⁴-Asp⁹⁹³ mainly contributes to maintaining a stable s2 open subconductance state. The Arg³⁵²-Asp⁹⁹³ salt bridge, in contrast, is involved in stabilizing both the s2 and full (f) open conductance states, with the main contribution being to the f state. The s1 subconductance state does not require either salt bridge. In confirmation of the role of Arg³⁵² and Asp⁹⁹³, channels bearing cysteines at these sites could be latched into a full open state using the bifunctional cross-linker 1,2-ethanediyl bismethanethiosulfonate, but only when applied in the open state. Channels remained latched open even after washout of ATP. The results suggest that these interacting residues contribute differently to stabilizing the open pore in different phases of the gating cycle.     

Mixed modes in opening of KcsA potassium channel from a targeted molecular dynamics simulation

Potassium channels conduct K⁺ flow selectively across the membrane through a central pore. During a process called gating, the potassium channels undergo a conformational change that opens or closes the ion-conducting pore. The potassium channel KcsA has been structurally determined in its closed state. However, the dynamic mechanism of the gating transition of the KcsA channel is still being investigated. Here, a targeted molecular dynamics simulation up to 150 ns is performed to investigate the detailed opening process of the KcsA channel with an open Kv1.2 structure serving as the target. The channel arrived at a self-determined quasi-stable state within 60 ns. The rigid-body and hinge-bending modes are observed mixed together in the remaining 90 ns long quasi-stable state. The mixed-mode movement seems come from the competition between the helix rigidity and the biased-applied gating force.     

A Revised Mechanism for the Activation of Complement C3 to C3b: A MOLECULAR EXPLANATION OF A DISEASE-ASSOCIATED POLYMORPHISM*

The solution structure of complement C3b is crucial for the understanding of complement activation and regulation. C3b is generated by the removal of C3a from C3. Hydrolysis of the C3 thioester produces C3u, an analog of C3b. C3b cleavage results in C3c and C3d (thioester-containing domain; TED). To resolve functional questions in relation to C3b and C3u, analytical ultracentrifugation and x-ray and neutron scattering studies were used with C3, C3b, C3u, C3c, and C3d, using the wild-type allotype with Arg¹⁰². In 50 mm NaCl buffer, atomistic scattering modeling showed that both C3b and C3u adopted a compact structure, similar to the C3b crystal structure in which its TED and macroglobulin 1 (MG1) domains were connected through the Arg¹⁰²–Glu¹⁰³² salt bridge. In physiological 137 mm NaCl, scattering modeling showed that C3b and C3u were both extended in structure, with the TED and MG1 domains now separated by up to 6 nm. The importance of the Arg¹⁰²–Glu¹⁰³² salt bridge was determined using surface plasmon resonance to monitor the binding of wild-type C3d(E1032) and mutant C3d(A1032) to immobilized C3c. The mutant did not bind, whereas the wild-type form did. The high conformational variability of TED in C3b in physiological buffer showed that C3b is more reactive than previously thought. Because the Arg¹⁰²-Glu¹⁰³² salt bridge is essential for the C3b-Factor H complex during the regulatory control of C3b, the known clinical associations of the major C3S (Arg¹⁰²) and disease-linked C3F (Gly¹⁰²) allotypes of C3b were experimentally explained for the first time.     

Characterization of the overall and local dynamics of a protein with intermediate rotational anisotropy: Differentiating between conformational exchange and anisotropic diffusion in the B3 domain of protein G

Because the overall tumbling provides a major contribution to protein spectral densities measured in solution, the choice of a proper model for this motion is critical for accurate analysis of protein dynamics. Here we study the overall and backbone dynamics of the B3 domain of protein G using ¹⁵N relaxation measurements and show that the picture of local motions is markedly dependent on the model of overall tumbling. The main difference is in the interpretation of the elevated R ₂ values in the -helix: the isotropic model results in conformational exchange throughout the entire helix, whereas no exchange is predicted by anisotropic models that place the longitudinal axis of diffusion tensor almost parallel to the helix axis. Due to small size (fast tumbling) of the protein, the T ₁ values have low sensitivity to NH bond orientation. The diffusion tensor derived from orientation dependence of R ₂/R ₁ is anisotropic (D _par/D _perp=1.4), with a small rhombic component. In order to distinguish the correct picture of motion, we apply model-independent methods that are sensitive to conformational exchange and do not require knowledge of protein structure or assumptions about its dynamics. A comparison of the CSA/dipolar cross-correlation rate constants with ¹⁵N relaxation rates and the estimation of R _ex terms from relaxation data at 9.4 and 14.1 T indicate no conformational exchange in the helix, in support of the anisotropic models. The experimentally derived diffusion tensor is in excellent agreement with theoretical predictions from hydrodynamic calculations; a detailed comparison with various hydrodynamic models revealed optimal parameters for hydrodynamic calculations.     

Cooperativity network of Trp‐cage miniproteins: probing salt‐bridges

Trp‐cage miniprotein was used to investigate the role of a salt‐bridge (Asp⁹–Arg¹⁶) in protein formation, by mutating residues at both sides, we mapped its contribution to overall stability and its role in folding mechanism. We found that both of the above side‐chains are also part of a dense interaction network composed of electrostatic, H‐bonding, hydrophobic, etc. components. To elucidate the fold stabilizing effects, we compared and contrasted electronic circular dichroism and NMR data of miniproteins equipped with a salt‐bridge with those of the salt‐bridge deleted mutants. Data were acquired both in neutral and in acidic aqueous solutions to decipher the pH dependency of both fully and partially charged partners. Our results indicate that the folding of Trp‐cage miniproteins is more complex than a simple two‐state process as we detected an intermediate state that differs significantly from the native fold. The intermediate formation is related to the salt‐bridge stabilization; in the miniprotein variants equipped with salt‐bridge the population of the intermediate state at acidic pH is significantly higher than it is for the salt‐bridge deleted mutants. In this molecular framework Arg¹⁶ stabilizes more than Asp⁹ does, because of its higher degree of 3D‐fold cooperation. In conclusion, the Xxx$^{9} \leftrightarrow$ Yyy¹⁶ salt‐bridge is not an isolated entity of this fold; rather it is an integrated part of a complex interaction network. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Redox-controlled backbone dynamics of human cytochrome c revealed by N NMR relaxation measurements

Redox-controlled backbone dynamics in cytochrome c (Cyt c) were revealed by 2D ¹⁵N NMR relaxation experiments. ¹⁵N T₁ and T₂ values and ¹H-¹⁵N NOEs of uniformly ¹⁵N-labeled reduced and oxidized Cyt c were measured, and the generalized order parameters (S²), the effective correlation time for internal motion (τ_e), the ¹⁵N exchange broadening contributions (R_ex) for each residue, and the overall correlation time (τ_m) were estimated by model-free dynamics formalism. These dynamic parameters clearly showed that the backbone dynamics of Cyt c are highly restricted due to the covalently bound heme that functions as the stable hydrophobic core. Upon oxidation of the heme iron in Cyt c, the average S² value was increased from 0.88 ± 0.01 to 0.92 ± 0.01, demonstrating that the mobility of the backbone is further restricted in the oxidized form. Such increases in the S² values were more prominent in the loop regions, including amino acid residues near the thioether bonds to the heme moiety and positively charged region around Lys87. Both of the regions are supposed to form the interaction site for cytochrome c oxidase (CcO) and the electron pathway from Cyt c to CcO. The redox-dependent mobility of the backbone in the interaction site for the electron transfer to CcO suggests an electron transfer mechanism regulated by the backbone dynamics in the Cyt c-CcO system.     

Functional roles of arginine residues in mung bean vacuolar H-pyrophosphatase

Plant vacuolar H⁺-translocating inorganic pyrophosphatase (V-PPase EC 3.6.1.1) utilizes inorganic pyrophosphate (PP_i) as an energy source to generate a H⁺ gradient potential for the secondary transport of ions and metabolites across the vacuole membrane. In this study, functional roles of arginine residues in mung bean V-PPase were determined by site-directed mutagenesis. Alignment of amino-acid sequence of K⁺-dependent V-PPases from several organisms showed that 11 of all 15 arginine residues were highly conserved. Arginine residues were individually substituted by alanine residues to produce R → A-substituted V-PPases, which were then heterologously expressed in yeast. The characteristics of mutant variants were subsequently scrutinized. As a result, most R → A-substituted V-PPases exhibited similar enzymatic activities to the wild-type with exception that R242A, R523A, and R609A mutants markedly lost their abilities of PP_i hydrolysis and associated H⁺-translocation. Moreover, mutation on these three arginines altered the optimal pH and significantly reduced K⁺-stimulation for enzymatic activities, implying a conformational change or a modification in enzymatic reaction upon substitution. In particular, R242A performed striking resistance to specific arginine-modifiers, 2,3-butanedione and phenylglyoxal, revealing that Arg²⁴² is most likely the primary target residue for these two reagents. The mutation at Arg²⁴² also removed F⁻ inhibition that is presumably derived from the interfering in the formation of substrate complex Mg²⁺-PP_i. Our results suggest accordingly that active pocket of V-PPase probably contains the essential Arg²⁴² which is embedded in a more hydrophobic environment.     

Sequence-dependent dynamical instability of the human prion protein: a comparative simulation study

The present study aimed to explore the most probable regions of the human prion protein backbone for which the initial steps of conformational transitions as a result of intrinsic and extrinsic perturbing factors on the protein structure can be assigned. A total of 0.3-μs molecular dynamics simulations on several analog structures of the protein have been performed. To mimic the impact of the extrinsic and intrinsic destructive parameters on the dynamical characteristics of the protein, mild acidic conditions and R208H mutation have been simulated. The findings indicated that distribution of conformational flexibilities along the protein chain was almost independent of the induced perturbing factors, and was mostly centralized on certain distinct parts of the structure comprising residues 132–145 and 187–203. Analyses also revealed that the segment comprising residues 187–203 may be considered as a peptide sequence, possessing high potential to start the initial steps of conformational rearrangements due to the induced physicochemical alterations. Sequence alignment and molecular dynamics data also revealed that segment 178–203 prefers to accommodate in extended structures rather than α-helices. Region 178–203 may be considered as a peptide switch capable of initiating the conformational transitions due to the introduced modifications and perturbing parameters.     

High-Field NMR and Circular Dichroism Solvent-Dependent Conformational Studies of the Bradykinin C-Terminal Tetrapeptide Ser-Pro-Phe-Arg

Abstract

The conformational properties of the tetrapeptide Ser¹-Pro²-Phe³-Arg⁴, the C-terminal fragment of the nonapeptide hormone bradykinin, have been studied by circular dichroism and two-dimensional NMR techniques. Measurements of coupling constants, NH temperature dependence rates and nuclear Overhauser effects (performed with rotating frame nuclear Overhauser spectroscopy, ROESY) in H₂O and CD₃OH/D₂O (80/20, v/v) reveal different conformations in the corresponding solvent. In aqueous solution the molecule exists in a random conformation or as an average of several conformations in rapid exchange. In CD₃OH/D₂O, however, the conformation is well-defined. The backbone of the peptide is extended, and the side-chains of Phe³ and Arg⁴ exhibit unusual rigidity for a peptide of this size. Evidently, the secondary structure is stabilized by a charge interaction between the guanidino group of Arg⁴ and the terminal carboxyl group, since experiments at various pH's show clearly that the definition of conformation decreases strongly upon protonation of the carboxyl function. A NH₃ ⁺(Ser¹)-COO^?(Arg⁴) salt bridge, as well as any form of turn stabilized by hydrogen bonds can be ruled out with certainty.     

Structure and flexibility of the thermophilic cold-shock protein of Thermus aquaticus

The thermophilic bacterium Thermus aquaticus is a well-known source of Taq polymerase. Here, we studied the structure and dynamics of the T. aquaticus cold-shock protein (Ta-Csp) to better understand its thermostability using NMR spectroscopy. We found that Ta-Csp has a five-stranded β-barrel structure with five salt bridges which are important for more rigid structure and a higher melting temperature (76 °C) of Ta-Csp compared to mesophilic and psychrophilic Csps. Microsecond to millisecond time scale exchange processes occur only at the β1–β2 surface region of the nucleic acid binding site with an average conformational exchange rate constant of 674 s⁻¹. The results imply that thermophilic Ta-Csp has a more rigid structure and may not need high structural flexibility to accommodate nucleic acids upon cold shock compared to its mesophile and psychrophile counterparts.     

A combined nuclear magnetic resonance and molecular dynamics study of the two structural motifs for mixed-linkage β-glucans: methyl β-cellobioside and methyl β-laminarabioside

The conformational and hydration properties of the two disaccharides methyl β-cellobioside and methyl β-laminarabioside were investigated by NMR spectroscopy and explicit solvation molecular dynamics simulations using the carbohydrate solution force field (CSFF). Adiabatic maps produced with this force field displayed 4 minima A: (Φ = 300°, Ψ = 280°), B: (Φ = 280°, Ψ = 210°), C: (Φ = 260°, Ψ = 60°), and D: (Φ = 60°, Ψ = 260°) for methyl β-cellobioside and 3 minima A: (Φ = 290°, Ψ = 130°), B: (Φ = 270°, Ψ = 290°), and C: (Φ = 60°, Ψ = 120°) for methyl β-laminarabioside. Molecular dynamics simulations were initiated from all minima. For each disaccharide, the simulation started from the A minimum was conducted for 50 ns, while the other minima were explored for 10 ns. The simulations revealed two stable minima for both compounds. For methyl β-cellobioside, the simulation minima in aqueous solution were shifted from their adiabatic map counterparts, while the simulation minima for methyl β-laminarabioside coincided with the corresponding adiabatic map minima. To validate the simulation results, NMR-derived NOEs and coupling constants across the glycoside linkage, ³J_HC and ³J_CH, were compared with values calculated from the MD trajectories. For each disaccharide, the best agreement was obtained for the simulations started at the A minimum. For both compounds, inter-ring water bridges in combination with the direct hydrogen bonds between the same groups were found to be determining factors for the overall solution structure of the disaccharides which differed from solid-state structures. Comparison with helical parameters showed that the preferred glycosidic dihedral configurations in the methyl β-cellobioside simulation were not highly compatible with the structure of cellulose, but that curdlan helix structures agreed relatively well with the methyl β-laminarabioside simulation. Polymers generated using glycosidic dihedral angles from the simulations revealed secondary structure motifs that that may help to elucidate polymer associations and small-molecule binding.     

Dynamic and Thermodynamic Response of the Ras Protein Cdc42Hs upon Association with the Effector Domain of PAK3

Human cell division cycle protein 42 (Cdc42Hs) is a small, Rho-type guanosine triphosphatase involved in multiple cellular processes through its interactions with downstream effectors. The binding domain of one such effector, the actin cytoskeleton-regulating p21-activated kinase 3, is known as PBD46. Nitrogen-15 backbone and carbon-13 methyl NMR relaxation was measured to investigate the dynamical changes in activated GMPPCP·Cdc42Hs upon PBD46 binding. Changes in internal motion of the Cdc42Hs, as revealed by methyl axis order parameters, were observed not only near the Cdc42Hs–PBD46 interface but also in remote sites on the Cdc42Hs molecule. The binding-induced changes in side-chain dynamics propagate along the long axis of Cdc42Hs away from the site of PBD46 binding with sharp distance dependence. Overall, the binding of the PBD46 effector domain on the dynamics of methyl-bearing side chains of Cdc42Hs results in a modest rigidification, which is estimated to correspond to an unfavorable change in conformational entropy of approximately − 10 kcal mol^− 1 at 298 K. A cluster of methyl probes closest to the nucleotide-binding pocket of Cdc42Hs becomes more rigid upon binding of PBD46 and is proposed to slow the catalytic hydrolysis of the γ phosphate moiety. An additional cluster of methyl probes surrounding the guanine ring becomes more flexible on binding of PBD46, presumably facilitating nucleotide exchange mediated by a guanosine exchange factor. In addition, the Rho insert helix, which is located at a site remote from the PBD46 binding interface, shows a significant dynamic response to PBD46 binding.     

Allosterism in human complement component 5a (hC5a): a damper of C5a receptor (C5aR) signaling

The phenomena of allosterism continues to advance the field of drug discovery, by illuminating gainful insights for many key processes, related to the structure–function relationships in proteins and enzymes, including the transmembrane G-protein coupled receptors (GPCRs), both in normal as well as in the disease states. However, allosterism is completely unexplored in the native protein ligands, especially when a small covalent change significantly modulates the pharmacology of the protein ligands toward the signaling axes of the GPCRs. One such example is the human C5a (^hC5a), the potent cationic anaphylatoxin that engages C5aR and C5L2 to elicit numerous immunological and non-immunological responses in humans. From the recently available structure–function data, it is clear that unlike the mouse C5a (^mC5a), the ^hC5a displays conformational heterogeneity. However, the molecular basis of such conformational heterogeneity, otherwise allosterism in ^hC5a and its precise contribution toward the overall C5aR signaling is not known. This study attempts to decipher the functional role of allosterism in ^hC5a, by exploring the inherent conformational dynamics in ^mC5a, ^hC5a and in its point mutants, including the proteolytic mutant des-Arg⁷⁴-^hC5a. Prima facie, the comparative molecular dynamics study, over total 500 ns, identifies Arg⁷⁴-Tyr²³ and Arg³⁷-Phe⁵¹ “cation-π” pairs as the molecular “allosteric switches” on ^hC5a that potentially functions as a damper of C5aR signaling.     

A Salt Bridge Linking the First Intracellular Loop with the C Terminus Facilitates the Folding of the Serotonin Transporter

The folding trajectory of solute carrier 6 (SLC6) family members is of interest because point mutations result in misfolding and thus cause clinically relevant phenotypes in people. Here we examined the contribution of the C terminus in supporting folding of the serotonin transporter (SERT; SLC6A4). Our working hypothesis posited that the amphipathic nature of the C-terminal α-helix (Thr⁶⁰³–Thr⁶¹³) was important for folding of SERT. Accordingly, we disrupted the hydrophobic moment of the α-helix by replacing Phe⁶⁰⁴, Ile⁶⁰⁸, or Ile⁶¹² by Gln. The bulk of the resulting mutants SERT-F604Q, SERT-I608Q, and SERT-I612Q were retained in the endoplasmic reticulum, but their residual delivery to the cell surface still depended on SEC24C. This indicates that the amphipathic nature of the C-terminal α-helix was dispensable to endoplasmic reticulum export. The folding trajectory of SERT is thought to proceed through the inward facing conformation. Consistent with this conjecture, cell surface expression of the misfolded mutants was restored by (i) introducing second site suppressor mutations, which trap SERT in the inward facing state, or (ii) by the pharmacochaperone noribogaine, which binds to the inward facing conformation. Finally, mutation of Glu⁶¹⁵ at the end of the C-terminal α-helix to Lys reduced surface expression of SERT-E615K. A charge reversal mutation in intracellular loop 1 restored surface expression of SERT-R152E/E615K to wild type levels. These observations support a mechanistic model where the C-terminal amphipathic helix is stabilized by an intramolecular salt bridge between residues Glu⁶¹⁵ and Arg¹⁵². This interaction acts as a pivot in the conformational search associated with folding of SERT.