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.     

The C2A-C2B linker defines the high affinity Ca(2+) binding mode of rabphilin-3A

The Ca(2+) binding properties of C2 domains are essential for the function of their host proteins. We present here the first crystal structures showing an unexpected Ca(2+) binding mode of the C2B domain of rabphilin-3A in atomic detail. Acidic residues from the linker region between the C2A and C2B domains of rabphilin-3A interact with the Ca(2+)-binding region of the C2B domain. Because of these interactions, the coordination sphere of the two bound Ca(2+) ions is almost complete. Mutation of these acidic residues to alanine resulted in a 10-fold decrease in the intrinsic Ca(2+) binding affinity of the C2B domain. Using NMR spectroscopy, we show that this interaction occurred only in the Ca(2+)-bound state of the C2B domain. In addition, this Ca(2+) binding mode was maintained in the C2 domain tandem fragment. In NMR-based liposome binding assays, the linker was not released upon phospholipid binding. Therefore, this unprecedented Ca(2+) binding mode not only shows how a C2 domain increases its intrinsic Ca(2+) affinity, but also provides the structural base for an atypical protein-Ca(2+)-phospholipid binding mode of rabphilin-3A.     

Conformational heterogeneity and dynamics in a βγ-crystallin from Hahella chejuensis

Most of the βγ-crystallins are structural proteins with high intrinsic stability, which gets enhanced by Ca(2+)-binding in microbial members. Functions of most of these proteins are yet to be known. However, a few of them were reported to be involved in Ca(2+)-dependent and stress-related functions. Hahellin, a microbial homolog, is a natively unfolded protein that acquires a well-folded structure upon Ca(2+) binding. Although the structure of βγ-crystallin domains is well understood, the dynamical features are yet to be explored. We have investigated for the first time the equilibrium dynamics, conformational heterogeneity and associated low-lying free-energy states of hahellin in its Ca(2+)-bound form using NMR spectroscopy to understand the dynamics of a βγ-crystallin domain. Hahellin shows large conformational heterogeneity with nearly 40% of the residues, some of which are part of Ca(2+)-binding loops, accessing alternative states. Further, out of the two Greek key motifs, which together constitute the βγ-crystallin domain, the second Greek key motif is floppy as compared to its relatively rigid counterpart. Taken together, we believe that these characteristics might be of importance to understand the stability and functions of βγ-crystallin domains.     

Low-frequency normal modes in horse liver alcohol dehydrogenase and motions of residues involved in the enzymatic reaction

Normal mode analysis using the elastic network model has provided characteristics and directions of the low-frequency large domain motions of horse liver alcohol dehydrogenase. Three normal modes (mode 1, mode 7, and mode 8) were identified as representative domain motions that may promote the onset of Near Attack Conformers or facilitate the product to be released. The pattern of the atomic displacement for some key residues (such as Val292 and Val203) revealed in this study is in line with experimental structural and kinetic studies and theoretical simulations.     

Conformational basis of the phospholipid requirement for the activity of SR Ca(2+)-ATPase

The delipidated sarcoplasmic reticulum (SR) Ca(2+)-ATPase was reconstituted into proteoliposomes containing different phospholipids. The result demonstrated the necessity of phosphatidylcholine (PC) for optimal ATPase activity and phosphatidylethanolamine (PE) for the optimal calcium transport activity. Fluorescence intensity of Fluorescein 5-isothiocyanate (FITC)-labeled enzyme at Lys515 as well as the measurement of the distance between 5-((2-[(iodoacetyl) amino] ethyl) amino)naphthalene-1-sulphonic acid (IAEDANS) label sites (Cys674/670) and Pr3+ demonstrated a conformational change of cytoplasmic domain, consequently, leading to the variation of the enzyme function with the proteoliposomes composition. Both the intrinsic fluorescence of Trp and its dynamic quenching by HB decreased with increasing PE content, revealing the conformational change of transmembrane domain. Time-resolved fluorescence study characterized three classes of Trp residues, which showed distinctive variation with the change in phospholipid composition. The phospholipid headgroup size caused the conformational change of SR Ca(2+)-ATPase, subsequent the ATPase activity and Ca2+ uptake.     

Transconformations of the SERCA1 Ca-ATPase: a normal mode study

