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.     

Orientational dynamics of indane dione spin-labeled myosin heads in relaxed and contracting skeletal muscle fibers.

We have used electron paramagnetic resonance (EPR) spectroscopy to study the orientation and rotational motions of spin-labeled myosin heads during steady-state relaxation and contraction of skinned rabbit psoas muscle fibers. Using an indane-dione spin label, we obtained EPR spectra corresponding specifically to probes attached to Cys 707 (SH1) on the catalytic domain of myosin heads. The probe is rigidly immobilized, so that it reports the global rotation of the myosin head, and the probe's principal axis is aligned almost parallel with the fiber axis in rigor, making it directly sensitive to axial rotation of the head. Numerical simulations of EPR spectra showed that the labeled heads are highly oriented in rigor, but in relaxation they have at least 90 degrees (Gaussian full width) of axial disorder, centered at an angle approximately equal to that in rigor. Spectra obtained in isometric contraction are fit quite well by assuming that 79 +/- 2% of the myosin heads are disordered as in relaxation, whereas the remaining 21 +/- 2% have the same orientation as in rigor. Computer-simulated spectra confirm that there is no significant population (> 5%) of heads having a distinct orientation substantially different (> 10 degrees) from that in rigor, and even the large disordered population of heads has a mean orientation that is similar to that in rigor. Because this spin label reports axial head rotations directly, these results suggest strongly that the catalytic domain of myosin does not undergo a transition between two distinct axial orientations during force generation. Saturation transfer EPR shows that the rotational disorder is dynamic on the microsecond time scale in both relaxation and contraction. These results are consistent with models of contraction involving 1) a transition from a dynamically disordered preforce state to an ordered (rigorlike) force-generating state and/or 2) domain movements within the myosin head that do not change the axial orientation of the SH1-containing catalytic domain relative to actin.     

Exploration of the transition state for tertiary structure formation between an RNA helix and a large structured RNA

Docking of the P1 duplex into the pre-folded core of the Tetrahymena group I ribozyme exemplifies the formation of tertiary interactions in the context of a complex, structured RNA. We have applied Phi-analysis to P1 docking, which compares the effects of modifications on the rate constant for docking (k(dock)) with the effects on the docking equilibrium (K(dock)). To accomplish this we used a single molecule fluorescence resonance energy transfer assay that allows direct determination of the rate constants for formation of thermodynamically favorable, as well as unfavorable, states. Modification of the eight groups of the P1 duplex that make tertiary interactions with the core and changes in solution conditions decrease K(dock) up to 500-fold, whereas k(dock) changes by 相似文献   

Phase transition in force during ramp stretches of skeletal muscle.

Active glycerinated rabbit psoas fibers were stretched at constant velocity (0.1-3.0 lengths/s) under sarcomere length control. As observed by previous investigators, force rose in two phases: an initial rapid increase over a small stretch (phase I), and a slower, more modest rise over the remainder of the stretch (phase II). The transition between the two phases occurred at a critical stretch (LC) of 7.7 +/- 0.1 nm/half-sarcomere that is independent of velocity. The force at critical stretch (PC) increased with velocity up to 1 length/s, then was constant at 3.26 +/- 0.06 times isometric force. The decay of the force response to a small step stretch was much faster during stretch than in isometric fibers. The addition of 3 mM vanadate reduced isometric tension to 0.08 +/- 0.01 times control isometric tension (P0), but only reduced PC to 0.82 +/- 0.06 times P0, demonstrating that prepowerstroke states contribute to force rise during stretch. The data can be explained by a model in which actin-attached cross-bridges in a prepowerstroke state are stretched into regions of high force and detach very rapidly when stretched beyond this region. The prepowerstroke state acts as a mechanical rectifier, producing large forces during stretch but small forces during shortening.     

Molecular modeling of the GABAC receptor ligand-binding domain

