Proline cis-trans isomers in calbindin D9k observed by X-ray crystallography

Computational Investigation of Mechanisms for pH Modulation of Human Chloride Channels

Many transmembrane proteins are modulated by intracellular or extracellular pH. Investigation of pH dependence generally proceeds by mutagenesis of a wide set of amino acids, guided by properties such as amino-acid conservation and structure. Prediction of pKas can streamline this process, allowing rapid and effective identification of amino acids of interest with respect to pH dependence. Commencing with the calcium-activated chloride channel bestrophin 1, the carboxylate ligand structure around calcium sites relaxes in the absence of calcium, consistent with a measured lack of pH dependence. By contrast, less relaxation in the absence of calcium in TMEM16A, and maintenance of elevated carboxylate sidechain pKas, is suggested to give rise to pH-dependent chloride channel activity. This hypothesis, modulation of calcium/proton coupling and pH-dependent activity through the extent of structural relaxation, is shown to apply to the well-characterised cytosolic proteins calmodulin (pH-independent) and calbindin D9k (pH-dependent). Further application of destabilised, ionisable charge sites, or electrostatic frustration, is made to other human chloride channels (that are not calcium-activated), ClC-2, GABAA, and GlyR. Experimentally determined sites of pH modulation are readily identified. Structure-based tools for pKa prediction are freely available, allowing users to focus on mutagenesis studies, construct hypothetical proton pathways, and derive hypotheses such as the model for control of pH-dependent calcium activation through structural flexibility. Predicting altered pH dependence for mutations in ion channel disorders can support experimentation and, ultimately, clinical intervention.

Proline Isomerization: From the Chemistry and Biology to Therapeutic Opportunities

Proline isomerization, the process of interconversion between the cis- and trans-forms of proline, is an important and unique post-translational modification that can affect protein folding and conformations, and ultimately regulate protein functions and biological pathways. Although impactful, the importance and prevalence of proline isomerization as a regulation mechanism in biological systems have not been fully understood or recognized. Aiming to fill gaps and bring new awareness, we attempt to provide a wholistic review on proline isomerization that firstly covers what proline isomerization is and the basic chemistry behind it. In this section, we vividly show that the cause of the unique ability of proline to adopt both cis- and trans-conformations in significant abundance is rooted from the steric hindrance of these two forms being similar, which is different from that in linear residues. We then discuss how proline isomerization was discovered historically followed by an introduction to all three types of proline isomerases and how proline isomerization plays a role in various cellular responses, such as cell cycle regulation, DNA damage repair, T-cell activation, and ion channel gating. We then explore various human diseases that have been linked to the dysregulation of proline isomerization. Finally, we wrap up with the current stage of various inhibitors developed to target proline isomerases as a strategy for therapeutic development.

Benefits and constrains of covalency: the role of loop length in protein stability and ligand binding

Protein folding is governed by non-covalent interactions under the benefits and constraints of the covalent linkage of the backbone chain. In the current work we investigate the influence of loop length variation on the free energies of folding and ligand binding in a small globular single-domain protein containing two EF-hand subdomains—calbindin D_9k. We introduce a linker extension between the subdomains and vary its length between 1 to 16 glycine residues. We find a close to linear relationship between the linker length and the free energy of folding of the Ca²⁺-free protein. In contrast, the linker length has only a marginal effect on the Ca²⁺ affinity and cooperativity. The variant with a single-glycine extension displays slightly increased Ca²⁺ affinity, suggesting that the slightly extended linker allows optimized packing of the Ca²⁺-bound state. For the extreme case of disconnected subdomains, Ca²⁺ binding becomes coupled to folding and assembly. Still, a high affinity between the EF-hands causes the non-covalent pair to retain a relatively high apparent Ca²⁺ affinity. Our results imply that loop length variation could be an evolutionary option for modulating properties such as protein stability and turnover without compromising the energetics of the specific function of the protein.

Effect of monovalent ion binding on molecular dynamics of the S100-family calcium-binding protein calbindin D9k

Calbindin D_9k is a member of the S100 subfamily of EF-hand calcium binding proteins, and has served as an important model system for biophysical studies. The fast timescale dynamics of the calcium-free (apo) state is characterized using molecular dynamics simulations. Order parameters for the backbone NH bond vectors are determined from the simulations and compared with experimentally derived values, with a focus on the dynamics of calcium-binding site I. There is a significant discrepancy between simulated and experimental order parameters for site I residues in the case of no ion bound in site I. However, it was found in the simulations that a Na⁺ ion can bind in site I, and the resulting order parameters determined from the simulations are in excellent agreement with experiment. Comparisons are made to X-ray structures of other S100 family members in which Na⁺ ions were observed or suggested to be bound in site I. © 2019 Wiley Periodicals, Inc.

Structural and dynamic responses of calcium ion binding loop residues in metallo-proteins

Conformational changes in bio-molecular systems are fundamental to several biological processes. It is important to study changes in responses of underlying microscopic variables, like dihedral angles as conformational change takes place. We perform all-atom simulations and modelling via Langevin equation to illustrate the changes in structural and dynamic responses of dihedral angles of calcium ion binding residues of different proteins in metal ion free (apo) and bound (holo) states. The equilibrium distributions of dihedral angles in apo- and holo-states represent structural response. Our studies show the presence of dihedrals with multiple peaks (isomeric states) separated by barrier heights is more frequent in apo- than in holo-state. The relaxation time-scale of dihedral fluctuations is found to increase linearly with decreasing barrier height due to more frequent barrier re-crossing events. The slow kinetic response of the dihedrals also contributes to slowing down of macro-scale fluctuations, which may be useful to understand kinetics of various bio-molecular processes.

