Experimentally exploring the conformational space sampled by domain reorientation in calmodulin

An Approximate Method in Using Molecular Mechanics Simulations To Study Slow Protein Conformational Changes

The broad range of characteristic motions in proteins has limited the applicability of molecular dynamics simulations in studying large-scale conformational transitions. We present an approximate method, making use of standard MD simulations and using a much larger integration time step, to obtain the structural changes for slow systematic motions of large complex systems. We show the applicability of this method by simulating the open to closed Calmodulin calcium binding domain conformational changes. Starting with the Ca2+-bound X-ray structure, and after the removal of the Ca2+ ions, our calculation yielded intermediate conformations during the rearrangement of helices in each Ca2+ binding pocket, leading to a structure with a lowest rmsd of 1.56 A compared to the NMR apo-calmodulin structure.

Long-Range-Distance NMR Effects in a Protein Labeled with a Lanthanide–DOTA Chelate

A two-thiol reactive lanthanide-DOTA (1,4,7,10-tetraazacyclododecane-N,N',N'',N'''-tetraacetic acid) chelate, CLaNP-3 (CLaNP=caged lanthanide NMR probe), was synthesized for the rigid attachment to cysteine groups on a protein surface, and used to obtain long-range-distance information from the {15N,1H} HSQC spectra of the protein-lanthanide complex. The DOTA ring exhibits several isomers that are in exchange; however, single resonances were observed for most amide groups in the protein, allowing determination of a single, apparent magnetic-susceptibility tensor. Pseudocontact shifts caused by Yb-containing CLaNP-3 were observed for atoms at 15-35 A from the metal. By using Gd-containing CLaNP-3, relaxation effects were observed, allowing distances up to 30 A from the paramagnetic center to be determined accurately. Similar results were obtained with a Gd-DTPA (diethylene-triaminepentaacetic acid) chelate, CLaNP-1, bound in the same bidentate manner to the protein. This study demonstrates that bidentate attachment of a paramagnetic probe enables determination of long-range distances.

Observation of microsecond time-scale protein dynamics in the presence of Ln3+ ions: application to the N-terminal domain of cardiac troponin C

The microsecond time-scale motions in the N-terminal domain of cardiac troponin C (NcTnC) loaded with lanthanide ions have been investigated by means of a (1)H(N) off-resonance spin-lock experiment. The observed relaxation dispersion effects strongly increase along the series of NcTnC samples containing La(3+), Ce(3+), and Pr(3+) ions. This rise in dispersion effects is due to modulation of long-range pseudocontact shifts by micros time-scale dynamics. Specifically, the motion in the coordination sphere of the lanthanide ion (i.e. in the NcTnC EF-hand motif) causes modulation of the paramagnetic susceptibility tensor which, in turn, causes modulation of pseudocontact shifts. It is also probable that opening/closing dynamics, previously identified in Ca(2+)-NcTnC, contributes to some of the observed dispersions. On the other hand, it is unlikely that monomer-dimer exchange in the solution of NcTnC is directly responsible for the dispersion effects. Finally, on-off exchange of the lanthanide ion does not seem to play any significant role. The amplification of dispersion effects by Ln(3+) ions is a potentially useful tool for studies of micros-ms motions in proteins. This approach makes it possible to observe the dispersions even when the local environment of the reporting spin does not change. This happens, for example, when the motion involves a 'rigid' structural unit such as individual alpha-helix. Even more significantly, the dispersions based on pseudocontact shifts offer better chances for structural characterization of the dynamic species. This method can be generalized for a large class of applications via the use of specially designed lanthanide-binding tags.

Two new chiral EDTA-based metal chelates for weak alignment of proteins in solution

[structure: see text] A short synthesis of EDTA-based metal chelates that can be attached to the cysteine residue of a protein via a disulfide bond is described. The complexes were used after coordination of lanthanides to align trigger factor and apo-calmodulin in solution to yield residual dipolar couplings and pseudocontact shifts. Alignment tensors for the new tags are linearly independent compared to those of previously published tags.

Structural basis for the observed differential magnetic anisotropic tensorial values in calcium binding proteins

