Characterization of NMR relaxation-active motions of a partially folded A-state analogue of ubiquitin

The dominant dynamics of a partially folded A-state analogue of ubiquitin that give rise to NMR 15N spin relaxation have been investigated using molecular dynamics (MD) computer simulations and reorientational quasiharmonic analysis. Starting from the X-ray structure of native ubiquitin with a protonation state corresponding to a low pH, the A-state analogue was generated by a MD simulation of a total length of 33 ns in a 60%/40% methanol/water mixture using a variable temperature scheme to control and speed up the structural transformation. The N-terminal half of the A-state analogue consists of loosely coupled native-like secondary structural elements, while the C-terminal half is mostly irregular in structure. Analysis of dipolar N-H backbone correlation functions reveals reorientational amplitudes and time-scale distributions that are comparable to those observed experimentally. Thus, the trajectory provides a realistic picture of a partially folded protein that can be used for gaining a better understanding of the various types of reorientational motions that are manifested in spin-relaxation parameters of partially folded systems. For this purpose, a reorientational quasiharmonic reorientational analysis was performed on the final 5 ns of the trajectory of the A-state analogue, and for comparison on a 5 ns trajectory of native ubiquitin. The largest amplitude reorientational modes show a markedly distinct behavior for the two states. While for native ubiquitin, such motions have a more local character involving loops and the C-terminal end of the polypeptide chain, the A-state analogue shows highly collective motions in the nanosecond time-scale range corresponding to larger-scale movements between different segments. Changes in reorientational backbone entropy between the A-state analogue and the native state of ubiquitin, which were computed from the reorientational quasiharmonic analyses, are found to depend significantly on motional correlation effects.     

Comparison of mode analyses at different resolutions applied to nucleic acid systems

More than two decades of different types of mode analyses has shown that these techniques can be useful in describing large-scale motions in protein systems. A number of mode analyses are available and include quasiharmonics, classical normal mode, block normal mode, and the elastic network model. Each of these methods has been validated for protein systems and this variety allows researchers to choose the technique that gives the best compromise between computational cost and the level of detail in the calculation. These same techniques have not been systematically tested for nucleic acid systems, however. Given the differences in interactions and structural features between nucleic acid and protein systems, the validity of these techniques in the protein regime cannot be directly translated into validity in the nucleic acid realm. In this work, we investigate the usefulness of the above mode analyses as applied to two RNA systems, i.e., the hammerhead ribozyme and a guanine riboswitch. We show that classical normal-mode analysis can match the magnitude and direction of residue fluctuations from the more detailed, anharmonic technique, quasiharmonic analysis of a molecular dynamics trajectory. The block normal-mode approximation is shown to hold in the nucleic acid systems studied. Only the mode analysis at the lowest level of detail, the elastic network model, produced mixed results in our calculations. We present data that suggest that the elastic network model, with the popular parameterization, is not best suited for systems that do not have a close packed structure; this observation also hints at why the elastic network model has been found to be valid for many globular protein systems. The different behaviors of block normal-mode analysis and the elastic network model, which invoke similar degrees of coarse-graining to the dynamics but use different potentials, suggest the importance of applying a heterogeneous potential function in a robust analysis of the dynamics of biomolecules, especially those that are not closely packed. In addition to these comparisons, we briefly discuss insights into the conformational space available to the hammerhead ribozyme.     

The entropic cost of protein-protein association: a case study on acetylcholinesterase binding to fasciculin-2

Protein-protein association is accompanied by a large reduction in translational and rotational (external) entropy. Based on a 15 ns molecular dynamics simulation of acetylcholinesterase (AChE) in complex with fasciculin 2 (Fas2), we estimate the loss in external entropy using quasiharmonic analysis and histogram-based approximations of the probability distribution function. The external entropy loss of AChE-Fas2 binding, ~30 cal/mol K, is found to be significantly larger than most previously characterized protein-ligand systems. However, it is less than the entropy loss estimated in an earlier study by A. V. Finkelstein and J. Janin, which was based on atomic motions in crystals.     

