Differential stability of beta-sheets and alpha-helices in beta-lactamase: a high temperature molecular dynamics study of unfolding intermediates

Beta-Lactamase, which catalyzes beta-lactam antibiotics, is prototypical of large alpha/beta proteins with a scaffolding formed by strong noncovalent interactions. Experimentally, the enzyme is well characterized, and intermediates that are slightly less compact and having nearly the same content of secondary structure have been identified in the folding pathway. In the present study, high temperature molecular dynamics simulations have been carried out on the native enzyme in solution. Analysis of these results in terms of root mean square fluctuations in cartesian and [phi, psi] space, backbone dihedral angles and secondary structural hydrogen bonds forms the basis for an investigation of the topology of partially unfolded states of beta-lactamase. A differential stability has been observed for alpha-helices and beta-sheets upon thermal denaturation to putative unfolding intermediates. These observations contribute to an understanding of the folding/unfolding processes of beta-lactamases in particular, and other alpha/beta proteins in general.

Transverse plane wave analysis of short elliptical chamber mufflers: An analytical approach

Short elliptical chamber mufflers are used often in the modern day automotive exhaust systems. The acoustic analysis of such short chamber mufflers is facilitated by considering a transverse plane wave propagation model along the major axis up to the low frequency limit. The one dimensional differential equation governing the transverse plane wave propagation in such short chambers is solved using the segmentation approaches which are inherently numerical schemes, wherein the transfer matrix relating the upstream state variables to the downstream variables is obtained. Analytical solution of the transverse plane wave model used to analyze such short chambers has not been reported in the literature so far. This present work is thus an attempt to fill up this lacuna, whereby Frobenius solution of the differential equation governing the transverse plane wave propagation is obtained. By taking a sufficient number of terms of the infinite series, an approximate analytical solution so obtained shows good convergence up to about 1300 Hz and also covers most of the range of muffler dimensions used in practice. The transmission loss (TL) performance of the muffler configurations computed by this analytical approach agrees excellently with that computed by the Matrizant approach used earlier by the authors, thereby offering a faster and more elegant alternate method to analyze short elliptical muffler configurations. (C) 2010 Elsevier Ltd. All rights reserved.

Coupled electrostatic-elastic analysis for topology optimization using material interpolation

