Functional dynamics in the active site of the ribonuclease binase

Protein secondary structure assignment using pc-polyline and convolutional neural network

MOTIVATION

The assignment of protein secondary structure elements (SSEs) underpins structural analysis and prediction. The backbone of a protein could be adequately represented using a pc-polyline that passes through the centers of its peptide planes. One salient feature of pc-polyline representation is that the secondary structure of a protein becomes recognizable in a matrix whose elements are the pairwise distances between two peptide plane centers. Thus, a pc-polyline could in turn be used to assign SSEs.

RESULTS

Using convolutional neural network (CNN) here we confirm that a pc-polyline indeed contains enough information for it to be used for the accurate assignments of the six SSE types: α-helix, β-sheet, β-bulge, 310 -helix, turn and loop. The applications to three large data sets show that the assignments by our CNN-based p2psse program agree very well with those by dssp, stride and quite well with those by five other programs. The analyses of their SSE assignments raise some general questions about the characterizations of protein secondary structure. In particular the analyses illustrate the difficulty with giving a quantitative and consistent definition for each of the six SSE types especially for 310 -helix, β-bulge, turn or loop in terms of either backbone H-bond patterns, or backbone dihedral angles, or Cα -polyline or pc-polyline. The difficulty suggests that the SSE space though being dominated by the regions for the six SSE types is to a certain degree continuous.

AVAILABILITY

The program is available at https://github.com/wlincong/p2pSSE.

Protein Dynamics revealed by NMR Relaxation Methods

Structural biology often focuses primarily on three-dimensional structures of biological macromolecules, deposited in the Protein Data Bank (PDB). This resource is a remarkable entity for the world-wide scientific and medical communities, as well as the general public, as it is a growing translation into three-dimensional space of the vast information in genomic databases, e.g. GENBANK. There is, however, significantly more to understanding biological function than the three-dimensional coordinate space for ground-state structures of biomolecules. The vast array of biomolecules experiences natural dynamics, interconversion between multiple conformational states, and molecular recognition and allosteric events that play out on timescales ranging from picoseconds to seconds. This wide range of timescales demands ingenious and sophisticated experimental tools to sample and interpret these motions, thus enabling clearer insight into functional annotation of the PDB. NMR spectroscopy is unique in its ability to sample this range of timescales at atomic resolution and in physiologically relevant conditions using spin relaxation methods. The field is constantly expanding to provide new creative experiments, to yield more detailed coverage of timescales, and to broaden the power of interpretation and analysis methods. This review highlights the current state of the methodology and examines the extension of analysis tools for more complex experiments and dynamic models. The future for understanding protein dynamics is bright, and these extended tools bring greater compatibility with developments in computational molecular dynamics, all of which will further our understanding of biological molecular functions. These facets place NMR as a key component in integrated structural biology.

Partial Metal Ion Saturation of C2 Domains Primes Synaptotagmin 1-Membrane Interactions

Synaptotagmin 1 (Syt1) is an integral membrane protein whose phospholipid-binding tandem C2 domains, C2A and C2B, act as Ca²⁺ sensors of neurotransmitter release. Our objective was to understand the role of individual metal-ion binding sites of these domains in the membrane association process. We used Pb²⁺, a structural and functional surrogate of Ca²⁺, to generate the protein states with well-defined protein-metal ion stoichiometry. NMR experiments revealed that binding of one divalent metal ion per C2 domain results in loss of conformational plasticity of the loop regions, potentially pre-organizing them for additional metal-ion and membrane-binding events. In C2A, a divalent metal ion in site 1 is sufficient to drive its weak association with phosphatidylserine-containing membranes, whereas in C2B, it enhances the interactions with the signaling lipid phosphatidylinositol-4,5-bisphosphate. In full-length Syt1, both Pb²⁺-complexed C2 domains associate with phosphatidylserine-containing membranes. Electron paramagnetic resonance experiments show that the extent of membrane insertion correlates with the occupancy of the C2 metal ion sites. Together, our results indicate that upon partial metal ion saturation of the intra-loop region, Syt1 adopts a dynamic, partially membrane-bound state. The properties of this state, such as conformationally restricted loop regions and positioning of C2 domains in close proximity to anionic lipid headgroups, "prime" Syt1 for cooperative binding of a full complement of metal ions and deeper membrane insertion.

Why the Energy Landscape of Barnase Is Hierarchical

We have used NMR and computational methods to characterize the dynamics of the ribonuclease barnase over a wide range of timescales in free and inhibitor-bound states. Using temperature- and denaturant-dependent measurements of chemical shift, we show that barnase undergoes frequent and highly populated hinge bending. Using relaxation dispersion, we characterize a slower and less populated motion with a rate of 750 ± 200 s⁻¹, involving residues around the lip of the active site, which occurs in both free and bound states and therefore suggests conformational selection. Normal mode calculations characterize correlated hinge bending motions on a very rapid timescale. These three measurements are combined with previous measurements and molecular dynamics calculations on barnase to characterize its dynamic landscape on timescales from picoseconds to milliseconds and length scales from 0.1 to 2.5 nm. We show that barnase has two different large-scale fluctuations: one on a timescale of 10⁻⁹−10⁻⁶ s that has no free energy barrier and is a hinge bending that is determined by the architecture of the protein; and one on a timescale of milliseconds (i.e., 750 s⁻¹) that has a significant free energy barrier and starts from a partially hinge-bent conformation. These two motions can be described as hierarchical, in that the more highly populated faster motion provides a platform for the slower (less probable) motion. The implications are discussed. The use of temperature and denaturant is suggested as a simple and general way to characterize motions on the intermediate ns-μs timescale.

Solution NMR views of dynamical ordering of biomacromolecules

BACKGROUND