Normal modes of symmetric protein assemblies. Application to the tobacco mosaic virus protein disk.

We use group theoretical methods to study the molecular dynamics of symmetric protein multimers in the harmonic or quasiharmonic approximation. The method explicitly includes the long-range correlations between protein subunits. It can thus address collective dynamic effects, such as cooperativity between subunits. The n lowest-frequency normal modes of each individual subunit are combined into symmetry coordinates for the entire multimer. The Hessian of the potential energy is thereby reduced to a series of blocks of order n or 2n. In the quasiharmonic approximation, the covariance matrix of the atomic oscillations is reduced to the same block structure by an analogous set of symmetry coordinates. The method is applied to one layer of the tobacco mosaic virus protein disk in vacuo, to gain insight into the role of conformational fluctuations and electrostatics in tobacco mosaic virus assembly. The system has 78,000 classical, positional, degrees of freedom, yet the calculation is reduced by symmetry to a problem of order 4,600. Normal modes in the 0-100 cm-1 range were calculated. The calculated correlations extend mainly from each subunit to its nearest neighbors. The network of core helices has weak correlations with the rest of the structure. Similarly, the inner loops 90-108 are uncorrelated with the rest of the structure. Thus, the model predicts that the dielectric response in the RNA-binding region is mainly due to the loops alone.     

Calculation of absolute protein-ligand binding affinity using path and endpoint approaches

A comparative analysis is provided of rigorous and approximate methods for calculating absolute binding affinities of two protein-ligand complexes: the FKBP protein bound with small molecules 4-hydroxy-2-butanone and FK506. Our rigorous approach is an umbrella sampling technique where a potential of mean force is determined by pulling the ligand out of the protein active site over several simulation windows. The results of this approach agree well with experimentally observed binding affinities. Also assessed is a commonly used approximate endpoint approach, which separately estimates enthalpy, solvation free energy, and entropy. We show that this endpoint approach has numerous variations, all of which are prone to critical shortcomings. For example, conventional harmonic and quasiharmonic entropy estimation procedures produce disparate results for the relatively simple protein-ligand systems studied in this work.     

Functional domain motions in proteins on the ~1-100 ns timescale: comparison of neutron spin-echo spectroscopy of phosphoglycerate kinase with molecular-dynamics simulation

Protein function often requires large-scale domain motion. An exciting new development in the experimental characterization of domain motions in proteins is the application of neutron spin-echo spectroscopy (NSE). NSE directly probes coherent (i.e., pair correlated) scattering on the ~1-100 ns timescale. Here, we report on all-atom molecular-dynamics (MD) simulation of a protein, phosphoglycerate kinase, from which we calculate small-angle neutron scattering (SANS) and NSE scattering properties. The simulation-derived and experimental-solution SANS results are in excellent agreement. The contributions of translational and rotational whole-molecule diffusion to the simulation-derived NSE and potential problems in their estimation are examined. Principal component analysis identifies types of domain motion that dominate the internal motion's contribution to the NSE signal, with the largest being classic hinge bending. The associated free-energy profiles are quasiharmonic and the frictional properties correspond to highly overdamped motion. The amplitudes of the motions derived by MD are smaller than those derived from the experimental analysis, and possible reasons for this difference are discussed. The MD results confirm that a significant component of the NSE arises from internal dynamics. They also demonstrate that the combination of NSE with MD is potentially useful for determining the forms, potentials of mean force, and time dependence of functional domain motions in proteins.     

Dynamic optimization of distributed biological systems using robust and efficient numerical techniques

