Inherent Dynamics of the Acid-Sensing Ion Channel 1 Correlates with the Gating Mechanism

The acid-sensing ion channel 1 (ASIC1) is a key receptor for extracellular protons. Although numerous structural and functional studies have been performed on this channel, the structural dynamics underlying the gating mechanism remains unknown. We used normal mode analysis, mutagenesis, and electrophysiological methods to explore the relationship between the inherent dynamics of ASIC1 and its gating mechanism. Here we show that a series of collective motions among the domains and subdomains of ASIC1 correlate with its acid-sensing function. The normal mode analysis result reveals that the intrinsic rotation of the extracellular domain and the collective motions between the thumb and finger induced by proton binding drive the receptor to experience a deformation from the extracellular domain to the transmembrane domain, triggering the channel pore to undergo “twist-to-open” motions. The movements in the transmembrane domain indicate that the likely position of the channel gate is around Leu440. These motion modes are compatible with a wide body of our complementary mutations and electrophysiological data. This study provides the dynamic fundamentals of ASIC1 gating.     

The requirement for mechanical coupling between head and S2 domains in smooth muscle myosin ATPase regulation and its implications for dimeric motor function

A combination of experimental structural data, homology modelling and elastic network normal mode analysis is used to explore how coupled motions between the two myosin heads and the dimerization domain (S2) in smooth muscle myosin II determine the domain movements required to achieve the inhibited state of this ATP-dependent molecular motor. These physical models rationalize the empirical requirement for at least two heptads of non-coiled alpha-helix at the junction between the myosin heads and S2, and the dependence of regulation on S2 length. The results correlate well with biochemical data regarding altered conformational-dependent solubility and stability. Structural models of the conformational transition between putative active states and the inhibited state show that torsional flexibility of the S2 alpha-helices is a key mechanical requirement for myosin II regulation. These torsional motions of the myosin heads about their coiled coil alpha-helices affect the S2 domain structure, which reciprocally affects the motions of the myosin heads. This inter-relationship may explain a large body of data on function of molecular motors that form dimers through a coiled-coil domain.     

Critical hydrophobic interactions between phosphorylation and actuator domains of Ca2+-ATPase for hydrolysis of phosphorylated intermediate

Functional roles of seven hydrophobic residues on the interface between the actuator (A) and phosphorylation (P) domains of sarcoplasmic reticulum Ca2+-ATPase were explored by alanine and serine substitutions. The residues examined were Ile179/Leu180/Ile232 on the A domain, Val705/Val726 on the P domain, and Leu119/Tyr122 on the loop linking the A domain and M2 (the second transmembrane helix). These residues gather to form a hydrophobic cluster around Tyr122 in the crystal structures of Ca2+-ATPase in Ca2+-unbound E2 (unphosphorylated) and E2P (phosphorylated) states but are far apart in those of Ca2+-bound E1 (unphosphorylated) and E1P (phosphorylated) states. The substitution-effects were also compared with those of Ile235 on the A domain/M3 linker and those of T181GE of the A domain, since they are in the immediate vicinity of the Tyr122-cluster. All these substitutions almost completely inhibited ATPase activity without inhibiting Ca2+-activated E1P formation from ATP. Substitutions of Ile235 and T181GE blocked the E1P to E2P transition, whereas those in the Tyr122-cluster blocked the subsequent E2P hydrolysis. Substitutions of Ile235 and Glu183 also blocked EP hydrolysis. Results indicate that the Tyr122-cluster is formed during the E1P to E2P transition to configure the catalytic site and position Glu183 properly for hydrolyzing the acylphosphate. Ile235 on the A domain/M3 linker likely forms hydrophobic interactions with the A domain and thereby allowing the strain of this linker to be utilized for large motions of the A domain during these processes. The Tyr122-cluster, Ile235, and T181GE thus seem to have different roles and are critical in the successive events in processing phosphorylated intermediates to transport Ca2+.     

The human estrogen receptor has transcriptional activator and repressor functions in the absence of ligand