We have constructed a molecular model of the ligand-binding domain of the GABA(C) receptor, which is a member of the Cys-loop ligand-gated ion channel family. The extracellular domains of these receptors share similar sequence homology (20%) with Limnaea acetylcholine-binding protein for which an X-ray crystal structure is available. We used this structure as a template for homology modeling of the GABA(C) receptor extracellular domain using FUGUE and MODELLER software. FlexX was then used to dock GABA into the receptor ligand-binding site, resulting in three alternative energetically favorable orientations. Residues located no more than 5 A from the docked GABA were identified for each model; of these, three were found to be common to all models with 14 others present only in certain models. Using data from experimental studies, we propose that the most likely orientation of GABA is with its amine close to Y198, and its carboxylate close to R104. These studies have therefore provided a model of the ligand-binding domain, which will be useful for both GABA(C) and GABA(A) receptor studies, and have also yielded an experimentally testable hypothesis of the location of GABA in the binding pocket. [Figure: see text].     

Structural changes in the transition state of protein folding: alternative interpretations of curved chevron plots.

The interpretation of folding rates is often rationalized within the context of transition state theory. This means that the reaction rate is linked to an activation barrier, the height of which is determined by the free energy difference between a ground state (the starting point) and an apparent transition state. Changes in the folding kinetics are thus caused by effects on either the ground state, the transition state, or both. However, structural changes of the transition state are rarely discussed in connection with experimental data, and kinetic anomalies are commonly ascribed to ground state effects alone, e.g., depletion or accumulation of structural intermediates upon addition of denaturant. In this study, we present kinetic data which are best described by transition state changes. We also show that ground state effects and transition state effects are in general difficult to distinguish kinetically. The analysis is based on the structurally homologous proteins U1A and S6. Both proteins display two-state behavior, but there is a marked difference in their kinetics. S6 exhibits a classical V-shaped chevron plot (log observed rate constant vs denaturant concentration), whereas U1A's chevron plot is symmetrically curved, like an inverted bell curve. However, S6 is readily mutated to display U1A-like kinetics. The seemingly drastic effects of these mutations are readily ascribed to transition state movements where large kinetic differences result from relatively small alterations of a common free energy profile and broad activation barriers.     

A computational comparison of the atomic models of the actomyosin interface

Several atomic models of the actomyosin interface have been proposed based on the docking together of their component structures using electron microscopy and resonance energy-transfer measurements. Although these models are in approximate agreement in the location of the binding interfaces when myosin is tightly bound to actin, their relationships to molecular docking simulations based on computational free-energy calculations are investigated here. Both rigid-docking and flexible-docking conformational search strategies were used to identify free-energy minima at the interfaces between atomic models of myosin and actin. These results suggest that the docking model produced by resonance energy-transfer data is closer to a free-energy minimum at the interface than are the available atomic models based on electron microscopy. The conformational searches were performed using both scallop and chicken skeletal muscle myosins and identified similarly oriented actin-binding interfaces that serve to validate that these models are at the global minimum. These results indicate that the existing docking models are close to but not precisely at the lowest-energy initial contact site for strong binding between myosin and actin that should represent an initial contact between the two proteins; therefore, conformational changes are likely to be important during the transition to a strongly bound complex.     

Kissing G Domains of MnmE Monitored by X-Ray Crystallography and Pulse Electron Paramagnetic Resonance Spectroscopy

MnmE, which is involved in the modification of the wobble position of certain tRNAs, belongs to the expanding class of G proteins activated by nucleotide-dependent dimerization (GADs). Previous models suggested the protein to be a multidomain protein whose G domains contact each other in a nucleotide dependent manner. Here we employ a combined approach of X-ray crystallography and pulse electron paramagnetic resonance (EPR) spectroscopy to show that large domain movements are coupled to the G protein cycle of MnmE. The X-ray structures show MnmE to be a constitutive homodimer where the highly mobile G domains face each other in various orientations but are not in close contact as suggested by the GDP-AlF_x structure of the isolated domains. Distance measurements by pulse double electron-electron resonance (DEER) spectroscopy show that the G domains adopt an open conformation in the nucleotide free/GDP-bound and an open/closed two-state equilibrium in the GTP-bound state, with maximal distance variations of 18 Å. With GDP and AlF_x, which mimic the transition state of the phosphoryl transfer reaction, only the closed conformation is observed. Dimerization of the active sites with GDP-AlF_x requires the presence of specific monovalent cations, thus reflecting the requirements for the GTPase reaction of MnmE. Our results directly demonstrate the nature of the conformational changes MnmE was previously suggested to undergo during its GTPase cycle. They show the nucleotide-dependent dynamic movements of the G domains around two swivel positions relative to the rest of the protein, and they are of crucial importance for understanding the mechanistic principles of this GAD.     