The transport of Ca(2+) by Ca-ATPase across the sarcoplasmic reticulum membrane is accompanied by several transconformations of the protein. Relying on the already established functional importance of low-frequency modes in dynamics of proteins, we report here a normal mode analysis of the Ca(2+)-ATPase based on the crystallographic structures of the E1Ca(2) and E2TG forms. The lowest-frequency modes reveal that the N and A(+Nter) domains undergo the largest amplitude movements. The dynamical domain analysis performed with the DomainFinder program suggests that they behave as rigid bodies, unlike the highly flexible P domain. We highlight two types of movements of the transmembrane helices: i), a concerted movement around an axis perpendicular to the membrane which "twists open" the lumenal side of the protein and ii), an individual translational and rotational mobility which is of lower amplitude for the helices hosting the calcium binding sites. Among all modes calculated for E1Ca, only three are enough to describe the transition to E2TG; the associated movements involve almost exclusively the A and N domains, reflecting the closure of the cytoplasmic headpiece and high displacement of the L7-8 lumenal loop. Subsequently, we discuss the potential contribution of the remaining low-frequency normal modes to the transconformations occurring within the overall calcium transport cycle.     

Diving Into the Lipid Bilayer to Investigate the Transmembrane Organization and Conformational State Transitions of P-type Ion ATPases

Although membrane proteins constitute more than 20% of the total proteins, the structures of only a few are known in detail. An important group of integral membrane proteins are ion-transporting ATPases of the P-type family, which share the formation of an acid-stable phosphorylated intermediate as part of their reaction cycle. There are several crystal structures of the sarcoplasmic reticulum Ca(2+) pump (SERCA) revealing different conformations, and recently, crystal structures of the H(+)-ATPase and the Na(+)/K(+)-ATPase were reported as well. However, there are no atomic resolution structures for other P-type ATPases including the plasma membrane calcium pump (PMCA), which is integral to cellular Ca(2+) signaling. Crystallization of these proteins is challenging because there is often no natural source from which the protein can be obtained in large quantities, and the presence of multiple isoforms in the same tissue further complicates efforts to obtain homogeneous samples suitable for crystallization. Alternative techniques to study structural aspects and conformational transitions in the PMCAs (and other P-type ATPases) have therefore been developed. Specifically, information about the structure and assembly of the transmembrane domain of an integral membrane protein can be obtained from an analysis of the lipid-protein interactions. Here, we review recent efforts using different hydrophobic photo-labeling methods to study the non-covalent interactions between the PMCA and surrounding phospholipids under different experimental conditions, and discuss how the use of these lipid probes can reveal valuable information on the membrane organization and conformational state transitions in the PMCA, Na(+)/K(+)-ATPase, and other P-type ATPases.     

Dynamics of the transition between open and closed conformations in a calmodulin C-terminal domain mutant

BACKGROUND: Calmodulin is a ubiquitous Ca(2+)-activated regulator of cellular processes in eukaryotes. The structures of the Ca(2+)-free (apo) and Ca(2+)-loaded states of calmodulin have revealed that Ca(2+) binding is associated with a transition in each of the two domains from a closed to an open conformation that is central to target recognition. However, little is known about the dynamics of this conformational switch. RESULTS: The dynamics of the transition between closed and open conformations in the Ca(2+)-loaded state of the E140Q mutant of the calmodulin C-terminal domain were characterized under equilibrium conditions. The exchange time constants (tau(ex)) measured for 42 residues range from 13 to 46 micros, with a mean of 21 +/- 3 micros. The results suggest that tau(ex) varies significantly between different groups of residues and that residues with similar values exhibit spatial proximity in the structures of apo and/or Ca(2+)-saturated wild-type calmodulin. Using data for one of these groups, we obtained an open population of p(o) = 0.50 +/- 0.17 and a closed --> open rate constant of k(o) = x 10(4) s(-1). CONCLUSIONS: The conformational exchange dynamics appear to involve locally collective processes that depend on the structural topology. Comparisons with previous results indicate that similar processes occur in the wild-type protein. The measured rates match the estimated Ca(2+) off rate, suggesting that Ca(2+) release may be gated by the conformational dynamics. Structural interpretation of estimated chemical shifts suggests a mechanism for ion release.     

Conformational dynamics of yeast calmodulin in the Ca(2+)-bound state probed using NMR relaxation dispersion