Most studies on the cloned human estrogen receptor (hER) have been conducted with a mutant receptor in which Gly400 is changed to Val. Here we describe two novel regulatory functions of wild-type hER that are hormone independent: (i) a constitutive activator function and (ii) a repressor activity. Mutations in the hormone-binding domain, including the Val400 mutation, impair both of these functions. In addition, DNA binding is strongly reduced in the mutant receptors. The hormone-binding domain of the hER thus controls DNA binding (and thereby the repressor function) of the hER as well as its constitutive activator function. Moreover, we find that the antiestrogen tamoxifen restores the constitutive activator function, the DNA binding, and the repressor function of the Val400 mutant, but has no effect on the constitutive activator function or DNA binding of the wild-type hER.     

Conformational Dynamics and Ligand Binding in the Multi-Domain Protein PDC109

PDC109 is a modular multi-domain protein with two fibronectin type II (Fn2) repeats joined by a linker. It plays a major role in bull sperm binding to the oviductal epithelium through its interactions with phosphorylcholines (PhCs), a head group of sperm cell membrane lipids. The crystal structure of the PDC109-PhC complex shows that each PhC binds to the corresponding Fn2 domain, while the two domains are on the same face of the protein. Long timescale explicit solvent molecular dynamics (MD) simulations of PDC109, in the presence and absence of PhC, suggest that PhC binding strongly correlates with the relative orientation of choline-phospholipid binding sites of the two Fn2 domains; unless the two domains tightly bind PhCs, they tend to change their relative orientation by deforming the flexible linker. The effective PDC109-PhC association constant of 28 M, estimated from their potential of mean force is consistent with the experimental result. Principal component analysis of the long timescale MD simulations was compared to the significantly less expensive normal mode analysis of minimized structures. The comparison indicates that difference between relative domain motions of PDC109 with bound and unbound PhC is captured by the first principal component in the principal component analysis as well as the three lowest normal modes in the normal mode analysis. The present study illustrates the use of detailed MD simulations to clarify the energetics of specific ligand-domain interactions revealed by a static crystallographic model, as well as their influence on relative domain motions in a multi-domain protein.     

Conformational dynamics of a multidomain protein by neutron scattering and computational analysis

The flexible conformations of a multidomain protein are responsible for its biological functions. Although MurD, a 47-kDa protein that consists of three domains, sequentially changes its domain conformation from an open form to a closed form through a semiclosed form in its enzymatic reaction, the domain dynamics in each conformation remains unclear. In this study, we verify the conformational dynamics of MurD in the corresponding three states (apo and ATP- and inhibitor-bound states) with a combination of small-angle x-ray and neutron scattering (SAXS and SANS), dynamic light scattering (DLS), neutron backscattering (NBS), neutron spin echo (NSE) spectroscopy, and molecular dynamics (MD) simulations. Applying principal component analysis of the MD trajectories, twisting and open-closed domain modes are identified as the major collective coordinates. The deviations of the experimental SAXS profiles from the theoretical calculations based on the known crystal structures become smaller in the ATP-bound state than in the apo state, and a further decrease is evident upon inhibitor binding. These results suggest that domain motions of the protein are suppressed step by step of each ligand binding. The DLS and NBS data yield collective and self-translational diffusion constants, respectively, and we used them to extract collective domain motions in nanometer and nanosecond scales from the NSE data. In the apo state, MurD shows both twisting and open-closed domain modes, whereas an ATP binding suppresses twisting domain motions, and a further reduction of open-closed mode is seen in the inhibitor-binding state. These observations are consistent with the structure modifications measured by the small-angle scattering as well as the MD simulations. Such changes in the domain dynamics associated with the sequential enzymatic reactions should be related to the affinity and reaction efficiency with a ligand that binds specifically to each reaction state.     

Model-free methods of analyzing domain motions in proteins from simulation: A comparison of normal mode analysis and molecular dynamics simulation of lysozyme

Model-free methods are introduced to determine quantities pertaining to protein domain motions from normal mode analyses and molecular dynamics simulations. For the normal mode analysis, the methods are based on the assumption that in low frequency modes, domain motions can be well approximated by modes of motion external to the domains. To analyze the molecular dynamics trajectory, a principal component analysis tailored specifically to analyze interdomain motions is applied. A method based on the curl of the atomic displacements is described, which yields a sharp discrimination of domains, and which defines a unique interdomain screw-axis. Hinge axes are defined and classified as twist or closure axes depending on their direction. The methods have been tested on lysozyme. A remarkable correspondence was found between the first normal mode axis and the first principal mode axis, with both axes passing within 3 Å of the alpha-carbon atoms of residues 2, 39, and 56 of human lysozyme, and near the interdomain helix. The axes of the first modes are overwhelmingly closure axes. A lesser degree of correspondence is found for the second modes, but in both cases they are more twist axes than closure axes. Both analyses reveal that the interdomain connections allow only these two degrees of freedom, one more than provided by a pure mechanical hinge. Proteins 27:425–437, 1997. © 1997 Wiley-Liss, Inc.     

