Sampling unfolding intermediates in calmodulin by single-molecule spectroscopy

We used single-pair fluorescence resonance energy transfer (spFRET) measurements to characterize denatured and partially denatured states of the multidomain calcium signaling protein calmodulin (CaM) in both its apo and Ca(2+)-bound forms. The results demonstrate the existence of an unfolding intermediate. A CaM mutant (CaM-T34C-T110C) was doubly labeled with fluorescent probes AlexaFlour 488 and Texas Red at opposing globular domains. Single-molecule distributions of the distance between fluorophores were obtained by spFRET at varying levels of the denaturant urea. Multiple conformational states of CaM were observed, and the amplitude of each conformation was dependent on urea concentration, with the amplitude of an extended conformation increasing upon denaturation. The distributions at intermediate urea concentrations could not be adequately described as a combination of native and denatured conformations, showing that CaM does not denature via a two-state process and demonstrating that at least one intermediate is present. The intermediate conformations formed upon addition of urea were different for Ca(2+)-CaM and apoCaM. An increase in the amplitude of a compact conformation in CaM was observed for apoCaM but not for Ca(2+)-CAM upon the addition of urea. The changes in the single-molecule distributions of CaM upon denaturation can be described by either a range of intermediate structures or by the presence of a single unfolding intermediate that grows in amplitude upon denaturation. A model for stepwise unfolding of CaM is suggested in which the domains of CaM unfold sequentially.     

Revealing two-state protein-protein interactions of calmodulin by single-molecule spectroscopy

We report a single-molecule fluorescence resonance energy transfer (FRET) and polarization study of conformational dynamics of calmodulin (CaM) interacting with a target peptide, C28W of a 28 amino acid oligomer. The C28W peptide represents the essential binding sequence domain of the Ca-ATPase protein interacting with CaM, which is important in cellular signaling for the regulation of energy in metabolism. However, the mechanism of the CaM/C28W recognition complex formation is still unclear. The amino-terminal (N-terminal) domain of the CaM was labeled with a fluorescein-based arsenical hairpin binder (FlAsH) that enables our unambiguous probing of the CaM N-terminal target-binding domain motions on a millisecond time scale without convolution of the probe-dye random motions. By analyzing the distribution of FRET efficiency between FlAsH labeled CaM and Texas Red labeled C28W and the polarization fluctuation dynamics and distributions of the CaM N-terminal domain, we reveal binding-unbinding motions of the N-terminal domain of the CaM in CaM/C28W complexes, which is strong evidence of a two-state binding interaction of CaM-mediated cell signaling.     

Myoglobin-CO substate structures and dynamics: multidimensional vibrational echoes and molecular dynamics simulations

Spectrally resolved infrared stimulated vibrational echo data were obtained for sperm whale carbonmonoxymyoglobin (MbCO) at 300 K. The measured dephasing dynamics of the CO ligand are in agreement with dephasing dynamics calculated with molecular dynamics (MD) simulations for MbCO with the residue histidine-64 (His64) having its imidazole epsilon nitrogen protonated (N(epsilon)-H). The two conformational substate structures B(epsilon) and R(epsilon) observed in the MD simulations are assigned to the spectroscopic A(1) and A(3) conformational substates of MbCO, respectively, based on the agreement between the experimentally measured and calculated dephasing dynamics for these substates. In the A(1) substate, the N(epsilon)-H proton and N(delta) of His64 are approximately equidistant from the CO ligand, while in the A(3) substate, the N(epsilon)-H of His64 is oriented toward the CO, and the N(delta) is on the surface of the protein. The MD simulations show that dynamics of His64 represent the major source of vibrational dephasing of the CO ligand in the A(3) state on both femtosecond and picosecond time scales. Dephasing in the A(1) state is controlled by His64 on femtosecond time scales, and by the rest of the protein and the water solvent on longer time scales.     

The N6-methyl group of adenine further increases the BI stability of DNA compared to C5-methyl groups

Methylated DNA bases are natural modifications which play an important role in protein-DNA interactions. Recent experimental and theoretical results have shown an influence of the base modification on the conformational behavior of the DNA backbone. MD simulations of four different B-DNA dodecamers (d(GC)(6), d(AT)(6), d(G(5mCG)(5)C), and d(A(T6mA)(5)T)) have been performed with the aim to examine the influence of methyl groups on the B-DNA backbone behavior. An additional control simulation of d(AU)(6) has also been performed to examine the further influence of the C5-methyl group in thymine. Methyl groups in the major groove (as in C5-methylcytosine, thymine, or N6-methyladenine) decrease the BII substate population of RpY steps. Due to methylation a clearer distinction of the BI substate stability between YpR and RpY (CpG/GpC or TpA/ApT) steps arises. A positive correlation between the BII substate population and base stacking distances is seen only for poly(GC). A methyl group added into the major groove increases mean water residence times around the purine N7 atom, which may stabilize the BI substate by improving the hydration network between the DNA backbone and the major groove. The N6-methyl group also forms a water molecule bridge between the N6 and O4 atoms, and thus further stabilizes the BI substate.     