Classification and stability of global inhomogeneous solutions of a macroscopic model of cell motion

Many micro-organisms use chemotaxis for aggregation, resulting in stable patterns. In this paper, the amoeba Dictyostelium discoideum serves as a model organism for understanding the conditions for aggregation and classification of resulting patterns. To accomplish this, a 1D nonlinear diffusion equation with chemotaxis that models amoeba behavior is analyzed. A classification of the steady state solutions is presented, and a Lyapunov functional is used to determine conditions for stability of inhomogenous solutions. Changing the chemical sensitivity, production rate of the chemical attractant, or domain length can cause the system to transition from having an asymptotic steady state, to having asymptotically stable single-step solution and multi-stepped stable plateau solutions.     

Conformational selection during weak binding at the actin and myosin interface

The molecular mechanism of the powerstroke in muscle is examined by resonance energy transfer techniques. Recent models suggesting a pre-cocking of the myosin head involving an enormous rotation between the lever arm and the catalytic domain were tested by measuring separation distances among myosin subfragment-2, the nucleotide site, and the regulatory light chain in the presence of nucleotide transition state analogs. Only small changes (<0.5 nm) were detected that are consistent with internal conformational changes of the myosin molecule, but not with extreme differences in the average lever arm position suggested by some atomic models. These results were confirmed by stopped-flow resonance energy transfer measurements during single ATP turnovers on myosin. To examine the participation of actin in the powerstroke process, resonance energy transfer between the regulatory light chain on myosin subfragment-1 and the C-terminus of actin was measured in the presence of nucleotide transition state analogs. The efficiency of energy transfer was much greater in the presence of ADP-AlF(4), ADP-BeF(x), and ADP-vanadate than in the presence of ADP or no nucleotide. These data detect profound differences in the conformations of the weakly and strongly attached cross-bridges that appear to result from a conformational selection that occurs during the weak binding of the myosin head to actin.     

Oligomerization state and supramolecular structure of the HIV-1 Vpu protein transmembrane segment in phospholipid bilayers

HIV‐1 Vpu is an 81‐residue protein with a single N‐terminal transmembrane (TM) helical segment that is involved in the release of new virions from host cell membranes. Vpu and its TM segment form ion channels in phospholipid bilayers, presumably by oligomerization of TM helices into a pore‐like structure. We describe measurements that provide new constraints on the oligomerization state and supramolecular structure of residues 1–40 of Vpu (Vpu_1–40), including analytical ultracentrifugation measurements to investigate oligomerization in detergent micelles, photo‐induced crosslinking experiments to investigate oligomerization in bilayers, and solid‐state nuclear magnetic resonance measurements to obtain constraints on intermolecular contacts between and orientations of TM helices in bilayers. From these data, we develop molecular models for Vpu TM oligomers. The data indicate that a variety of oligomers coexist in phospholipid bilayers, so that a unique supramolecular structure can not be defined. Nonetheless, since oligomers of various sizes have similar intermolecular contacts and orientations, molecular models developed from our data are most likely representative of Vpu TM oligomers that exist in host cell membranes.     

The crystal structure of a Cys25 -> Ala mutant of human procathepsin S elucidates enzyme-prosequence interactions

The crystal structure of the active-site mutant Cys25 --> Ala of glycosylated human procathepsin S is reported. It was determined by molecular replacement and refined to 2.1 Angstrom resolution, with an R-factor of 0.198. The overall structure is very similar to other cathepsin L-like zymogens of the C1A clan. The peptidase unit comprises two globular domains, and a small third domain is formed by the N-terminal part of the prosequence. It is anchored to the prosegment binding loop of the enzyme. Prosegment residues beyond the prodomain dock to the substrate binding cleft in a nonproductive orientation. Structural comparison with published data for mature cathepsin S revealed that procathepsin S residues Phe146, Phe70, and Phe211 adopt different orientations. Being part of the S1' and S2 pockets, they may contribute to the selectivity of ligand binding. Regarding the prosequence, length, orientation and anchoring of helix alpha3p differ from related zymogens, thereby possibly contributing to the specificity of propeptide-enzyme interaction in the papain family. The discussion focuses on the functional importance of the most conserved residues in the prosequence for structural integrity, inhibition and folding assistance, considering scanning mutagenesis data published for procathepsin S and for its isolated propeptide.     