Steric interactions of valines 1, 5, and 7 in [valine 5, D-alanine 8] gramicidin A channels.

When the central valine residues 6, 7, and 8 of gramicidin A (gA) are shifted by one position, the resulting [Val(5), D-Ala(8)]gA forms right-handed channels with a single-channel conductance and average duration somewhat less than gA channels. The reduction in channel duration has been attributed to steric conflict between the side chains of Val(1) and Val(5) in opposing monomers (Koeppe, R. E. II, D. V. Greathouse, A. Jude, G. Saberwal, L. L. Providence, and O. S. Andersen. 1994. J. Biol. Chem. 269:12567-12576). To investigate the orientations and motions of valines in [Val(5), D-Ala(8)]gA, we have incorporated (2)H labels at Val 1, 5, or 7 and recorded (2)H-NMR spectra of oriented and nonoriented samples in hydrated dimyristoylphosphatidylcholine. Spectra of nonoriented samples at 4 degrees C reveal powder patterns that indicate rapid side chain "hopping" for Val(5), and an intermediate rate of hopping for Val(1) and Val(7) that is somewhat slower than in gA. Oriented samples of deuterated Val(1) and Val(7) show large changes in the methyl and C(beta)-(2)H quadrupolar splittings (Deltanu(q)) when Ala(5) in native gA is changed to Val(5). Three or more peaks for the Val(1) methyls with Deltanu(q) values that vary with the echo delay, together with an intermediate spectrum for nonoriented samples at 4 degrees C, suggest unusual side chain dynamics for Val(1) in [Val(5), D-Ala(8)]gA. These results are consistent with a steric conflict that has been introduced between the two opposing monomers. In contrast, the acylation of gA has little influence on the side chain dynamics of Val(1), regardless of the identity of residue 5.     

Cooperative nature of gating transitions in K+ channels as seen from dynamic importance sampling calculations

The growing dataset of K⁺ channel x‐ray structures provides an excellent opportunity to begin a detailed molecular understanding of voltage‐dependent gating. These structures, while differing in sequence, represent either a stable open or closed state. However, an understanding of the molecular details of gating will require models for the transitions and experimentally testable predictions for the gating transition. To explore these ideas, we apply dynamic importance sampling to a set of homology models for the molecular conformations of K⁺ channels for four different sets of sequences and eight different states. In our results, we highlight the importance of particular residues upstream from the Pro‐Val‐Pro (PVP) region to the gating transition. This supports growing evidence that the PVP region is important for influencing the flexibility of the S6 helix and thus the opening of the gating domain. The results further suggest how gating on the molecular level depends on intra‐subunit motions to influence the cooperative behavior of all four subunits of the K⁺ channel. We hypothesize that the gating process occurs in steps: first sidechain movement, then inter‐S5‐S6 subunit motions, and lastly the large‐scale domain rearrangements. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Mechanism of molecular interactions for tRNA(Val) recognition by valyl-tRNA synthetase

The molecular interactions between valyl-tRNA synthetase (ValRS) and tRNA(Val), with the C34-A35-C36 anticodon, from Thermus thermophilus were studied by crystallographic analysis and structure-based mutagenesis. In the ValRS-bound structure of tRNA(Val), the successive A35-C36 residues (the major identity elements) of tRNA(Val) are base-stacked upon each other, and fit into a pocket on the alpha-helix bundle domain of ValRS. Hydrogen bonds are formed between ValRS and A35-C36 of tRNA(Val) in a base-specific manner. The C-terminal coiled-coil domain of ValRS interacts electrostatically with A20 and hydrophobically with the G19*C56 tertiary base pair. The loss of these interactions by the deletion of the coiled-coil domain of ValRS increased the K(M) value for tRNA(Val) 28-fold and decreased the k(cat) value 19-fold in the aminoacylation. The tRNA(Val) K(M) and k(cat) values were increased 21-fold and decreased 32-fold, respectively, by the disruption of the G18*U55 and G19*C56 tertiary base pairs, which associate the D- and T-loops for the formation of the L-shaped tRNA structure. Therefore, the coiled-coil domain of ValRS is likely to stabilize the L-shaped tRNA structure during the aminoacylation reaction.     

