Probing the local dynamics of nucleotide-binding pocket coupled to the global dynamics: myosin versus kinesin

Based on the elastic network model, we develop a new analysis for protein complexes, which probes the local dynamics of a subsystem that is elastically coupled to a fluctuating environment. This method is applied to a comparative dynamical analysis of the nucleotide-binding pocket of two motor proteins-myosins and kinesins. In myosins, the observed structural changes in the nucleotide-pocket from the transition state to the rigorlike state are dominated by the lowest normal mode that involves significant movements in both switch I and switch II; in kinesins, the measured conformational changes in the nucleotide-pocket are also dominated by the lowest mode, which, however, only involves large movement in switch I. We then compute the global structural changes induced by the nucleotide-pocket deformations as described by the dominant pocket-mode, which yield encouraging results: in myosins, multiple hinge motions involving the opening/closing of the cleft between the upper and lower 50 -kDa subdomains and the swinging movement of the converter are induced, which are dominated by precisely the same global mode that has been recently identified by us as important to the dynamical correlations among the nucleotide-pocket, the actin-binding site, and the converter; in kinesins, the induced global conformational changes are well described by a highly collective global mode which hints for a dynamical pathway spanning from the nucleotide-pocket to the neck-linker via the H6 helix.     

Toward the mechanism of dynamical couplings and translocation in hepatitis C virus NS3 helicase using elastic network model

Hepatitis C virus NS3 helicase is an enzyme that unwinds double-stranded polynucleotides in an ATP-dependent reaction. It provides a promising target for small molecule therapeutic agents against hepatitis C. Design of such drugs requires a thorough understanding of the dynamical nature of the mechanochemical functioning of the helicase. Despite recent progress, the detailed mechanism of the coupling between ATPase activity and helicase activity remains unclear. Based on an elastic network model (ENM), we apply two computational analysis tools to probe the dynamical mechanism underlying the allosteric coupling between ATP binding and polynucleotide binding in this enzyme. The correlation analysis identifies a network of hot-spot residues that dynamically couple the ATP-binding site and the polynucleotide-binding site. Several of these key residues have been found by mutational experiments as functionally important, while our analysis also reveals previously unexplored hot-spot residues that are potential targets for future mutational studies. The conformational changes between different crystal structures of NS3 helicase are found to be dominated by the lowest frequency mode solved from the ENM. This mode corresponds to a hinge motion of the highly flexible domain 2. This motion simultaneously modulates the opening/closing of the domains 1-2 cleft where ATP binds, and the domains 2-3 cleft where the polynucleotide binds. Additionally, a small twisting motion of domain 1, observed in both mode 1 and the computed ATP binding induced conformational change, fine-tunes the binding affinity of the domains 1-3 interface for the polynucleotide. The combination of these motions facilitates the translocation of a single-stranded polynucleotide in an inchworm-like manner.     

Rapid activation of transducin by mutations distant from the nucleotide-binding site: evidence for a mechanistic model of receptor-catalyzed nucleotide exchange by G proteins

G proteins act as molecular switches in which information flow depends on whether the bound nucleotide is GDP ("off") or GTP ("on"). We studied the basal and receptor-catalyzed nucleotide exchange rates of site-directed mutants of the alpha subunit of transducin. We identified three amino acid residues (Thr-325, Val-328, and Phe-332) in which mutation resulted in dramatic increases (up to 165-fold) in basal nucleotide exchange rates in addition to enhanced receptor-catalyzed nucleotide exchange rates. These three residues are located on the inward facing surface of the alpha5 helix, which lies between the carboxyl-terminal tail and a loop contacting the nucleotide-binding pocket. Mutation of amino acid residues on the outward facing surface of the same alpha5 helix caused a decrease in receptor-catalyzed nucleotide exchange. We propose that the alpha5 helix comprises a functional microdomain in G proteins that affects basal nucleotide release rates and mediates receptor-catalyzed nucleotide exchange at a distance from the nucleotide-binding pocket.     

Subunit-selective contribution to channel gating of the M4 domain of the nicotinic receptor