ABSTRACT: BACKGROUND: Systems biology allows the analysis of biological systems behavior under different conditions through in silico experimentation. The possibility of perturbing biological systems in different manners calls for the design of perturbations to achieve particular goals. Examples would include, the design of a chemical stimulation to maximize the amplitude of a given cellular signal or to achieve a desired pattern in pattern formation systems, etc. Such design problems can be mathematically formulated as dynamic optimization problems which are particularly challenging when the system is described by partial differential equations. This work addresses the numerical solution of such dynamic optimization problems for spatially distributed biological systems. The usual nonlinear and large scale nature of the mathematical models related to this class of systems and the presence of constraints on the optimization problems, impose a number of difficulties, such as the presence of suboptimal solutions, which call for robust and efficient numerical techniques. RESULTS: Here, the use of a control vector parameterization approach combined with efficient and robust hybrid global optimization methods and a reduced order model methodology is proposed. The capabilities of this strategy are illustrated considering the solution of a two challenging problems: bacterial chemotaxis and the FitzHugh-Nagumo model. CONCLUSIONS: In the process of chemotaxis the objective was to efficiently compute the time-varying optimal concentration of chemotractant in one of the spatial boundaries in order to achieve predefined cell distribution profiles. Results are in agreement with those previously published in the literature. The FitzHugh-Nagumo problem is also efficiently solved and it illustrates very well how dynamic optimization may be used to force a system to evolve from an undesired to a desired pattern with a reduced number of actuators. The presented methodology can be used for the efficient dynamic optimization of generic distributed biological systems.     

无细胞蛋白表达体系研究进展及在生物制药领域中的应用

作为一种快速高效的体外蛋白合成手段,无细胞蛋白表达体系(Cell-free Protein Synthesis,CFPS)一直以来就被广泛应用于基础生物学领域的研究。与传统的基于细胞的体内表达体系相比,CFPS突破了细胞的生理限制,其可调控性强、对毒性蛋白的耐受力高,使得许多很难在体内合成的复杂蛋白在体外顺利表达。近年来随着研究人员不断对CFPS进行优化,通过简化制备工艺、开发价格低廉的能量再生系统、稳定底物供应、促进蛋白正确折叠等方式,成功研发出生产效率高、成本低廉、反应体积大的表达体系。凭借其高通量和大规模的蛋白表达优势,CFPS为解决生物制药领域中面临的难题提供了新的解决思路,并成功地应用于高通量药物筛选、大规模生产重组蛋白药物、个体化定制肿瘤疫苗等领域,显示出其在生物制药领域的重要应用潜力。     

Poly(ethylene glycol) interfaces: an approach for enhanced performance of microfluidic systems

Microfluidic systems are extensively used platform for analytical and therapeutic applications. One of the major problems encountered in these systems is the loss of material due to non-specific surface interactions. When biological solutions are flowed through microchannels, they tend to adsorb on the surface due to the negative charge of the surface. This results in a reduced efficiency of the system which can be critical in sensitive biological analysis. Poly(ethylene glycol) (PEG) is known to form non-fouling interfaces on silicon and glass which are common materials used in microfluidic systems. The most common approach for modifying silicon/glass with PEG involves a solution phase protocol. Since the micro/nanofluidic systems have channel sizes ranging from hundreds of microns to submicron with variety of complicated network, this surface modification approach is not sufficient in forming uniform, conformal, and ultrathin films on the surface. Due to the enclosed features in these systems, the properties of liquids such as viscosity and surface tension play an important role in the clogging and eventually biofouling of these microchannels. Hence, we have developed a solvent-free vapor deposition protocol for modifying silicon/glass surfaces with PEG. Various concentrations of protein solutions were flowed through unmodified and PEG-modified glass microcapillaries of different lengths at different flow rates. PEG surfaces formed on silicon have shown 80% reduction in protein adsorption in static conditions.     

Cancer-related mutations in BRCA1-BRCT cause long-range structural changes in protein-protein binding sites: a molecular dynamics study