To understand the mechanisms related to the 'dynamical ordering' of macromolecules and biological systems, it is crucial to monitor, in detail, molecular interactions and their dynamics across multiple timescales. Solution nuclear magnetic resonance (NMR) spectroscopy is an ideal tool that can investigate biophysical events at the atomic level, in near-physiological buffer solutions, or even inside cells.

SCOPE OF REVIEW

In the past several decades, progress in solution NMR has significantly contributed to the elucidation of three-dimensional structures, the understanding of conformational motions, and the underlying thermodynamic and kinetic properties of biomacromolecules. This review discusses recent methodological development of NMR, their applications and some of the remaining challenges.

MAJOR CONCLUSIONS

Although a major drawback of NMR is its difficulty in studying the dynamical ordering of larger biomolecular systems, current technologies have achieved considerable success in the structural analysis of substantially large proteins and biomolecular complexes over 1MDa and have characterised a wide range of timescales across which biomolecular motion exists. While NMR is well suited to obtain local structure information in detail, it contributes valuable and unique information within hybrid approaches that combine complementary methodologies, including solution scattering and microscopic techniques.

GENERAL SIGNIFICANCE

For living systems, the dynamic assembly and disassembly of macromolecular complexes is of utmost importance for cellular homeostasis and, if dysregulated, implied in human disease. It is thus instructive for the advancement of the study of the dynamical ordering to discuss the potential possibilities of solution NMR spectroscopy and its applications. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

Dynamic Response of the C2 Domain of Protein Kinase Cα to Ca2+ Binding

Ca²⁺-dependent conserved-region 2 (C2) domains target their host signaling proteins to anionic membranes. The Ca²⁺-binding event is a prerequisite for membrane association. Here, we investigate multiscale metal-ion-dependent dynamics of the C2 domain of protein kinase Cα (C2α) using NMR spectroscopy. Interactions with metal ions attenuate microsecond-timescale motions of the loop regions, indicating that preorganization of the metal-binding loops occurs before membrane insertion. Binding of a full complement of Ca²⁺ ions has a profound effect on the millisecond-timescale dynamics of the N- and C-terminal regions of C2α. We propose that Ca²⁺ binding allosterically destabilizes the terminal regions of C2α and thereby facilitates the conformational rearrangement necessary for full membrane insertion and activation of protein kinase Cα.

A critical assessment of methods to recover information from averaged data

Conformational heterogeneity is key to the function of many biomacromolecules, but only a few groups have tried to characterize it until recently. Now, thanks to the increased throughput of experimental data and the increased computational power, the problem of the characterization of protein structural variability has become more and more popular. Several groups have devoted their efforts in trying to create quantitative, reliable and accurate protocols for extracting such information from averaged data. We analyze here different approaches, discussing strengths and weaknesses of each. All approaches can roughly be clustered into two groups: those satisfying the maximum entropy principle and those recovering ensembles composed of a restricted number of molecular conformations. In the first case, the solution focuses on the features that are common to all the infinite solutions satisfying the experimental data; in the second case, the reconstructed ensemble shows the conformational regions where a large probability can be placed. The upper limits for conformational probabilities (MaxOcc) can also be calculated. We also give an overview of the mainstream experimental observables, with considerations on the assumptions underlying their usage.

Role of Proteome Physical Chemistry in Cell Behavior

We review how major cell behaviors, such as bacterial growth laws, are derived from the physical chemistry of the cell’s proteins. On one hand, cell actions depend on the individual biological functionalities of their many genes and proteins. On the other hand, the common physics among proteins can be as important as the unique biology that distinguishes them. For example, bacterial growth rates depend strongly on temperature. This dependence can be explained by the folding stabilities across a cell’s proteome. Such modeling explains how thermophilic and mesophilic organisms differ, and how oxidative damage of highly charged proteins can lead to unfolding and aggregation in aging cells. Cells have characteristic time scales. For example, E. coli can duplicate as fast as 2–3 times per hour. These time scales can be explained by protein dynamics (the rates of synthesis and degradation, folding, and diffusional transport). It rationalizes how bacterial growth is slowed down by added salt. In the same way that the behaviors of inanimate materials can be expressed in terms of the statistical distributions of atoms and molecules, some cell behaviors can be expressed in terms of distributions of protein properties, giving insights into the microscopic basis of growth laws in simple cells.

Using NMR spectroscopy to elucidate the role of molecular motions in enzyme function

Conformational motions play an essential role in enzyme function, often facilitating the formation of enzyme-substrate complexes and/or product release. Although considerable debate remains regarding the role of molecular motions in the conversion of enzymatic substrates to products, numerous examples have found motions to be crucial for optimization of enzyme scaffolds, effective substrate binding, and product dissociation. Conformational fluctuations are often rate-limiting to enzyme catalysis, primarily through product release, with the chemical reaction occurring much more quickly. As a result, the direct involvement of motions at various stages along the enzyme reaction coordinate remains largely unknown and untested. In the following review, we describe the use of solution NMR techniques designed to probe various timescales of molecular motions and detail examples in which motions play a role in propagating catalytic effects from the active site and directly participate in essential aspects of enzyme function.

[Identification of 2',3'-cGMP as an intermediate of RNA catalytic cleavage by binase and evaluation of its biological action]