Accurate Electron-Nucleus Distances from Paramagnetic Relaxation Enhancements

Measurements of paramagnetic relaxation enhancements (PREs) in ¹H NMR spectra are an important tool to obtain long-range distance information in proteins, but quantitative interpretation is easily compromised by nonspecific intermolecular PREs. Here we show that PREs generated by lanthanides with anisotropic magnetic susceptibilities offer a route to accurate calibration-free distance measurements. As these lanthanides change ¹H chemical shifts due to pseudocontact shifts, the relaxation rates in the paramagnetic and diamagnetic state can be measured with a single sample that simultaneously contains the protein labeled with a paramagnetic and a diamagnetic lanthanide ion. Nonspecific intermolecular PREs are thus automatically subtracted when calculating the PREs as the difference in nuclear relaxation rates between paramagnetic and diamagnetic protein. Although PREs from lanthanides with anisotropic magnetic susceptibilities are complicated by additional cross-correlation effects and residual dipolar couplings (RDCs) in the paramagnetic state, these effects can be controlled by the choice of lanthanide ion and experimental conditions. Using calbindin D_9k with erbium, we succeeded in measuring intramolecular PREs with unprecedented accuracy, resulting in distance predictions with a root-mean-square-deviation of <0.9 Å in the range 11-24 Å.

Allosteric Fine-Tuning of the Binding Pocket Dynamics in the ITK SH2 Domain by a Distal Molecular Switch: An Atomistic Perspective

Although the regulation of function of proteins by allosteric interactions has been identified in many subcellular processes, molecular switches are also known to induce long-range conformational changes in proteins. A less well understood molecular switch involving cis-trans isomerization of a peptidyl-prolyl bond could induce a conformational change directly to the backbone that is propagated to other parts of the protein. However, these switches are elusive and hard to identify because they are intrinsic to biomolecules that are inherently dynamic. Here, we explore the conformational dynamics and free energy landscape of the SH2 domain of interleukin-2-inducible T-cell or tyrosine kinase (ITK) to fully understand the conformational coupling between the distal cis-trans molecular switch and its binding pocket of the phosphotyrosine motif. We use multiple microsecond-long all-atom molecular dynamics simulations in explicit water for over a total of 60 μs. We show that cis-trans isomerization of the Asn286-Pro287 peptidyl-prolyl bond is directly coupled to the dynamics of the binding pocket of the phosphotyrosine motif, in agreement with previous NMR experiments. Unlike the cis state that is localized and less dynamic in a single free energy basin, the trans state samples two distinct conformations of the binding pocket-one that recognizes the phosphotyrosine motif and the other that is somewhat similar to that of the cis state. The results provide an atomic-level description of a less well understood allosteric regulation by a peptidyl-prolyl cis-trans molecular switch that could aid in the understanding of normal and aberrant subcellular processes and the identification of these elusive molecular switches in other proteins.

On the Ca2+ binding and conformational change in EF-hand domains: Experimental evidence of Ca2+-saturated intermediates of N-domain of calmodulin

Double mutation of Q41L and K75I in the N-domain of calmodulin (N-Cam) stabilizes the closed form of N-Cam such that binding of Ca²⁺ in solution no longer triggers a conformational change to the open form, and its Ca²⁺ binding affinity decreases dramatically. To further investigate the solvation effects on the structure, Ca²⁺ binding affinity and conformational dynamics of this N-Cam double mutant in the Ca²⁺ saturated state, we solved its X-ray structure. Surprisingly, the structure revealed an open conformation of the domain which contradicts its closed conformation in solution. Here we provide evidence that crystallization conditions were responsible for this Ca²⁺-saturated domain open conformation in the crystal. Importantly, we demonstrate that the presence of the crystallization co-precipitant and alcohols were able to induce a progressive opening of the closed form of this domain, in Ca²⁺ saturated state, in solution. However, in the Ca²⁺ depleted state, addition of alcohols was unable to induce any opening of this domain in solution. In addition, in the Ca²⁺ saturated state, the molecular dynamics simulations show that while N-Cam can populate the open and closed conformation, the N-Cam double mutant exclusively populates the closed conformation. Our results provide experimental evidence of intermediate conformations of Ca²⁺-saturated N-Cam in solution. We propose that conformational change of Ca²⁺ sensor EF-hand domains depends on solvation energetics, Ca²⁺ binding to promote the full open form, Ca²⁺ depleted state conformational dynamics, and the chemical properties of the molecules nearby key residues such as those at positions 41 and 75 in N-Cam.

Using Paramagnetism to Slow Down Nuclear Relaxation in Protein NMR

Paramagnetic metal ions accelerate nuclear spin relaxation; this effect is widely used for distance measurement and called paramagnetic relaxation enhancement (PRE). Theoretical predictions established that, under special circumstances, it is also possible to achieve a reduction in nuclear relaxation rates (negative PRE). This situation would occur if the mechanism of nuclear relaxation in the diamagnetic state is counterbalanced by a paramagnetic relaxation mechanism caused by the metal ion. Here we report the first experimental evidence for such a cross-correlation effect. Using a uniformly ¹⁵N-labeled mutant of calbindin D_9k loaded with either Tm³⁺ or Tb³⁺, reduced R₁ and R₂ relaxation rates of backbone ¹⁵N spins were observed compared with the diamagnetic reference (the same protein loaded with Y³⁺). The effect arises from the compensation of the chemical shift anisotropy tensor by the anisotropic dipolar shielding generated by the unpaired electron spin.