Cancer-associated mutations in the BRCT domain of BRCA1 (BRCA1-BRCT) abolish its tumor suppressor function by disrupting interactions with other proteins such as BACH1. Many cancer-related mutations do not cause sufficient destabilization to lead to global unfolding under physiological conditions, and thus abrogation of function probably is due to localized structural changes. To explore the reasons for mutation-induced loss of function, the authors performed molecular dynamics simulations on three cancer-associated mutants, A1708E, M1775R, and Y1853ter, and on the wild type and benign M1652I mutant, and compared the structures and fluctuations. Only the cancer-associated mutants exhibited significant backbone structure differences from the wild-type crystal structure in BACH1-binding regions, some of which are far from the mutation sites. Backbone differences of the A1708E mutant from the liganded wild type structure in these regions are much larger than those of the unliganded wild type X-ray or molecular dynamics structures. These BACH1-binding regions of the cancer-associated mutants also exhibited increases in their fluctuation magnitudes compared with the same regions in the wild type and M1562I mutant, as quantified by quasiharmonic analysis. Several of the regions of increased fluctuation magnitude correspond to correlated motions of residues in contact that provide a continuous path of fluctuating amino acids in contact from the A1708E and Y1853ter mutation sites to the BACH1-binding sites with altered structure and dynamics. The increased fluctuations in the disease-related mutants suggest an increase in vibrational entropy in the unliganded state that could result in a larger entropy loss in the disease-related mutants upon binding BACH1 than in the wild type. To investigate this possibility, vibrational entropies of the A1708E and wild type in the free state and bound to a BACH1-derived phosphopeptide were calculated using quasiharmonic analysis, to determine the binding entropy difference DeltaDeltaS between the A1708E mutant and the wild type. DeltaDeltaS was determined to be -4.0 cal mol(-1) K(-1), with an uncertainty of 2 cal mol(-1) K(-1); that is, the entropy loss upon binding the peptide is 4.0 cal mol(-1) K(-1) greater for the A1708E mutant, corresponding to an entropic contribution to the DeltaDeltaG of binding (-TDeltaDeltaS) 1.1 kcal mol(-1) more positive for the mutant. The observed differences in structure, flexibility, and entropy of binding likely are responsible for abolition of BACH1 binding, and illustrate that many disease- related mutations could have very long-range effects. The methods described here have potential for identifying correlated motions responsible for other long-range effects of deleterious mutations.     

Mathematical modeling of differentiation in dictyostelium discoideum

Methods for the dynamic analysis of biochemical differentiation are presented. These are demonstrated in the analysis of biochemical differentiation of the carbohydrate system in D. discoideum. Procedures for simplification which are presented are projection and contraction of the system trajectory in state space and the generation of reduced equivalent dynamic metabolic networks. The importance of the hierarchical structure of differentiating systems is discussed and the concept of a dynamic embedding diagram is introduced. It is shown that complex systems must be analyzed on an epoch by epoch basis, each epoch being a period of time characterized by a constant dynamic embedding diagram, and that widely different time scales and state space scales may be necessary in different epochs. In particular there is no a priori lower limit to the time scale which may be necessary during the analysis. Some problems in mathematically defining differentiation are discussed.     

Induced fit and the entropy of structural adaptation in the complexation of CAP and lambda-repressor with cognate DNA sequences

Molecular dynamics (MD) simulations of 5 ns on protein-DNA complexes of catabolite-activator protein (CAP), lambda-repressor, and their corresponding uncomplexed protein and DNA, are reported. These cases represent two extremes of DNA bending, with CAP DNA bent severely and the lambda-operator nearly straight when complexed with protein. The calculations were performed using the AMBER suite of programs and the parm94 force field, validated for these studies by good agreement with experimental nuclear magnetic resonance data on DNA. An explicit computational model of structural adaptation and computation of the quasiharmonic entropy of association were obtained from the MD. The results indicate that, with respect to canonical B-form DNA, the extreme bending of the DNA in the complex with CAP is approximately 60% protein-induced and 40% intrinsic to the sequence-dependent structure of the free oligomer. The DNA in the complex is an energetically strained form, and the MD results are consistent with a conformational-capture mechanism. The calculated quasiharmonic entropy change accounts for the entropy difference between the two cases. The calculated entropy was decomposed into contributions from protein adaptation, DNA adaptation, and protein-DNA structural correlations. The origin of the entropy difference between CAP and lambda-repressor complexation arises more from the additional protein adaptation in the case of lambda, than to DNA bending and entropy contribution from DNA bending. The entropy arising from protein DNA cross-correlations, a contribution not previously discussed, is surprisingly large.     