Binase--Bacillus pumilus RNase--endonuclease cleaves the phosphodiester bond between the 3'-guanylic residue and the 5'-OH residue of adjacent nucleotides with the formation of corresponding intermediate 2',3'-cGMP. Subsequent hydrolysis of 2',3'-cGMP into 3'-phosphate is highly specific and occurs slowly So the question arises about the existing time of that positional isomer during RNA catalytic cleavage by binase and about 2',3'-cGMP role in antitumor activity of the enzyme: In present work by means of enzyme-linked immunosorbent assay we established that during catalytic cleavage of RNA by binase 2',3'-cGMP is preserved in reaction mixture for an hour, at the same time phosphodiesterases activation doesn't lead to the total elimination of 2',3'-cGMP. The highest amount of 2',3'-cGMP was observed under the pH 8.5, it reaches nanomolar concentration at initial RNA concentration of 100-1000 μg/mL. Exogenous 2',3'-cGMP, like its positional isomer 3',5'-cGMP, doesn't trigger an apoptosis of human lung adenocarcinoma A549 cells, which are sensitive to binase apoptogenic action. Taking into account data about binase internalization and activation of mitochondrial pores opening by 2',3'-cyclic guanosine phosphates we may consider that 2',3'-cGMP can contribute to the apoptosis initiated by binase only when 2',3'-cGMP is generated intracellularly.

Changing the topology of protein backbone: the effect of backbone cyclization on the structure and dynamics of a SH3 domain

Understanding of the effects of the backbone cyclization on the structure and dynamics of a protein is essential for using protein topology engineering to alter protein stability and function. Here we have determined, for the first time, the structure and dynamics of the linear and various circular constructs of the N-SH3 domain from protein c-Crk. These constructs differ in the length and amino acid composition of the cyclization region. The backbone cyclization was carried out using intein-mediated intramolecular chemical ligation between the juxtaposed N- and the C-termini. The structure and backbone dynamics studies were performed using solution NMR. Our data suggest that the backbone cyclization has little effect on the overall three-dimensional structure of the SH3 domain: besides the termini, only minor structural changes were found in the proximity of the cyclization region. In contrast to the structure, backbone dynamics are significantly affected by the cyclization. On the subnanosecond time scale, the backbone of all circular constructs on average appears more rigid than that of the linear SH3 domain; this effect is observed over the entire backbone and is not limited to the cyclization site. The backbone mobility of the circular constructs becomes less restricted with increasing length of the circularization loop. In addition, significant conformational exchange motions (on the sub-millisecond time scale) were found in the N-Src loop and in the adjacent β-strands in all circular constructs studied in this work. These effects of backbone cyclization on protein dynamics have potential implications for the stability of the protein fold and for ligand binding.

Chemical exchange in biomacromolecules: past, present, and future

The perspective reviews quantitative investigations of chemical exchange phenomena in proteins and other biological macromolecules using NMR spectroscopy, particularly relaxation dispersion methods. The emphasis is on techniques and applications that quantify the populations, interconversion kinetics, and structural features of sparsely populated conformational states in equilibrium with a highly populated ground state. Applications to folding, molecular recognition, catalysis, and allostery by proteins and nucleic acids are highlighted.

Measuring dynamic and kinetic information in the previously inaccessible supra-τ(c) window of nanoseconds to microseconds by solution NMR spectroscopy

Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful tool that has enabled experimentalists to characterize molecular dynamics and kinetics spanning a wide range of time-scales from picoseconds to days. This review focuses on addressing the previously inaccessible supra-tc window (defined as τ(c) < supra-τ(c) < 40 μs; in which tc is the overall tumbling time of a molecule) from the perspective of local inter-nuclear vector dynamics extracted from residual dipolar couplings (RDCs) and from the perspective of conformational exchange captured by relaxation dispersion measurements (RD). The goal of the first section is to present a detailed analysis of how to extract protein dynamics encoded in RDCs and how to relate this information to protein functionality within the previously inaccessible supra-τ(c) window. In the second section, the current state of the art for RD is analyzed, as well as the considerable progress toward pushing the sensitivity of RD further into the supra-τ(c) scale by up to a factor of two (motion up to 25 μs). From the data obtained with these techniques and methodology, the importance of the supra-τ(c) scale for protein function and molecular recognition is becoming increasingly clearer as the connection between motion on the supra-τ(c) scale and protein functionality from the experimental side is further strengthened with results from molecular dynamics simulations.

Protein dynamics by ¹⁵N nuclear magnetic relaxation

Nitrogen-15 relaxation is the most ubiquitous source of information about protein (backbone) dynamics used by NMR spectroscopists. It provides the general characteristics of hydrodynamics as well as internal motions on subnanosecond, micro- and millisecond timescales of a biomolecule. Here, we present a full protocol to perform and analyze a series of experiments to measure the (15)N longitudinal relaxation rate, the (15)N transverse relaxation rate under an echo train or a single echo, the (15)N-(1)H dipolar cross-relaxation rate, as well as the longitudinal and transverse cross-relaxation rates due to the cross-correlation of the nitrogen-15 chemical shift anisotropy and the dipolar coupling with the adjacent proton. These rates can be employed to carry out model-free analyses and can be used to quantify accurately the contribution of chemical exchange to transverse relaxation.

Temperature and solvent dependence of the dynamical landscape of tau protein conformations

We report the variation with temperature of the ensemble distribution of conformations spanned by the tau protein in its dynamical states measured by small-angle X-ray scattering (SAXS) using synchrotron radiation. The SAXS data show a clear temperature variation of the distribution of occupied protein conformations from 293 to 318 K. More conformations with a smaller radius of gyration are occupied at higher temperature. The protein-solvent interactions are shown by computer simulation to be essential for controlling the dynamics of protein conformations, providing evidence for the key role of water solvent in the protein dynamics, as proposed by Giorgio Careri.

Comparison of NMR and crystal structures highlights conformational isomerism in protein active sites