Certain metal ions are inhibitors of cytochrome b6f complex 'Rieske' iron-sulfur protein domain movements

Many current models of the Q cycle for the cytochrome (cyt) b6f and the cyt bc1 complexes incorporate 'Rieske' iron-sulfur protein (ISP) domain movements to gate electron transfer and to ensure high yields of proton shuttling. It was previously proposed that copper ions, which bind at a site distant from the quinol oxidase (Q(o)) site, inhibit plastoquinol (PQH2) binding by restraining the hydrophilic head domain of the ISP [Rao B. K., S., Tyryshkin, A. M., Roberts, A. G., Bowman, M. K., and Kramer, D. M. (1999) Biochemistry 38, 3285-3296]. The present work presents evidence that this is indeed the case for both copper ions and Zn2+, which appear to inhibit by similar mechanisms. Electron paramagnetic resonance (EPR) spectra show that Cu2+ and Zn2+ binding to the cyt b6f complex displaces the Q(o) site inhibitor 2,5-dibromo-3-methyl-6-isopropylbenzoquinone (DBMIB). At high concentrations, both DBMIB and Cu2+ or Zn2+ can bind simultaneously, altering the Rieske 2Fe2S cluster and Cu2+ EPR spectra, suggesting perturbations in their respective binding sites. Both Zn2+ and Cu1+ altered the orientations of the Rieske 2Fe2S cluster with respect to the membrane plane, but had no effect on that of the cyt b6 hemes. Cu2+ was found to change the orientation of the cyt f heme plane, consistent with binding on the cyt f protein. Within conservative constraints, the data suggest that the ISP is shifted into a position intermediate between the ISP(C) position, when the Q(o) site is unoccupied, and the ISP(B) position, when the Q(o) site is occupied by inhibitors such as DBMIB or stigmatellin. These results support the role of ISP domain movements in Q(o) site catalysis.     

Characterization of the hydrodynamic properties of the folding transition state of an SH3 domain by magnetization transfer NMR spectroscopy

Protein folding kinetic data have been obtained for the marginally stable N-terminal Src homology 3 domain of the Drosophila protein drk (drkN SH3) in an investigation of the hydrodynamic properties of its folding transition state. Due to the presence of NMR resonances of both folded and unfolded states at equilibrium, kinetic data can be derived from NMR magnetization transfer techniques under equilibrium conditions. Kinetic analysis as a function of urea (less than approximately 1 M) and glycerol enables determination of alpha values, measures of the energetic sensitivity of the transition state to the perturbation relative to the end states of the protein folding reaction (the folded and unfolded states). Both end states have previously been studied experimentally by NMR spectroscopic and other biophysical methods in great detail and under nondenaturing conditions. Combining these results with the kinetic folding data obtained here, we can characterize the folding transition state without requiring empirical models for the unfolded state structure. We are thus able to give a reliable measure of the solvent-accessible surface area of the transition state of the drkN SH3 domain (4730 +/- 360 A(2)) based on urea titration data. Glycerol titration data give similar results and additionally demonstrate that folding of this SH3 domain is dependent on solvent viscosity, which is indicative of at least partial hydration of the transition state. Because SH3 domains appear to fold by a common folding mechanism, the data presented here provide valuable insight into the transition states of the drkN and other SH3 domains.     

On the dynamic nature of the transition state for protein-protein association as determined by double-mutant cycle analysis and simulation