Nuclear overhauser spectroscopy of chiral CHD methylene groups

Nuclear magnetic resonance spectroscopy (NMR) can provide a great deal of information about structure and dynamics of biomolecules. The quality of an NMR structure strongly depends on the number of experimental observables and on their accurate conversion into geometric restraints. When distance restraints are derived from nuclear Overhauser effect spectroscopy (NOESY), stereo-specific assignments of prochiral atoms can contribute significantly to the accuracy of NMR structures of proteins and nucleic acids. Here we introduce a series of NOESY-based pulse sequences that can assist in the assignment of chiral CHD methylene protons in random fractionally deuterated proteins. Partial deuteration suppresses spin-diffusion between the two protons of CH2 groups that normally impedes the distinction of cross-relaxation networks for these two protons in NOESY spectra. Three and four-dimensional spectra allow one to distinguish cross-relaxation pathways involving either of the two methylene protons so that one can obtain stereospecific assignments. In addition, the analysis provides a large number of stereospecific distance restraints. Non-uniform sampling was used to ensure optimal signal resolution in 4D spectra and reduce ambiguities of the assignments. Automatic assignment procedures were modified for efficient and accurate stereospecific assignments during automated structure calculations based on 3D spectra. The protocol was applied to calcium-loaded calbindin D9k. A large number of stereospecific assignments lead to a significant improvement of the accuracy of the structure.

A microfluidic platform for quantitative measurements of effective protein charges and single ion binding in solution

Microfluidic electrophoresis enables the comparison of dry sequence and solvated protein charges, and the detection of protein–ion binding.

Structures and energies of the transition between two conformations of the alternate frame folding calbindin-D9k protein: a theoretical study

We carried out molecular dynamics simulations and energy calculations for the two states of the alternate frame folding (AFF) calbindin-D_9k protein and their conformational transition in Ca²⁺-free form to address their dynamical transition mechanism.

The structure of cardiac troponin C regulatory domain with bound Cd2+ reveals a closed conformation and unique ion coordination