Most calmodulin (CaM) in apo and Ca(2+)-bound states show a dumb-bell-like structure, involving the N- and C-terminal domains, connected with a flexible linker. However, Ca(2+)-bound yeast calmodulin (yCaM) takes on a unique globular structure; the target-binding site of this protein is autoinhibited. We applied NMR relaxation dispersion experiments to yCaM in the Ca(2+)-bound state. The amide (15)N and (1)H(N) relaxation dispersion profiles indicated the presence of conformational dynamics for specific residues at the interface between the N- and C-terminal domains. We conclude that these conformational dynamics were derived from the mobility of the C-terminal domain.     

Locating the thapsigargin-binding site on Ca(2+)-ATPase by cryoelectron microscopy

Thapsigargin (TG) is a potent inhibitor of Ca(2+)-ATPase from sarcoplasmic and endoplasmic reticula. Previous enzymatic studies have concluded that Ca(2+)-ATPase is locked in a dead-end complex upon binding TG with an affinity of <1 nM and that this complex closely resembles the E(2) enzymatic state. We have studied the structural effects of TG binding by cryoelectron microscopy of tubular crystals, which have previously been shown to comprise Ca(2+)-ATPase molecules in the E(2) conformation. In particular, we have compared 3D reconstructions of Ca(2+)-ATPase in the absence and presence of either TG or its dansylated derivative. The overall molecular shape of Ca(2+)-ATPase in the reconstructions is very similar, demonstrating that the TG/Ca(2+)-ATPase complex does indeed physically resemble the E(2) conformation, in contrast to massive domain movements that appear to be induced by Ca(2+) binding. Difference maps reveal a consistent difference on the lumenal side of the membrane, which we conclude corresponds to the thapsigargin-binding site. Modeling the atomic structure for Ca(2+)-ATPase into our density maps reveals that this binding site is composed of the loops between transmembrane segments M3/M4 and M7/M8. Indirect effects are proposed to explain the effects of the S3 stalk segment on thapsigargin affinity as well as thapsigargin-induced changes in ATP affinity. Indeed, a second difference density was observed at the decavanadate-binding site within the three cytoplasmic domains, which we believe reflects an altered affinity as a result of the long-range conformational coupling that drives the reaction cycle of this family of ATP-dependent ion pumps.     

Mechanism of action of the novel plasma membrane Ca(2+)-pump inhibitor caloxin

Caloxin 2A1 is a novel inhibitor of the plasma membrane (PM) Ca(2+)-pump [Am. J. Physiol. Cell Physiol. 280 (2001) C1027]. The PM Ca(2+)-pump is a Ca(2+)-Mg(2+)-ATPase that expels Ca(2+) from cells to help them maintain low concentrations of cytosolic Ca(2+). Caloxin 2A1 inhibits Ca(2+)-Mg(2+)-ATPase in human erythrocyte leaky ghosts. Here we report that this inhibition is non-competitive with respect to the substrates Ca(2+) and ATP and the activator calmodulin. This was anticipated since the high affinity binding site for Ca(2+) and sites for ATP and calmodulin are intracellular whereas caloxin 2A1 is a peptide selected for binding to the second extracellular domain of the pump. Caloxin 2A1 also inhibited the Ca(2+)-dependent formation of the acid stable 140 kDa acylphosphate intermediate from 32P-gamma-ATP. However, it did not inhibit the formation of the acylphosphate intermediate in the reverse direction-from 32P-orthophosphate. Consistent with results on mutagenesis of transmembrane residues in the pump protein, we suggest that caloxin 2A1 inhibits conformational changes required during the reaction cycle of the pump.     

Conformational changes of (Ca2+-Mg2+)-ATPase of erythrocyte plasma membrane caused by calmodulin and phosphatidylserine as revealed by circular dichroism and fluorescence studies