The process of protein-protein association starts with their random collision, which may develop into an encounter complex followed by a transition state and final complex formation. Here we aim to experimentally characterize the nature of the transition state of protein-protein association for three different protein-protein interactions; Barnase-Barstar, TEM1-BLIP and IFNalpha2-IFNAR2, and use the data to model the transition state structures. To model the transition state, we determined inter-protein distance-constraints of the activated complex by using double mutant cycles (DMC) assuming that interacting residues are spatially close. Significant DeltaDeltaG(double dagger)(int) values were obtained only between residues on Barnase and Barstar. However, introducing specific mutations that optimize the charge complementarity between BLIP and TEM1 resulted in the introduction of significant DeltaDeltaG(double dagger)(int) values also between residues of these two proteins. While electrostatic interactions make major contributions towards stabilizing the transition state, we show two examples where steric hindrance exerts an effect on the transition state as well. To model the transition-state structures from the experimental DeltaDeltaG(double dagger)(int) values, we introduced a method for structure perturbation, searching for those inter-protein orientations that best support the experimental DeltaDeltaG(double dagger)(int) values. Two types of transition states were found, specific as observed for Barnase-Barstar and the electrostatically optimized TEM1-BLIP mutants, and diffusive as shown for wild-type TEM1-BLIP and IFNalpha2-IFNAR2. The specific transition states are characterized by defined inter-protein orientations, which cannot be modeled for the diffusive transition states. Mutations introduced through rational design can change the transition state from diffusive to specific. Together, these data provide a structural view of the mechanism allowing rates of association to differ by five orders of magnitude between different protein complexes.     

Structure of the base of the L7/L12 stalk of the Haloarcula marismortui large ribosomal subunit: analysis of L11 movements

Initiation factors, elongation factors, and release factors all interact with the L7/L12 stalk of the large ribosomal subunit during their respective GTP-dependent cycles on the ribosome. Electron density corresponding to the stalk is not present in previous crystal structures of either 50 S subunits or 70 S ribosomes. We have now discovered conditions that result in a more ordered factor-binding center in the Haloarcula marismortui (H.ma) large ribosomal subunit crystals and consequently allows the visualization of the full-length L11, the N-terminal domain (NTD) of L10 and helices 43 and 44 of 23 S rRNA. The resulting model is currently the most complete reported structure of a L7/L12 stalk in the context of a ribosome. This region contains a series of intermolecular interfaces that are smaller than those typically seen in other ribonucleoprotein interactions within the 50 S subunit. Comparisons of the L11 NTD position between the current structure, which is has an NTD splayed out with respect to previous structures, and other structures of ribosomes in different functional states demonstrates a dynamic range of L11 NTD movements. We propose that the L11 NTD moves through three different relative positions during the translational cycle: apo-ribosome, factor-bound pre-GTP hydrolysis and post-GTP hydrolysis. These positions outline a pathway for L11 NTD movements that are dependent on the specific nucleotide state of the bound ligand. These three states are represented by the orientations of the L11 NTD relative to the ribosome and suggest that L11 may play a more specialized role in the factor binding cycle than previously appreciated.     

A metastable form of the large envelope protein of duck hepatitis B virus: low-pH release results in a transition to a hydrophobic, potentially fusogenic conformation

We have examined the structure and fusion potential of the duck hepatitis B virus (DHBV) envelope proteins by treating subviral particles with deforming agents known to release envelope proteins of viruses from a metastable to a fusion-active state. Exposure of DHBV particles to low pH triggered a major structural change in the large envelope protein (L), resulting in exposure of trypsin sites within its S domain but without affecting the same region in the small surface protein (S) subunits. This conformational change was associated with increased hydrophobicity of the particle surface, most likely arising from surface exposure of the hydrophobic first transmembrane domain (TM1). In the hydrophobic conformation, DHBV particles were able to bind to liposomes and intact cells, while in their absence these particles aggregated, resulting in viral inactivation. These results suggests that some L molecules are in a spring-loaded metastable state which, when released, exposes a previously hidden hydrophobic domain, a transition potentially representing the fusion-active state of the envelope.     

Transition state analysis of the complex between coagulation factor VIIa and tissue factor: suggesting a sequential domain-binding pathway