The amino-terminal domain of cardiac troponin C (cNTnC) is an essential Ca(2+) sensor found in cardiomyocytes. It undergoes a conformational change upon Ca(2+) binding and transduces the signal to the rest of the troponin complex to initiate cardiac muscle contraction. Two classical EF-hand motifs (EF1 and EF2) are present in cNTnC. Under physiological conditions, only EF2 binds Ca(2+); EF1 is a vestigial site that has lost its function in binding Ca(2+) owing to amino-acid sequence changes during evolution. Proteins with EF-hand motifs are capable of binding divalent cations other than calcium. Here, the crystal structure of wild-type (WT) human cNTnC in complex with Cd(2+) is presented. The structure of Cd(2+)-bound cNTnC with the disease-related mutation L29Q, as well as a structure with the residue differences D2N, V28I, L29Q and G30D (NIQD), which have been shown to have functional importance in Ca(2+) sensing at lower temperatures in ectothermic species, have also been determined. The structures resemble the overall conformation of NMR structures of Ca(2+)-bound cNTnC, but differ significantly from a previous crystal structure of Cd(2+)-bound cNTnC in complex with deoxycholic acid. The subtle structural changes observed in the region near the mutations may play a role in the increased Ca(2+) affinity. The 1.4 Å resolution WT cNTnC structure, which is the highest resolution structure yet obtained for cardiac troponin C, reveals a Cd(2+) ion coordinated in the canonical pentagonal bipyramidal geometry in EF2 despite three residues in the loop being disordered. A Cd(2+) ion found in the vestigial ion-binding site of EF1 is coordinated in a noncanonical `distorted' octahedral geometry. A comparison of the ion coordination observed within EF-hand-containing proteins for which structures have been solved in the presence of Cd(2+) is presented. A refolded WT cNTnC structure is also presented.

Protein surface and core dynamics show concerted hydration-dependent activation

By specifically labeling leucine/valine methyl groups and lysine side chains "inside" and "outside" dynamics of proteins on the nanosecond timescale are compared using neutron scattering. Surprisingly, both groups display similar dynamics as a function of temperature, and the buried hydrophobic core is sensitive to hydration and undergoes a dynamical transition.

Calculation of Free Energy Landscape in Multi-Dimensions with Hamiltonian-Exchange Umbrella Sampling on Petascale Supercomputer

An extremely scalable computational strategy is described for calculations of the potential of mean force (PMF) in multidimensions on massively distributed supercomputers. The approach involves coupling thousands of umbrella sampling (US) simulation windows distributed to cover the space of order parameters with a Hamiltonian molecular dynamics replica-exchange (H-REMD) algorithm to enhance the sampling of each simulation. In the present application, US/H-REMD is carried out in a two-dimensional (2D) space and exchanges are attempted alternatively along the two axes corresponding to the two order parameters. The US/H-REMD strategy is implemented on the basis of parallel/parallel multiple copy protocol at the MPI level, and therefore can fully exploit computing power of large-scale supercomputers. Here the novel technique is illustrated using the leadership supercomputer IBM Blue Gene/P with an application to a typical biomolecular calculation of general interest, namely the binding of calcium ions to the small protein Calbindin D9k. The free energy landscape associated with two order parameters, the distance between the ion and its binding pocket and the root-mean-square deviation (rmsd) of the binding pocket relative the crystal structure, was calculated using the US/H-REMD method. The results are then used to estimate the absolute binding free energy of calcium ion to Calbindin D9k. The tests demonstrate that the 2D US/H-REMD scheme greatly accelerates the configurational sampling of the binding pocket, thereby improving the convergence of the potential of mean force calculation.

Purification, crystallization and preliminary crystallographic analysis of the SH2 domain of IL-2-inducible T-cell kinase

Proline is a unique amino acid owing to the relatively small energy difference between the cis and trans conformations of its peptide bond. The X-Pro imide bond readily undergoes cis-trans isomerization in the context of short peptides as well as some proteins. However, the direct detection of cis-trans proline isomerization in folded proteins is technically challenging. NMR spectroscopy is well suited to the direct detection of proline isomerization in folded proteins. It is less clear how well X-ray crystallography can reveal this conformational exchange event in folded proteins. Conformational heterogeneity owing to cis-trans proline isomerization in the Src homology 2 (SH2) domain of the IL-2-inducible T-cell kinase (ITK) has been extensively characterized by NMR. Using the ITK SH2 domain as a test system, an attempt was made to determine whether proline isomerization could be detected in a crystal structure of the ITK SH2 domain. As a first step towards this goal, the purification, crystallization and preliminary characterization of the ITK SH2 domain are described.

Insights into modulation of calcium signaling by magnesium in calmodulin, troponin C and related EF-hand proteins

The Ca(2+)-binding helix-loop-helix structural motif called "EF-hand" is a common building block of a large family of proteins that function as intracellular Ca(2+)-receptors. These proteins respond specifically to micromolar concentrations of Ca(2+) in the presence of ~1000-fold excess of the chemically similar divalent cation Mg(2+). The intracellular free Mg(2+) concentration is tightly controlled in a narrow range of 0.5-1.0mM, which at the resting Ca(2+) levels is sufficient to fully or partially saturate the Ca(2+)-binding sites of many EF-hand proteins. Thus, to convey Ca(2+) signals, EF-hand proteins must respond differently to Ca(2+) than to Mg(2+). In this review the structural aspects of Mg(2+) binding to EF-hand proteins are considered and interpreted in light of the recently proposed two-step Ca(2+)-binding mechanism (Grabarek, Z., J. Mol. Biol., 2005, 346, 1351). It is proposed that, due to stereochemical constraints imposed by the two-EF-hand domain structure, the smaller Mg(2+) ion cannot engage the ligands of an EF-hand in the same way as Ca(2+) and defaults to stabilizing the apo-like conformation of the EF-hand. It is proposed that Mg(2+) plays an active role in the Ca(2+)-dependent regulation of cellular processes by stabilizing the "off state" of some EF-hand proteins, thereby facilitating switching off their respective target enzymes at the resting Ca(2+) levels. Therefore, some pathological conditions attributed to Mg(2+) deficiency might be related to excessive activation of underlying Ca(2+)-regulated cellular processes. This article is part of a Special Issue entitled: 11th European Symposium on Calcium.

Force field-dependent structural divergence revealed during long time simulations of Calbindin d9k

The structural and the dynamic features of the Calbindin (CaB) protein in its holo and apo states are compared using molecular dynamics simulations under nine different force fields (FFs) (G43a1, G53a6, Opls-AA, Amber94, Amber99, Amber99p, AmberGS, AmberGSs, and Amber99sb). The results show that most FFs reproduce reasonably well the majority of the experimentally derived features of the CaB protein. However, in several cases, there are significant differences in secondary structure properties, root mean square deviations (RMSDs), root mean square fluctuations (RMSFs), and S(2) order parameters among the various FFs. What is more, in certain cases, these parameters differed from the experimentally derived values. Some of these deviations became noticeable only after 50 ns. A comparison with experimental data indicates that, for CaB, the Amber94 shows overall best agreement with the measured values, whereas several others seem to deviate from both crystal and nuclear magnetic resonance data.

Comprehensive determination of (3)J (HNHalpha) for unfolded proteins using (13)C'-resolved spin-echo difference spectroscopy

An experiment is presented to determine (3)J(HNHalpha) coupling constants, with significant advantages for applications to unfolded proteins. The determination of coupling constants for the peptide chain using 1D (1)H, or 2D and 3D (1)H-(15)N correlation spectroscopy is often hampered by extensive resonance overlap when dealing with flexible, disordered proteins. In the experiment detailed here, the overlap problem is largely circumvented by recording (1)H-(13)C' correlation spectra, which demonstrate superior resolution for unfolded proteins. J-coupling constants are extracted from the peak intensities in a pair of 2D spin-echo difference experiments, affording rapid acquisition of the coupling data. In an application to the cytoplasmic domain of human neuroligin-3 (hNlg3cyt) data were obtained for 78 residues, compared to 54 coupling constants obtained from a 3D HNHA experiment. The coupling constants suggest that hNlg3cyt is intrinsically disordered, with little propensity for structure.

Multiscale characterization of protein conformational ensembles

We propose a multiscale exploration method to characterize the conformational space populated by a protein at equilibrium. The method efficiently obtains a large set of equilibrium conformations in two stages: first exploring the entire space at a coarse-grained level of detail, then narrowing a refined exploration to selected low-energy regions. The coarse-grained exploration periodically adds all-atom detail to selected conformations to ensure that the search leads to regions which maintain low energies in all-atom detail. The second stage reconstructs selected low-energy coarse-grained conformations in all-atom detail. A low-dimensional energy landscape associated with all-atom conformations allows focusing the exploration to energy minima and their conformational ensembles. The lowest energy ensembles are enriched with additional all-atom conformations through further multiscale exploration. The lowest energy ensembles obtained from the application of the method to three different proteins correctly capture the known functional states of the considered systems.

De novo protein structure generation from incomplete chemical shift assignments

NMR chemical shifts provide important local structural information for proteins. Consistent structure generation from NMR chemical shift data has recently become feasible for proteins with sizes of up to 130 residues, and such structures are of a quality comparable to those obtained with the standard NMR protocol. This study investigates the influence of the completeness of chemical shift assignments on structures generated from chemical shifts. The Chemical-Shift-Rosetta (CS-Rosetta) protocol was used for de novo protein structure generation with various degrees of completeness of the chemical shift assignment, simulated by omission of entries in the experimental chemical shift data previously used for the initial demonstration of the CS-Rosetta approach. In addition, a new CS-Rosetta protocol is described that improves robustness of the method for proteins with missing or erroneous NMR chemical shift input data. This strategy, which uses traditional Rosetta for pre-filtering of the fragment selection process, is demonstrated for two paramagnetic proteins and also for two proteins with solid-state NMR chemical shift assignments.

Overview of protein structural and functional folds

This overview provides an illustrated, comprehensive survey of some commonly observed protein‐fold families and structural motifs, chosen for their functional significance. It opens with descriptions and definitions of the various elements of protein structure and associated terminology. Following is an introduction into web‐based structural bioinformatics that includes surveys of interactive web servers for protein fold or domain annotation, protein‐structure databases, protein‐structure‐classification databases, structural alignments of proteins, and molecular graphics programs available for personal computers. The rest of the overview describes selected families of protein folds in terms of their secondary, tertiary, and quaternary structural arrangements, including ribbon‐diagram examples, tables of representative structures with references, and brief explanations pointing out their respective biological and functional significance.

Structure of Ca2+-bound S100A4 and its interaction with peptides derived from nonmuscle myosin-IIA

S100A4, also known as mts1, is a member of the S100 family of Ca2+-binding proteins that is directly involved in tumor invasion and metastasis via interactions with specific protein targets, including nonmuscle myosin-IIA (MIIA). Human S100A4 binds two Ca2+ ions with the typical EF-hand exhibiting an affinity that is nearly 1 order of magnitude tighter than that of the pseudo-EF-hand. To examine how Ca2+ modifies the overall organization and structure of the protein, we determined the 1.7 A crystal structure of the human Ca2+-S100A4. Ca2+ binding induces a large reorientation of helix 3 in the typical EF-hand. This reorganization exposes a hydrophobic cleft that is comprised of residues from the hinge region,helix 3, and helix 4, which afford specific target recognition and binding. The Ca2+-dependent conformational change is required for S100A4 to bind peptide sequences derived from the C-terminal portion of the MIIA rod with submicromolar affinity. In addition, the level of binding of Ca2+ to both EF-hands increases by 1 order of magnitude in the presence of MIIA. NMR spectroscopy studies demonstrate that following titration with a MIIA peptide, the largest chemical shift perturbations and exchange broadening effects occur for residues in the hydrophobic pocket of Ca2+-S100A4. Most of these residues are not exposed in apo-S100A4 and explain the Ca2+ dependence of formation of theS100A4-MIIA complex. These studies provide the foundation for understanding S100A4 target recognition and may support the development of reagents that interfere with S100A4 function.

Networks of coupled rotators: relationship between structures and internal dynamics in metal-binding proteins. Applications to apo- and holo-calbindin

This article presents an analysis of the internal dynamics of the Ca2+-binding protein calbindin, based on the Networks of Coupled Rotators (NCRs) introduced recently. Several fundamental and practical issues raised by this approach are investigated. The roles of various parameters of the model are examined. The NCR model is shown to account for the modifications of the internal dynamics upon Ca2+ binding by calbindin. Two alternative strategies to estimate local internal effective correlation times of the protein are proposed, which offer good agreement between predictions and experiment.

"Four-dimensional" protein structures: examples from metalloproteins

The fact that an object, for example, a protein, possesses a three-dimensional structure seems an obvious concept. However, when the object is flexible, the concept is less obvious. Growing experimental data over several decades show that proteins are not rigid objects, but they may sample more or less wide ranges of different conformations. To stress this concept, we propose to call the range of sampled conformations the "fourth dimension" of the protein structure. Nuclear magnetic resonance is a precious technique to define this fourth dimension. Examples of conformational heterogeneity taken from the realm of metalloproteins and their functional implications are discussed.

Chemically accurate protein structures: validation of protein NMR structures by comparison of measured and predicted pKa values

A new method is presented for evaluating the quality of protein structures obtained by NMR. This method exploits the dependence between measurable chemical properties of a protein, namely pKa values of acidic residues, and protein structure. The accurate and fast empirical computational method employed by the PROPKA program ( http://www.propka.chem.uiowa.edu) allows the user to test the ability of a given structure to reproduce known pKa values, which in turn can be used as a criterion for the selection of more accurate structures. We demonstrate the feasibility of this novel idea for a series of proteins for which both NMR and X-ray structures, as well as pKa values of all ionizable residues, have been determined. For the 17 NMR ensembles used in this study, this criterion is shown effective in the elimination of a large number of NMR structure ensemble members.

Structural basis for diversity of the EF-hand calcium-binding proteins

The calcium binding proteins of the EF-hand super-family are involved in the regulation of all aspects of cell function. These proteins exhibit a great diversity of composition, structure, Ca2+-binding and target interaction properties. Here, our current understanding of the Ca2+-binding mechanism is assessed. The structures of the EF-hand motifs containing 11-14 amino acid residues in the Ca2+-binding loop are analyzed within the framework of the recently proposed two-step Ca2+-binding mechanism. A hypothesis is put forward that in all EF-hand proteins the Ca2+-binding and the resultant conformational responses are governed by the central structure connecting the Ca2+-binding loops in the two-EF-hand domain. This structure, named EFbeta-scaffold, defines the position of the bound Ca2+, and coordinates the function of the N-terminal (variable and flexible) with the C-terminal (invariable and rigid) parts of the Ca2+-binding loop. It is proposed that the nature of the first ligand of the Ca2+-binding loop is an important determinant of the conformational change. Additional factors, including the interhelical contacts, the length, structure and flexibility of the linker connecting the EF-hand motifs, and the overall energy balance provide the fine-tuning of the Ca2+-induced conformational change in the EF-hand proteins.

The crystal structure of the cis-proline to glycine variant (P114G) of ribonuclease A

Replacement of a cis-proline by glycine at position 114 in ribonuclease A leads to a large decrease in thermal stability and simplifies the refolding kinetics. A crystallographic approach was used to determine whether the decrease in thermal stability results from the presence of a cis glycine peptide bond, or from a localized structural rearrangement caused by the isomerization of the mutated cis 114 peptide bond. The structure was solved at 2.0 A resolution and refined to an R-factor of 19.5% and an R(free) of 21.9%. The overall conformation of the protein was similar to that of wild-type ribonuclease A; however, there was a large localized rearrangement of the mutated loop (residues 110-117-a 9.3 A shift of the Calpha atom of residue 114). The peptide bond before Gly114 is in the trans configuration. Interestingly, a large anomalous difference density was found near residue 114, and was attributed to a bound cesium ion present in the crystallization experiment. The trans isomeric configuration of the peptide bond in the folded state of this mutant is consistent with the refolding kinetics previously reported, and the associated protein conformational change provides an explanation for the decreased thermal stability.

Calcium-modulated S100 protein-phospholipid interactions. An NMR study of calbindin D9k and DPC

The cellular functions of several S100 proteins involve specific interactions with phospholipids and the cell membrane. The interactions between calbindin D(9k) (S100D) and the detergent dodecyl phosphocholine (DPC) were studied using NMR spectroscopy. In the absence of Ca(2+), the protein associates with DPC micelles. The micelle-associated state has intact helical secondary structures but no apparent tertiary fold. At neutral pH, Ca(2+)-loaded calbindin D(9k) does not associate with DPC micelles. However, a specific interaction is observed with individual DPC molecules at a site close to the linker between the two EF-hands. Binding to this site occurs only when Ca(2+) is bound to the protein. A reduction in pH in the absence of Ca(2+) increases the stability of the micelle-associated state. This along with the corresponding reduction in Ca(2+) affinity causes a transition to the micelle-associated state also in the presence of Ca(2+) when the pH is lowered. Site-specific analysis of the data indicates that calbindin D(9k) has a core of three tightly packed helices (A, B, and D), with a dynamic fourth helix (C) more loosely associated. Evidence is presented that the Ca(2+)-binding characteristics of the two EF-hands are distinctly different in a micelle environment. The role of calbindin D(9k) in the cell is discussed, along with the broader implications for the function of the S100 protein family.

13C direct detected experiments: optimization for paramagnetic signals

To optimize 13C direct detected experiments for the observation of signals close to a paramagnetic center, we have assessed the sensitivity of different sequences based on CO-Cali coherence transfer. Features of CACO experiments were tested for Calbindin D9k, in which one of the two native Ca2+ ions is replaced by the paramagnetic Ce3+ ion. We have studied the comparison of single vs multiple quantum coherence transfer evolution as well as the influence of in-phase vs anti-phase detection of 13CO signals and finally the comparison of a coherence transfer step based on a CyO in plane with respect to a Cy ali in plane. The acquisition of the anti-phase component of the signal, accomplished by the removal of the last refocusing steps, allowed the identification of some signals unobserved with other pathways. The structural dependency of paramagnetism-induced nuclear relaxation is such that the identification of the most suitable coherence transfer pathway is not known "a priori" but it is driven by the relative proximity of Cali and CO to the paramagnetic center.

Solution structure and internal dynamics of NSCP, a compact calcium-binding protein

The solution structure of Nereis diversicolor sarcoplasmic calcium-binding protein (NSCP) in the calcium-bound form was determined by NMR spectroscopy, distance geometry and simulated annealing. Based on 1859 NOE restraints and 262 angular restraints, 17 structures were generated with a rmsd of 0.87 A from the mean structure. The solution structure, which is highly similar to the structure obtained by X-ray crystallography, includes two open EF-hand domains, which are in close contact through their hydrophobic surfaces. The internal dynamics of the protein backbone were determined by studying amide hydrogen/deuterium exchange rates and 15N nuclear relaxation. The two methods revealed a highly compact and rigid structure, with greatly restricted mobility at the two termini. For most of the amide protons, the free energy of exchange-compatible structural opening is similar to the free energy of structural stability, suggesting that isotope exchange of these protons takes place through global unfolding of the protein. Enhanced conformational flexibility was noted in the unoccupied Ca2+-binding site II, as well as the neighbouring helices. Analysis of the experimental nuclear relaxation and the molecular dynamics simulations give very similar profiles for the backbone generalized order parameter (S2), a parameter related to the amplitude of fast (picosecond to nanosecond) movements of N(H)-H vectors. We also noted a significant correlation between this parameter, the exchange rate, and the crystallographic B factor along the sequence.

A united residue force-field for calcium-protein interactions

United-residue potentials are derived for interactions of the calcium cation with polypeptide chains in energy-based prediction of protein structure with a united-residue (UNRES) force-field. Specific potentials were derived for the interaction of the calcium cation with the Asp, Glu, Asn, and Gln side chains and the peptide group. The analytical expressions for the interaction energies for each of these amino acids were obtained by averaging the electrostatic interaction energy, expressed by a multipole series over the dihedral angles not considered in the united-residue model, that is, the side-chain dihedral angles chi and the dihedral angles lambda for the rotation of peptide groups about the C(alpha)...C(alpha) virtual-bond axes. For the side-chains that do not interact favorably with calcium, simple excluded-volume potentials were introduced. The parameters of the potentials were obtained from ab initio quantum mechanical calculations of model systems at the Restricted Hartree-Fock (RHF) level with the 6-31G(d,p) basis set. The energy surfaces of pairs consisting of Ca(2+)-acetate, Ca(2+)-propionate, Ca(2+)-acetamide, Ca(2+)-propionamide, and Ca(2+)-N-methylacetamide systems (modeling the Ca(2+)-Asp(-), Ca(2+)-Glu(-), Ca(2+)-Asn, Ca(2+)-Gln, and Ca(2+)-peptide group interactions) at different distances and orientations were calculated. For each pair, the restricted free energy (RFE) surfaces were calculated by numerical integration over the degrees of freedom lost when switching from the all-atom model to the united-residue model. Finally, the analytical expressions for each pair were fitted to the RFE surfaces. This force-field was able to distinguish the EF-hand motif from all potential binding sites in the crystal structures of bovine alpha-lactalbumin, whiting parvalbumin, calbindin D9K, and apo-calbindin D9K.

Ab initio simulations of protein-folding pathways by molecular dynamics with the united-residue model of polypeptide chains

We report the application of Langevin dynamics to the physics-based united-residue (UNRES) force field developed in our laboratory. Ten trajectories were run on seven proteins [PDB ID codes 1BDD (alpha; 46 residues), 1GAB (alpha; 47 residues), 1LQ7 (alpha; 67 residues), 1CLB (alpha; 75 residues), 1E0L (beta; 28 residues), and 1E0G (alpha+beta; 48 residues), and 1IGD (alpha+beta; 61 residues)] with the UNRES force field parameterized by using our recently developed method for obtaining a hierarchical structure of the energy landscape. All alpha-helical proteins and 1E0G folded to the native-like structures, whereas 1IGD and 1E0L yielded mostly nonnative alpha-helical folds although the native-like structures are lowest in energy for these two proteins, which can be attributed to neglecting the entropy factor in the current parameterization of UNRES. Average folding times for successful folding simulations were of the order of nanoseconds, whereas even the ultrafast-folding proteins fold only in microseconds, which implies that the UNRES time scale is approximately three orders of magnitude larger than the experimental time scale because the fast motions of the secondary degrees of freedom are averaged out. Folding with Langevin dynamics required 2-10 h of CPU time on average with a single AMD Athlon MP 2800+ processor depending on the size of the protein. With the advantage of parallel processing, this process leads to the possibility to explore thousands of folding pathways and to predict not only the native structure but also the folding scenario of a protein together with its quantitative kinetic and thermodynamic characteristics.

Electrostatic contributions to the kinetics and thermodynamics of protein assembly

The role of electrostatic interactions in the assembly of a native protein structure was studied using fragment complementation. Contributions of salt, pH, or surface charges to the kinetics and equilibrium of calbindin D(9k) reconstitution was measured in the presence of Ca(2+) using surface plasmon resonance and isothermal titration calorimetry. Whereas surface charge substitutions primarily affect the dissociation rate constant, the association rates are correlated with subdomain net charge in a way expected for Coulomb interactions. The affinity is reduced in all mutants, with the largest effect (260-fold) observed for the double mutant K25E+K29E. At low net charge, detailed charge distribution is important, and charges remote from the partner EF-hand have less influence than close ones. The effects of salt and pH on the reconstitution are smaller than mutational effects. The interaction between the wild-type EF-hands occurs with high affinity (K(A) = 1.3 x 10(10) M(-1); K(D) = 80 pM). The enthalpy of association is overall favorable and there appears to be a very large favorable entropic contribution from the desolvation of hydrophobic surfaces that become buried in the complex. Electrostatic interactions contribute significantly to the affinity between the subdomains, but other factors, such as hydrophobic interactions, dominate.

Monitoring the Early Steps of Unfolding of Dicalcium and Mono-Ce3+-Substituted Forms of P43M Calbindin D9k†

Early steps of unfolding of P43M Calbindin D(9k) have been evaluated by NMR spectroscopy on the native dicalcium and on the paramagnetic monocerium-substituted derivative. Although at 2 M GdmHCl the protein core maintains its overall folding and structure, amide (15)N R(2) measurements and cross correlation rates between N-H dipole-dipole relaxation and (15)N CSA relaxation reveal a closer and stronger packing of the hydrophobic interactions in the protein as a response to the presence of denaturing agents in solution. A complete reorientation of the Met43 side chain toward the hydrophobic core is accomplished by the disappearance of the millisecond dynamics observed on the native form of Calbindin D(9k), while cross correlation rates provide evidence that the two-way hydrogen bond between Leu23 and Val61 is broken or substantially weakened. The substitution of the calcium ion in site II with the paramagnetic Ce(3+) ion allowed us to obtain a number of long-range nonconventional constraints, namely, pseudocontact shifts, which were used, together with the NOEs collected on the native state, to monitor subtle structural variations occurring in the non-native state of the protein. Although the average rmsd between the structures of native and non-native states is small (0.48 A), structural rearrangements could be reliably identified. Our results provide unprecedented information about the behavior of Calbindin D(9k) during the early steps of unfolding. Furthermore, they constitute strong evidence of the efficiency of paramagnetism-based constraints in monitoring subtle structural changes that are beyond the sensitivity of an approach based only on NOE.

Correlation between normal modes in the 20-200 cm-1 frequency range and localized torsion motions related to certain collective motions in proteins

In certain biologically relevant collective motions, such as protein domain motions and sub-domain motions, large amplitude movements are localized in one or a few flexible regions consisting of a small number of residues. This paper explores the possible use of normal mode analysis in probing localized vibrational torsion motions in these flexible regions that may be related to certain collective motions. The normal modes of 10 structures of five proteins in different conformation (TRP repressor, calmodulin, calbindin D(9k), HIV-1 protease and troponin C), known to have shear or hinge domain or sub-domain motion, respectively, are analyzed. Our study identifies, for each structure, unique normal modes in the 20-200 cm-1 frequency range, whose corresponding motions are primarily concentrated in the region where large amplitude torsion movements of a known domain or sub-domain motion occur. This suggests possible correlation between normal modes at 20-200 cm-1 frequency range and initial fluctuational motions leading to localized collective motions in proteins, and thus the potential application of normal mode analysis in facilitating the study of biologically important localized motions in biomolecules.

Native states of adenylate kinase are two active sub-ensembles

There are two kinds of conformational forms of adenylate kinase (AK) in equilibrium in solution with different ANS-binding properties. Furthermore, the nature of AP(5)A inhibition suggests also that the native forms of AK for binding with different substrates pre-exist in the absence of substrates. In the present study, a kinetics approach was used to explore the native forms distinguished by ANS-binding properties and by the nature of AP(5)A inhibition. The results revealed that the native forms distinguished by ANS probe are two conformational sub-ensembles. Both sub-ensembles are active and consist of a series of forms, which pre-exist in solution and can bind with different substrates. The K(m) values of N(1) for AMP, ADP and MgATP are larger than that of N(2), indicating that the N(2) sub-ensemble is more specific for binding substrates. This is consistent with the previous observation that the activity of N(2) is about 1.8-fold of that of N(1).

Distance-scaled, finite ideal-gas reference state improves structure-derived potentials of mean force for structure selection and stability prediction

The distance-dependent structure-derived potentials developed so far all employed a reference state that can be characterized as a residue (atom)-averaged state. Here, we establish a new reference state called the distance-scaled, finite ideal-gas reference (DFIRE) state. The reference state is used to construct a residue-specific all-atom potential of mean force from a database of 1011 nonhomologous (less than 30% homology) protein structures with resolution less than 2 A. The new all-atom potential recognizes more native proteins from 32 multiple decoy sets, and raises an average Z-score by 1.4 units more than two previously developed, residue-specific, all-atom knowledge-based potentials. When only backbone and C(beta) atoms are used in scoring, the performance of the DFIRE-based potential, although is worse than that of the all-atom version, is comparable to those of the previously developed potentials on the all-atom level. In addition, the DFIRE-based all-atom potential provides the most accurate prediction of the stabilities of 895 mutants among three knowledge-based all-atom potentials. Comparison with several physical-based potentials is made.

Measurement of Ca2+-binding constants of proteins and presentation of the CaLigator software

The complexity of Ca2+ cell signaling is dependent on a plethoria of Ca2+-binding proteins that respond to signals in different ranges of Ca2+ concentrations. Since the function of these proteins is directly coupled to their Ca2+-binding properties, there is a need for accurately determined equilibrium Ca2+-binding constants. In this work we outline the experimental techniques available to determine Ca2+-binding constants in proteins, derive the models used to describe the binding, and present CaLigator, software for least-square fitting directly to the measured quantity. The use of the software is illustrated for Ca2+-binding data obtained for two deamidated forms of calbindin D(9k), either an isospartate-56 (beta form) or a normal Asp-56 (alpha form). Here, the Ca2+-binding properties of the two isoforms have been studied using the chelator method. The alpha form shows similar Ca2+-binding properties to the wild type while the beta form has lost both cooperativety and affinity.