Comparison of negative and positive ion electrospray ionization mass spectra of calmodulin and its complex with trifluoperazine

The protein calmodulin (apoCaM) undergoes a conformational change when it binds calcium. This structure of the protein (Ca4CaM) is a dumbbell-shaped molecule that undergoes a further profound conformational change on binding of the antipsychotic drug trifluoperazine (TFP). Experimental conditions were developed to prepare samples of apoCaM, Ca4CaM and Ca4CaM/TFP that were substantially free of sodium. The effects of the conformational changes of calmodulin on the charge-state distributions observed in positive ion and negative ion electrospray ionization (ESI) mass spectra were examined. Conversion of apoCaM into Ca4CaM was concomitant with a change in the negative ion ESI mass spectrum whereby the 16- ion was the most abundant ion observed for the apo form and the 8- ion was the most abundant for the complex. In contrast, in the positive ion ESI mass spectra of apoCaM and Ca4CaM, the most abundant species in each case was the 8+ ion. When a complex of Ca4CaMwith TFP was prepared, the most abundant species was the 5+ ion. This is consistent with a conformational change of Ca4CaM that rendered some basic sites inaccessible to ionization in the ESI process. Using the same Ca4CaM/TFP mixture, no complex with TFP was observed in negative ion ESI mass spectra. These observations are discussed in the context of the structural changes that are known to occur in calmodulin, and suggestions are made to explain the apparently conflicting data. The results reported here reflect on the validity of using differences in charge-state distributions observed in ESI mass spectra to assess conformational changes in proteins.     

Distance measurements on spin-labelled biomacromolecules by pulsed electron paramagnetic resonance

The biological function of protein, DNA, and RNA molecules often depends on relative movements of domains with dimensions of a few nanometers. This length scale can be accessed by distance measurements between spin labels if pulsed electron paramagnetic resonance (EPR) techniques such as electron-electron double resonance (ELDOR) and double-quantum EPR are used. The approach does not require crystalline samples and is well suited to biomacromolecules with an intrinsic flexibility as distributions of distances can be measured. Furthermore, oligomerization or complexation of biomacromolecules can also be studied, even if it is incomplete. The sensitivity of the technique and the reliability of the measured distance distribution depend on careful optimization of the experimental conditions and procedures for data analysis. Interpretation of spin-to-spin distance distributions in terms of the structure of the biomacromolecules furthermore requires a model for the conformational distribution of the spin labels.     

Subtle pH differences trigger single residue motions for moderating conformations of calmodulin

This study reveals the essence of ligand recognition mechanisms by which calmodulin (CaM) controls a variety of Ca(2+) signaling processes. We study eight forms of calcium-loaded CaM each with distinct conformational states. Reducing the structure to two degrees of freedom conveniently describes main features of the conformational changes of CaM via simultaneous twist-bend motions of the two lobes. We utilize perturbation-response scanning (PRS) technique, coupled with molecular dynamics simulations. PRS is based on linear response theory, comprising sequential application of directed forces on selected residues followed by recording the resulting protein coordinates. We analyze directional preferences of the perturbations and resulting conformational changes. Manipulation of a single residue reproduces the structural change more effectively than that of single/pairs/triplets of collective modes of motion. Our findings also give information on how the flexible linker acts as a transducer of binding information to distant parts of the protein. Furthermore, by perturbing residue E31 located in one of the EF hand motifs in a specific direction, it is possible to induce conformational change relevant to five target structures. Independently, using four different pK(a) calculation strategies, we find this particular residue to be the charged residue (out of a total of 52), whose ionization state is most sensitive to subtle pH variations in the physiological range. It is plausible that at relatively low pH, CaM structure is less flexible. By gaining charged states at specific sites at a pH value around 7, such as E31 found in the present study, local conformational changes in the protein will lead to shifts in the energy landscape, paving the way to other conformational states. These findings are in accordance with Fluorescence Resonance Energy Transfer (FRET) measured shifts in conformational distributions towards more compact forms with decreased pH. They also corroborate mutational studies and proteolysis results which point to the significant role of E31 in CaM dynamics.     

