Robust refocusing of 13C magnetization in multidimensional NMR experiments by adiabatic fast passage pulses

We show that adiabatic fast passage (AFP) pulses are robust refocusing elements of transverse ¹³C magnetization in multidimensional NMR experiments. A pair of identical AFP pulses can refocus selected parts or a complete ¹³ C chemical shift range in ¹³C spectra. In the constant time ¹³C-¹H HSQC, replacement of attenuated rectangular pulses by selective AFP pulses results in a sensitivity enhancement of up to a factor of 1.8. In the 3D CBCA(CO)NH the signal-to-noise ratio is increased by a factor of up to 1.6.     

Determination of the rotational dynamics and pH dependence of the hydrogen exchange rates of the arginine guanidino group using NMR spectroscopy

Summary An improved version of the constant-time HSQC experiment is presented that gives uniform sensitivity over the complete ¹³C bandwidth in ¹³C−¹H correlation experiments without creating artifacts in the methyl and aromatic regions of the spectra. The improvement is achieved by replacing the refocussing ¹³C 180° pulse in the evolution time by a combination of a full-power (22 kHz) hyperbolic secant 180° pulse that inverts and refocusses the entire ¹³C window, immediately followed by a selective 180° pulse on the CO region. Further improvement in signal-to-noise in the aromatic and methyl regions, although less spectacular, is obtained by replacing the other two 180° ¹³C pulses in the INEPT parts of the pulse sequence by full-power hyperbolic secant pulses. Results of simulations and experimental data are presented that demonstrate the excellent performance of the hyperbolic secant pulse for broadband inversion and show that refocussing of transverse magnetization occurs over the same bandwidth, albeit with a ¹³C signal phase that depends quadratically on offset. A further modification, in which one of the selective pulses on the CO region is omitted, is also presented. Implications for other 2D and 3D experiments performed at high fields, where uniform ¹³C inversion and refocussing is desirable, are discussed.     

Tailoring broadband inversion pulses for MAS Solid state NMR

A simple approach is demonstrated for designing optimised broadband inversion pulses for MAS solid state NMR studies of biological systems. The method involves a two step numerical optimisation procedure and takes into account experimental requirements such as the pulse length, resonance offset range and extent of H₁ inhomogeneity compensation needed. A simulated annealing protocol is used initially to find appropriate values for the parameters that define the well known tanh/tan adiabatic pulse such that a satisfactory spin inversion is achieved with minimum RF field strength. This information is then used in the subsequent stage of refinement where the RF pulse characteristics are further tailored via a local optimisation procedure without imposing any restrictions on the amplitude and frequency modulation profiles. We demonstrate that this approach constitutes a generally applicable tool for obtaining pulses with good inversion characteristics. At moderate MAS frequencies the efficacy of the method is experimentally demonstrated for generating double-quantum NMR spectra via the zero-quantum dipolar recoupling scheme RFDR.     

Multiscale modeling of cellular actin filaments: from atomistic molecular to coarse-grained dynamics

In this article, we present a computational multiscale model for the characterization of subcellular proteins. The model is encoded inside a simulation tool that builds coarse-grained (CG) force fields from atomistic simulations. Equilibrium molecular dynamics simulations on an all-atom model of the actin filament are performed. Then, using the statistical distribution of the distances between pairs of selected groups of atoms at the output of the MD simulations, the force field is parameterized using the Boltzmann inversion approach. This CG force field is further used to characterize the dynamics of the protein via Brownian dynamics simulations. This combination of methods into a single computational tool flow enables the simulation of actin filaments with length up to 400 nm, extending the time and length scales compared to state-of-the-art approaches. Moreover, the proposed multiscale modeling approach allows to investigate the relationship between atomistic structure and changes on the overall dynamics and mechanics of the filament and can be easily (i) extended to the characterization of other subcellular structures and (ii) used to investigate the cellular effects of molecular alterations due to pathological conditions.     

Unsteady processes during stimulated emission from a relativistic electron beam in a quasi-longitudinal electrostatic pump field