Lanthanide ions (Ln(3+)), which have ionic radii similar to those of Ca(2+), can displace the latter in a calcium binding protein, without affecting its tertiary structure. The paramagnetic Ln(3+) possesses large anisotropic magnetic susceptibilities and produce pseudocontact shifts (PCSs), which have r(-3) dependence. The PCS can be seen for spins as far as 45 A from the paramagnetic ion. They aid in structure refinement of proteins by providing long-range distance constraints. Besides, they can be used to determine the interdomain orientation in multidomain proteins. This is particularly important in the context of a calcium binding protein from Entamoeba histolytica (EhCaBP), which consists of two globular domains connected by a flexible linker region containing 8 residues. As a first step to obtain the interdomain orientation in EhCaBP, a suite of 2D and 3D heteronuclear experiments were recorded on EhCaBP by displacing calcium with Ce(3+), Ho(3+), Er(3+), Tm(3+), Dy(3+), and Yb(3+) ions in separate experiments, and the PCS of (1)H(N) and (15)N spins were measured. Such data have been used in the refinement of the individual domain structures of the protein in parallel with the calculation of the respective magnetic anisotropy tensorial values, which differ substantially (2.1-2.8 times) from what is found in other Ca(2+) binding loops. This study provides a structural basis for such variations in the magnetic anisotropy tensorial values.

EF-hand protein dynamics and evolution of calcium signal transduction: an NMR view

Calcium signaling, one of the most widespread signaling mechanisms in cells, is generally carried out by EF-hand proteins, characterized by a helix-loop-helix motif paired in functional domains. EF-hand proteins may be viewed as molecular switches activated by calcium concentration transients. The EF-hand structural database has grown to a point where meaningful inferences on the functional conformational rearrangements upon calcium binding can be made by comparing a fair number of pairs of end points, i.e., the structures of the apo and calcium-bound forms. More compact descriptors of the movement associated with calcium binding, in terms of principal component analysis of the six interhelical angles, have also become available. Dynamic information obtained by NMR, also with the aid of calcium substitution with paramagnetic lanthanides, is shedding light on the intrinsic amplitude of the conformational degrees of freedom sampled by the various members of the EF-hand superfamily, as well as on the time scales of the motions. Particularly, NMR of lanthanide derivatives helps in capturing long time scale motions. Both static and dynamic pictures reveal a large variety of behaviors. It is increasingly recognized that the EF-hand machinery has differentiated its behavior during evolution in several ways, e.g., by modifying one of the loops, by undergoing a further duplication after the initial motif duplication that originated the functional domain, or by acquiring the ability to dimerize.

The Ig Doublet Z1Z2: A Model System for the Hybrid Analysis of Conformational Dynamics in Ig Tandems from Titin

Titin is a gigantic elastic filament that determines sarcomere ultrastructure and stretch response in vertebrate muscle. It folds into numerous Ig and FnIII domains connected in tandem. Data on interdomain arrangements and dynamics are essential for understanding the function of this filament. Here, we report a mechanistic analysis of the conformational dynamics of two Ig domains from the N terminus of titin, Z1Z2, by using X-ray crystallography, SAXS, NMR relaxation data, and residual dipolar couplings in combination. Z1Z2 preferentially adopts semiextended conformations in solution, with close-hinge arrangements representing low-probability states. Although interdomain contacts are not observed, the linker appears to acquire moderate rigidity via small, local hydrophobic interactions. Thus, Z1Z2 constitutes an adaptable modular system with restricted dynamics. We speculate that its preexistent conformation contributes to the selective recruitment of the binding partner telethonin onto the repetitive surface of the filament. The structural interconversion of four Z1Z2 conformers is analyzed.

A structural and dynamic characterization of the EF-hand protein CLSP

The structure and dynamics of human calmodulin-like skin protein (CLSP) have been characterized by NMR spectroscopy. The mobility of CLSP has been found to be different for the N-terminal and C-terminal domains. The isolated domains were also expressed and analyzed. The structure of the isolated C-terminal domain is presented. The N-terminal domain is characterized by four stable helices, which experience large fluctuations. This is shown to be due to mutations in the hydrophobic core. The overall N-terminal domain behavior is similar both in the full-length protein and in the isolated domain. By exploiting the capability of Tb3+ bound to CLSP to induce partial orientation of the molecule in a magnetic field, restricted motion of one domain with respect to the other was proved. By using NMR, ITC, and ESI-MS, the calcium and magnesium binding properties were investigated. Finally, CLSP is framed into the evolutionary scheme of the calmodulin-like family.

