Molecular statistics of cytochrome <Emphasis Type="Italic">c</Emphasis>: structural plasticity and molecular environment

Nuclear magnetic resonance experiments performed on yeast mitochondrial cytochrome c (Cytc), a paradigmatic electron transfer protein, reveal that the two oxidation states have similar structures, but different mobility: despite the few structural differences compared with the reduced form, the oxidized form displays a larger unfolding propensity. Molecular dynamics simulations performed on both NMR reduced and NMR oxidized forms show that the reduced form has a larger solvent-accessible surface area (SASA). Starting from this observation, a molecular statistical approach was then applied in order to correlate the molecular surface to molecular mobility. Simulations started from biased initial conditions corresponding to different molecular sizes were combined with the maximal constrained entropy method. The NMR structure of oxidized Cytc is more suited to expose a smaller SASA than the NMR structure of the reduced form, but the accessible conformational landscape at 300 K around the NMR oxidized structure is flatter than for the NMR reduced structure. Protein configurations of smaller SASA and size display larger plasticity when they resemble the NMR oxidized structure, whereas they are more rigid when they resemble the NMR reduced structure. Implications of the results for the protein properties during its functional process are discussed. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Dehydrogenases reducing NAD+ in the presence of D (-) 3- phosphoglycerate in mycobacteria (BCG, H37Ra, M. phlei)]

Referring to the elution volume on a Sephadex G-150 column only one specific peak is obtained, the same for the BCG, H37Ra and Mycobacterium phlei strains grown on Sauton synthetic medium. Some properties of these partially purified dehydrogenases are studied (conservation and dialysis in media of different salt concentrations, equilibrium constant, Km, heat stability). All enzyme preparations from tubercle bacilli (BCG, H37Ra) are readily inactivated by heat and are very unstable in solutions of low ionic strength. In contrast, under the same experimental conditions, all enzyme samples from M. phlei are, comparatively, much more stable towards these factors [heat, salt (potassium phosphate, NaCl) concentration].     

Microscale NMR

NMR spectroscopy is increasingly being used to characterize microliter and smaller-volume samples. Substances at picomole levels have been identified using NMR spectrometers equipped with microcoil-based probes. NMR probes that incorporate multiple sample chambers enable higher-throughput NMR experiments. Hyphenation of capillary-scale separations and microcoil NMR has also decreased analysis time of mixtures. For example, capillary isotachophoresis/NMR allows the highest mass sensitivity nanoliter-volume flow cells to be used with low microliter volume samples because isotachophoresis concentrates the microliter volume sample into the nanoliter volume NMR detection probe. In addition, the diagnostic capabilities of NMR spectroscopy allow the physico-chemical aspects of a capillary separation process to be characterized on-line. Because of such advances, the application of NMR to smaller samples continues to grow.     

Sensitivity enhancement for membrane proteins reconstituted in parallel and perpendicular oriented bicelles obtained by using repetitive cross-polarization and membrane-incorporated free radicals

Multidimensional separated local-field and spin-exchange experiments employed by oriented-sample solid-state NMR are essential for structure determination and spectroscopic assignment of membrane proteins reconstituted in macroscopically aligned lipid bilayers. However, these experiments typically require a large number of scans in order to establish interspin correlations. Here we have shown that a combination of optimized repetitive cross polarization (REP-CP) and membrane-embedded free radicals allows one to enhance the signal-to-noise ratio by factors 2.4-3.0 in the case of Pf1 coat protein reconstituted in magnetically aligned bicelles with their normals being either parallel or perpendicular to the main magnetic field. Notably, spectral resolution is not affected at the 2:1 radical-to-protein ratio. Spectroscopic assignment of Pf1 coat protein in the parallel bicelles has been established as an illustration of the method. The proposed methodology will advance applications of oriented-sample NMR technique when applied to samples containing smaller quantities of proteins and three-dimensional experiments.     

High-throughput tissue extraction protocol for NMR- and MS-based metabolomics