Insights into Protein Aggregation by NMR Characterization of Insoluble SH3 Mutants Solubilized in Salt-Free Water

Protein aggregation in vivo has been extensively associated with a large spectrum of human diseases. On the other hand, mechanistic insights into protein aggregation in vitro were incomplete due to the inability in solubilizing insoluble proteins for high-resolution biophysical investigations. However, a new avenue may be opened up by our recent discovery that previously-thought insoluble proteins can in fact be solubilized in salt-free water. Here we use this approach to study the NMR structural and dynamic properties of an insoluble SH3 mutant with a naturally-occurring insertion of Val22 at the tip of the diverging turn. The obtained results reveal: 1) regardless of whether the residue is Val, Ala, Asp or Arg, the insertion will render the first hNck2 SH3 domain to be insoluble in buffers. Nevertheless, all four mutants could be solubilized in salt-free water and appear to be largely unfolded as evident from their CD and NMR HSQC spectra. 2) Comparison of the chemical shift deviations reveals that while in V22-SH3 the second helical region is similarly populated as in the wild-type SH3 at pH 2.0, the first helical region is largely unformed. 3) In V22-SH3, many non-native medium-range NOEs manifest to define non-native helical conformations. In the meanwhile a small group of native-like long-range NOEs still persists, indicating the existence of a rudimentary native-like tertiary topology. 4) Although overall, V22-SH3 has significantly increased backbone motions on the ps-ns time scale, some regions still own restricted backbone motions as revealed by analyzing ¹⁵N relaxation data. Our study not only leads to the establishment of the first high-resolution structural and dynamic picture for an insoluble protein, but also shed more light on the molecular events for the nonhierarchical folding mechanism. Furthermore, a general mechanism is also proposed for in vivo protein aggregation triggered by the genetic mutation and posttranslational modification.     

Filtering molecular dynamics trajectories to reveal low-frequency collective motions: phospholipase A2

A novel method for analysing molecular dynamics trajectories has been developed, which filters out high frequencies using digital signal processing techniques and facilitates focusing on the low-frequency collective motions of proteins. These motions involve low energy slow motions, which lead to important biological phenomena such as domain closure and allosteric effects in enzymes. The filtering method treats each of the atomic trajectories obtained from the molecular dynamics simulation as a "signal". The trajectories of each of the atoms in the system (or any subset of interest) are Fourier transformed to the frequency domain, a filtering function is applied and then an inverse transformation back to the time domain yields the filtered trajectory. The filtering method has been used to study the dynamics of the enzyme phospholipase A2. In the filtered trajectory, all the high frequency bond and valence angle vibrations were eliminated, leaving only low-frequency motion, mainly fluctuations in torsions and conformational transitions. Analysis of this trajectory revealed interesting motions of the protein, including concerted movements of helices, and changes in shape of the active site cavity. Unlike normal mode analysis, which has been used to study the motion of proteins, this method does not require converged minimizations or diagonalization of a matrix of second derivatives. In addition, anharmonicity, multiple minima and conformational transitions are treated explicitly. Thus, the filtering method avoids most of the approximations implicit in other investigations of the dynamic behaviour of large systems.     

Extracellular-Regulated Kinase 2 Is Activated by the Enhancement of Hinge Flexibility

Protein motions underlie conformational and entropic contributions to enzyme catalysis; however, relatively little is known about the ways in which this occurs. Studies of the mitogen-activated protein kinase ERK2 (extracellular-regulated protein kinase 2) by hydrogen-exchange mass spectrometry suggest that activation enhances backbone flexibility at the linker between N- and C-terminal domains while altering nucleotide binding mode. Here, we address the hypothesis that enhanced backbone flexibility within the hinge region facilitates kinase activation. We show that hinge mutations enhancing flexibility promote changes in the nucleotide binding mode consistent with domain movement, without requiring phosphorylation. They also lead to the activation of monophosphorylated ERK2, a form that is normally inactive. The hinge mutations bypass the need for pTyr but not pThr, suggesting that Tyr phosphorylation controls hinge motions. In agreement, monophosphorylation of pTyr enhances both hinge flexibility and nucleotide binding mode, measured by hydrogen-exchange mass spectrometry. Our findings demonstrate that regulated protein motions underlie kinase activation. Our working model is that constraints to domain movement in ERK2 are overcome by phosphorylation at pTyr, which increases hinge dynamics to promote the active conformation of the catalytic site.     

Domain motions in bacteriophage T4 lysozyme: A comparison between molecular dynamics and crystallographic data

A comparison of a series of extended molecular dynamics (MD) simulations of bacteriophage T4 lysozyme in solvent with X-ray data is presented. Essential dynamics analyses were used to derive collective fluctuations from both the simulated trajectories and a distribution of crystallographic conformations. In both cases the main collective fluctuations describe domain motions. The protein consists of an N- and C-terminal domain connected by a long helix. The analysis of the distribution of crystallographic conformations reveals that the N-terminal helix rotates together with either of these two domains. The main domain fluctuation describes a closure mode of the two domains in which the N-terminal helix rotates concertedly with the C-terminal domain, while the domain fluctuation with second largest amplitude corresponds to a twisting mode of the two domains, with the N-terminal helix rotating concertedly with the N-terminal domain. For the closure mode, the difference in hinge-bending angle between the most open and most closed X-ray structure along this mode is 49 degrees. In the MD simulation that shows the largest fluctuation along this mode, a rotation of 45 degrees was observed. Although the twisting mode has much less freedom than the closure mode in the distribution of crystallographic conformations, experimental results suggest that it might be functionally important. Interestingly, the twisting mode is sampled more extensively in all MD simulations than it is in the distribution of X-ray conformations. Proteins 31:116–127, 1998. © 1998 Wiley-Liss, Inc.     

Normal‐mode‐based modeling of allosteric couplings that underlie cyclic conformational transition in F1 ATPase

F₁ ATPase, a rotary motor comprised of a central stalk ( γ subunit) enclosed by three α and β subunits alternately arranged in a hexamer, features highly cooperative binding and hydrolysis of ATP. Despite steady progress in biophysical, biochemical, and computational studies of this fascinating motor, the structural basis for cooperative ATPase involving its three catalytic sites remains not fully understood. To illuminate this key mechanistic puzzle, we have employed a coarse‐grained elastic network model to probe the allosteric couplings underlying the cyclic conformational transition in F₁ ATPase at a residue level of detail. We will elucidate how ATP binding and product (ADP and phosphate) release at two catalytic sites are coupled with the rotation of γ subunit via various domain motions in α ₃ β ₃ hexamer (including intrasubunit hinge‐bending motions in β subunits and intersubunit rigid‐body rotations between adjacent α and β subunits). To this end, we have used a normal‐mode‐based correlation analysis to quantify the allosteric couplings of these domain motions to local motions at catalytic sites and the rotation of γ subunit. We have then identified key amino acid residues involved in the above couplings, some of which have been validated against past studies of mutated and γ ‐truncated F₁ ATPase. Our finding strongly supports a binding change mechanism where ATP binding to the empty catalytic site triggers a series of intra‐ and intersubunit domain motions leading to ATP hydrolysis and product release at the other two closed catalytic sites. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Envisioning the Loop Movements and Rotation of the Two Subdomains of Dihydrofolate Reductase by Elastic Normal Mode Analysis

Abstract

Normal mode analysis, using the elastic network model, has been employed to envision the low frequency normal mode motion trends in the structures of five intermediates and a transition state in the kinetic pathway of E. coli dihydrofolate reductase (DHFR). Five of the reaction pathway analog structures and a crystal structure resembling the transition state, using X-ray analyses determined by Kraut et al., have been adapted as structural models. The motions that poise pathways of the M20 loop transitions from closed to occluded conformations and sub domain rotation to close the substrate cleft, have been predicted and envisioned for the first time by this study. Pathway entries to the movement of the substrate binding cleft helices are also envisioned. These motions play roles in transition structure stabilization and in regulating the release of the product tetrahydrofolate (THF). The motions observed push the ground state conformation of each intermediate towards a higher energy sub state conformation. A set of conserved residues involved in the catalytic reactions and conformational changes, previously studied by kinetic, theoretical and NMR, have been analyzed. The importance of these motions in terms of protein dynamics are revealed and envisioned by the normal mode analysis. Additional residues are proposed as candidates for further study of their potential promotional function.     

Structural basis for double-sieve discrimination of L-valine from L-isoleucine and L-threonine by the complex of tRNA(Val) and valyl-tRNA synthetase

Valyl-tRNA synthetase (ValRS) strictly discriminates the cognate L-valine from the larger L-isoleucine and the isosteric L-threonine by the tRNA-dependent "double sieve" mechanism. In this study, we determined the 2.9 A crystal structure of a complex of Thermus thermophilus ValRS, tRNA(Val), and an analog of the Val-adenylate intermediate. The analog is bound in a pocket, where Pro(41) allows accommodation of the Val and Thr moieties but precludes the Ile moiety (the first sieve), on the aminoacylation domain. The editing domain, which hydrolyzes incorrectly synthesized Thr-tRNA(Val), is bound to the 3' adenosine of tRNA(Val). A contiguous pocket was found to accommodate the Thr moiety, but not the Val moiety (the second sieve). Furthermore, another Thr binding pocket for Thr-adenylate hydrolysis was suggested on the editing domain.     

Microsecond-to-millisecond conformational dynamics demarcate the GluR2 glutamate receptor bound to agonists glutamate, quisqualate, and AMPA

Chemical shift changes and internal motions on microsecond-to-millisecond time scales of the S1S2 ligand-binding domain of the GluR2 ionotropic glutamate receptor have been studied by NMR spectroscopy in the presence of the agonists glutamic acid (glutamate), quisqualic acid (quisqualate), and alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid (AMPA). Although the crystal structures of the three agonist-bound forms of GluR2 S1S2 ligand-binding domain are very similar, chemical shift changes imply that AMPA-bound GluR2 S1S2 is conformationally distinct from glutamate- and quisqualate-bound forms of GluR2 S1S2. NMR spin relaxation measurements for backbone amide (15)N nuclei reveal that GluR2 S1S2 exhibits reduced chemical exchange line broadening, resulting from microsecond-to-millisecond conformational dynamics, in AMPA-bound compared to glutamate- and quisqualate-bound states. The largest changes in line broadening are observed for two regions of GluR2 S1S2: Val683 and the segment around Lys716-Cys718. The differences in binding affinity of these agonists do not explain the differences in microsecond-to-millisecond conformational dynamics because quisqualate and AMPA bind with similar affinities that are 10-fold greater than the affinity of glutamate. Differences in conformational mobility may reflect differences in the binding mode of AMPA in the GluR2 S1S2 active site compared to the other two ligands. The sites of conformational mobility in GluR2 S1S2 imply that subtle differences exist between the agonists glutamate, quisqualate, and AMPA in modulating glutamate receptor function.     

Characterization of the {alpha}-helix region in domain 3 of the haemolytic lectin CEL-III: implications for self-oligomerization and haemolytic processes

CEL-III is a haemolytic lectin, which has two beta-trefoil domains (domains 1 and 2) and a beta-sheet-rich domain (domain 3). In domain 3 (residues 284-432), there is a hydrophobic region containing two alpha-helices (H8 and H9, residues 317-357) and a loop between them, in which alternate hydrophobic residues, especially Val residues, are present. To elucidate the role of the alpha-helix region in the haemolytic process, peptides corresponding to different parts of this region were synthesized and characterized. The peptides containing the sequence that corresponded to the loop and second alpha-helix (H9) showed the strongest antibacterial activity for Staphylococcus aureus and Bacillus subtilis through a marked permeabilization of the bacterial cell membrane. The recombinant glutathione S-transferase (GST)-fusion proteins containing domain 3 or the alpha-helix region peptide formed self-oligomers, whereas mutations in the alternate Val residues in the alpha-helix region lead to decreased oligomerization ability of the fusion proteins. These results suggest that the alpha-helix region, particularly its alternate Val residues are important for oligomerization of CEL-III in target cell membranes, which is also required for a subsequent haemolytic action.