The muscle nicotinic receptor (AChR) is a pentamer of four different subunits, each of which contains four transmembrane domains (M1-M4). We recently showed that channel opening and closing rates of the AChR depend on a hydrogen bond involving a threonine at position 14' of the M4 domain in the alpha-subunit. To determine whether residues in equivalent positions in non-alpha-subunits contribute to channel gating, we mutated deltaT14', betaT14', and epsilonS14' and evaluated changes in the kinetics of acetylcholine-activated currents. The mutation epsilonS14'A profoundly slows the rate of channel closing, an effect opposite to that produced by mutation of alphaT14'. Unlike mutations of alphaT14', epsilonS14'A does not affect the rate of channel opening. Mutations in deltaT14' and betaT14' do not affect channel opening or closing kinetics, showing that conserved residues are not functionally equivalent in all subunits. Whereas alphaT14'A and epsilonS14'A subunits contribute additively to the closing rate, they contribute nonadditively to the opening rate. Substitution of residues preserving the hydrogen bonding ability at position 14' produce nearly normal gating kinetics. Thus, we identify subunit-specific contributions to channel gating of equivalent residues in M4 and elucidate the underlying mechanistic and structural bases.     

Fluorescence labeling and computational analysis of the strut of myosin's 50kDa cleft

A new fluorescent labeling procedure specific for the strut sequence of myosin subfragment-1's 50kDa cleft was developed using CY3 N-hydroxy succinimidyl ester as a hydrophobic tag and hydrophobic interaction chromatography to purify the major labeled species which retained actin-activated ATPase activity. Stern-Volmer analysis suggests that the CY3 is in close proximity to basic residues, consistent with inspection of the mapped labeling site in the atomic model. Fluorescence polarization indicates that the CY3 becomes more mobile upon actin binding, supporting a location near the actomyosin interface. In contrast, nucleotide binding to myosin had little impact on the CY3. Molecular mechanics and stochastic dynamics simulations suggest that this labeling site is sensitive to forced cleft opening and closure, but the upper 50kDa cleft does not move easily. In addition, there appear to be some long-range effects of forced cleft opening and closing that could impact the lever arm position.     

Mechanisms of isoform-specific residue influence on GTP-bound HRas,KRas, and NRas

HRas, KRas, and NRas are GTPases with a common set of effectors that control many cell-signaling pathways, including proliferation through Raf kinase. Their G-domains are nearly identical in sequence, with a few isoform-specific residues that have an effect on dynamics and biochemical properties. Here, we use accelerated molecular dynamics (aMD) simulations consistent with solution x-ray scattering experiments to elucidate mechanisms through which isoform-specific residues associated with each Ras isoform affects functionally important regions connected to the active site. HRas-specific residues cluster in loop 8 to stabilize the nucleotide-binding pocket, while NRas-specific residues on helix 3 directly affect the conformations of switch I and switch II. KRas, the most globally flexible of the isoforms, shows greatest fluctuations in the switch regions enhanced by a KRas-specific residue in loop 7 and a highly dynamic loop 8 region. The analysis of isoform-specific residue effects on Ras proteins is supported by NMR experiments and is consistent with previously published biochemical data.     

Conserved gating hinge in ligand- and voltage-dependent K+ channels

Ion channels open and close their pore in a process called gating. On the basis of crystal structures of two voltage-independent K(+) channels, KcsA and MthK, a conformational change for gating has been proposed whereby the inner helix bends at a glycine hinge point (gating hinge) to open the pore and straightens to close it. Here we ask if a similar gating hinge conformational change underlies the mechanics of pore opening of two eukaryotic voltage-dependent K(+) channels, Shaker and BK channels. In the Shaker channel, substitution of the gating hinge glycine with alanine and several other amino acids prevents pore opening, but the ability to open is recovered if a secondary glycine is introduced at an adjacent position. A proline at the gating hinge favors the open state of the Shaker channel as if by preventing inner helix straightening. In BK channels, which have two adjacent glycine residues, opening is significantly hindered in a graded manner with single and double mutations to alanine. These results suggest that K(+) channels, whether ligand- or voltage-dependent, open when the inner helix bends at a conserved glycine gating hinge.     

A structural analysis of asymmetry required for catalytic activity of an ABC-ATPase domain dimer