In metabolomics, tissues typically are extracted by grinding in liquid nitrogen followed by the stepwise addition of solvents. This is time-consuming and difficult to automate, and the multiple steps can introduce variability. Here we optimize tissue extraction methods compatible with high-throughput, reproducible nuclear magnetic resonance (NMR) spectroscopy- and mass spectrometry (MS)-based metabolomics. Previously, we concluded that methanol/chloroform/water extraction is preferable for metabolomics, and we further optimized this here using fish liver and an automated Precellys 24 bead-based homogenizer, allowing rapid extraction of multiple samples without carryover. We compared three solvent addition strategies: stepwise, two-step, and all solvents simultaneously. Then we evaluated strategies for improved partitioning of metabolites between solvent phases, including the addition of extra water and different partition times. Polar extracts were analyzed by NMR and principal components analysis, and the two-step approach was preferable based on lipid partitioning, reproducibility, yield, and throughput. Longer partitioning or extra water increased yield and decreased lipids in the polar phase but caused metabolic decay in these extracts. Overall, we conclude that the two-step method with extra water provides good quality data but that the two-step method with 10 min partitioning provides a more accurate snapshot of the metabolome. Finally, when validating the two-step strategy using NMR and MS metabolomics, we showed that technical variability was considerably smaller than biological variability.     

Signal Intensities Derived from Different NMR Probes and Parameters Contribute to Variations in Quantification of Metabolites

We discovered that serious issues could arise that may complicate interpretation of metabolomic data when identical samples are analyzed at more than one NMR facility, or using slightly different NMR parameters on the same instrument. This is important because cross-center validation metabolomics studies are essential for the reliable application of metabolomics to clinical biomarker discovery. To test the reproducibility of quantified metabolite data at multiple sites, technical replicates of urine samples were assayed by 1D-¹H-NMR at the University of Alberta and the University of Michigan. Urine samples were obtained from healthy controls under a standard operating procedure for collection and processing. Subsequent analysis using standard statistical techniques revealed that quantitative data across sites can be achieved, but also that previously unrecognized NMR parameter differences can dramatically and widely perturb results. We present here a confirmed validation of NMR analysis at two sites, and report the range and magnitude that common NMR parameters involved in solvent suppression can have on quantitated metabolomics data. Specifically, saturation power levels greatly influenced peak height intensities in a frequency-dependent manner for a number of metabolites, which markedly impacted the quantification of metabolites. We also investigated other NMR parameters to determine their effects on further quantitative accuracy and precision. Collectively, these findings highlight the importance of and need for consistent use of NMR parameter settings within and across centers in order to generate reliable, reproducible quantified NMR metabolomics data.     

Vibrational averaging of chemical shift anisotropies in model peptides

The effects of chemical shift anisotropy (CSA) are evident in line-shapes or side-band analysis in solid-state NMR, in the observed line positions in partially oriented samples, and in relaxation effects in liquid-state studies. In all of these cases, the effective shielding tensor is influenced by fast vibrational averaging in addition to larger-amplitude internal motions and to overall libration or rotation. Here we compute the contributions of vibrational averaging (including zero-point motions) to the CSA relaxation strengths for the nitrogen and carbonyl carbon in two simple peptide models, and for snapshots taken from a path-integral simulation of a small protein. Because the ¹⁵N shielding tensor is determined by all the atoms of the peptide group, it is less influenced by vibrational motion than (for example) the N–H dipolar interaction, which is more sensitive to the motion of the light hydrogen atom. Computed order parameters for CSA averaging are hence much closer to unity than are N–H dipolar order parameters. This leads to a reduction by about 9% in the magnitude of the amide nitrogen CSA that is needed to fit liquid-state relaxation data. Similar considerations apply to the carbonyl carbon shielding tensor, but in this case the differences between dipolar and CSA averaging are smaller. These considerations will be important for making comparisons between CSA tensors extracted from various NMR experiments, and for comparisons to quantum chemical calculations carried out on static conformers.     

Simple algorithms for the estimation of the initial number of individuals in commingled skeletal remains

The determination of the number of individuals represented within commingled remains is based on two types of estimators, those assessing the minimum number of individuals and those assessing the most likely number of individuals. Much as the latter produce improved results, they still exhibit significant drawbacks, which are related to the misidentification of the number of pairs between the existing bilateral elements. This article addresses these problems through the use of two computer algorithms. One algorithm produces a number of potential pairs between bilateral elements and the other estimates the number of individuals in a commingled sample by incorporating the percentages of lost and altered bones into the analysis. These algorithms were validated using hypothetical and actual skeletal samples, and are more effective in comparison to any conventional estimators, particularly in cases, where the elements are poorly preserved.     

High-resolution proton and carbon-13 NMR of membranes: why sonicate?