The problem of stimulated emission from a relativistic electron beam in an external electrostatic pump field is studied. A set of nonlinear time-dependent equations for the spatiotemporal dynamics of the undulator radiation amplitude and the amplitude of the beam space charge field is derived. The beam electrons are described by a modified version of the macroparticle method. The regimes of the single-particle and collective Cherenkov effects during convective and absolute instabilities are considered. The nonlinear dynamics of radiation pulses emitted during the instabilities of the beam in its interaction with the forward and backward electromagnetic waves is investigated.     

Amplification of radiation damping in a 600-MHz NMR spectrometer: Application to the study of water-protein interactions

Summary A new application of a recently developed electronic radiation-damping (RD) control system is presented. It is possible to amplify radiation damping so as to make the water magnetization return back to its equilibrium direction in a time shorter than the characteristic RD time. Certain types of experiments involving radiation damping as a selective inversion pulse can be significantly improved by this new method. Moreover, amplification of RD is shown to improve water suppression and consequently the dynamics of 2D NOESY experiments on proteins.     

The biphasic behavior of incoherent feed-forward loops in biomolecular regulatory networks

An incoherent feed-forward loop (FFL) is one of the most-frequently observed motifs in biomolecular regulatory networks. It has been thought that the incoherent FFL is designed simply to induce a transient response shaped by a 'fast activation and delayed inhibition'. We find that the dynamics of various incoherent FFLs can be further classified into two types: time-dependent biphasic responses and dose-dependent biphasic responses. Why do the structurally identical incoherent FFLs play such different dynamical roles? Through computational studies, we show that the dynamics of the two types of incoherent FFLs are mutually exclusive. Following from further computational results and experimental observations, we hypothesize that incoherent FFLs have been optimally designed to achieve distinct biological function arising from different cellular contexts. Additional Supporting Information may be found in the online version of the article.     

The early effects of chemical carcinogens on adult rat hepatocytes in primary culture. II. Effects on unscheduled DNA synthesis, cell division and alpha-fetoprotein production.

The effects of exposure of rat hepatocytes in primary maintenance culture to chemical carcinogens has been studied with respect cytotoxicity and alterations in mitotic index, unscheduled DNA synthesis and alpha-fetoprotein (AFP) production. All compounds tested produced cytotoxicity. Increases in mitotic index and unscheduled DNA synthesis and the production of AFP were observed after treatment of the cells with the carcinogens but not after treatment with the non-carcinogenic isomers. These increases were dose-dependent and depended on the time of exposure and the time incubated postexposure. The patterns of the increase in mitotic index and AFP production after cessation of carcinogen exposure were very similar, with the increase in mitotic index occurring slightly before that for the AFP production and it is suggested from this and other data that the production of AFP is dependent on the generation of a cell species functionally distinct from the non-dividing hepatocytes. It is also suggested that measurement of unscheduled DNA synthesis in conjunction with that of AFP production in cultured hepatocytes may be useful as part of a screening programme for chemical carcinogens.     

Computation and analysis of time-dependent sensitivities in Generalized Mass Action systems

Understanding biochemical system dynamics is becoming increasingly important for insights into the functioning of organisms and for biotechnological manipulations, and additional techniques and methods are needed to facilitate investigations of dynamical properties of systems. Extensions to the method of Ingalls and Sauro, addressing time-dependent sensitivity analysis, provide a new tool for executing such investigations. We present here the results of sample analyses using time-dependent sensitivities for three model systems taken from the literature, namely an anaerobic fermentation pathway in yeast, a negative feedback oscillator modeling cell-cycle phenomena, and the Mitogen Activated Protein (MAP) kinase cascade. The power of time-dependent sensitivities is particularly evident in the case of the MAPK cascade. In this example it is possible to identify the emergence of a concentration of MAPKK that provides the best response with respect to rapid and efficient activation of the cascade, while over- and under-expression of MAPKK relative to this concentration have qualitatively different effects on the transient response of the cascade. Also of interest is the quite general observation that phase-plane representations of sensitivities in oscillating systems provide insights into the manner with which perturbations in the envelope of the oscillation result from small changes in initial concentrations of components of the oscillator. In addition to these applied analyses, we present an algorithm for the efficient computation of time-dependent sensitivities for Generalized Mass Action (GMA) systems, the most general of the canonical system representations of Biochemical Systems Theory (BST). The algorithm is shown to be comparable to, or better than, other methods of solution, as exemplified with three biochemical systems taken from the literature.     

Effects of warming rate,temperature, and antifreeze proteins on the survival of mouse spermatozoa frozen at an optimal rate

We have recently reported that the survival of mouse spermatozoa is decreased when they are warmed at a suboptimal rate after being frozen at an optimal rate. We proposed that this drop in survival is caused by physical damage derived from the recrystallization of extracellular ice during slow warming. The first purpose of the present study was to determine the temperatures over which the decline in survival occurs during slow warming and the kinetics of the decline at fixed subzero temperatures. The second purpose was to examine the effects of antifreeze proteins (AFP) on the survival of slowly warmed mouse spermatozoa, the rationale being that AFP have the property of inhibiting ice recrystallization. With respect to the first point, a substantial loss in motility occurred when slow warming was continued to higher than -50 degrees C and the survival of the sperm decreased with an increase in the temperature at which slow warming was terminated. In contrast, the motility of sperm that were warmed rapidly to these temperatures remained high initially but dropped with increased holding time. At -30 degrees C, most of the drop occurred in 5 min. These results are consistent with the hypothesis that damage develops as a consequence of the recrystallization of the external ice. AFP ought to inhibit such recrystallization, but we found that the addition of AFP-I, AFP-III, and an antifreeze glycoprotein at concentrations of 1-100 microg/ml did not protect the frozen-thawed cells; rather it led to a decrease in survival that was proportional to the concentration. There was no decrease in survival from exposure to the AFP in the absence of freezing. AFP are known to produce changes in the structure and habit of ice crystals, and some have reported deleterious consequences associated with those structural changes. We suggest that such changes may be the basis of the adverse effects of AFP on the survival of the sperm, especially since mouse sperm are exquisitely sensitive to a variety of mechanical stresses.     

Analysis of ice-binding sites in fish type II antifreeze protein by quantum mechanics

Many organisms living in cold environments can survive subzero temperatures by producing antifreeze proteins (AFPs) or antifreeze glycoproteins. In this paper we investigate the ice-binding surface of type II AFP by quantum mechanical methods, which, to the best of our knowledge, represents the first time that molecular orbital computational approaches have been applied to AFPs. Molecular mechanical approaches, including molecular docking, energy minimization, and molecular dynamics simulation, were used to obtain optimal systems for subsequent quantum mechanical analysis. We selected 17 surface patches covering the entire surface of the type II AFP and evaluated the interaction energy between each of these patches and two different ice planes using semi-empirical quantum mechanical methods. We have demonstrated the weak orbital overlay phenomenon and the change of bond orders in ice. These results consistently indicate that a surface patch containing 19 residues (K37, L38, Y20, E22, Y21, I19, L57, T56, F53, M127, T128, F129, R17, C7, N6, P5, G10, Q1, and W11) is the most favorable ice-binding site for both a regular ice plane and an ice plane where water O atoms are randomly positioned. Furthermore, for the first time the computation results provide new insights into the weakening of the ice lattice upon AFP binding, which may well be a primary factor leading to AFP-induced ice growth inhibition.     

Computational investigation of the dynamic control of cAMP signaling by PDE4 isoform types

