Search PubMed⌕ Search

Biomedical subjects

Martin Billeter

Publications and source records attributed to Martin Billeter.

16 recordsLinked to original sources

Structural characterization of an unusually stable cyclic peptide, kalata B2 from Oldenlandia affinis.

Kalata peptides are isolated from an African medicinal plant, Oldenlandia affinis, an aqueous decoction of which can be ingested to accelerate uterine contraction during childbirth. The closely packed disulfide core of kalata peptides confers unusual stability against thermal, chemical, and enzymatic degradation. The molecular arrangement may hamper NMR-assisted disulfide connectivity assignment. We have combined NMR with high-resolution mass spectrometry (MS) and MS/MS of native and chemically derivatized kalata B2 to determine its amino acid sequence and disulfide connectivity. Infrared multiphoton dissociation establishes the disulfide bond linkages in kalata B2 as I-IV, II-V and III-VI.

Amino Acid Sequence↗

Robust and versatile interpretation of spectra with coupled evolution periods using multi-way decomposition.

Coupling of evolution periods in NMR experiments is a very promising approach that has recently attracted much attention for its substantial savings in measurement time. Another novel concept, which has already proved useful in many types of NMR applications is multi-way decomposition. The PRODECOMP tool described here addresses the combination of these two modern tools in protein NMR. The following properties of this approach are described and illustrated. Highly similar information to what would be found in a corresponding full-dimensional spectrum is extracted from a set of projected spectra. All experimental spectra are used simultaneously, avoiding sensitivity loss associated with individual examination of the spectra. Aliasing caused by the linear combinations of individual shifts in the projected spectra is automatically resolved, allowing for better resolution due to smaller spectral widths. Reconstructions of various high-dimensional spectra, including the corresponding full-dimensional spectrum become straightforward. Spectral overlap is efficiently resolved. The capabilities of PRODECOMP are illustrated for a 14 kD protein, for which 12 projections of a 5-dimensional spectrum with the nuclei N, HN, CO, Calpha and Halpha are analysed.

Algorithms↗

Multiway decomposition of NMR spectra with coupled evolution periods.

Coupling evolution periods in NMR experiments on proteins has recently attracted much attention for its substantial savings in measurement time. Using the concept of multiway decomposition, which already proved useful in many types of NMR applications, the novel tool PRODECOMP decomposes sets of spectra with coupled evolution periods and attempts to provide all information that would be found in a corresponding full-dimensional spectrum. It offers the following advantages: (1) all experimental spectra are used simultaneously, avoiding sensitivity loss associated with individual examination of the spectra, (2) folding (aliasing) caused by the linear combinations of individual shifts in the projected spectra is automatically resolved, allowing for better resolution due to smaller spectral widths, (3) various types of schemes for the coupling of evolution periods are possible, (4) reconstructions of various high-dimensional spectra, including the corresponding full-dimensional spectrum, are straightforward, (5) overlap in both the projected and the directly detected dimensions is efficiently resolved. The capabilities of PRODECOMP are illustrated for a 14 kD protein, for which 12 projections of a five-dimensional spectrum with the nuclei N, HN, CO, Calpha, and Halpha are analyzed.

Journal Article↗

Accuracy and robustness of three-way decomposition applied to NMR data.

Three-way decomposition is a very versatile analysis tool with applications in a variety of protein NMR fields. It has been used to extract structural data from 3D NOESYs, to determine relaxation rates in large proteins, to identify ligand binding in screening for lead compounds, and to complement non-uniformly recorded (sparse) spectra. All applications so far concerned experimental data sets; it thus remains to address questions of accuracy and robustness of the method using simulated data where the correct answer is known. Systematic tests are presented for relaxation and NOESY data sets. Mixtures of real and synthetic data are used to allow control of various parameters and comparisons with correct reference data, while working with input that is as realistic as possible. The influence of the following parameters is evaluated: signal-to-noise, overlap of signals and the use of a regularization procedure within the algorithm. The main criteria used for the evaluation are accuracy and precision. It is shown that deterioration of accuracy is indicated by internal checks such as decrease of precision. Both with relaxation data and when interpreting NOESY spectra, three-way decomposition exhibits a robust behavior in situations with severe signal overlap and/or poor signal-to-noise, e.g., by avoiding false positives in the NOE shapes of NOESY decompositions. As a complement to this study, three-way decomposition is compared to other methods that achieve the same type of results.

Algorithms↗

Signal identification in NMR spectra with coupled evolution periods.