We have obtained high-field (11.7-T) proton and carbon-13 Fourier transform (FT) nuclear magnetic resonance (NMR) spectra of egg lecithin and egg lecithin-cholesterol (1:1) multibilayers, using "magic-angle" sample spinning (MASS) techniques, and sonicated egg lecithin and egg lecithin-cholesterol (1:1) vesicles, using conventional FT NMR methods. Resolution of the proton and carbon-13 MASS NMR spectra of the pure egg lecithin samples is essentially identical with that of sonicated samples, but spectra of the unsonicated lipid, using MASS, can be obtained very much faster than with the more dilute, sonicated systems. With the 1:1 lecithin-cholesterol systems, proton MASS NMR spectra are virtually identical with conventional FT spectra of sonicated samples, while with 13C NMR, we demonstrate that most 13C nuclei in the cholesterol moiety can be monitored, even though these same nuclei are essentially invisible, i.e., are severely broadened, in the corresponding sonicated systems. In addition, 13C MASS NMR, spectra can again be recorded much faster than with sonicated samples, due to concentration effects. Taken together, these results strongly suggest there will seldom be need in the future to resort to ultrasonic disruption of lipid bilayer membranes in order to obtain high-resolution proton or carbon-13 NMR spectra.     

The application of DOSY NMR and molecular dynamics simulations to explore the mechanism(s) of micelle binding of antimicrobial peptides containing unnatural amino acids

Anionic and zwitterionic micelles are often used as simple models for the lipids found in bacterial and mammalian cell membranes to investigate antimicrobial peptide‐lipid interactions. In our laboratory we have employed a variety of 1D, 2D, and diffusion ordered (DOSY) NMR experiments to investigate the interactions of antimicrobial peptides containing unnatural amino acids with SDS and DPC micelles. Complete assignment of the proton spectra of these peptides is prohibited by the incorporation of a high percentage of unnatural amino acids which don't contain amide protons into the backbone. However preliminary assignment of the TOCSY spectra of compound 23 in the presence of both micelles indicated multiple conformers are present as a result of binding to these micelles. Chemical Shift Indexing agreed with previously collected CD spectra that indicated on binding to SDS micelles compound 23 adopts a mixture of α‐helical structures and on binding to DPC micelles this peptide adopts a mixture of helical and β‐turn/sheet like structures. DOSY NMR experiments also indicated that the total positive charge and the relative placement of that charge at the N‐terminus or C‐terminus are important in determining the mole fraction of the peptide that will bind to the different micelles. DOSY and ¹H‐NMR experiments indicated that the length of Spacer #1 plays a major role in defining the binding conformation of these analogs with SDS micelles. Results obtained from molecular simulations studies of the binding of compounds 23 and 36 with SDS micelles were consistent with the observed NMR results. © 2013 Wiley Periodicals, Inc. Biopolymers 99: 548–561, 2013.     

NMR relaxation parameters of methyl groups as a tool to map the interfaces of helix–helix interactions in membrane proteins

A number of recent advances in the field of magic-angle-spinning (MAS) solid-state NMR have enabled its application to a range of biological systems of ever increasing complexity. To retain biological relevance, these samples are increasingly studied in a hydrated state. At the same time, experimental feasibility requires the sample preparation process to attain a high sample concentration within the final MAS rotor. We discuss these considerations, and how they have led to a number of different approaches to MAS NMR sample preparation. We describe our experience of how custom-made (or commercially available) ultracentrifugal devices can facilitate a simple, fast and reliable sample preparation process. A number of groups have since adopted such tools, in some cases to prepare samples for sedimentation-style MAS NMR experiments. Here we argue for a more widespread adoption of their use for routine MAS NMR sample preparation.     

What do we mean by biological complexity?

The purpose of the present paper is to offer a precise definition of the concepts of integration, emergence and complexity in biological networks through the use of the information theory. If two distinct properties of a network are expressed by two discrete variables, the classical subadditivity principle of Shannon's information theory applies when all the nodes of the network are associated with these properties. If not, the subadditivity principle may not apply. This situation is often to be encountered with enzyme and metabolic networks, for some nodes may well not be associated with these two properties. This is precisely what is occurring with an enzyme that binds randomly its two substrates. This situation implies that an enzyme, or a metabolic network, may display a joint entropy equal, smaller, or larger than the corresponding sum of individual entropies of component sub-systems. In the first case, the collective properties of the network can be reduced to the individual properties of its components. Moreover, the network is devoid of any information. In the second case, the system displays integration effects, behaves as a coherent whole, and has positive information. But if the joint entropy of the network is smaller than the sum of the individual entropies of its components, then the system has emergent collective properties and can be considered complex. Moreover, under these conditions, its information is negative. The extent of negative information is enhanced if the enzyme, or the metabolic network, is far away from equilibrium.     

Integrative biology: linking levels of organization