In this paper, we present a novel analytical formulation for the coupled partial differential equations governing electrostatically actuated constrained elastic structures of inhomogeneous material composition. We also present a computationally efficient numerical framework for solving the coupled equations over a reference domain with a fixed finite-element mesh. This serves two purposes: (i) a series of problems with varying geometries and piece-wise homogeneous and/or inhomogeneous material distribution can be solved with a single pre-processing step, (ii) topology optimization methods can be easily implemented by interpolating the material at each point in the reference domain from a void to a dielectric or a conductor. This is attained by considering the steady-state electrical current conduction equation with a `leaky capacitor' model instead of the usual electrostatic equation. This formulation is amenable for both static and transient problems in the elastic domain coupled with the quasi-electrostatic electric field. The procedure is numerically implemented on the COMSOL Multiphysics (R) platform using the weak variational form of the governing equations. Examples have been presented to show the accuracy and versatility of the scheme. The accuracy of the scheme is validated for the special case of piece-wise homogeneous material in the limit of the leaky-capacitor model approaching the ideal case.

Surface-ray tracing on hybrid surfaces of revolution for UTD mutual coupling analysis

An angle invariance property based on Hertz's principle of particle dynamics is employed to facilitate the surface-ray tracing on nondevelopable hybrid quadric surfaces of revolution (h-QUASOR's). This property, when used in conjunction with a Geodesic Constant Method, yields analytical expressions for all the ray-parameters required in the UTD formulation. Differential geometrical considerations require that some of the ray-parameters (defined heuristically in the UTD for the canonical convex surfaces) be modified before the UTD can be applied to such hybrid surfaces. Mutual coupling results for finite-dimensional slots have been presented as an example on a satellite launch vehicle modeled by general paraboloid of revolution and right circular cylinder.

Analysis of a conformation-specific epitope of the alpha subunit of human chorionic gonadotropin: Study using monoclonal antibody probes

Monoclonal antibodies (MAbs) have been used extensively for identification of sequence-specific epitopes using either the ELISA or/and IRMA methods, However, attempts to use MAbs for identification of conformation-specific epitopes have been very few as they are considered very labile. We have investigated the stability of conformation-specific epitopes of human chorionic gonadotropin (hCG) using a quantitative solid-phase radioimmnunoassay (SPRIA) technique. Several epitopes are stable to mild modification (chemical and proteolytic) conditions, and epitopes show differential stability for these modifications. Based on these observations, a monoclonal antibody (MAb 16) for an a-subunit-specific epitope of hCG has been used to monitor changes at the epitopic site (identified as epitope 16) on modification of hCG, using SPRIA with immobilized MAb 16. Modifications of amino groups, hydroxyl group of tyrosine as well as carboxyl group of Asp/Glu all bring about sufficient changes in the epitope integrity. Peptide bond hydrolysis at lysine residues damages the epitope, but not at arginine residues, Hydrolysis at tyrosine does not affect the epitope, though modification of the side-chain of tyrosine inactivates the epitope. Destruction of the epitope occurs on reduction of the disulphide bonds. Partial retention of the epitope activity is seen on modification of carboxyl or the epsilon-amino groups of lysine. Based on these results four to six amino acids have been identified to be at the epitopic site, and the data suggest that two peptide segments are brought together by the disulphide bond Cys10-Cys60 to form the epitope.

Analysis of a long thick orthotropic circular cylindrical shell panel

An exact three-dimensional elasticity solution has been obtained for an infinitely long, thick transversely isotropic circular cylindrical shell panel, simply supported along the longitudinal edges and subjected to a radial patch load. Using a set of three displacement functions, the boundary value problem is reduced to Bessel's differential equation. Numerical results are presented for different thickness to mean radius ratios and semicentral angles of the shell panel. Classical and first-order shear deformation orthotropic shell theories have been examined in comparison with the present elasticity solution.

A note on a class of singular integro-differential equations

A simplified analysis is employed to handle a class of singular integro-differential equations for their solutions

Stability analysis of stochastically parametered nonconservative columns

Nonconservatively loaded columns. which have stochastically distributed material property values and stochastic loadings in space are considered. Young's modulus and mass density are treated to constitute random fields. The support stiffness coefficient and tip follower load are considered to be random variables. The fluctuations of external and distributed loadings are considered to constitute a random field. The variational formulation is adopted to get the differential equation and boundary conditions. The non self-adjoint operators are used at the boundary of the regularity domain. The statistics of vibration frequencies and modes are obtained using the standard perturbation method, by treating the fluctuations to be stochastic perturbations. Linear dependence of vibration and stability parameters over property value fluctuations and loading fluctuations are assumed. Bounds for the statistics of vibration frequencies are obtained. The critical load is first evaluated for the averaged problem and the corresponding eigenvalue statistics are sought. Then, the frequency equation is employed to transform the eigenvalue statistics to critical load statistics. Specialization of the general procedure to Beck, Leipholz and Pfluger columns is carried out. For Pfluger column, nonlinear transformations are avoided by directly expressing the critical load statistics in terms of input variable statistics.

Violin string shape functions for finite element analysis of rotating Timoshenko beams

Violin strings are relatively short and stiff and are well modeled by Timoshenko beam theory. We use the static part of the homogeneous differential equation of violin strings to obtain new shape functions for the finite element analysis of rotating Timoshenko beams. For deriving the shape functions, the rotating beam is considered as a sequence of violin strings. The violin string shape functions depend on rotation speed and element position along the beam length and account for centrifugal stiffening effects as well as rotary inertia and shear deformation on dynamic characteristics of rotating Timoshenko beams. Numerical results show that the violin string basis functions perform much better than the conventional polynomials at high rotation speeds and are thus useful for turbo machine applications. (C) 2011 Elsevier B.V. All rights reserved.

Physics based basis function for vibration analysis of high speed rotating beams

The natural frequencies of continuous systems depend on the governing partial differential equation and can be numerically estimated using the finite element method. The accuracy and convergence of the finite element method depends on the choice of basis functions. A basis function will generally perform better if it is closely linked to the problem physics. The stiffness matrix is the same for either static or dynamic loading, hence the basis function can be chosen such that it satisfies the static part of the governing differential equation. However, in the case of a rotating beam, an exact closed form solution for the static part of the governing differential equation is not known. In this paper, we try to find an approximate solution for the static part of the governing differential equation for an uniform rotating beam. The error resulting from the approximation is minimized to generate relations between the constants assumed in the solution. This new function is used as a basis function which gives rise to shape functions which depend on position of the element in the beam, material, geometric properties and rotational speed of the beam. The results of finite element analysis with the new basis functions are verified with published literature for uniform and tapered rotating beams under different boundary conditions. Numerical results clearly show the advantage of the current approach at high rotation speeds with a reduction of 10 to 33% in the degrees of freedom required for convergence of the first five modes to four decimal places for an uniform rotating cantilever beam.

Performance analysis of the random early detection algorithm

In this article we consider a finite queue with its arrivals controlled by the random early detection algorithm. This is one of the most prominent congestion avoidance schemes in the Internet routers. The aggregate arrival stream from the population of transmission control protocol sources is locally considered stationary renewal or Markov modulated Poisson process with general packet length distribution. We study the exact dynamics of this queue and provide the stability and the rates of convergence to the stationary distribution and obtain the packet loss probability and the waiting time distribution. Then we extend these results to a two traffic class case with each arrival stream renewal. However, computing the performance indices for this system becomes computationally prohibitive. Thus, in the latter half of the article, we approximate the dynamics of the average queue length process asymptotically via an ordinary differential equation. We estimate the error term via a diffusion approximation. We use these results to obtain approximate transient and stationary performance of the system. Finally, we provide some computational examples to show the accuracy of these approximations.

A new beam finite element for the analysis of functionally graded materials

A new beam element is developed to study the thermoelastic behavior of functionally graded beam structures. The element is based on the first-order shear deformation theory and it accounts for varying elastic and thermal properties along its thickness. The exact solution of static part of the governing differential equations is used to construct interpolating polynomials for the element formulation. Consequently, the stiffness matrix has super-convergent property and the element is free of shear locking. Both exponential and power-law variations of material property distribution are used to examine different stress variations. Static, free vibration and wave propagation problems are considered to highlight the behavioral difference of functionally graded material beam with pure metal or pure ceramic beams. (C) 2003 Elsevier Science Ltd. All rights reserved.

Analysis of an LMS linear equalizer for fading channels in decision directed mode

We consider a time varying wireless fading channel, equalized by an LMS linear equalizer in decision directed mode (DD-LMS-LE). We study how well this equalizer tracks the optimal Wiener equalizer. Initially we study a fixed channel.For a fixed channel, we obtain the existence of DD attractors near the Wiener filter at high SNRs using an ODE (Ordinary Differential Equation) approximating the DD-LMS-LE. We also show, via examples, that the DD attractors may not be close to the Wiener filters at low SNRs. Next we study a time varying fading channel modeled by an Auto-regressive (AR) process of order 2. The DD-LMS equalizer and the AR process are jointly approximated by the solution of a system of ODEs. We show via examples that the LMS equalizer ODE show tracks the ODE corresponding to the instantaneous Wiener filter when the SNR is high. This may not happen at low SNRs.

Tracking analysis of an LMS decision feedback equalizer for a wireless channel

We consider a time varying wireless fading channel, equalized by an LMS Decision Feedback equalizer (DFE). We study how well this equalizer tracks the optimal MMSEDFE (Wiener) equalizer. We model the channel by an Autoregressive (AR) process. Then the LMS equalizer and the AR process are jointly approximated by the solution of a system of ODEs (ordinary differential equations). Using these ODEs, we show via some examples that the LMS equalizer moves close to the instantaneous Wiener filter after initial transience. We also compare the LMS equalizer with the instantaneous optimal DFE (the commonly used Wiener filter) designed assuming perfect previous decisions and computed using perfect channel estimate (we will call it as IDFE). We show that the LMS equalizer outperforms the IDFE almost all the time after initial transience.

On Convergence of Differential Evolution Over a Class of Continuous Functions With Unique Global Optimum

Differential evolution (DE) is arguably one of the most powerful stochastic real-parameter optimization algorithms of current interest. Since its inception in the mid 1990s, DE has been finding many successful applications in real-world optimization problems from diverse domains of science and engineering. This paper takes a first significant step toward the convergence analysis of a canonical DE (DE/rand/1/bin) algorithm. It first deduces a time-recursive relationship for the probability density function (PDF) of the trial solutions, taking into consideration the DE-type mutation, crossover, and selection mechanisms. Then, by applying the concepts of Lyapunov stability theorems, it shows that as time approaches infinity, the PDF of the trial solutions concentrates narrowly around the global optimum of the objective function, assuming the shape of a Dirac delta distribution. Asymptotic convergence behavior of the population PDF is established by constructing a Lyapunov functional based on the PDF and showing that it monotonically decreases with time. The analysis is applicable to a class of continuous and real-valued objective functions that possesses a unique global optimum (but may have multiple local optima). Theoretical results have been substantiated with relevant computer simulations.