Cyclic adenosine monophosphate (cAMP) is a generic signaling molecule that, through precise control of its signaling dynamics, exerts distinct cellular effects. Consequently, aberrant cAMP signaling can have detrimental effects. Phosphodiesterase 4 (PDE4) enzymes profoundly control cAMP signaling and comprise different isoform types wherein enzymatic activity is modulated by differential feedback mechanisms. Because these feedback dynamics are non-linear and occur coincidentally, their effects are difficult to examine experimentally but can be well simulated computationally. Through understanding the role of PDE4 isoform types in regulating cAMP signaling, PDE4-targeted therapeutic strategies can be better specified. Here, we established a computational model to study how feedback mechanisms on different PDE4 isoform types lead to dynamic, isoform-specific control of cAMP signaling. Ordinary differential equations describing cAMP dynamics were implemented in the VirtualCell environment. Simulations indicated that long PDE4 isoforms exert the most profound control on oscillatory cAMP signaling, as opposed to the PDE4-mediated control of single cAMP input pulses. Moreover, elevating cAMP levels or decreasing PDE4 levels revealed different effects on downstream signaling. Together these results underline that cAMP signaling is distinctly regulated by different PDE4 isoform types and that this isoform specificity should be considered in both computational and experimental follow-up studies to better define PDE4 enzymes as therapeutic targets in diseases in which cAMP signaling is aberrant.     

Peptide backbone circularization enhances antifreeze protein thermostability

Antifreeze proteins (AFPs) are a class of ice‐binding proteins that promote survival of a variety of cold‐adapted organisms by decreasing the freezing temperature of bodily fluids. A growing number of biomedical, agricultural, and commercial products, such as organs, foods, and industrial fluids, have benefited from the ability of AFPs to control ice crystal growth and prevent ice recrystallization at subzero temperatures. One limitation of AFP use in these latter contexts is their tendency to denature and irreversibly lose activity at the elevated temperatures of certain industrial processing or large‐scale AFP production. Using the small, thermolabile type III AFP as a model system, we demonstrate that AFP thermostability is dramatically enhanced via split intein‐mediated N‐ and C‐terminal end ligation. To engineer this circular protein, computational modeling and molecular dynamics simulations were applied to identify an extein sequence that would fill the 20‐Å gap separating the free ends of the AFP, yet impose little impact on the structure and entropic properties of its ice‐binding surface. The top candidate was then expressed in bacteria, and the circularized protein was isolated from the intein domains by ice‐affinity purification. This circularized AFP induced bipyramidal ice crystals during ice growth in the hysteresis gap and retained 40% of this activity even after incubation at 100°C for 30 min. NMR analysis implicated enhanced thermostability or refolding capacity of this protein compared to the noncyclized wild‐type AFP. These studies support protein backbone circularization as a means to expand the thermostability and practical applications of AFPs.     

Chemical shift correlation via RFDR: elimination of resonance offset effects

It is shown that it is possible to effectively execute RFDR experiments with adiabatic inversion pulses and obtain resonance offset compensation that is superior to what can be achieved by conventional rectangular pulses. Employing 40-s tanh/tan adiabatic pulses at a power level of 38 kHz and a spinning speed of 12 kHz it is demonstrated that the range of resonance offset compensation achieved is sufficient to generate, via a single experiment, homonuclear chemical shift correlation spectra in the entire ¹³C chemical shift range in peptides/proteins at the currently available field strengths.     

Robustness Analysis of Stochastic Biochemical Systems

We propose a new framework for rigorous robustness analysis of stochastic biochemical systems that is based on probabilistic model checking techniques. We adapt the general definition of robustness introduced by Kitano to the class of stochastic systems modelled as continuous time Markov Chains in order to extensively analyse and compare robustness of biological models with uncertain parameters. The framework utilises novel computational methods that enable to effectively evaluate the robustness of models with respect to quantitative temporal properties and parameters such as reaction rate constants and initial conditions. We have applied the framework to gene regulation as an example of a central biological mechanism where intrinsic and extrinsic stochasticity plays crucial role due to low numbers of DNA and RNA molecules. Using our methods we have obtained a comprehensive and precise analysis of stochastic dynamics under parameter uncertainty. Furthermore, we apply our framework to compare several variants of two-component signalling networks from the perspective of robustness with respect to intrinsic noise caused by low populations of signalling components. We have successfully extended previous studies performed on deterministic models (ODE) and showed that stochasticity may significantly affect obtained predictions. Our case studies demonstrate that the framework can provide deeper insight into the role of key parameters in maintaining the system functionality and thus it significantly contributes to formal methods in computational systems biology.     