The ATP-binding cassette (ABC)-transporter haemolysin (Hly)B, a central element of a Type I secretion machinery, acts in concert with two additional proteins in Escherichia coli to translocate the toxin HlyA directly from the cytoplasm to the exterior. The basic set of crystal structures necessary to describe the catalytic cycle of the isolated HlyB-NBD (nucleotide-binding domain) has now been completed. This allowed a detailed analysis with respect to hinge regions, functionally important key residues and potential energy storage devices that revealed many novel features. These include a structural asymmetry within the ATP dimer that was significantly enhanced in the presence of Mg2+, indicating a possible functional asymmetry in the form of one open and one closed phosphate exit tunnel. Guided by the structural analysis, we identified two amino acids, closing one tunnel by an apparent salt bridge. Mutation of these residues abolished ATP-dependent cooperativity of the NBDs. The implications of these new findings for the coupling of ATP binding and hydrolysis to functional activity are discussed.     

Coupling between Normal Modes Drives Protein Conformational Dynamics: Illustrations Using Allosteric Transitions in Myosin II

Structure-based elastic network models (ENMs) have been remarkably successful in describing conformational transitions in a variety of biological systems. Low-frequency normal modes are usually calculated from the ENM that characterizes elastic interactions between residues in contact in a given protein structure with a uniform force constant. To explore the dynamical effects of nonuniform elastic interactions, we calculate the robustness and coupling of the low-frequency modes in the presence of nonuniform variations in the ENM force constant. The variations in the elastic interactions, approximated here by Gaussian noise, approximately account for perturbation effects of heterogeneous residue-residue interactions or evolutionary sequence changes within a protein family. First-order perturbation theory provides an efficient and qualitatively correct estimate of the mode robustness and mode coupling for finite perturbations to the ENM force constant. The mode coupling analysis and the mode robustness analysis identify groups of strongly coupled modes that encode for protein functional motions. We illustrate the new concepts using myosin II motor protein as an example. The biological implications of mode coupling in tuning the allosteric couplings among the actin-binding site, the nucleotide-binding site, and the force-generating converter and lever arm in myosin isoforms are discussed. We evaluate the robustness of the correlation functions that quantify the allosteric couplings among these three key structural motifs.     

Segment TM7 from the cytoplasmic hemi-channel from VO-H+-V-ATPase includes a flexible region that has a potential role in proton translocation

A 900-MHz NMR study is reported of peptide sMTM7 that mimics the cytoplasmic proton hemi-channel domain of the seventh transmembrane segment (TM7) from subunit a of H(+)-V-ATPase from Saccharomyces cerevisiae. The peptide encompasses the amino acid residues known to actively participate in proton translocation. In addition, peptide sMTM7 contains the amino acid residues that upon mutation cause V-ATPase to become resistant against the inhibitor bafilomycin. 2D TOCSY and NOESY (1)H-(1)H NMR spectra are obtained of sMTM7 dissolved in d(6)-DMSO and are used to calculate the three-dimensional structure of the peptide. The NMR-based structures and corresponding dynamical features of peptide sMTM7 show that sMTM7 is composed of two alpha-helical regions. These regions are separated by a flexible hinge of two residues. The hinge acts as a ball-and-joint socket and both helical segments move independently with respect to one another. This movement in TM7 is suggested to cause the opening and closing of the cytoplasmic proton hemi-channel and enables proton translocation.     

A concerted mechanism for opening the GDP binding pocket and release of the nucleotide in hetero-trimeric G-proteins

G-protein hetero-trimers play a fundamental role in cell function. Their dynamic behavior at the atomic level remains to be understood. We have studied the Gi hetero-trimer through a combination of molecular dynamics simulations and normal mode analyses. We showed that these big proteins could undergo large-amplitude conformational changes, without any energy penalty and with an intrinsic dynamics centered on their GDP binding pocket. Among the computed collective motions, one of the modes (mode 17) was particularly able to significantly open both the base and the phosphate sides of the GDP binding pocket. This mode describing mainly a motion between the Ras-like and the helical domains of G_α was in close agreement with some available X-ray data and with many other biochemical/biophysical observations including the kink of helix α5. The use of a new protocol, which allows extraction of the GDP ligand along the computed normal modes, supported that the exit of GDP was largely coupled to an opening motion along mode 17. We propose for the first time a “concerted mechanism” model in which the opening of the GDP pocket and the kink of the α5 helix occur concomitantly and favor GDP release from G_αβγ complexes. This model is discussed in the context of the G-protein-coupled receptor/G-protein interaction close to the cell membrane.     

How the headpiece hinge angle is opened: New insights into the dynamics of integrin activation