Novel multidimensional NMR experiments rely on modified time-domain sampling schemes to provide significant savings of experimental time. Several approaches are based on the coupling of evolution times resulting in a reduction of the dimensionality of the recorded spectra, and a concomitant saving of experimental time. We present a consistent and general tool, called EVOCOUP, for the analysis of these reduced dimensionality spectra. The approach is flexible in the sense that the input can consist of various forms of reduced dimensionality spectra, that any piece of information can be removed (provided enough information is left), e.g., signals undetectable due to poor signal-to-noise or covered by artifacts, and that it can be applied to spectra involving any number of nuclei. The use of a general optimization procedure and an appropriate target function provides for a robust approach with well-defined results and ensures optimal use of redundant information normally present in the input. Spectral overlap in the directly detected dimension is resolved in a fully automated manner, avoiding the assessment of signal quality and its use in combinatorial trials. The positions of all peaks in a corresponding full-dimensional spectrum are obtained without need for reconstruction of this spectrum. In a systematic analysis of a complete spectrum recorded for the 14 kDa protein azurin and involving five different nuclei, only four spin systems were missed and no false spins systems were detected.

Artifacts↗

Specific DNA recognition by the Antp homeodomain: MD simulations of specific and nonspecific complexes.

Four molecular dynamics simulation trajectories of complexes between the wild-type or a mutant Antennapedia homeodomain and 2 DNA sequences were generated in order to probe the mechanisms governing the specificity of DNA recognition. The starting point was published affinity measurements showing that a single protein mutation combined with a replacement of 2 base pairs yields a new high-affinity complex, whereas the other combinations, with changes on only 1 macromolecule, exhibited lower affinity. The simulations of the 4 complexes yielded fluctuating networks of interaction. On average, these networks differ significantly, explaining the switch of affinity caused by the alterations in the macromolecules. The network of mostly hydrogen-bonding interactions involving several water molecules, which was suggested both by X-ray and NMR structures of the wild-type homeodomain and its DNA operator sequence, could be reproduced in the trajectory. More interestingly, the high-affinity complex with alterations in both the protein and the DNA yielded again a dynamic but very tight network of intermolecular interactions, however, attributing a significantly stronger role to direct hydrophobic interactions at the expense of water bridges. The other 2 homeodomain-DNA complexes, with only 1 molecule altered, show on average over the trajectories a clearly reduced number of protein-DNA interactions. The observations from these simulations suggest specific experiments and thus close the circle formed by biochemical, structural, and computational studies. The shift from a water-dominated to a more "dry" interface may prove important in the design of proteins binding DNA in a specific manner.

Antennapedia Homeodomain Protein↗

Gated electron transfers and electron pathways in azurin: a NMR dynamic study at multiple fields and temperatures.

Dynamic properties of electron transfer pathways in a small blue copper cupredoxin are explored using an extensive 15N NMR relaxation study of reduced Pseudomonas aeruginosa azurin at four magnetic fields (500-900 MHz) and at two temperatures chosen well below the melting point of the protein. Following a careful model-free analysis, several protein regions with different dynamic regimes are identified. Nanosecond time-scale mobility characterizes various residues of the hydrophobic surface patch believed to mark the natural entry point for electrons, notably the surface-exposed copper-ligand His117. These findings are consistent with a gated electron transfer process according to the "dynamic docking" model. Residues 47-49 along intramolecular pathways of electrons show rigidity that is remarkably conserved when increasing the temperature. Three different conformational exchange processes were observed in the millisecond range, one near the only disulfide bridge in the molecule and two near the copper ion. The latter two processes are consistent with previous data such as crystal structures at various pH values and NMR relaxation dispersion experiments; they may indicate an additional gated electron transfer mechanism at slower time-scales.

Azurin↗

DNA adopts normal B-form upon incorporation of highly fluorescent DNA base analogue tC: NMR structure and UV-Vis spectroscopy characterization.

The influence of the highly fluorescent tricyclic cytosine base analogue (tC) on duplex DNA conformation is investigated. The duplex properties are characterized by absorbance and circular dichroism (CD) for all combinations of neighbouring bases to tC, and an NMR structure is determined for one tC-containing sequence. For the oligonucleotides with one tC incorporated instead of cytosine, the melting temperature is increased on average by 2.7 degrees C above that for the unmodified ones. CD spectra are practically identical for modified and unmodified sequences, indicating an unperturbed B-DNA conformation. The NMR structure determination of the self-complementary sequence 5'-CTC(tC)ACGTGGAG shows a DNA conformation consistent with B-form for the whole duplex. The root-mean-square distance for the nucleotides of the eight central base pairs between the 10 structures with lowest CYANA target functions and a mean structure is 0.45 +/- 0.17 A. The NMR data confirm correct base pairing for tC by the observation of both intrastrand and interstrand imino proton NOEs. Altogether, this suggests that tC works well as a cytosine analogue, i.e. it is situated in the base stack, forming hydrogen bonds with G in the complementary strand, without distorting the DNA backbone conformation. This first example of an artificial, highly fluorescent DNA base that does not perturb the DNA conformation could have valuable applications for the study of the structure and dynamics of nucleic acid systems.

