Dataset Information

Finite element analysis of the time-dependent Smoluchowski equation for acetylcholinesterase reaction rate calculations.

ABSTRACT: This article describes the numerical solution of the time-dependent Smoluchowski equation to study diffusion in biomolecular systems. Specifically, finite element methods have been developed to calculate ligand binding rate constants for large biomolecules. The resulting software has been validated and applied to the mouse acetylcholinesterase (mAChE) monomer and several tetramers. Rates for inhibitor binding to mAChE were calculated at various ionic strengths with several different time steps. Calculated rates show very good agreement with experimental and theoretical steady-state studies. Furthermore, these finite element methods require significantly fewer computational resources than existing particle-based Brownian dynamics methods and are robust for complicated geometries. The key finding of biological importance is that the rate accelerations of the monomeric and tetrameric mAChE that result from electrostatic steering are preserved under the non-steady-state conditions that are expected to occur in physiological circumstances.

SUBMITTER: Cheng Y

PROVIDER: S-EPMC1853150 | biostudies-literature | 2007 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Finite element analysis of the time-dependent Smoluchowski equation for acetylcholinesterase reaction rate calculations.

Cheng Yuhui Y Suen Jason K JK Zhang Deqiang D Bond Stephen D SD Zhang Yongjie Y Song Yuhua Y Baker Nathan A NA Bajaj Chandrajit L CL Holst Michael J MJ McCammon J Andrew JA

Biophysical journal 20070216 10

This article describes the numerical solution of the time-dependent Smoluchowski equation to study diffusion in biomolecular systems. Specifically, finite element methods have been developed to calculate ligand binding rate constants for large biomolecules. The resulting software has been validated and applied to the mouse acetylcholinesterase (mAChE) monomer and several tetramers. Rates for inhibitor binding to mAChE were calculated at various ionic strengths with several different time steps. ...[more]

PMID: 17307827

Similar Datasets

Project description:We consider the design of an effective and reliable adaptive finite element method (AFEM) for the nonlinear Poisson-Boltzmann equation (PBE). We first examine the two-term regularization technique for the continuous problem recently proposed by Chen, Holst, and Xu based on the removal of the singular electrostatic potential inside biomolecules; this technique made possible the development of the first complete solution and approximation theory for the Poisson-Boltzmann equation, the first provably convergent discretization, and also allowed for the development of a provably convergent AFEM. However, in practical implementation, this two-term regularization exhibits numerical instability. Therefore, we examine a variation of this regularization technique which can be shown to be less susceptible to such instability. We establish a priori estimates and other basic results for the continuous regularized problem, as well as for Galerkin finite element approximations. We show that the new approach produces regularized continuous and discrete problems with the same mathematical advantages of the original regularization. We then design an AFEM scheme for the new regularized problem, and show that the resulting AFEM scheme is accurate and reliable, by proving a contraction result for the error. This result, which is one of the first results of this type for nonlinear elliptic problems, is based on using continuous and discrete a priori L(∞) estimates to establish quasi-orthogonality. To provide a high-quality geometric model as input to the AFEM algorithm, we also describe a class of feature-preserving adaptive mesh generation algorithms designed specifically for constructing meshes of biomolecular structures, based on the intrinsic local structure tensor of the molecular surface. All of the algorithms described in the article are implemented in the Finite Element Toolkit (FETK), developed and maintained at UCSD. The stability advantages of the new regularization scheme are demonstrated with FETK through comparisons with the original regularization approach for a model problem. The convergence and accuracy of the overall AFEM algorithm is also illustrated by numerical approximation of electrostatic solvation energy for an insulin protein.

Project description:Radiofrequency ablation is the most common minimally invasive therapy used in the United States to treat hepatocellular carcinoma and liver metastases. The ability to perform real-time temperature imaging while a patient is undergoing ablation therapy may help reduce the high recurrence rates following ablation therapy. Ultrasound echo signals undergo time shifts with increasing temperature due to sound speed and thermal expansion, which are tracked using both 1D cross correlation and 2D block matching based speckle tracking methods. In this paper, we present a quantitative evaluation of the accuracy and precision of temperature estimation using the above algorithms on both simulated and experimental data. A finite element analysis simulation of radiofrequency ablation of hepatic tissue was developed. Finite element analysis provides a method to obtain the exact temperature distribution along with a mapping of the tissue displacement due to thermal expansion. These local displacement maps were combined with the displacement due to speed of sound changes and utilized to generate ultrasound radiofrequency frames at specified time increments over the entire ablation procedure. These echo signals provide an ideal test-bed to evaluate the performance of both speckle tracking methods, since the estimated temperature results can be compared directly to the exact finite element solution. Our results indicate that the 1D cross-correlation (CC) method underestimates the cumulative displacement by 0.20mm, while the underestimation with 2D block matching (BM) is about 0.14 mm after 360 s of ablation. The 1D method also overestimates the size of the ablated region by 5.4% when compared to 2.4% with the 2D method after 720 s of ablation. Hence 2D block matching provides better tracking of temperature variations when compared to the 1D cross-correlation method over the entire duration of the ablation procedure. In addition, results obtained using 1D cross-correlation diverge from the ideal finite element results after 7 min of ablation and for temperatures greater than 65 degrees C. In a similar manner, experimental results presented using a tissue-mimicking phantom also demonstrate that the maximum percent difference with 2D block matching was 5%, when compared to 31% with the 1D method over the 700 s heating duration on the phantom.