How the integrin head transitions to the high-affinity conformation is debated. Although experiments link activation with the opening of the hinge angle between the betaA and hybrid domains in the ligand-binding headpiece, this hinge is closed in the liganded alpha(v)beta3 integrin crystal structure. We replaced the RGD peptide ligand of this structure with the 10th type III fibronectin module (FnIII10) and discovered through molecular dynamics (MD) equilibrations that when the conformational constraints of the leg domains are lifted, the betaA/hybrid hinge opens spontaneously. Together with additional equilibrations on the same nanosecond timescale in which small structural variations impeded hinge-angle opening, these simulations allowed us to identify the allosteric pathway along which ligand-induced strain propagates via elastic distortions of the alpha1 helix to the betaA/hybrid domain hinge. Finally, we show with steered MD how force accelerates hinge-angle opening along the same allosteric pathway. Together with available experimental data, these predictions provide a novel framework for understanding integrin activation.     

Crystal structure and mutational analysis of the Saccharomyces cerevisiae cell cycle regulatory protein Cks1: implications for domain swapping, anion binding and protein interactions

BACKGROUND: The Saccharomyces cerevisiae protein Cks1 (cyclin-dependent kinase subunit 1) is essential for cell-cycle progression. The biological function of Cks1 can be modulated by a switch between two distinct molecular assemblies: the single domain fold, which results from the closing of a beta-hinge motif, and the intersubunit beta-strand interchanged dimer, which arises from the opening of the beta-hinge motif. The crystal structure of a cyclin-dependent kinase (Cdk) in complex with the human Cks homolog CksHs1 single-domain fold revealed the importance of conserved hydrophobic residues and charged residues within the beta-hinge motif. RESULTS: The 3.0 A resolution Cks1 structure reveals the strict structural conservation of the Cks alpha/beta-core fold and the beta-hinge motif. The beta hinge identified in the Cks1 structure includes a novel pivot and exposes a cluster of conserved tyrosine residues that are involved in Cdk binding but are sequestered in the beta-interchanged Cks homolog suc1 dimer structure. This Cks1 structure confirms the conservation of the Cks anion-binding site, which interacts with sidechain residues from the C-terminal alpha helix of another subunit in the crystal. CONCLUSIONS: The Cks1 structure exemplifies the conservation of the beta-interchanged dimer and the anion-binding site in evolutionarily distant yeast and human Cks homologs. Mutational analyses including in vivo rescue of CKS1 disruption support the dual functional roles of the beta-hinge residue Glu94, which participates in Cdk binding, and of the anion-binding pocket that is located 22 A away and on an opposite face to Glu94. The Cks1 structure suggests a biological role for the beta-interchanged dimer and the anion-binding site in targeting Cdks to specific phosphoproteins during cell-cycle progression.     

Dynamic and coordinating domain motions in the active subunits of the F1-ATPase molecular motor

F1-ATPase is a rotary molecular motor crucial for various cellular functions. In F1-ATPase, the rotation of the gammadeltaepsilon subunits against the hexameric alpha(3)beta(3) subunits is highly coordinative, driven by ATP hydrolysis and structural changes at three beta subunits. However, the dynamical and coordinating structural transitions in the beta subunits are not fully understood at the molecular level. Here we examine structural transitions and domain motions in the active subunits of F1-ATPase via dynamical domain analysis of the alpha(3)beta(3)gammadeltaepsilon complex. The domain movement and hinge axes and bending residues have been identified and determined for various conformational changes of the beta-subunits. P-loop and the ATP-binding pocket are for the first time found to play essential mechanical functions additional to the catalytic roles. The cooperative conformational changes pertaining to the rotary mechanism of F1-ATPase appears to be more complex than Boyer's 'bi-site' activity. These findings provide unique molecular insights into dynamic and cooperative domain motions in F1-ATPase.     

Modulating Hinge Flexibility in the APP Transmembrane Domain Alters γ-Secretase Cleavage

