Potential energy landscape and robustness of a gene regulatory network: Toggle switch

Finding a multidimensional potential landscape is the key for addressing important global issues, such as the robustness of cellular networks. We have uncovered the underlying potential energy landscape of a simple gene regulatory network: a toggle switch. This was realized by explicitly constructing the steady state probability of the gene switch in the protein concentration space in the presence of the intrinsic statistical fluctuations due to the small number of proteins in the cell. We explored the global phase space for the system. We found that the protein synthesis rate and the unbinding rate of proteins to the gene were small relative to the protein degradation rate; the gene switch is monostable with only one stable basin of attraction. When both the protein synthesis rate and the unbinding rate of proteins to the gene are large compared with the protein degradation rate, two global basins of attraction emerge for a toggle switch. These basins correspond to the biologically stable functional states. The potential energy barrier between the two basins determines the time scale of conversion from one to the other. We found as the protein synthesis rate and protein unbinding rate to the gene relative to the protein degradation rate became larger, the potential energy barrier became larger. This also corresponded to systems with less noise or the fluctuations on the protein numbers.

Dynamics and intrinsic statistical fluctuations of a gene switch

We studied a simple gene regulatory network, the toggle switch. Specifically, we examined the means and statistical fluctuations in numbers of proteins. We found that when omega, the ratio of rates of protein-gene unbinding to protein degradation, was between similar to 10(-3) and similar to 10, the fluctuations were much larger than those we would have expected from Poisson statistics. In addition, we examined characteristic time values for system relaxation and found both that they increased with omega and that they have significant phase transition effects, with a secondary time scale appearing near the boundary between bistable and other phases. Last, we discuss the bistability of the toggle switch.

Context and conformation dictate function of a transcription antitermination switch

In bacteriophage, transcription elongation is regulated by the N protein, which binds a nascent mRNA hairpin ( termed boxB) and enables RNA polymerase to read through distal terminators. We have examined the structure, energetics and in vivo function of a number of N boxB complexes derived from in vitro protein selection. Trp18 fully stacks on the RNA loop in the wild-type structure, and can become partially or completely unstacked when the sequence context is changed three or four residues away, resulting in a recognition interface in which the best binding residues depend on the sequence context. Notably, in vivo antitermination activity correlates with the presence of a stacked aromatic residue at position 18, but not with N boxB binding affinity. Our work demonstrates that RNA polymerase responds to subtle conformational changes in cis-acting regulatory complexes and that approximation of components is not sufficient to generate a fully functional transcription switch.

Electromechanical Model of RF MEMS Switches

With the recent rapid growth of Radio Frequency Micro-Electro-Mechanical Systems (RF MEMS) switches, there has developed an emergent requirement for more accurate theoretical models to predict their electromechanical behaviors. Many parameters exist in the analysis of the behavior of the switch, and it is inconvenient for further study. In this paper, an improved model is introduced, considering simultaneously axial stress, residual stress, and fringing-field effect of the fixed-fixed bridge structure. To avoid any unnecessary repetitive model tests and numerical simulation for RF MEMS switches, some dimensionless numbers are derived by making governing equation dimensionless. The electromechanical behavior of the fixed-fixed bridge structure of RF MEMS switches is totally determined by these dimensionless numbers.

Influence of van der Waals and Casimir Forces on Electrostatic Torsional Actuators

The influence of van der Waals (vdW) and Casimir forces on the stability of the electrostatic torsional nanoelectromechanical systems (NEMS) actuators is analyzed in the paper. With the consideration of vdW and Casimir effects, the dependence of the critical tilting angle and pull-in voltage on the sizes of structure is investigated. And the influence of vdW torque is compared with that of Casimir torque. The modified coefficients of vdW and Casimir torques on the pull-in voltage are, respectively, calculated. When the gap is sufficiently small, pull-in can still take place with arbitrary small angle perturbation because of the action of vdW and Casimir torques even if there is not electrostatic torque. And the critical pull-in gaps for two cases are, respectively, derived.

Casimir Effect on the Pull-in Parameters of Nanometer Switches

Casimir effect on the critical pull-in gap and pull-in voltage of nanoelectromechanical switches is studied. An approximate analytical expression of the critical pull-in gap with Casimir force is presented by the perturbation theory. The corresponding pull-in parameters are computed numerically, from which one can notice the nonlinear effect of Casimir force on the pull-in parameters. The detachment length has been presented, which increases with increasing thickness of the beam.