Adiabatic TOCSY for C,C and H,H J-transfer

Adiabatic pulses have been widely used for broadband decoupling and spin inversion at high magnetic fields. In this paper we propose adiabatic pulses and supercycles that can be used at high magnetic fields like 800 or 900 MHz to obtain broadband TOCSY sequences with C,C or H,H J-transfer. The new mixing sequences are equal or even superior to the well known DIPSI-2,3 experiments with respect to bandwidth. They prove robust against pulse miscalibration and B₁ inhomogeneity and are therefore attractive for fully automated spectrometer environments. These adiabatic mixing sequences have been incorporated in a novel z-filter HCCH-TOCSY experiment.     

Interaction Mechanisms of Cavitation Bubbles Induced by Spatially and Temporally Separated fs-Laser Pulses

The emerging use of femtosecond lasers with high repetition rates in the MHz regime together with limited scan speed implies possible mutual optical and dynamical interaction effects of the individual cutting spots. In order to get more insight into the dynamics a time-resolved photographic analysis of the interaction of cavitation bubbles is presented. Particularly, we investigated the influence of fs-laser pulses and their resulting bubble dynamics with various spatial as well as temporal separations. Different time courses of characteristic interaction effects between the cavitation bubbles were observed depending on pulse energy and spatio-temporal pulse separation. These ranged from merely no interaction to the phenomena of strong water jet formation. Afterwards, the mechanisms are discussed regarding their impact on the medical application of effective tissue cutting lateral to the laser beam direction with best possible axial precision: the mechanical forces of photodisruption as well as the occurring water jet should have low axial extend and a preferably lateral priority. Furthermore, the overall efficiency of energy conversion into controlled mechanical impact should be maximized compared to the transmitted pulse energy and unwanted long range mechanical side effects, e.g. shock waves, axial jet components. In conclusion, these experimental results are of great importance for the prospective optimization of the ophthalmic surgical process with high-repetition rate fs-lasers.     

Bursting dynamics remarkably improve the performance of neural networks on liquid computing

Burst firings are functionally important behaviors displayed by neural circuits, which plays a primary role in reliable transmission of electrical signals for neuronal communication. However, with respect to the computational capability of neural networks, most of relevant studies are based on the spiking dynamics of individual neurons, while burst firing is seldom considered. In this paper, we carry out a comprehensive study to compare the performance of spiking and bursting dynamics on the capability of liquid computing, which is an effective approach for intelligent computation of neural networks. The results show that neural networks with bursting dynamic have much better computational performance than those with spiking dynamics, especially for complex computational tasks. Further analysis demonstrate that the fast firing pattern of bursting dynamics can obviously enhance the efficiency of synaptic integration from pre-neurons both temporally and spatially. This indicates that bursting dynamic can significantly enhance the complexity of network activity, implying its high efficiency in information processing.     

Chemical exchange effects during refocusing pulses in constant-time CPMG relaxation dispersion experiments

In the analysis of the constant-time Carr-Purcell-Meiboom-Gill (CT-CPMG) relaxation dispersion experiment, chemical exchange parameters, such as rate of exchange and population of the exchanging species, are typically optimized using equations that predict experimental relaxation rates recorded as a function of effective field strength. In this process, the effect of chemical exchange during the CPMG pulses is typically assumed to be the same as during the free-precession. This approximation may introduce systematic errors into the analysis of data because the number of CPMG pulses is incremented during the constant-time relaxation period, and the total pulse duration therefore varies as a function of the effective field strength. In order to estimate the size of such errors, we simulate the time-dependence of magnetization during the entire constant time period, explicitly taking into account the effect of the CPMG pulses on the spin relaxation rate. We show that in general the difference in the relaxation dispersion profile calculated using a practical pulse width from that calculated using an extremely short pulse width is small, but under certain circumstances can exceed 1 s^?1. The difference increases significantly when CPMG pulses are miscalibrated.