Towards a molecular understanding of human diseases using Dictyostelium discoideum

The social amoeba Dictyostelium discoideum is increasingly being used as a simple model for the investigation of problems that are relevant to human health. This article focuses on several recent examples of Dictyostelium-based biomedical research, including the analysis of immune-cell disease and chemotaxis, centrosomal abnormalities and lissencephaly, bacterial intracellular pathogenesis, and mechanisms of neuroprotective and anti-cancer drug action. The combination of cellular, genetic and molecular biology techniques that are available in Dictyostelium often makes the analysis of these problems more amenable to study in this system than in mammalian cell culture. Findings that have been made in these areas using Dictyostelium have driven research in mammalian systems and have established Dictyostelium as a powerful model for human-disease analysis.     

Predicting the Effects of Woody Encroachment on Mammal Communities,Grazing Biomass and Fire Frequency in African Savannas

With grasslands and savannas covering 20% of the world’s land surface, accounting for 30–35% of worldwide Net Primary Productivity and supporting hundreds of millions of people, predicting changes in tree/grass systems is priority. Inappropriate land management and rising atmospheric CO₂ levels result in increased woody cover in savannas. Although woody encroachment occurs world-wide, Africa’s tourism and livestock grazing industries may be particularly vulnerable. Forecasts of responses of African wildlife and available grazing biomass to increases in woody cover are thus urgently needed. These predictions are hard to make due to non-linear responses and poorly understood feedback mechanisms between woody cover and other ecological responders, problems further amplified by the lack of long-term and large-scale datasets. We propose that a space-for-time analysis along an existing woody cover gradient overcomes some of these forecasting problems. Here we show, using an existing woody cover gradient (0–65%) across the Kruger National Park, South Africa, that increased woody cover is associated with (i) changed herbivore assemblage composition, (ii) reduced grass biomass, and (iii) reduced fire frequency. Furthermore, although increased woody cover is associated with reduced livestock production, we found indigenous herbivore biomass (excluding elephants) remains unchanged between 20–65% woody cover. This is due to a significant reorganization in the herbivore assemblage composition, mostly as a result of meso-grazers being substituted by browsers at increasing woody cover. Our results suggest that woody encroachment will have cascading consequences for Africa’s grazing systems, fire regimes and iconic wildlife. These effects will pose challenges and require adaptation of livelihoods and industries dependent on conditions currently prevailing.     

Flexibility and conformational entropy in protein-protein binding

To better understand the interplay between protein-protein binding and protein dynamics, we analyzed molecular dynamics simulations of 17 protein-protein complexes and their unbound components. Complex formation does not restrict the conformational freedom of the partner proteins as a whole, but, rather, it leads to a redistribution of dynamics. We calculate the change in conformational entropy for seven complexes with quasiharmonic analysis. We see significant loss, but also increased or unchanged conformational entropy. Where comparison is possible, the results are consistent with experimental data. However, stringent error estimates based on multiple independent simulations reveal large uncertainties that are usually overlooked. We observe substantial gains of pseudo entropy in individual partner proteins, and we observe that all complexes retain residual stabilizing intermolecular motions. Consequently, protein flexibility has an important influence on the thermodynamics of binding and may disfavor as well as favor association. These results support a recently proposed unified model for flexible protein-protein association.     

Inert filter media for the biofiltration of waste gases – characteristics and biomass control