Drug detection based on the conformational changes of calmodulin and the fluorescence of its enhanced green fluorescent protein fusion partner

A homogeneous phase protein-based assay for the high throughput screening of drugs was developed using enhanced green fluorescent protein (EGFP) as the reporter. For that, a fusion protein between calmodulin (CaM) and EGFP was constructed in order to monitor the conformational changes induced in CaM upon binding to tricyclic anti-depressant drugs. In the presence of Ca²⁺, CaM undergoes a conformational change exposing a hydrophobic pocket that interacts with CaM-binding proteins, peptides, and drugs. Further, the conformational changes induced in CaM upon binding to Ca²⁺ and the target analyte drug, leads to a change in the microenvironment of EGFP concomitant with a change in its fluorescence intensity. The observed change in fluorescence intensity of EGFP can be correlated to the concentration of the analyte present in the sample. Further, the response of CaM–EGFP fusion protein in the presence of Ca²⁺ to increasing concentrations of phenothiazines and structurally related tricyclic anti-depressants was investigated. Dose-response curves for various tricyclic anti-depressants were prepared. Moreover, this assay can serve as a model system for other homogeneous binding assays for pharmaceuticals employing genetically fused binding proteins with reporter proteins and may find applications in the high throughput screening of tricyclic anti-depressants.     

Freezing a single distal motion in dihydrofolate reductase

Constraining a single motion between distal residues separated by approximately 28 A in hybrid quantum/classical molecular dynamics simulations is found to increase the free energy barrier for hydride transfer in dihydrofolate reductase by approximately 3 kcal/mol. Our analysis indicates that a single distal constraint alters equilibrium motions throughout the enzyme on a wide range of time scales. This alteration of the conformational sampling of the entire system is sufficient to significantly increase the free energy barrier and decrease the rate of hydride transfer. Despite the changes in conformational sampling introduced by the constraint, the system assumes a similar transition state conformation with a donor-acceptor distance of approximately 2.72 A to enable the hydride transfer reaction. The modified thermal sampling leads to a substantial increase in the average donor-acceptor distance for the reactant state, however, thereby decreasing the probability of sampling the transition state conformations with the shorter distances required for hydride transfer. These simulations indicate that fast thermal fluctuations of the enzyme, substrate, and cofactor lead to conformational sampling of configurations that facilitate hydride transfer. The fast thermal motions are in equilibrium as the reaction progresses along the collective reaction coordinate, and the overall average equilibrium conformational changes occur on the slower time scale measured experimentally. Recent single molecule experiments suggest that at least some of these thermally averaged equilibrium conformational changes occur on the millisecond time scale of the hydride transfer reaction. Thus, introducing a constraint that modifies the conformational sampling of an enzyme could significantly impact its catalytic activity.     

Mechanism of Nitric Oxide Synthase Regulation: Electron Transfer and Interdomain Interactions

Nitric oxide synthase (NOS), a flavo-hemoprotein, tightly regulates nitric oxide (NO) synthesis and thereby its dual biological activities as a key signaling molecule for vasodilatation and neurotransmission at low concentrations, and also as a defensive cytotoxin at higher concentrations. Three NOS isoforms, iNOS, eNOS and nNOS (inducible, endothelial, and neuronal NOS), achieve their key biological functions by tight regulation of interdomain electron transfer (IET) process via interdomain interactions. In particular, the FMN-heme IET is essential in coupling electron transfer in the reductase domain with NO synthesis in the heme domain by delivery of electrons required for O(2) activation at the catalytic heme site. Compelling evidence indicates that calmodulin (CaM) activates NO synthesis in eNOS and nNOS through a conformational change of the FMN domain from its shielded electron-accepting (input) state to a new electron-donating (output) state, and that CaM is also required for proper alignment of the domains. Another exciting recent development in NOS enzymology is the discovery of importance of the the FMN domain motions in modulating reactivity and structure of the catalytic heme active site (in addition to the primary role of controlling the IET processes). In the absence of a structure of full-length NOS, an integrated approach of spectroscopic (e.g. pulsed EPR, MCD, resonance Raman), rapid kinetics (laser flash photolysis and stopped flow) and mutagenesis methods is critical to unravel the molecular details of the interdomain FMN/heme interactions. This is to investigate the roles of dynamic conformational changes of the FMN domain and the docking between the primary functional FMN and heme domains in regulating NOS activity. The recent developments in understanding of mechanisms of the NOS regulation that are driven by the combined approach are the focuses of this review. An improved understanding of the role of interdomain FMN/heme interaction and CaM binding may serve as the basis for the design of new selective inhibitors of NOS isoforms.     