Biological systems are composed of different levels of organization. Usually, one considers the atomic, molecular, cellular, individual, population, community and ecosystem levels. These levels of organization also correspond to different levels of observation of the system, from microscopic to macroscopic, i.e., to different time and space scales. The more microscopic the level is, the faster the time scale and the smaller the space scale are. The dynamics of the complete system is the result of the coupled dynamical processes that take place in each of its levels of organization at different time scales. Variables aggregation methods take advantage of these different time scales to reduce the dimension of mathematical models such as a system of ordinary differential equations. We are going to study the dynamics of a system which is hierarchically organized in the sense that it is composed of groups of elements that can be themselves divided into further smaller sub-groups and so on. The hierarchical structure of the system results from the fact that the intra-group interactions are assumed to be larger than inter-group ones. We present aggregation methods that allow one to build a reduced model that governs a few global variables at the slow time scale.     

Automated projection spectroscopy in solid-state NMR

Given that solid-state NMR is being used for protein samples of increasing molecular weight and complexity, higher-dimensionality methods are likely to be more and more indispensable for unambiguous chemical shift assignments in the near future. In addition, solid-state NMR spectral properties are increasingly comparable with solution NMR, allowing adaptation of more sophisticated solution NMR strategies for the solid state in addition to the conventional methodology. Assessing first principles, here we demonstrate the application of automated projection spectroscopy for a micro-crystalline protein in the solid state.     

Thermotolerance and regeneration in the brittle star species complex Ophioderma longicauda: A preliminary study comparing lineages and Mediterranean basins

Global warming is expected to change marine species distributions; it is thus critical to understand species current thermotolerance. The brittle star species complex Ophioderma longicauda comprises a broadcast spawning lineage L1 and a brooding lineage L3. We collected L1 specimens from Marseilles and Crete, and L3 specimens from Crete. We monitored survival, autotomy and arm regeneration at 17, 26 and 30 °C during 14 weeks. Globally O. longicauda showed good resistance to elevated temperatures compared to other published studies on ophiuroids. The L3 sample displayed a better thermotolerance than L1 samples. Yet, more research is needed to establish whether these differences are due to lineages, geographic origin, or random effects. We provided for the first time individual regeneration trajectories, and showed that regeneration followed a growth curve and was highly influenced by temperature in both lineages. Our results highlight the importance of taking into account the presence of cryptic species when studying the potential effects of global warming.     

Stochastic models for circadian rhythms: effect of molecular noise on periodic and chaotic behaviour

Circadian rhythms are endogenous oscillations that occur with a period close to 24 h in nearly all living organisms. These rhythms originate from the negative autoregulation of gene expression. Deterministic models based on such genetic regulatory processes account for the occurrence of circadian rhythms in constant environmental conditions (e.g., constant darkness), for entrainment of these rhythms by light-dark cycles, and for their phase-shifting by light pulses. When the numbers of protein and mRNA molecules involved in the oscillations are small, as may occur in cellular conditions, it becomes necessary to resort to stochastic simulations to assess the influence of molecular noise on circadian oscillations. We address the effect of molecular noise by considering the stochastic version of a deterministic model previously proposed for circadian oscillations of the PER and TIM proteins and their mRNAs in Drosophila. The model is based on repression of the per and tim genes by a complex between the PER and TIM proteins. Numerical simulations of the stochastic version of the model are performed by means of the Gillespie method. The predictions of the stochastic approach compare well with those of the deterministic model with respect both to sustained oscillations of the limit cycle type and to the influence of the proximity from a bifurcation point beyond which the system evolves to stable steady state. Stochastic simulations indicate that robust circadian oscillations can emerge at the cellular level even when the maximum numbers of mRNA and protein molecules involved in the oscillations are of the order of only a few tens or hundreds. The stochastic model also reproduces the evolution to a strange attractor in conditions where the deterministic PER-TIM model admits chaotic behaviour. The difference between periodic oscillations of the limit cycle type and aperiodic oscillations (i.e. chaos) persists in the presence of molecular noise, as shown by means of Poincaré sections. The progressive obliteration of periodicity observed as the number of molecules decreases can thus be distinguished from the aperiodicity originating from chaotic dynamics. As long as the numbers of molecules involved in the oscillations remain sufficiently large (of the order of a few tens or hundreds, or more), stochastic models therefore provide good agreement with the predictions of the deterministic model for circadian rhythms.     

Nuclear Magnetic Resonance Spectroscopy for the Study of B-Cell Epitopes