Paramagnetic tagging of diamagnetic proteins for solution NMR

In this article, approaches towards the paramagnetic tagging of diamagnetic proteins are reviewed. Alignment can be achieved by adding paramagnetic fusion proteins or peptides to the C- or the N-terminus or by attaching paramagnetic tags to Cystein residues. Applications for the study of homodimer structures and protein/ligand interactions, as well as protein domain dynamics, are reviewed.

Analysis of interdomain dynamics in a two-domain protein using residual dipolar couplings together with 15N relaxation data

In this paper, we propose the idea that simultaneous analysis of NMR relaxation data and residual dipolar couplings (RDCs) can provide information about interdomain dynamics in a multidomain protein, which cannot be derived from each data set separately. Specifically, such an approach can be useful when the interdomain motions occur on a timescale comparable to or slower than the overall tumbling in solution. We analyze residual dipolar couplings together with 15N relaxation data for Lys48-linked di-ubiquitin (Ub2), in which interdomain dynamics are described as interconversion between two distinct conformational states of the protein. Our results show that 15N relaxation and residual dipolar coupling data can be used as two complementary experimental data sets for consistent characterization of interdomain conformations and dynamics in this dual-domain protein.

Inter-domain orientation and motions in VAT-N explored by residual dipolar couplings and 15N backbone relaxation

The N-terminal domain of VAT (Valosine-containing protein-like ATPase of Thermoplasma acidophilum), VAT-N (20.5 kDa), is considered to be the primary substrate-recognition site of the complex. The solution structure of VAT-N derived in our laboratory using conventionally obtained NMR restraints shows the existence of two equally sized sub-domains, VAT-Nn and VAT-Nc, together forming a kidney-shaped overall structure. The putative substrate-binding site of VAT-N involves free loops and a highly charged groove located on the surface of the protein. Alternatively, the opening of the cleft between the domains to accommodate substrate has been proposed to be part of the functional mechanism. We have used the residual dipolar couplings (RDCs) obtained in a bicelle medium to refine the structure of VAT-N. The long-range information available from RDCs both defines the sub-domain orientation and probes possible inter-domain motions. In addition, 15N backbone relaxation data were obtained and analysed within the model-free framework. Together, the data provides a refined structure with improved local geometry, but with the overall kidney shape intact. Further, the protein is rigid overall, with no evidence of inter-domain motions.

Ensemble-based convergence analysis of biomolecular trajectories

Assessing the convergence of a biomolecular simulation is an essential part of any careful computational investigation, because many fundamental aspects of molecular behavior depend on the relative populations of different conformers. Here we present a physically intuitive method to self-consistently assess the convergence of trajectories generated by molecular dynamics and related methods. Our approach reports directly and systematically on the structural diversity of a simulation trajectory. Straightforward clustering and classification steps are the key ingredients, allowing the approach to be trivially applied to systems of any size. Our initial study on met-enkephalin strongly suggests that even fairly long trajectories (approximately 50 ns) may not be converged for this small--but highly flexible--system.

ARTIST: An activated method in internal coordinate space for sampling protein energy landscapes

We present the first applications of an activated method in internal coordinate space for sampling all-atom protein conformations, the activation-relaxation technique for internal coordinate space trajectories (ARTIST). This method differs from all previous internal coordinate-based studies aimed at folding or refining protein structures in that conformational changes result from identifying and crossing well-defined saddle points connecting energy minima. Our simulations of four model proteins containing between 4 and 47 amino acids indicate that this method is efficient for exploring conformational space in both sparsely and densely packed environments, and offers new perspectives for applications ranging from computer-aided drug design to supramolecular assembly.

Weak alignment NMR: a hawk-eyed view of biomolecular structure