Two spectroscopic techniques, circular dichroism and steady-state fluorescence, were employed in order to study conformational changes of the purified, detergent-solubilized (Ca2+-Mg2+)-ATPase of porcine erythrocyte ghost membranes. Circular dichroism (CD) spectra in the peptide region were obtained from the purified (Ca2+-Mg2+)-ATPase of porcine erythrocyte ghost membranes with the aim to investigate the secondary structure of the enzyme in the presence of calmodulin (CaM) or phosphatidylserine (PS), as well as in the E1 and E2 states. The E1 conformation was stabilized by 10 microM free Ca2+, while the E2 conformation was stabilized by 0.1 mM ethylene glycol bis(2-aminoethyl ether)-N,N,N',N'-tetraacetic acid (EGTA). It was found that the E1 and E2 states of the enzyme strikingly differed in their secondary structure (66% and 46% of calculated alpha-helix content, respectively). In the presence of Ca2+, PS decreased the helical content of the ATPase to 61%, while CaM to 55%. Quenching of intrinsic fluorescence of (Ca2+-Mg2+)-ATPase by acrylamide, performed in the presence of Ca2+, gave evidence for a single class of tryptophan residues with Stern-Volmer constant (KSV) of 10 M-1. Accessibility of tryptophan residues varied depending on the conformational status of the enzyme. Addition of PS and CaM decreased the KSV value to 7.6 M-1 and 8.5 M-1, respectively. In the absence of Ca2+, KSV was 7.0 M-1. KI and CsCl were less effective as quenchers. The fluorescence energy transfer between (Ca2+-Mg2+)-ATPase tryptophan residues and dansyl derivative of covalently labeled CaM occurred in the presence of EGTA, but was further promoted by Ca2+. It is concluded that the interaction of CaM and PS with (Ca2+-Mg2+)-ATPase results in different conformational states of the enzyme. CD and fluorescence spectroscopy allowed to distinguish these states from the E1 and E2 conformational forms of the ATPase.     

Functional consequences of alterations to Thr247, Pro248, Glu340, Asp813, Arg819, and Arg822 at the interfaces between domain P, M3, and L6-7 of sarcoplasmic reticulum Ca2+-ATPase. Roles in Ca2+ interaction and phosphoenzyme processing

Point mutants with alterations to amino acid residues Thr(247), Pro(248), Glu(340), Asp(813), Arg(819), and Arg(822) of sarcoplasmic reticulum Ca(2+)-ATPase were analyzed by transient kinetic measurements. In the Ca(2+)-ATPase crystal structures, most of these residues participate in a hydrogen-bonding network between the phosphorylation domain (domain P), the third transmembrane helix (M3), and the cytoplasmic loop connecting the sixth and the seventh transmembrane helices (L6-7). In several of the mutants, a pronounced phosphorylation "overshoot" was observed upon reaction of the Ca(2+)-bound enzyme with ATP, because of accumulation of dephosphoenzyme at steady state. Mutations of Glu(340) and its partners, Thr(247) and Arg(822), in the bonding network markedly slowed the Ca(2+) binding transition (E2 --> E1 --> Ca(2)E1) as well as Ca(2+) dissociation from Ca(2+) site II back toward the cytosol but did not affect the apparent affinity for vanadate. These mutations may have caused a slowing, in both directions, of the conformational change associated directly with Ca(2+) interaction at Ca(2+) site II. Because mutation of Asp(813) inhibited the Ca(2+) binding transition, but not Ca(2+) dissociation, and increased the apparent affinity for vanadate, the effect on the Ca(2+) binding transition seems in this case to be exerted by slowing the E2 --> E1 conformational change. Because the rate was not significantly enhanced by a 10-fold increase of the Ca(2+) concentration, the slowing is not the consequence of reduced affinity of any pre-binding site for Ca(2+). Furthermore, the mutations interfered in specific ways with the phosphoenzyme processing steps of the transport cycle; the transition from ADP-sensitive phosphoenzyme to ADP-insensitive phosphoenzyme (Ca(2)E1P --> E2P) was accelerated by mutations perturbing the interactions mediated by Glu(340) and Asp(813) and inhibited by mutation of Pro(248), and mutations of Thr(247) induced charge-specific changes of the rate of dephosphorylation of E2P.     

Immunological relatedness of the sarcoplasmic reticulum Ca(2+)-ATPase and the Na+,K(+)-ATPase.

The effect of anti-ATPase antibodies with epitopes near Asp-351 (PR-8), Lys-515 (PR-11) and the ATP binding domain (D12) of the Ca(2+)-ATPase of sarcoplasmic reticulum (EC 3.6.1.38) was analyzed. The PR-8 and D12 antibodies reacted freely with the Ca(2+)-ATPase in the native membrane, indicating that their epitopes are exposed on the cytoplasmic surface. Both PR-8 and D12 interfered with the crystallization of the Ca(2+)-ATPase, suggesting that their binding sites are at interfaces between ATPase molecules. PR-11 had no effect on ATPase-ATPase interactions or on the ATPase activity of sarcoplasmic reticulum. The epitope of PR-11 is suggested to be the VIDRC sequence at residues 520-525, while that of D12 at residues 670-720 of the Ca(2+)-ATPase. The use of predictive algorithms of antigenicity for identification of potential antigenic determinants in the Ca(2+)-ATPase is analyzed.     