High-resolution nuclear magnetic resonance (NMR) spectroscopy is a structural technique that is finding increasing use in the study of antibody–antigen interactions. In this review we describe how the dynamic structural parameters obtained from NMR spectroscopy can further our understanding of B-cell epitopes and their function. Specific applications of NMR spectroscopy to examine the residues on peptides and proteins that contact the antibody combining site are also described. These include “footprinting” techniques using H–D exchange–COSY NMR spectroscopy, which are particularly useful for epitope mapping of protein antigens. For smaller systems, such as Fab–or Fv–peptide complexes, nuclear magnetization transfer difference NMR spectroscopy, transferred nuclear Overhauser effect spectroscopy, double-quantum-filtered NOE spectroscopy, and isotope editing techniques have been applied. The interpretation and limitations of the data obtained from these procedures and anticipated improvements in these applications in the future are discussed.     

Dynamics of the Hck-SH3 domain: comparison of experiment with multiple molecular dynamics simulations

Molecular dynamics calculations provide a method by which the dynamic properties of molecules can be explored over timescales and at a level of detail that cannot be obtained experimentally from NMR or X-ray analyses. Recent work (Philippopoulos M, Mandel AM, Palmer AG III, Lim C, 1997, Proteins 28:481-493) has indicated that the accuracy of these simulations is high, as measured by the correspondence of parameters extracted from these calculations to those determined through experimental means. Here, we investigate the dynamic behavior of the Src homology 3 (SH3) domain of hematopoietic cell kinase (Hck) via 5N backbone relaxation NMR studies and a set of four independent 4 ns solvated molecular dynamics calculations. We also find that molecular dynamics simulations accurately reproduce fast motion dynamics as estimated from generalized order parameter (S2) analysis for regions of the protein that have experimentally well-defined coordinates (i.e., stable secondary structural elements). However, for regions where the coordinates are not well defined, as indicated by high local root-mean-square deviations among NMR-determined structural family members or high B-factors/low electron density in X-ray crystallography determined structures, the parameters calculated from a short to moderate length (less than 5-10 ns) molecular dynamics trajectory are dependent on the particular coordinates chosen as a starting point for the simulation.     

Residual dipolar coupling constants and structure determination of large DNA duplexes

Several NMR works have shown that long-range information provided by residual dipolar couplings (RDCs) significantly improve the global structure definition of RNAs and DNAs. Most of these are based on the use of a large set of RDCs, the collect of which requires samples labeled with ¹³C, ¹⁵N, and sometimes, ²H. Here, we carried out torsion-angle dynamics simulations on a non-self complementary DNA fragment of 17 base-pairs, d(GGAAAATATCTAGCAGT).(ACTGCTAGAGATTTTCC). This reproduces the U5 LTR distal end of the HIV-1 cDNA that contains the enzyme integrase binding site. Simulations aimed at evaluating the impact of RDCs on the structure definition of long oligonucleotides, were performed in incorporating (i) nOe-distances at both < 4.5 Å and < 5 Å; (ii) a small set of ¹³C-¹H RDCs, easily detectable at the natural abundance, and (iii) a larger set of RDCs only accessible through the ¹³C labeling of DNAs. Agreement between a target structure and a simulated structure was measured in terms of precision and accuracy. Results allowed to define conditions in which accurate DNA structures can be determined. We confirmed the strong impact of RDCs on the structure determination, and, above all, we found that a small set of RDC constraints (ca. 50) detectable at the natural abundance is sufficient to accurately derive the global and local DNA duplex structures when used in conjunction with nOe-distances < 5 Å.     

Small angle X-ray scattering as a complementary tool for high-throughput structural studies

Structural crystallography and nuclear magnetic resonance (NMR) spectroscopy are the predominant techniques for understanding the biological world on a molecular level. Crystallography is constrained by the ability to form a crystal that diffracts well and NMR is constrained to smaller proteins. Although powerful techniques, they leave many soluble, purified structurally uncharacterized protein samples. Small angle X-ray scattering (SAXS) is a solution technique that provides data on the size and multiple conformations of a sample, and can be used to reconstruct a low-resolution molecular envelope of a macromolecule. In this study, SAXS has been used in a high-throughput manner on a subset of 28 proteins, where structural information is available from crystallographic and/or NMR techniques. These crystallographic and NMR structures were used to validate the accuracy of molecular envelopes reconstructed from SAXS data on a statistical level, to compare and highlight complementary structural information that SAXS provides, and to leverage biological information derived by crystallographers and spectroscopists from their structures. All the ab initio molecular envelopes calculated from the SAXS data agree well with the available structural information. SAXS is a powerful albeit low-resolution technique that can provide additional structural information in a high-throughput and complementary manner to improve the functional interpretation of high-resolution structures.