Imposing a very slight deviation from the isotropic random distribution of macromolecules in solution in an NMR sample tube permits the measurement of residual internuclear dipolar couplings (RDCs). Such interactions are very sensitive functions of the time-averaged orientation of the corresponding internuclear vectors and thereby offer highly precise structural information. In recent years, advances have been made both in the technology to measure RDCs and in the computational procedures that integrate this information in the structure determination process. The exceptional precision with which RDCs can be measured under weakly aligned conditions is also starting to reveal the mostly, but not universally, subtle effects of internal protein dynamics. Importantly, RDCs potentially can reveal motions taking place on a timescale slower than rotational diffusion and analysis is uniquely sensitive to the direction of motion, not just its amplitude.

Unwinding the helical linker of calcium-loaded calmodulin: A molecular dynamics study

The fold of calmodulin (CaM) consists of two globular domains connected by a helical segment (the linker), whose conformational properties play a crucial role for the protein's molecular recognition processes. Here we investigate the structural properties of the linker by performing a 11.5 ns molecular dynamics (MD) simulation of calcium-loaded human CaM in aqueous solution. The calculations are based on the AMBER force field. The calculated S2 order parameters are in good accord with NMR data: The structure of the linker in our simulations is much more flexible than that emerging from the Homo sapiens X-ray structure, consistently with the helix unwinding observed experimentally in solution. This process occurs spontaneously in a nanosecond timescale, as observed also in a very recent simulation based on the GROMOS force field. A detailed description of the mechanism that determines the linker unwinding is provided, in which electrostatic contacts between the two globular domains play a critical role. The orientation of the domains emerging from our MD calculations is consistent both with former X-ray scattering data and a recent NMR work. Based on our findings, a rationale for the experimentally measured entropy cost associated to binding to the protein's cellular partners is also given.

NMR Spectroscopy of Paramagnetic Metalloproteins

This article deals with the solution structure determination of paramagnetic metalloproteins by NMR spectroscopy. These proteins were believed not to be suitable for NMR investigations for structure determination until a decade ago, but eventually novel experiments and software protocols were developed, with the aim of making the approach suitable for the goal and as user-friendly and safe as possible. In the article, we also give hints for the optimization of experiments with respect to each particular metal ion, with the aim of also providing a handy tool for nonspecialists. Finally, a section is dedicated to the significant progress made on 13C direct detection, which reduces the negative effects of paramagnetism and may constitute a new chapter in the whole field of NMR spectroscopy.

Probing metal-protein interactions using a de novo design approach

De novo design of metalloproteins provides a valuable tool for understanding the structural constraints and functional attributes of natural biological systems using first principles. This review focuses on recent research aimed primarily at probing the subtle interactions between metals and proteins in designed systems. Considerable attention has focussed on redefining novel design methods used in mimicking natural hemeproteins, mononuclear and dinuclear metallopeptides and functional biological electron-transfer proteins. The present results indicate that the field of metalloprotein design is contributing significantly to the understanding of metals in biology.

Biophysical characterization of ERp29. Evidence for a key structural role of cysteine 125

ERp29 is a major resident of the endoplasmic reticulum (ER) that seemingly plays an important role in most animal cells. Although a protein-folding association is widely supported, ERp29's specific molecular function remains unknown. A chaperone activity was postulated from evidence that ERp29 forms multimers like the classical ER chaperones, but conflicting results have emerged from our recent studies. Here a biophysical approach was used to clarify this issue and also reveal a key structural role for ERp29's characteristic cysteine, Cys-125. Applying hydrodynamic parameters derived from sedimentation and dynamic light-scattering analyses, a model of ERp29's quaternary structure was assembled from existing tertiary substructures. Comparison with Windbeutel, an ERp29-like protein from fruit fly with specialized chaperone activity, revealed similar tri-lobar gross structures but some finer differences consistent with functional divergence. Solubility and hydrophobic probe assays revealed moderate surface hydrophobicity, which was reduced in mutant ERp29 in which serine replaced Cys-125. This mutant was also relatively labile to proteolytic degradation, providing two reasons for the strict conservation of Cys-125. No multimerization was observed with untagged ERp29, which existed as tight homodimers (K(d) < 50 nm), whereas His-tagged ERp29 artifactually formed 670-kDa oligomers. These findings distinguish ERp29 biophysically from its peers in the ER including Windbeutel, endorsing our postulate that ERp29 adds a distinct type of folding activity to the ER machinery. By invoking novel functional associations for Cys-125 and the adjoining linker, new clues about how ERp29 might work have also arisen.