The JCSG has recently developed a protocol for systematic comparisons of high-quality crystal and NMR structures of proteins. In this paper, the extent to which this approach can provide function-related information on the two functionally annotated proteins TM1081, a Thermotoga maritima anti-σ factor antagonist, and A2LD1 (gi:13879369), a mouse γ-glutamylamine cyclotransferase, is explored. The NMR structures of the two proteins have been determined in solution at 313 and 298 K, respectively, using the current JCSG protocol based on the software package UNIO for extensive automation. The corresponding crystal structures were solved by the JCSG at 100 K and 1.6 Å resolution and at 100 K and 1.9 Å resolution, respectively. The NMR and crystal structures of the two proteins share the same overall molecular architectures. However, the precision of the structure determination along the amino-acid sequence varies over a significantly wider range in the NMR structures than in the crystal structures. Thereby, in each of the two NMR structures about 65% of the residues have displacements below the average and in both proteins the less well ordered residues include large parts of the active sites, in addition to some highly solvent-exposed surface areas. Whereas the latter show increased disorder in the crystal and in solution, the active-site regions display increased displacements only in the NMR structures, where they undergo local conformational exchange on the millisecond time scale that appears to be frozen in the crystals. These observations suggest that a search for molecular regions showing increased structural disorder and slow dynamic processes in solution while being well ordered in the corresponding crystal structure might be a valid initial step in the challenge of identifying putative active sites in functionally unannotated proteins with known three-dimensional structure.

The EVB as a quantitative tool for formulating simulations and analyzing biological and chemical reactions

Recent years have seen dramatic improvements in computer power, allowing ever more challenging problems to be approached. In light of this, it is imperative to have a quantitative model for examining chemical reactivity, both in the condensed phase and in solution, as well as to accurately quantify physical organic chemistry (particularly as experimental approaches can often be inconclusive). Similarly, computational approaches allow for great progress in studying enzyme catalysis, as they allow for the separation of the relevant energy contributions to catalysis. Due to the complexity of the problems that need addressing, there is a need for an approach that can combine reliability with an ability to capture complex systems in order to resolve long-standing controversies in a unique way. Herein, we will demonstrate that the empirical valence bond (EVB) approach provides a powerful way to connect the classical concepts of physical organic chemistry to the actual energies of enzymatic reactions by means of computation. Additionally, we will discuss the proliferation of this approach, as well as attempts to capture its basic chemistry and repackage it under different names. We believe that the EVB approach is the most powerful tool that is currently available for studies of chemical processes in the condensed phase in general and enzymes in particular, particularly when trying to explore the different proposals about the origin of the catalytic power of enzymes.

Selective characterization of microsecond motions in proteins by NMR relaxation

The three-dimensional structures of macromolecules fluctuate over a wide range of time-scales. Separating the individual dynamic processes according to frequency is of importance in relating protein motions to biological function and stability. We present here a general NMR method for the specific characterization of microsecond motions at backbone positions in proteins even in the presence of other dynamics such as large-amplitude nanosecond motions and millisecond chemical exchange processes. The method is based on measurement of relaxation rates of four bilinear coherences and relies on the ability of strong continuous radio frequency fields to quench millisecond chemical exchange. The utility of the methodology is demonstrated and validated through two specific examples focusing on the thermo-stable proteins, ubiquitin and protein L, where it is found that small-amplitude microsecond dynamics are more pervasive than previously thought. Specifically, these motions are localized to alpha helices, loop regions, and regions along the rim of beta sheets in both of the proteins examined. A third example focuses on a 28 kDa ternary complex of the chaperone Chz1 and the histones H2A.Z/H2B, where it is established that pervasive microsecond motions are localized to a region of the chaperone that is important for stabilizing the complex. It is further shown that these motions can be well separated from extensive millisecond dynamics that are also present and that derive from exchange of Chz1 between bound and free states. The methodology is straightforward to implement, and data recorded at only a single static magnetic field are required.

Global and local mobility of apocalmodulin monitored through fast-field cycling relaxometry

Calmodulin (CaM) is a ubiquitous eukaryotic protein with two conformationally independent domains that can bind up to two calcium ions each. In the calcium-bound state, CaM is able to regulate a vast number of cellular activities by binding to a multiplicity of target proteins in different modes. Its versatility has been ascribed to its anomalously high flexibility. The calcium-free form (apoCaM), which is the resting state of CaM in cells, is also able to functionally bind a number of protein targets, but its dynamics has received less attention. At variance with the calcium-bound form, the crystal structure of apoCaM shows a compact organization of the two domains, but NMR measurements could not detect any contact between them, thus indicating the presence of mobility in solution. The mobility of apoCaM is here investigated through protein proton relaxation rate measurements performed with a high-sensitivity fast-field cycling relaxometer. Such measurements provide direct access to the spectral density function and show that 1), the reorientation time is in agreement with a closed form of the protein; but 2), the collective order parameter is much smaller than for other well folded compact proteins, indicating that a remarkably large side-chain mobility must be present.

Protein proton-proton dynamics from amide proton spin flip rates

Residue-specific amide proton spin-flip rates K were measured for peptide-free and peptide-bound calmodulin. K approximates the sum of NOE build-up rates between the amide proton and all other protons. This work outlines the theory of multi-proton relaxation, cross relaxation and cross correlation, and how to approximate it with a simple model based on a variable number of equidistant protons. This model is used to extract the sums of K-rates from the experimental data. Error in K is estimated using bootstrap methodology. We define a parameter Q as the ratio of experimental K-rates to theoretical K-rates, where the theoretical K-rates are computed from atomic coordinates. Q is 1 in the case of no local motion, but decreases to values as low as 0.5 with increasing domination of sidechain protons of the same residue to the amide proton flips. This establishes Q as a monotonous measure of local dynamics of the proton network surrounding the amide protons. The method is applied to the study of proton dynamics in Ca(2+)-saturated calmodulin, both free in solution and bound to smMLCK peptide. The mean Q is 0.81 +/- 0.02 for free calmodulin and 0.88 +/- 0.02 for peptide-bound calmodulin. This novel methodology thus reveals the presence of significant interproton disorder in this protein, while the increase in Q indicates rigidification of the proton network upon peptide binding, confirming the known high entropic cost of this process.