Formation Of Dendritic Nanostructures In Pyrex Glass Anodically Bonded To Silicon Coated With An Aluminum Thin Film

Anodic bonding with thin films of metal or alloy as an intermediate layer, finds increasing applications in micro/nanoelectromechanical systems. At the bonding temperature of 350 degrees C, voltage of 400 V, and 30 min duration, the anodic bonding is completed between Pyrex glass and crystalline silicon coated with an aluminum thin film with a thickness comprised between 50 and 230 nm. Sodium-depleted layers and dendritic nanostructures were observed in Pyrex 7740 glass adjacent to the bonding interface. The sodium depletion width does not increase remarkably with the thickness of aluminum film. The dendritic nanostructures result from aluminum diffusion into the Pyrex glass. This experimental research is expected to enhance the understanding of how the depletion layer and dendritic nanostructures affect the quality of anodic bonding. (C) 2007 Elsevier B.V. All rights reserved.

Fractal Pattern Formation In Anodic Bonding Of Pyrex Glass/Al/Si

Anodic bonding of Pyrex glass/Al/Si is an important bonding technique in micro/nanoelectromechanical systems (MEMS/NEMS) industry. The anodic bonding of Pyrex 7740 glass/Aluminum film/Silicon is completed at the temperature from 300 degrees C to 375 degrees C with a bonding voltage between 150 V and 450 V. The fractal patterns are formed in the intermediate Al thin film. This pattern has the fractal dimension of the typical two-dimensional diffusion-limited aggregation (2D DLA) process, and the fractal dimension is around 1.7. The fractal patterns consist of Al and Si crystalline grains, and their occurrences are due to the limited diffusion, aggregation, and crystallization of Si and Al atoms in the intermediate Al layers. The formation of the fractal pattern is helpful to enhance the bonding strength between the Pyrex 7740 glass and the aluminum thin film coated on the crystal silicon substrates.

Energy Principle And Nonlinear Electric-Mechanical Behavior Of Ferroelectric Ceramics