Circular Dichroism↗

Accurate relaxation parameters for large proteins.

The practical applicability, performance, and robustness of three-way decomposition (TWD) for the extraction of relaxation parameters are demonstrated for a large protein with 370 residues, the maltose binding protein. An ordinary set of seven relaxation-modulated (15)N HSQC spectra, recorded at another site, is systematically analyzed. For all 341 assigned backbone amide groups, including 21 pairs and one group of three overlapped peaks, T1 decay values were determined. On isolated peaks, TWD extracts T1 values with systematically lower error bounds compared to conventional tools, although for these simple cases the improvements remain limited. However, in the presence of spectral artifacts, the decrease in errors can become significant, demonstrating the higher robustness of TWD. For about half of the peaks in overlapped regions, the decomposition allowed separation of the signals, yielding significantly different T1 values between overlapping signals. For the rest, similarity of the decay times for the two or three overlapping signals could be confirmed within usually low error bounds. The use of TWD thus leads to a significant increase in the number of accessible relaxation probes in large proteins. With a newly implemented graphical user interface, the application of TWD requires merely a peak list, and thus no additional effort compared to conventional approaches is needed.

Algorithms↗

Fully automated sequence-specific resonance assignments of hetero- nuclear protein spectra.

Full automation of the analysis of spectra is a prerequisite for high-throughput NMR studies in structural or functional genomics. Sequence-specific assignments often form the major bottleneck. Here, we present a procedure that yields nearly complete backbone and side chain resonance assignments starting from a set of heteronuclear three-dimensional spectra. Neither manual intervention, e.g., to correct lists obtained from peak picking before feeding these to an assignment program, nor protein-specific information, e.g., structures of homologous proteins, were required. By combining two earlier published procedures, AUTOPSY [Koradi et al. (1998) J. Magn. Reson., 135, 288-297] and GARANT [Bartels et al. (1996) J. Biomol. NMR, 7, 207-213], with a new program, PICS, all necessary steps from spectra analyses to sequence-specific assignments were performed fully automatically. Characteristic features of the present approach are a flexible design allowing as input almost any combination of NMR spectra, applicability to side chains, robustness with respect to parameter choices (such as noise levels) and reproducibility. In this study, automated resonance assignments were obtained for the 14 kD blue copper protein azurin from P. aeruginosa using five spectra: HNCACB, HNHA, HCCH-TOCSY, 15N-NOESY-HSQC and 13C-NOESY-HSQC. Peaks from these three-dimensional spectra were filtered and calibrated with the help of two two-dimensional spectra: 15N-HSQC and 13C-HSQC. The rate of incorrect assignments is less than 1.5% for backbone nuclei and about 3.5% when side chain protons are also considered.

Azurin↗

Optimizing resolution in multidimensional NMR by three-way decomposition.

Resolution depends on the number of points sampled in a FID; in indirectly detected dimensions it is an important determinant of the total experiment time. Based on the high redundancy present in NMR data, we propose the following timesaving scheme for three-dimensional spectra. An extensive grid of discrete t1- and t2-values is used, which increases resolution while preserving the spectral width. Total experiment time is reduced by avoiding the recording of t3-FIDs for selected pairs of t1 and t2; typically the recording is omitted for about 75% of the (t1,t2) combinations. These data sets are referred to as sparse, and post-experimental processing making optimal use of spectral redundancy provides the missing, non-recorded data. We have previously shown that three-way decomposition (TWD) within the MUNIN approach provides a practical way to process dense NMR data sets. Here, a novel TWD algorithm [Ibraghimov, (2002) Numer. Linear Algebra Appl. 9, 551-565] is used to complement a sparsely recorded time-domain data set by providing the missing FIDs for all (t1,t2) combinations omitted in the experiment. A necessary condition is that for each t1-value at least a few FIDs are recorded, and similar for each t2-value. The method is demonstrated on non-uniformly sampled (15)N-NOESY-HSQC data sets recorded for the 14 kD protein azurin. The spectra obtained by TWD, reconstruction and ordinary transform to frequency-domain are, in spite of the large number of signals and the high dynamic range typical for NOESYs, highly similar to a corresponding reference spectrum, for which all (t1,t2) combinations were recorded.

Magnetic Resonance Spectroscopy↗

NMR detection of multiple transitions to low-populated states in azurin.