Injury of a blood vessel exposes membrane-bound tissue factor (TF) to blood, which allows binding of coagulation factor VIIa (FVIIa). This initiation of the coagulation cascade is dictated by a specific multi-domain interaction between FVIIa and TF. To examine the energies involved in the transition state of the FVIIa:TF complex, various residues in the extracellular part of TF (sTF) that are known to interact with FVIIa were replaced with a smaller cysteine residue. Determination of Phi values in each of the positions using surface plasmon resonance measurements enabled us to characterize the transition state complex between the resulting sTF variants and FVIIa. We found that the interactions in the transition state seemed to be most pronounced between the protease domain of FVIIa and sTF while detailed specific interactions between the Gla-domain and sTF were missing. Thus, the transition state energy data indicate a sequential binding event between these two macromolecules.     

Monomeric state of S100P protein: Experimental and molecular dynamics study

S100 proteins constitute a large subfamily of the EF-hand superfamily of calcium binding proteins. They possess one classical EF-hand Ca²⁺-binding domain and an atypical EF-hand domain. Most of the S100 proteins form stable symmetric homodimers. An analysis of literature data on S100 proteins showed that their physiological concentrations could be much lower than dissociation constants of their dimeric forms. It means that just monomeric forms of these proteins are important for their functioning. In the present work, thermal denaturation of apo-S100P protein monitored by intrinsic tyrosine fluorescence has been studied at various protein concentrations within the region from 0.04–10 μM. A transition from the dimeric to monomeric form results in a decrease in protein thermal stability shifting the mid-transition temperature from 85 to 75 °C. Monomeric S100P immobilized on the surface of a sensor chip of a surface plasmon resonance instrument forms calcium dependent 1 to 1 complexes with human interleukin-11 (equilibrium dissociation constant 1.2 nM). In contrast, immobilized interleukin-11 binds two molecules of dimeric S100P with dissociation constants of 32 nM and 288 nM. Since effective dissociation constant of dimeric S100P protein is very low (0.5 μM as evaluated from our data) the sensitivity of the existing physical methods does not allow carrying out a detailed study of S100P monomer properties. For this reason, we have used molecular dynamics methods to evaluate structural changes in S100P upon its transition from the dimeric to monomeric state. 80-ns molecular dynamics simulations of kinetics of formation of S100P, S100B and S100A11 monomers from the corresponding dimers have been carried out. It was found that during the transition from the homo-dimer to monomer form, the three S100 monomer structures undergo the following changes: (1) the helices in the four-helix bundles within each monomer rotate in order to shield the exposed non-polar residues; (2) almost all lost contacts at the dimer interface are substituted with equivalent and newly formed interactions inside each monomer, and new stabilizing interactions are formed; and (3) all monomers recreate functional hydrophobic cores. The results of the present study show that both dimeric and monomeric forms of S100 proteins can be functional.     

Model of glycoprotein hormone receptor ligand binding and signaling

Studies described here were initiated to develop a model of glycoprotein hormone receptor structure and function. We found that the region that links the lutropin receptor leucine-rich repeat domain (LRD) to its transmembrane domain (TMD) has substantial roles in ligand binding and signaling, hence we term it the signaling specificity domain (SSD). Theoretical considerations indicated the short SSDs in marmoset lutropin and salmon follitropin receptors have KH domain folds. We assembled models of lutropin, follitropin, and thyrotropin receptors by aligning models of their LRD, TMD, and shortened SSD in a manner that explains how substitutions in follitropin and thyrotropin receptors distant from their apparent ligand binding sites enable them to recognize lutropins. In these models, the SSD is parallel to the concave surface of the LRD and makes extensive contacts with TMD outer loops 1 and 2. The LRD appears to contact TMD outer loop 3 and a few residues in helices 1, 5, 6, and 7. We propose that signaling results from contacts of the ligands with the SSD and LRD that alter the LRD, which then moves TMD helices 6 and 7. The positions of the LRD and SSD support the notion that the receptor can be activated by hormones that dock with these domains in either of two different orientations. This would account for the abilities of some ligands and ligand chimeras to bind multiple receptors and for some receptors to bind multiple ligands. This property of the receptor may have contributed significantly to ligand-receptor co-evolution.