Binding studies of nNOS-active amphibian peptides and Ca2+ calmodulin, using negative ion electrospray ionisation mass spectrometry

Amphibian peptides which inhibit the formation of nitric oxide by neuronal nitric oxide synthase (nNOS) do so by binding to the protein cofactor, Ca2+calmodulin (Ca2+CaM). Complex formation between active peptides and Ca2+CaM has been demonstrated by negative ion electrospray ionisation mass spectrometry using an aqueous ammonium acetate buffer system. In all cases studied, the assemblies are formed with a 1:1:4 calmodulin/peptide/Ca2+ stoichiometry. In contrast, the complex involving the 20-residue binding domain of the plasma Ca2+ pump C20W (LRRGQILWFRGLNRIQTQIK-OH) with CaM has been shown by previous two-dimensional nuclear magnetic resonance (2D NMR) studies to involve complexation of the C-terminal end of CaM. Under identical conditions to those used for the amphibian peptide study, the ESI complex between C20W and CaM shows specific 1:1:2 stoichiometry. Since complex formation with the studied amphibian peptides requires Ca2+CaM to contain its full complement of four Ca2+ ions, this indicates that the amphibian peptides require both ends of the CaM to effect complex formation. Charge-state analysis and an H/D exchange experiment (with caerin 1.8) suggest that complexation involves Ca2+CaM undergoing a conformational change to a more compact structure.     

Fourier-transform infrared study of the switching process in a ferroelectric liquid crystalline polymer

The implementation of the step-scan technique on our FT-IR spectrometer enabled us to follow the electric-field induced reorientation dynamics of different molecular segments of a ferroelectric liquid crystalline polymer on a sub-millisecond time scale. It was detected that not only the mesogen but also the spacer and at least part of the backbone take part in the reorientation process.     

Investigating and characterizing the binding activity of the immobilized calmodulin to calmodulin-dependent protein kinase I binding domain with atomic force microscopy

Protein–protein interactions are responsible for many biological processes, and the study of how proteins undergo a conformational change induced by other proteins in the immobilized state can help us to understand a protein’s function and behavior, empower the current knowledge on molecular etiology of disease, as well as the discovery of putative protein targets of therapeutic interest. In this study, a bottom-up approach was utilized to fabricate micro/nanometer-scale protein patterns. One cysteine mutated calmodulin (CaM), as a model protein, was immobilized on thiol-terminated pattern surfaces. Atomic Force Microscopy (AFM) was then employed as a tool to investigate the interactions between CaM and CaM kinase I binding domain, and show that the immobilized CaM retains its activity to interact with its target protein. Our work demonstrate the potential of employing AFM to the research and assay works evolving surface-based protein–protein interactions biosensors, bioelectronics or drug screening.     

Rate processes with non-Markovian dynamical disorder

Rate processes with dynamical disorder are investigated within a simple framework provided by unidirectional electron transfer (ET) with fluctuating transfer rate. The rate fluctuations are assumed to be described by a non-Markovian stochastic jump process which reflects conformational dynamics of an electron transferring donor-acceptor molecular complex. A tractable analytical expression is obtained for the relaxation of the donor population (in the Laplace-transformed time domain) averaged over the stationary conformational fluctuations. The corresponding mean transfer time is also obtained in an analytical form. The case of two-state fluctuations is studied in detail for a model incorporating substate diffusion within one of the conformations. It is shown that an increase of the conformational diffusion time results in a gradual transition from the regime of fast modulation characterized by the averaged ET rate to the regime of quasistatic disorder. This transition occurs at the conformational mean residence time intervals fixed much less than the inverse of the corresponding ET rates. An explanation of this paradoxical effect is provided. Moreover, its presence is also manifested for the simplest, exactly solvable non-Markovian model with a biexponential distribution of the residence times in one of the conformations. The nontrivial conditions for this phenomenon to occur are found.     

Electron transfer between cofactors in protein domains linked by a flexible tether

Some key electron-transfer (ET) proteins have domains containing redox cofactors connected by a flexible tether. The relative motion of the domains is essential for reaction because of the strong ET rate dependence on distance. The constrained conformational flexibility produces a kinetic regime intermediate between “unimolecular” and “bimolecular”. We used a simple model for ET coupled to conformational diffusion to explore the structure dependence of the ET kinetics. The model describes the evolution of an initially prepared reduced donor state, and recovers the diffusion and electron tunneling-limited regimes. The restriction of the conformational space imposed by the tether introduces an entropic component to the effective donor–acceptor interaction potential. As such, the tether length may control the transition between the electron tunneling and diffusion-limited regimes.     