Intramembrane cleavage of the β-amyloid precursor protein C99 substrate by γ-secretase is implicated in Alzheimer’s disease pathogenesis. Biophysical data have suggested that the N-terminal part of the C99 transmembrane domain (TMD) is separated from the C-terminal cleavage domain by a di-glycine hinge. Because the flexibility of this hinge might be critical for γ-secretase cleavage, we mutated one of the glycine residues, G38, to a helix-stabilizing leucine and to a helix-distorting proline. Both mutants impaired γ-secretase cleavage and also altered its cleavage specificity. Circular dichroism, NMR, and backbone amide hydrogen/deuterium exchange measurements as well as molecular dynamics simulations showed that the mutations distinctly altered the intrinsic structural and dynamical properties of the substrate TMD. Although helix destabilization and/or unfolding was not observed at the initial ε-cleavage sites of C99, subtle changes in hinge flexibility were identified that substantially affected helix bending and twisting motions in the entire TMD. These resulted in altered orientation of the distal cleavage domain relative to the N-terminal TMD part. Our data suggest that both enhancing and reducing local helix flexibility of the di-glycine hinge may decrease the occurrence of enzyme-substrate complex conformations required for normal catalysis and that hinge mobility can thus be conducive for productive substrate-enzyme interactions.     

Crystal structure of the Apo,unactivated insulin-like growth factor-1 receptor kinase. Implication for inhibitor specificity

The x-ray structure of the unactivated kinase domain of insulin-like growth factor-1 receptor (IGFRK-0P) is reported here at 2.7 A resolution. IGFRK-0P is composed of two lobes connected by a hinge region. The N-terminal lobe of the kinase is a twisted beta-sheet flanked by a single helix, and the C-terminal lobe comprises eight alpha-helices and four short beta-strands. The ATP binding pocket and the catalytic center reside at the interface of the two lobes. Despite the overall similarity to other receptor tyrosine kinases, three notable conformational modifications are observed: 1) this kinase adopts a more closed structure, with its two lobes rotated further toward each other; 2) the conformation of the proximal end of the activation loop (residues 1121-1129) is different; 3) the orientation of the nucleotide-binding loop is altered. Collectively, these alterations lead to a different ATP-binding pocket that might impact on inhibitor designs for IGFRK-0P. Two molecules of IGFRK-0P are seen in the asymmetric unit; they are associated as a dimer with their ATP binding clefts facing each other. The ordered N terminus of one monomer approaches the active site of the other, suggesting that the juxtamembrane region of one molecule could come into close proximity to the active site of the other.     

A gating hinge in Na+ channels; a molecular switch for electrical signaling

Voltage-gated sodium channels are members of a large family with similar pore structures. The mechanism of opening and closing is unknown, but structural studies suggest gating via bending of the inner pore helix at a glycine hinge. Here we provide functional evidence for this gating model for the bacterial sodium channel NaChBac. Mutation of glycine 219 to proline, which would strongly favor bending of the alpha helix, greatly enhances activation by shifting its voltage dependence -51 mV and slowing deactivation by 2000-fold. The mutation also slows voltage-dependent inactivation by 1200-fold. The effects are specific because substitutions of proline at neighboring positions and substitutions of other amino acids at position 219 have much smaller functional effects. Our results fit a model in which concerted bending at glycine 219 in the S6 segments of NaChBac serves as a gating hinge. This gating motion may be conserved in most members of this large ion channel protein family.     

Characterization of the ATP-binding domain of the sarco(endo)plasmic reticulum Ca(2+)-ATPase: probing nucleotide binding by multidimensional NMR.

The skeletal muscle sarco(endo)plasmic reticulum Ca(2+)-ATPase (SERCA1a) mediates muscle relaxation by pumping Ca(2+) from the cytosol to the ER/SR lumen. In efforts aimed at understanding the structural basis for the conformational changes accompanying the reaction cycle catalyzed by SERCA1a, we have studied the ATP-binding domain of SERCA1a in both nucleotide-bound and -free forms by NMR. Limited proteolysis analyses guided us to express a 28 kDa stably folded fragment containing the nucleotide-binding domain of SERCA1a spanning residues Thr357-Leu600. ATP binding activity was demonstrated for this fragment by a FITC competition assay. A nearly complete backbone resonance assignment of this 28 kDa ATP-binding fragment, in both the AMP-PNP-bound and -free forms, was obtained by means of heteronuclear multidimensional NMR techniques. NMR titration experiments with AMP-PNP revealed a confined nucleotide-binding site which coincides with a cytoplasmic pocket region identified in the crystal structure of apo-SERCA1a. These results are consistent with previous site-directed mutagenesis studies of SERCA1a.