Transitions to conformational states with very low populations were detected for the reduced blue copper protein azurin from Pseudomonas aeruginosa by applying constant relaxation time CPMG measurements to the backbone (15)N nuclei at three magnetic fields (11.7, 14.1, and 18.8 T) and three temperatures (25.7, 35.4, and 44.8 degrees C). Two exchange processes with different rate constants could be discriminated despite populations of the excited states below 1% and spatial neighborhood of the two processes. The group of (15)N nuclei involved in the faster process exhibits at 44.8 degrees C a forward rate constant of 11.7+/-2.4 s(-1) and a population of the exited state of 0.39+/-0.07%. They surround the aromatic ring of histidine 35 whose protonation state is coupled to the flipping of a neighboring peptide plane. For the slower process, the forward rate constant and population of the exited state at 44.8 degrees C are 4.1+/-0.1 s(-1) and 0.45+/-0.02%, respectively. The residues involved cluster nearby the copper ion, which is separated from the protonation site of histidine 35 by about 8 A, indicating conformational rearrangements involving the copper coordinating loops. The dependence of the equilibrium constant on the temperature is consistent with an enthalpy-dominated transition around the copper, but an entropy-controlled transition near histidine 35. The detection by nuclear magnetic resonance of millisecond to second conformational transitions near the copper ion suggests a low energy-cost rearrangement of the copper-binding site that may be necessary for efficient electron transfer.

Azurin↗

Automated analysis of large sets of heteronuclear correlation spectra in NMR-based drug discovery.

Drug discovery procedures based on NMR typically require the analysis of thousands of NMR spectra. For example, in "SAR by NMR", two-dimensional NMR spectra are recorded for a target protein mixed with ligand candidates from a comprehensive library of small molecules and are compared to the corresponding spectrum for the protein alone. We present an automated procedure for the comparative analysis of large sets of heteronuclear single quantum coherence spectra, which is based on three-way decomposition and implemented as the software package MUNIN. In a single step, spectra with differences in the peak positions (indicating ligand binding) and the affected peaks are identified. By omission of peak picking, ad hoc scoring of the quality of doubtful peaks is avoided. The procedure has been tested on the bacterial ribonuclease barnase, with a protein concentration of only 50 microM, using several small molecules including the substrate analogue 3'-GMP. Sets of 51 spectra were processed simultaneously, and it is concluded that spectra with binding ligands can be unambiguously identified from much larger sets of spectra.

Bacterial Proteins↗

Dynamics-modulated biological activity of transforming growth factor beta3.

Transforming growth factor beta3 (TGF-beta3) is an important mediator of growth, maintenance, and repair processes in human cells. Internal dynamic properties have been derived from (15)N NMR relaxation data and mapped onto the spatial structure of TGF-beta3. The pattern of internal dynamics in the structure identifies potential "hot spots" of binding free energy and reveals the importance of conformational entropy in the interaction of TGF-beta3 with the receptors. The observed internal dynamics set TGF-beta3 apart from other TGF-beta isoforms, with which it shares the same fold. These findings may explain functional differences among the various TGF-beta isoforms and thus prove essential in the search for related therapeutic agents.

Amino Acid Sequence↗

Essential domain motions in barnase revealed by MD simulations.

The wealth of data accumulated on the bacterial ribonuclease barnase is complemented by molecular dynamics trajectories starting from four different experimental structures and covering a total of >10 ns. Using principal component analysis, the simulations are interpreted in view of dynamic domains and hinges promoting relative motions of these domains. Two domains with residues 7-22 and 52-108 for the first domain and residues 25-51 for the second domain were consistently observed. Hinge regions consist primarily of Tyr24, Ser50, Ile51, and Gly52. Earlier mutation studies have demonstrated that the residues of the hinge regions play essential roles for the stability and activity of barnase. The domain motions are correlated to inter-domain interactions involving functionally important active site residues, such as Lys27 and Glu73. A model is presented that combines the observation of dynamic domains and their motions with the extensive mutation data from the literature. Enthalpic energy contributions originating from specific inter-domain interactions as well as entropic energy contributions due to the domain motions are discussed in the frame of this model and compared with destabilization energies measured for corresponding mutants.

Bacterial Proteins↗

Three-way decomposition of a complete 3D 15N-NOESY-HSQC.

Three-way decomposition is applied for the structural analysis of a complete three-dimensional (15)N-NOESY-HSQC of the 128 residues long protein azurin. The procedure presented includes decomposition using the software MUNIN, providing an initial characterization of the complete spectrum by 355 components. This is followed by post-processing yielding a final list of 149 components, 123 of which characterize 1859 NOE peaks from backbone N-H groups. Components from three-way decomposition are defined as direct products of one-dimensional shapes along the three dimensions. Thus, a complete set of distance constraints from this spectrum is obtained by one-dimensional peak picking of the shapes along the NOE dimension. Correctness and completeness of this set of NOEs are tested for all backbone amide groups against both an independent peak picking algorithm and the three-dimensional crystal structure of azurin, and a coincidence of about 95% is observed. Automated 'demixing' of components that are 'mixed' in a complex manner due to overlap of the HN and/or (15)N frequencies is illustrated.

Azurin↗