A whole-cell assay for the high throughput screening of calmodulin antagonists

Cell-based screening systems for pharmaceuticals are desired over molecular biosensing systems because of the information they provide on toxicity and bioavailability. However, the majority of sensing systems developed are molecular biosensing type screening systems and cannot be easily adapted to cell-based screening. In this study, we demonstrate that protein-based molecular sensing systems that employ a fluorescent protein as a signal transducer are amenable to cell-based sensing by expressing the protein molecular sensing system in the cell and employing these cells for screening of desired molecules. To achieve this, we expressed a molecular sensing system based on the fusion protein of calmodulin (CaM) and enhanced green fluorescent protein (EGFP) in bacterial cells, and utilized these cells for the screening of CaM antagonists. In the presence of Ca²⁺, CaM undergoes a conformational change exposing a hydrophobic pocket that interacts with CaM-binding proteins, peptides, and drugs. This conformational change induced in CaM leads to a change in the microenvironment of EGFP, resulting in a change in its fluorescence intensity. The observed change in fluorescence intensity of EGFP can be correlated to the concentration of the analyte present in the sample. Dose-response curves for various tricyclic antidepressants were generated using cells containing CaM-EGFP fusion protein. Additionally, we demonstrate the versatility of our system for studying protein-protein interactions by using cells to study the binding of a peptide to CaM. The study showed that the CaM-EGFP fusion protein within the intact cells responds similarly to that of the isolated fusion protein, hence eliminating the need for any isolation and purification steps. We have demonstrated that this system can be used for the rapid screening of various CaM antagonists that are potential antipsychotic drugs.     

Direct measurement by laser flash photolysis of intraprotein electron transfer in a rat neuronal nitric oxide synthase

Intraprotein interdomain electron transfer (IET) from flavin mononucleotide (FMN) to heme is essential in nitric oxide (NO) synthesis by NO synthase (NOS). Our previous laser flash photolysis studies have provided a direct determination of the kinetics of IET between the FMN and heme domains in truncated oxyFMN constructs of rat neuronal NOS (nNOS) and murine inducible NOS (iNOS), in which only the oxygenase and FMN domains along with the calmodulin (CaM) binding site are present [Feng, C. J.; Tollin, G.; Holliday, M. A.; Thomas, C.; Salerno, J. C.; Enemark, J. H.; Ghosh, D. K. Biochemistry 2006, 45, 6354-6362. Feng, C. J.; Thomas, C.; Holliday, M. A.; Tollin, G.; Salerno, J. C.; Ghosh, D. K.; Enemark, J. H. J. Am. Chem. Soc. 2006, 128, 3808-3811]. Here, we report the kinetics of IET between the FMN and heme domains in a rat nNOS holoenzyme in the presence and absence of added CaM using laser flash photolysis of CO dissociation in comparative studies on partially reduced NOS and a single domain NOS oxygenase construct. The IET rate constant in the presence of CaM is 36 s-1, whereas no IET was observed in the absence of CaM. The kinetics reported here are about an order of magnitude slower than the kinetics in a rat nNOS oxyFMN construct with added CaM (262 s-1). We attribute the slower IET between FMN and heme in the holoenzyme to the additional step of dissociation of the FMN domain from the reductase complex before reassociation with the oxygenase domain to form the electron-transfer competent output state complex. This work provides the first direct measurement of CaM-controlled electron transfer between catalytically significant redox couples of FMN and heme in a nNOS holoenzyme.     

Tracking Conformational Changes in Calmodulin in vitro,in Cell Extract,and in Cells by Electron Paramagnetic Resonance Distance Measurements

It is an open question whether the conformations of proteins sampled in dilute solutions are the same as in the cellular environment. Here we address this question by double electron-electron resonance (DEER) distance measurements with Gd(III) spin labels to probe the conformations of calmodulin (CaM) in vitro, in cell extract, and in human HeLa cells. Using the CaM mutants N53C/T110C and T34C/T117C labeled with maleimide-DOTA-Gd(III) in the N- and C-terminal domains, we observed broad and varied interdomain distance distributions. The in vitro distance distributions of apo-CaM and holo-CaM in the presence and absence of the IQ target peptide can be described by combinations of closed, open, and collapsed conformations. In cell extract, apo- and holo-CaM bind to target proteins in a similar way as apo- and holo-CaM bind to IQ peptide in vitro. In HeLa cells, however, in the presence or absence of elevated in-cell Ca²⁺ levels CaM unexpectedly produced more open conformations and very broad distance distributions indicative of many different interactions with in-cell components. These results show-case the importance of in-cell analyses of protein structures.