Many experimental observations have shown that a single domain in a ferroelectric material switches by progressive movement of domain walls, driven by a combination of electric field and stress. The mechanism of the domain switch involves the following steps: initially, the domain has a uniform spontaneous polarization; new domains with the reverse polarization direction nucleate, mainly at the surface, and grow though the crystal thickness; the new domain expands sideways as a new domain continues to form; finally, the domain switch coalesces to complete the polarization reversal. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of the ferroelectric material and used to study the nonlinear constitutive behavior of a ferroelectric body in this paper. The principle of stationary total potential energy is put forward in which the basic unknown quantities are the displacement u(i), electric displacement D-i and volume fraction rho(I) of the domain switching for the variant I. The mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total potential energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion established by Hwang et al. [ 1]. Based on the domain switching criterion, a set of linear algebraic equations for determining the volume fraction rho(I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. If the volume fraction rho(I) of domain switching for each domain is prescribed, the unknown displacement and electric potential can be obtained based on the conventional finite element procedure. It is assumed that a domain switches if the reduction in potential energy exceeds a critical energy barrier. According to the experimental results, the energy barrier will strengthen when the volume fraction of the domain switching increases. The external mechanical and electric loads are increased step by step. The volume fraction rho(I) of domain switching for each element obtained from the last loading step is used as input to the constitutive equations. Then the strain and electric fields are calculated based on the conventional finite element procedure. The finite element analysis is carried out on the specimens subjected to uniaxial coupling stress and electric field. Numerical results and available experimental data are compared and discussed. The present theoretic prediction agrees reasonably with the experimental results.

Long-range automaton models of earthquakes: Power-law accelerations, correlation evolution, and mode-switching

We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r(-p) dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.

Numerical exploration of new features on three-dimensional transition of a cylinder wake

The three-dimensional transition of the wake flow behind a circular cylinder is studied in detail by direct numerical simulations using 3D incompressible N-S equations for Reynolds number ranging from 200 to 300. New features and vortex dynamics of the 3D transition of the wake are found and investigated. At Re = 200, the flow pattern is characterized by mode A instability. However, the spanwise characteristic length of the cylinder determines the transition features. Particularly for the specific spanwise characteristic length linear stable mode may dominate the wake in place of mode A and determine the spanwise phase difference of the primary vortices shedding. At Re = 250 and 300 it is found that the streamwise vortices evolve into a new type of mode - "dual vortex pair mode" downstream. The streamwise vortex structures switch among mode A, mode B and dual vortex pair mode from near wake to downstream wake. At Re = 250, an independent low frequency f(m) in addition to the vortex shedding frequency f(s) is identified. Frequency coupling between f(m) and f(s) occurs. These result in the irregularity of the temporal signals and become a key feature in the transition of the wake. Based on the formation analysis of the streamwise vorticity in the vicinity of cylinder, it is suggested that mode A is caused by the emergence of the spanwise velocity due to three dimensionality of the incoming flow past the cylinder. Energy distribution on various wave numbers and the frequency variation in the wake are also described.

Electric characteristics of pseudo-spark during breakdown phase

An analysis of the time-dependent resistive voltage and power deposition during the breakdown phase of pseudo-spark is presented. The voltage and current were measured by specially designed low-inductance capacitive voltage divider and current measuring resistor. The measured waveforms of voltage and current are digitized and processed by a computer program to remove the inductive component, so as to obtain resistive voltage and power deposition. The influence of pressure, cathode geometry and charging voltage of storage capacitors on the electrical properties in the breakdown phase are investigated. The results suggest that the breakdown phase of pseudo-spark consists of three stages. The first stage is mainly hollow cathode discharge. In the second stage, field-enhanced thermionic emission takes place, resulting in a fast voltage drop and sharp rise of discharge current. The third stage of discharge depends simply on the parameters of the discharge circuit.

The permanence of density waves of galaxy and the numerical investigation of waser mechanism

This paper deals with in detail the permanence of the spiral structure of galaxies andthe characters of waser mechanism. A simplified model of galaxy is adopted. Variousdynamical characters of density waves are studied using numerical calculation method. Theresults verify very well the switch character f waser and the tunnel effect of density wavesat the potential barrier of corotation circle as is shown in a previous work of the author.

A Proteomic Study on Cell Cycle Progression of Endothelium Exposed to Tumor Conditioned Medium and the Possible Role of Cyclin D-1/E

This study was designed to comprehensively analyze the differential expression of proteins from human umbilical vein endothelial cells (HUVECs) exposed to tumor conditioned medium (TCM) and to identify the key regulator in the cell cycle progression. The HUVECs were exposed to TCM from breast carcinoma cell line MDA-MB-231, then their cell cycle distribution was measured by flow cytometer (FCM). The role of protein in cell cycle progression was detected via two-dimensional polyacrylamide gel electrophoresis (2-DE) and western blotting. Following the stimulation of TCM, HUVECs showed a more cells in the S phase than did the negative control group (ECGF-free medium with 20% FBS), but the HUVECs' level was similar to the positive control group (medium with 25 mug/ml ECGF and 20% FBS). Increased expression of cyclin D-1/E and some changes in other related proteins occurred after incubation with TCM. From our results, we can conclude that breast carcinoma cell line MDA-MB-231 may secrete soluble pro-angiogenic factors that induce the HUVEC angiogenic switch, during which the expression of cell cycle regulator cyclin D-1/E increases and related proteins play an important role in this process.

Research on Q-Switched Nd : YAG laser engraving

IN this paper, the engraving process with Q-Switched Nd:YAG laser is investigated. High power density is the pre- requisition to vapor materials, and high repetition rate makes the engraving process highly efficient. An acousto- optic Q-Switch is applied in the cavity of CW 200 W Nd:YAG laser to achieve the high peak power density and the high pulse repetition rate. Different shape craters are formed in a patterned structure on the material surface when the laser beam irradiates on it by controlling power density, pulse repetition rate, pulse quantity and pulse interval. In addition, assisting oxygen gas is used for not only improving combustion to deepen the craters but also removing the plasma that generated on the top of craters. Off-focus length classified as negative and positive has a substantial effect on crater diameters. According to the message of rotating angle positions from material to be engraved and the information of graph pixels from computer, a special graph is imparted to the material by integrating the Q- Switched Nd:YAG laser with the computer graph manipulation and the numerically controlled worktable. The crater diameter depends on laser beam divergence and laser focal length. The crater diameter changes from 50 micrometers to 300 micrometers , and the maximum of crater depth reaches one millimeter.