Accessing ns-micros side chain dynamics in ubiquitin with methyl RDCs

This study presents the first application of the model-free analysis (MFA) (Meiler in J Am Chem Soc 123:6098-6107, 2001; Lakomek in J Biomol NMR 34:101-115, 2006) to methyl group RDCs measured in 13 different alignment media in order to describe their supra-tau (c) dynamics in ubiquitin. Our results indicate that methyl groups vary from rigid to very mobile with good correlation to residue type, distance to backbone and solvent exposure, and that considerable additional dynamics are effective at rates slower than the correlation time tau (c). In fact, the average amplitude of motion expressed in terms of order parameters S (2) associated with the supra-tau (c) window brings evidence to the existence of fluctuations contributing as much additional mobility as those already present in the faster ps-ns time scale measured from relaxation data. Comparison to previous results on ubiquitin demonstrates that the RDC-derived order parameters are dominated both by rotameric interconversions and faster libration-type motions around equilibrium positions. They match best with those derived from a combined J-coupling and residual dipolar coupling approach (Chou in J Am Chem Soc 125:8959-8966, 2003) taking backbone motion into account. In order to appreciate the dynamic scale of side chains over the entire protein, the methyl group order parameters are compared to existing dynamic ensembles of ubiquitin. Of those recently published, the broadest one, namely the EROS ensemble (Lange in Science 320:1471-1475, 2008), fits the collection of methyl group order parameters presented here best. Last, we used the MFA-derived averaged spherical harmonics to perform highly-parameterized rotameric searches of the side chains conformation and find expanded rotamer distributions with excellent fit to our data. These rotamer distributions suggest the presence of concerted motions along the side chains.

Approximate reconstruction of continuous spatially complex domain motions by multialignment NMR residual dipolar couplings

NMR spectroscopy is one of the most powerful techniques for studying the internal dynamics of biomolecules. Current formalisms approximate the dynamics using simple continuous motional models or models involving discrete jumps between a small number of states. However, no approach currently exists for interpreting NMR data in terms of continuous spatially complex motional paths that may feature more than one distinct maneuver. Here, we present an approach for approximately reconstructing spatially complex continuous motions of chiral domains using NMR anisotropic interactions. The key is to express Wigner matrix elements, which can be determined experimentally using residual dipolar couplings, as a line integral over a curve in configuration space containing an ensemble of conformations and to approximate the curve using a series of geodesic segments. Using this approach and five sets of synthetic residual dipolar couplings computed for five linearly independent alignment conditions, we show that it is theoretically possible to reconstruct salient features of a multisegment interhelical motional trajectory obtained from a 65 ns molecular dynamics simulation of a stem-loop RNA. Our study shows that the 3-D atomic reconstruction of complex motions in biomolecules is within experimental reach.

Role of loop-loop interactions in coordinating motions and enzymatic function in triosephosphate isomerase

The enzyme triosephosphate isomerase (TIM) has been used as a model system for understanding the relationship between protein sequence, structure, and biological function. The sequence of the active site loop (loop 6) in TIM is directly correlated with a conserved motif in loop 7. Replacement of loop 7 of chicken TIM with the corresponding loop 7 sequence from an archaeal homologue caused a 10(2)-fold loss in enzymatic activity, a decrease in substrate binding affinity, and a decrease in thermal stability. Isotope exchange studies performed by one-dimensional (1)H NMR showed that the substrate-derived proton in the enzyme is more susceptible to solvent exchange for DHAP formation in the loop 7 mutant than for WT TIM. TROSY-Hahn Echo and TROSY-selected R(1rho) experiments indicate that upon mutation of loop 7, the chemical exchange rate for active site loop motion is nearly doubled and that the coordinated motion of loop 6 is reduced relative to that of the WT. Temperature dependent NMR experiments show differing activation energies for the N- and C-terminal hinges in this mutant enzyme. Together, these data suggest that interactions between loop 6 and loop 7 are necessary to provide the proper chemical context for the enzymatic reaction to occur and that the interactions play a significant role in modulating the chemical dynamics near the active site.

Residual dipolar couplings as a tool to study molecular recognition of ubiquitin

RDCs (residual dipolar couplings) in NMR spectroscopy provide information about protein dynamics complementary to NMR relaxation methods, especially in the previously inaccessible time window between the protein correlation time tau(c) and 50 micros. For ubiquitin, new modes of motion of the protein backbone could be detected using RDC-based techniques. An ensemble of ubiquitin based on these RDC values is found to comprise all different conformations that ubiquitin adopts upon binding to different recognition proteins. These conformations in protein-protein complexes had been derived from 46 X-ray structures. Thus, for ubiquitin recognition by other proteins, conformational selection rather than induced fit seems to be the dominant mechanism.

Resonant response to temperature modulation for enzymatic dynamics characterization

We consider enzymes involved in a three-state Michaelis-Menten kinetics and submitted to well-chosen temperature modulations of small amplitude. From the first-order amplitudes of concentration oscillations, we design a response function that is maximum for targeted values of the chemical relaxation times. This resonant function can be used to screen a large set of enzymes and identify the one governed by the desired kinetics. The method gives access to all the dynamical parameters of the targeted enzyme without resorting to a fit. We show how to estimate the precision of this parameter determination and give some hints for experimental validation.

Structural biology by NMR: structure, dynamics, and interactions

The function of bio-macromolecules is determined by both their 3D structure and conformational dynamics. These molecules are inherently flexible systems displaying a broad range of dynamics on time-scales from picoseconds to seconds. Nuclear Magnetic Resonance (NMR) spectroscopy has emerged as the method of choice for studying both protein structure and dynamics in solution. Typically, NMR experiments are sensitive both to structural features and to dynamics, and hence the measured data contain information on both. Despite major progress in both experimental approaches and computational methods, obtaining a consistent view of structure and dynamics from experimental NMR data remains a challenge. Molecular dynamics simulations have emerged as an indispensable tool in the analysis of NMR data.