Soil biofilters and related systems based onthe use of natural filter beds have been usedfor several years for solving specific airpollution problems. Over the past decade,significant improvements have been brought tothese original bioprocesses, among which thedevelopment and use of new inert packingmaterials. The present paper overviews the mostcommon inert packings used in biofiltration ofwaste gases and their major characteristics. Apotential problem recently encountered whenusing inert filter beds is the heterogenousdistribution of biomass on the packingmaterial, and the excessive growth andaccumulation of biomass when treating highorganic loads, eventually leading to cloggingof the biofilter and reduced efficiency.Several strategies that have been proposed forsolving such problems are described in thispaper. Technologies for controlling excessbiomass accumulation can be grouped into fourcategories based on the use of mechanicalforces, the use of specific chemicals, thereduction of microbial growth, and predation.     

Endogenous regulatory oligopeptides: their structure, functions and localization

A number of problems connected with the study of several hundreds of known endogenous peptide molecules have been considered on the basis of data obtained from EROP-Moscow database. A large number of peptide structures can be reduced to a relatively small number of peptide families formed on the basis of homology of amino acid sequence. The ability of most peptides to participate in different regulatory systems of an organism has been demonstrated and the physical and chemical basis for specificity and ambiguity has been discussed. Most problems of structure, function, and localization of endogenous regulatory peptides are similar both in higher and in lower organisms.     

Standardisation of DNA quantitation by image analysis: quality control of instrumentation.

BACKGROUND: DNA image analysis is frequently performed in clinical practice as a prognostic tool and to improve diagnosis. The precision of prognosis and diagnosis depends on the accuracy of analysis and particularly on the quality of image analysis systems. It has been reported that image analysis systems used for DNA quantification differ widely in their characteristics (Thunissen et al.: Cytometry 27: 21-25, 1997). This induces inter-laboratory variations when the same sample is analysed in different laboratories. In microscopic image analysis, the principal instrumentation errors arise from the optical and electronic parts of systems. They bring about problems of instability, non-linearity, and shading and glare phenomena. METHODS: The aim of this study is to establish tools and standardised quality control procedures for microscopic image analysis systems. Specific reference standard slides have been developed to control instability, non-linearity, shading and glare phenomena and segmentation efficiency. RESULTS: Some systems have been controlled with these tools and these quality control procedures. Interpretation criteria and accuracy limits of these quality control procedures are proposed according to the conclusions of a European project called PRESS project (Prototype Reference Standard Slide). Beyond these limits, tested image analysis systems are not qualified to realise precise DNA analysis. CONCLUSIONS: The different procedures presented in this work determine if an image analysis system is qualified to deliver sufficiently precise DNA measurements for cancer case analysis. If the controlled systems are beyond the defined limits, some recommendations are given to find a solution to the problem.     

Electrochemical arrays coupled with magnetic separators for immunochemistry

Electrochemical methods are increasingly applied to immunoassays, because they overcome problems associated with other modes of detection. In particular, with respect to conventional immunoassays, electrochemical immunosensors show versatility, reliability, and fast analysis time. In immunosensor strategy, the antigen or antibody can be immobilized directly onto the surface of the electrochemical transducer that will finally be used to reveal the amount of the affinity reaction. However, the use of the electrode surface as a solid phase as well as an electrochemical transducer presents some problems: a shielding of the surface by biospecifically bound antibody molecules can cause hindrance in the electron transfer, resulting in a reduced voltammetric signal. Thus, as an alternative solid phase, magnetic beads because of their low toxicity and high biocompatibility have gained much attention in chemistry, associated with various analytical techniques, due to their suitability for immobilization of biomolecules. Magnetic micro- or nanobeads can be separated easily and quickly by magnetic forces and will be used together with bioaffine ligands, e.g., antibodies or proteins with a high affinity to the target. The special advantages of magnetic separation techniques are the fast and simple handling of a sample vial and the opportunity to deal with large sample volumes without the need for time-consuming centrifugation steps. This also makes biomagnetic separation ideal for automated assay/analysis systems which will play a very important role in the near future. This review presents some examples of immunochemical assay developed using magnetic beads as a solid phase coupled with electrochemical detection techniques, in particular, using electrochemical arrays as transducers. Applications related to static measurements, together with in-flow detection systems are presented.