Conformational fluctuations of the Ca2+-ATPase in the native membrane environment. Effects of pH, temperature, catalytic substrates, and thapsigargin

Digestion with proteinase K or trypsin yields complementary information on conformational transitions of the Ca(2+)-ATPase (SERCA) in the native membrane environment. Distinct digestion patterns are obtained with proteinase K, revealing interconversion of E1 and E2 or E1 approximately P and E2-P states. The pH dependence of digestion patterns shows that, in the presence of Mg(2+), conversion of E2 to E1 pattern occurs (even when Ca(2+) is absent) as H(+) dissociates from acidic residues. Mutational analysis demonstrates that the Glu(309) and Glu(771) acidic residues (empty Ca(2+)-binding sites I and II) are required for stabilization of E2. Glu(309) ionization is most important to yield E1. However, a further transition produced by Ca(2+) binding to E1 (i.e. E1.2Ca(2+)) is still needed for catalytic activation. Following ATP utilization, H(+)/Ca(2+) exchange is involved in the transition from the E1 approximately P.2Ca(2+) to the E2-P pattern, whereby alkaline pH will limit this conformational transition. Complementary experiments on digestion with trypsin exhibit high temperature dependence, indicating that, in the E1 and E2 ground states, the ATPase conformation undergoes strong fluctuations related to internal protein dynamics. The fluctuations are tightly constrained by ATP binding and phosphoenzyme formation, and this constraint must be overcome by thermal activation and substrate-free energy to allow enzyme turnover. In fact, a substantial portion of ATP free energy is utilized for conformational work related to the E1 approximately P.2Ca(2+) to E2-P transition, thereby disrupting high affinity binding and allowing luminal diffusion of Ca(2+). The E2 state and luminal path closure follow removal of conformational constraint by phosphate.     

Phosphorylated Ca2+-ATPase stable enough for structural studies

The atomic structure of sarcoplasmic reticulum Ca(2+)-ATPase, in a Ca(2+)-bound conformation, has recently been elucidated (Toyoshima, C., Nakasako, M., Nomura, H. & Ogawa, H. (2000) Nature 405, 647-655). Important steps for further understanding the mechanism of ion pumps will be the atomic structural characterization of different key conformational intermediates of the transport cycle, including phosphorylated intermediates. Following our previous report (Champeil, P., Henao, F., Lacapère, J.-J. & McIntosh, D. B. (2000) J. Biol. Chem. 276, 5795-5803), we show here that it is possible to prepare a phosphorylated form of sarcoplasmic reticulum Ca(2+)-ATPase (labeled with fluorescein isothiocyanate) with a week-long stability both in membranes and in mixed lipid-detergent micelles. We show that this phosphorylated fluorescein isothiocyanate-ATPase can form two-dimensional arrays in membranes, similar to those that were used previously to reconstruct from cryoelectron microscopy images the three-dimensional structure of Ca(2+)-free unphosphorylated ATPase. The results also provide hope that crystals of phosphorylated Ca(2+)-ATPase suitable for x-ray crystallography will be achieved.     

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.     

Interaction of chlorpromazine with low molecular mass ion-transporting ATPase modulator proteins from rat brain cytosol.

Two low molecular mass proteins (13 kDa which inhibits Na+,K(+)-ATPase and 12 kDa which modulates Ca2+, Mg(2+)- and Ca(2+)-ATPases), purified from rat brain cytosol form complexes with chlorpromazine (CPZ) on incubation. The conformational characteristics of the proteins and their complex have been studied by comparing the fluorescence and CD spectra. The tryptophan fluorescence data show that the inhibitor-CPZ complex does not quench the fluorescence of NA+,K(+)-ATPase significantly. CD spectra indicate that the structure of the inhibitor is changed on formation of the complex. The inhibitor-CPZ complex significantly changes the conformation of Na+,K(+)-ATPase. The regulator protein-CPZ complex does not have any appreciable effect on Ca2+, Mg(2+)- and Ca(2+)-ATPase activities. The Trp-fluorescence of Ca2+,Mg(2+)- and Ca(2+)-ATPase are not significantly affected in presence of the complex. CD spectra indicate that the structure of the regulator is abruptly affected on formation of the complex. The conformations of Ca2+,Mg(2+)- and Ca(2+)-ATPases are found to be altered in presence of the complex.