Eta(z)/kappa: a transverse relaxation optimized spectroscopy NMR experiment measuring longitudinal relaxation interference

NMR spin relaxation experiments provide a powerful tool for the measurement of global and local biomolecular rotational dynamics at subnanosecond time scales. Technical limitations restrict most spin relaxation studies to biomolecules weighing less than 10 kDa, considerably smaller than the average protein molecular weight of 30 kDa. In particular, experiments measuring eta(z), the longitudinal (1)H(N)-(15)N dipole-dipole (DD)/(15)N chemical shift anisotropy (CSA) cross-correlated relaxation rate, are among those least suitable for use with larger biosystems. This is unfortunate because these experiments yield valuable insight into the variability of the (15)N CSA tensor over the polypeptide backbone, and this knowledge is critical to the correct interpretation of most (15)N-NMR backbone relaxation experiments, including R(2) and R(1). In order to remedy this situation, we present a new (1)H(N)-(15)N transverse relaxation optimized spectroscopy experiment measuring eta(z) suitable for applications with larger proteins (up to at least 30 kDa). The presented experiment also yields kappa, the site-specific rate of longitudinal (1)H(N)-(1)H(') DD cross relaxation. We describe the eta(z)/kappa experiment's performance in protonated human ubiquitin at 30.0 degrees C and in protonated calcium-saturated calmodulin/peptide complex at 20.0 degrees C, and demonstrate preliminary experimental results for a deuterated E. coli DnaK ATPase domain construct at 34 degrees C.

Barnase as a new therapeutic agent triggering apoptosis in human cancer cells

Background

RNases are currently studied as non-mutagenic alternatives to the harmful DNA-damaging anticancer drugs commonly used in clinical practice. Many mammalian RNases are not potent toxins due to the strong inhibition by ribonuclease inhibitor (RI) presented in the cytoplasm of mammalian cells.

Methodology/Principal Findings

In search of new effective anticancer RNases we studied the effects of barnase, a ribonuclease from Bacillus amyloliquefaciens, on human cancer cells. We found that barnase is resistant to RI. In MTT cell viability assay, barnase was cytotoxic to human carcinoma cell lines with half-inhibitory concentrations (IC₅₀) ranging from 0.2 to 13 µM and to leukemia cell lines with IC₅₀ values ranging from 2.4 to 82 µM. Also, we characterized the cytotoxic effects of barnase-based immunoRNase scFv 4D5-dibarnase, which consists of two barnase molecules serially fused to the single-chain variable fragment (scFv) of humanized antibody 4D5 that recognizes the extracellular domain of cancer marker HER2. The scFv 4D5-dibarnase specifically bound to HER2-positive cells and was internalized via receptor-mediated endocytosis. The intracellular localization of internalized scFv 4D5-dibarnase was determined by electronic microscopy. The cytotoxic effect of scFv 4D5-dibarnase on HER2-positive human ovarian carcinoma SKOV-3 cells (IC₅₀ = 1.8 nM) was three orders of magnitude greater than that of barnase alone. Both barnase and scFv 4D5-dibarnase induced apoptosis in SKOV-3 cells accompanied by internucleosomal chromatin fragmentation, membrane blebbing, the appearance of phosphatidylserine on the outer leaflet of the plasma membrane, and the activation of caspase-3.

Conclusions/Significance

These results demonstrate that barnase is a potent toxic agent for targeting to cancer cells.

Gel-based oligonucleotide microarray approach to analyze protein-ssDNA binding specificity

Gel-based oligonucleotide microarray approach was developed for quantitative profiling of binding affinity of a protein to single-stranded DNA (ssDNA). To demonstrate additional capabilities of this method, we analyzed the binding specificity of ribonuclease (RNase) binase from Bacillus intermedius (EC 3.1.27.3) to ssDNA using generic hexamer oligodeoxyribonucleotide microchip. Single-stranded octamer oligonucleotides were immobilized within 3D hemispherical gel pads. The octanucleotides in individual pads 5'-{N}N(1)N(2)N(3)N(4)N(5)N(6){N}-3' consisted of a fixed hexamer motif N(1)N(2)N(3)N(4)N(5)N(6) in the middle and variable parts {N} at the ends, where {N} represent A, C, G and T in equal proportions. The chip has 4096 pads with a complete set of hexamer sequences. The affinity was determined by measuring dissociation of the RNase-ssDNA complexes with the temperature increasing from 0 degrees C to 50 degrees C in quasi-equilibrium conditions. RNase binase showed the highest sequence-specificity of binding to motifs 5'-NNG(A/T/C)GNN-3' with the order of preference: GAG > GTG > GCG. High specificity towards G(A/T/C)G triplets was also confirmed by measuring fluorescent anisotropy of complexes of binase with selected oligodeoxyribonucleotides in solution. The affinity of RNase binase to other 3-nt sequences was also ranked. These results demonstrate the applicability of the method and provide the ground for further investigations of nonenzymatic functions of RNases.

Structure and dynamics of pin1 during catalysis by NMR

The link between internal enzyme motions and catalysis is poorly understood. Correlated motions in the microsecond-to-millisecond timescale may be critical for enzyme function. We have characterized the backbone dynamics of the peptidylprolyl isomerase (Pin1) catalytic domain in the free state and during catalysis. Pin1 is a prolyl isomerase of the parvulin family and specifically catalyzes the isomerization of phosphorylated Ser/Thr-Pro peptide bonds. Pin1 has been shown to be essential for cell-cycle progression and to interact with the neuronal tau protein inhibiting its aggregation into fibrillar tangles as found in Alzheimer's disease. (15)N relaxation dispersion measurements performed on Pin1 during catalysis reveal conformational exchange processes in the microsecond timescale. A subset of active site residues undergo kinetically similar exchange processes even in the absence of a substrate, suggesting that this area is already "primed" for catalysis. Furthermore, structural data of the turning-over enzyme were obtained through inter- and intramolecular nuclear Overhauser enhancements. This analysis together with a characterization of the substrate concentration dependence of the conformational exchange allowed the distinguishing of regions of the enzyme active site that are affected primarily by substrate binding versus substrate isomerization. Together these data suggest a model for the reaction trajectory of Pin1 catalysis.

Transition state theory can be used in studies of enzyme catalysis: lessons from simulations of tunnelling and dynamical effects in lipoxygenase and other systems

The idea that enzyme catalysis involves special factors such as coherent fluctuations, quantum mechanical tunnelling and non-equilibrium solvation (NES) effects has gained popularity in recent years. It has also been suggested that transition state theory (TST) cannot be used in studies of enzyme catalysis. The present work uses reliable state of the art simulation approaches to examine the above ideas. We start by demonstrating that we are able to simulate any of the present catalytic proposals using the empirical valence bond (EVB) potential energy surfaces, the dispersed polaron model and the quantized classical path (QCP) approach, as well as the approximate vibronic method. These approaches do not treat the catalytic effects by phenomenological treatments and thus can be considered as first principles approaches (at least their ability to compare enzymatic reaction to the corresponding solution reactions). This work will consider the lipoxygenase reaction, and to lesser extent other enzymes, for specific demonstration. It will be pointed out that our study of the lipoxygenase reaction reproduces the very large observed isotope effect and the observed rate constant while obtaining no catalytic contribution from nuclear quantum mechanical (NQM) effects. Furthermore, it will be clarified that our studies established that the NQM effect decreases rather than increases when the donor-acceptor distance is compressed. The consequences of these findings in terms of the temperature dependence of the kinetic isotope effect and in terms of different catalytic proposals will be discussed. This paper will also consider briefly the dynamical effects and conclude that such effects do not contribute in a significant way to enzyme catalysis. Furthermore, it will be pointed out that, in contrast to recent suggestions, NES effects are not dynamical effects and should therefore be part of the activation free energy rather than the transmission factor. In view of findings of the present work and our earlier works, it seems that TST provides a quantitative tool for studies of enzyme catalysis and that the key open questions are related to the nature of the factors that lead to transition state stabilization.

Faithful estimation of dynamics parameters from CPMG relaxation dispersion measurements

This work examines the robustness of fitting of parameters describing conformational exchange (k(ex), p(a/b), and Deltaomega) processes from CPMG relaxation dispersion data. We have analyzed the equations describing conformational exchange processes for the intrinsic inter-dependence of their parameters that leads to the existence of multiple equivalent solutions, which equally satisfy the experimental data. We have used Monte-Carlo simulations and fitting to the synthetic data sets as well as the direct 3-D mapping of the parameter space of k(ex), p(a/b), and Deltaomega to quantitatively assess the degree of the parameter inter-dependence. The demonstrated high correlation between parameters can preclude accurate dynamics parameter estimation from NMR spin-relaxation data obtained at a single static magnetic field. The strong parameter inter-dependence can readily be overcome through acquisition of spin-relaxation data at more than one static magnetic field thereby allowing accurate assessment of conformational exchange properties.

Characterizing a Partially Folded Intermediate of the Villin Headpiece Domain Under Non-denaturing Conditions: Contribution of His41 to the pH-dependent Stability of the N-terminal Subdomain

The contribution of interactions involving the imidazole ring of His41 to the pH-dependent stability of the villin headpiece (HP67) N-terminal subdomain has been investigated by nuclear magnetic resonance (NMR) spin relaxation. NMR-derived backbone N-H order parameters (S2) for wild-type (WT) HP67 and H41Y HP67 indicate that reduced conformational flexibility of the N-terminal subdomain in WT HP67 is due to intramolecular interactions with the His41 imidazole ring. These interactions, together with desolvation effects, contribute to significantly depress the pKa of the buried imidazole ring in the native state. 15N R1rho relaxation dispersion data indicate that WT HP67 populates a partially folded intermediate state that is 10.9 kJ mol(-1) higher in free energy than the native state under non-denaturing conditions at neutral pH. The partially folded intermediate is characterized as having an unfolded N-terminal subdomain while the C-terminal subdomain retains a native-like fold. Although the majority of the residues in the N-terminal subdomain sample a random-coil distribution of conformations, deviations of backbone amide 1H and 15N chemical shifts from canonical random-coil values for residues within 5A of the His41 imidazole ring indicate that a significant degree of residual structure is maintained in the partially folded ensemble. The pH-dependence of exchange broadening is consistent with a linear three-state exchange model whereby unfolding of the N-terminal subdomain is coupled to titration of His41 in the partially folded intermediate with a pKa,I=5.69+/-0.07. Although maintenance of residual interactions with the imidazole ring in the unfolded N-terminal subdomain appears to reduce pKa,I compared to model histidine compounds, protonation of His41 disrupts these interactions and reduces the difference in free energy between the native state and partially folded intermediate under acidic conditions. In addition, chemical shift changes for residues Lys70-Phe76 in the C-terminal subdomain suggest that the HP67 actin binding site is disrupted upon unfolding of the N-terminal subdomain, providing a potential mechanism for regulating the villin-dependent bundling of actin filaments.

PhosphoThr peptide binding globally rigidifies much of the FHA domain from Arabidopsis receptor kinase-associated protein phosphatase

A net increase in the backbone rigidity of the kinase-interacting FHA domain (KI-FHA) from the Arabidopsis receptor kinase-associated protein phosphatase (KAPP) accompanies the binding of a phosphoThr peptide from its CLV1 receptor-like kinase partner, according to (15)N NMR relaxation at 11.7 and 14.1 T. All of the loops of free KI-FHA display evidence of nanosecond-scale motions. Many of these same residues have residual dipolar couplings that deviate from structural predictions. Binding of the CLV1 pT868 peptide seems to reduce nanosecond-scale fluctuations of all loops, including half of the residues of recognition loops. Residues important for affinity are found to be rigid, i.e., conserved residues and residues of the subsite for the key pT+3 peptide position. This behavior parallels SH2 and PTB domain recognition of pTyr peptides. PhosphoThr peptide binding increases KI-FHA backbone rigidity (S(2)) of three recognition loops, a loop nearby, seven strands from the beta-sandwich, and a distal loop. Compensating the trend of increased rigidity, binding enhances fast mobility at a few sites in four loops on the periphery of the recognition surface and in two loops on the far side of the beta-sandwich. Line broadening evidence of microsecond- to millisecond-scale fluctuations occurs across the six-stranded beta-sheet and nearby edges of the beta-sandwich; this forms a network connected by packing of interior side chains and H-bonding. A patch of the slowly fluctuating residues coincides with the site of segment-swapped dimerization in crystals of the FHA domain of human Chfr. Phosphopeptide binding introduces microsecond- to millisecond-scale fluctuations to more residues of the long 8/9 recognition loop of KI-FHA. The rigidity of this FHA domain appears to couple as a whole to pThr peptide binding.

Staphylococcus aureus Sortase A transpeptidase. Calcium promotes sorting signal binding by altering the mobility and structure of an active site loop

Many virulence factors in gram-positive bacteria are covalently anchored to the cell-wall peptidoglycan by sortase enzymes, a group of widely distributed cysteine transpeptidases. The Staphylococcus aureus Sortase A protein (SrtA) is the archetypal member of the Sortase family and is activated by Ca2+, an adaptation that may facilitate host colonization as elevated concentrations of this ion are encountered in human tissue. Here we show that a single Ca2+ ion bound to an ordered pocket on SrtA allosterically activates catalysis by modulating both the structure and dynamics of a large active site loop. Detailed nitrogen-15 relaxation measurements indicate that Ca2+ may facilitate the adaptive recognition of the substrate by inducing slow micro- to millisecond time-scale dynamics in the active site. Interestingly, relaxation compensated Carr-Purcell-Meiboom-Gill experiments suggest that the time scale of these motions is directly correlated with ion binding. The results of site-directed mutagenesis indicate that this motional coupling is mediated by the side chain of Glu-171, which is positioned within the beta6/beta7 loop and shown to contribute to Ca2+ binding. The available structural and dynamics data are compatible with a loop closure model of Ca2+ activation, in which the beta6/beta7 loop fluctuates between a binding competent closed form that is stabilized by Ca2+, and an open, highly flexible state that removes key substrate contacting residues from the active site.

Estimating the accuracy of protein structures using residual dipolar couplings

It has been commonly recognized that residual dipolar coupling data provide a measure of quality for protein structures. To quantify this observation, a database of 100 single-domain proteins has been compiled where each protein was represented by two independently solved structures. Backbone 1H-15N dipolar couplings were simulated for the target structures and then fitted to the model structures. The fits were characterized by an R-factor which was corrected for the effects of non-uniform distribution of dipolar vectors on a unit sphere. The analyses show that favorable R values virtually guarantee high accuracy of the model structure (where accuracy is defined as the backbone coordinate rms deviation). On the other hand, unfavorable R values do not necessarily suggest low accuracy. Based on the simulated data, a simple empirical formula is proposed to estimate the accuracy of protein structures. The method is illustrated with a number of examples, including PDZ2 domain of human phosphatase hPTP1E.

Improvement of duty-cycle heating compensation in NMR spin relaxation experiments

To reliably measure NMR relaxation properties of macromolecules is a prerequisite for precise experiments that identify subtle variations in relaxation rates, as required for the determination of rotational diffusion anisotropy, CSA tensor determination, advanced motional modeling or entropy difference estimations. An underlying problem with current NMR relaxation measurement protocols is maintaining constant sample temperature throughout the execution of the relaxation series especially when rapid data acquisition is required. Here, it is proposed to use a combination of a heating compensation and a proton saturation sequence at the beginning of the NMR relaxation pulse scheme. This simple extension allows reproducible, robust and rapid acquisition of NMR spin relaxation data sets. The method is verified with (15)N spin relaxation measurements for human ubiquitin.

A new amide proton R1rho experiment permits accurate characterization of microsecond time-scale conformational exchange

A new off-resonance spin-lock experiment to record relaxation dispersion profiles of amide protons is presented. The sensitivity-enhanced HSQC-type sequence is designed to minimize the interference from cross-relaxation effects and ensure that the dispersion profiles in the absence of micros-ms time-scale dynamics are flat. Toward this end (i) the proton background is eliminated by sample deuteration (Ishima et al., 1998), (ii) 1H spin lock is applied to two-spin modes 2(H(x)Sin theta + H(z)Cos theta) N(z), and (iii) the tilt angle theta approximately 35 degrees is maintained throughout the series of measurements (Desvaux et al. Mol. Phys., 86 (1995) 1059). The relaxation dispersion profiles recorded in this manner sample a wide range of effective rf field strengths (up to and in excess of 20 kHz) which makes them particularly suitable for studies of motions on the time scale < or = 100 micros. The new experiment has been tested on the Ca2+-loaded regulatory domain of cardiac troponin C. Many residues show pronounced dispersions with remarkably similar correlation times of approximately 30 micros. Furthermore, these residues are localized in the regions that have been previously implicated in conformational changes (Spyracopoulos et al. Biochemistry, 36 (1997) 12138).