Project description:PurposeThe study presents an averaged anterior eye geometry model combined with a localised material model that is straightforward, appropriate and amenable for implementation in finite element (FE) modelling.MethodsBoth right and left eye profile data of 118 subjects (63 females and 55 males) aged 22-67 years (38.5 ± 7.6) were used to build an averaged geometry model. Parametric representation of the averaged geometry model was achieved through two polynomials dividing the eye into three smoothly connected volumes. This study utilised the collagen microstructure x-ray data of 6 ex-vivo healthy human eyes, 3 right eyes and 3 left eyes in pairs from 3 donors, 1 male and 2 females aged between 60 and 80 years, to build a localised element-specific material model for the eye.ResultsFitting the cornea and the posterior sclera sections to a 5th-order Zernike polynomial resulted in 21 coefficients. The averaged anterior eye geometry model recorded a limbus tangent angle of 37° at a radius of 6.6 mm from the corneal apex. In terms of material models, the difference between the stresses generated in the inflation simulation up to 15 mmHg in the ring-segmented material model and localised element-specific material model were significantly different (p < 0.001) with the ring-segmented material model recording average Von-Mises stress 0.0168 ± 0.0046 MPa and the localised element-specific material model recording average Von-Mises stress 0.0144 ± 0.0025 MPa.ConclusionsThe study illustrates an averaged geometry model of the anterior human eye that is easy to generate through two parametric equations. This model is combined with a localised material model that can be used either parametrically through a Zernike fitted polynomial or non-parametrically as a function of the azimuth angle and the elevation angle of the eye globe. Both averaged geometry and localised material models were built in a way that makes them easy to implement in FE analysis without additional computation cost compared to the limbal discontinuity so-called idealised eye geometry model or ring-segmented material model.

Project description:UnlabelledAmong the most fundamental questions in viral evolutionary biology are how fast viruses evolve and how evolutionary rates differ among viruses and fluctuate through time. Traditionally, viruses are loosely classed into two groups: slow-evolving DNA viruses and fast-evolving RNA viruses. As viral evolutionary rate estimates become more available, it appears that the rates are negatively correlated with the measurement timescales and that the boundary between the rates of DNA and RNA viruses might not be as clear as previously thought. In this study, we collected 396 viral evolutionary rate estimates across almost all viral genome types and replication strategies, and we examined their rate dynamics. We showed that the time-dependent rate phenomenon exists across multiple levels of viral taxonomy, from the Baltimore classification viral groups to genera. We also showed that, by taking the rate decay dynamics into account, a clear division between the rates of DNA and RNA viruses as well as reverse-transcribing viruses could be recovered. Surprisingly, despite large differences in their biology, our analyses suggested that the rate decay speed is independent of viral types and thus might be useful for better estimation of the evolutionary time scale of any virus. To illustrate this, we used our model to reestimate the evolutionary timescales of extant lentiviruses, which were previously suggested to be very young by standard phylogenetic analyses. Our analyses suggested that these viruses are millions of years old, in agreement with paleovirological evidence, and therefore, for the first time, reconciled molecular analyses of ancient and extant viruses.ImportanceThis work provides direct evidence that viral evolutionary rate estimates decay with their measurement timescales and that the rate decay speeds do not differ significantly among viruses despite the vast differences in their molecular features. After adjustment for the rate decay dynamics, the division between the rates of double-stranded DNA (dsDNA), single-stranded RNA (ssRNA), and ssDNA/reverse-transcribing viruses could be seen more clearly than before. Our results provide a guideline for further improvement of the molecular clock. As a demonstration of this, we used our model to reestimate the timescales of modern lentiviruses, which were previously thought to be very young, and concluded that they are millions of years old. This result matches the estimate from paleovirological analyses, thus bridging the gap between ancient and extant viral evolutionary studies.

Dataset Information

Finite element analysis of the time-dependent Smoluchowski equation for acetylcholinesterase reaction rate calculations.

Publications

Finite element analysis of the time-dependent Smoluchowski equation for acetylcholinesterase reaction rate calculations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets