Browse
Submit Data
Databases
API
Help

Dataset Information

16 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Finite element analysis of drug electrostatic diffusion: inhibition rate studies in N1 neuraminidase.

ABSTRACT: This article describes a numerical solution of the steady-state Poisson-Boltzmann-Smoluchowski (PBS) and Poisson-Nernst-Planck (PNP) equations to study diffusion in biomolecular systems. Specifically, finite element methods have been developed to calculate electrostatic interactions and ligand binding rate constants for large biomolecules. The resulting software has been validated and applied to the wild-type and several mutated avian influenza neurominidase crystal structures. The calculated rates show very good agreement with recent experimental 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 electrostatic steering plays the important role in the drug binding process of the neurominidase.

SUBMITTER: Cheng Y

PROVIDER: S-EPMC3107071 | biostudies-literature | 2009

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Finite element analysis of drug electrostatic diffusion: inhibition rate studies in N1 neuraminidase.

Cheng Yuhui Y Holst Michael J MJ McCammon J A JA

Pacific Symposium on Biocomputing. Pacific Symposium on Biocomputing 20090101

This article describes a numerical solution of the steady-state Poisson-Boltzmann-Smoluchowski (PBS) and Poisson-Nernst-Planck (PNP) equations to study diffusion in biomolecular systems. Specifically, finite element methods have been developed to calculate electrostatic interactions and ligand binding rate constants for large biomolecules. The resulting software has been validated and applied to the wild-type and several mutated avian influenza neurominidase crystal structures. The calculated ra ...[more]

PMID: 19209708

Similar Datasets

Tetrameric mouse acetylcholinesterase: continuum diffusion rate calculations by solving the steady-state Smoluchowski equation using finite element methods.

Project description:The tetramer is the most important form for acetylcholinesterase in physiological conditions, i.e., in the neuromuscular junction and the nervous system. It is important to study the diffusion of acetylcholine to the active sites of the tetrameric enzyme to understand the overall signal transduction process in these cellular components. Crystallographic studies revealed two different forms of tetramers, suggesting a flexible tetramer model for acetylcholinesterase. Using a recently developed finite element solver for the steady-state Smoluchowski equation, we have calculated the reaction rate for three mouse acetylcholinesterase tetramers using these two crystal structures and an intermediate structure as templates. Our results show that the reaction rates differ for different individual active sites in the compact tetramer crystal structure, and the rates are similar for different individual active sites in the other crystal structure and the intermediate structure. In the limit of zero salt, the reaction rates per active site for the tetramers are the same as that for the monomer, whereas at higher ionic strength, the rates per active site for the tetramers are approximately 67%-75% of the rate for the monomer. By analyzing the effect of electrostatic forces on ACh diffusion, we find that electrostatic forces play an even more important role for the tetramers than for the monomer. This study also shows that the finite element solver is well suited for solving the diffusion problem within complicated geometries.

| S-EPMC1305222 | biostudies-literature

Considerations for reporting finite element analysis studies in biomechanics.

Project description:Simulation-based medicine and the development of complex computer models of biological structures is becoming ubiquitous for advancing biomedical engineering and clinical research. Finite element analysis (FEA) has been widely used in the last few decades to understand and predict biomechanical phenomena. Modeling and simulation approaches in biomechanics are highly interdisciplinary, involving novice and skilled developers in all areas of biomedical engineering and biology. While recent advances in model development and simulation platforms offer a wide range of tools to investigators, the decision making process during modeling and simulation has become more opaque. Hence, reliability of such models used for medical decision making and for driving multiscale analysis comes into question. Establishing guidelines for model development and dissemination is a daunting task, particularly with the complex and convoluted models used in FEA. Nonetheless, if better reporting can be established, researchers will have a better understanding of a model's value and the potential for reusability through sharing will be bolstered. Thus, the goal of this document is to identify resources and considerate reporting parameters for FEA studies in biomechanics. These entail various levels of reporting parameters for model identification, model structure, simulation structure, verification, validation, and availability. While we recognize that it may not be possible to provide and detail all of the reporting considerations presented, it is possible to establish a level of confidence with selective use of these parameters. More detailed reporting, however, can establish an explicit outline of the decision-making process in simulation-based analysis for enhanced reproducibility, reusability, and sharing.

| S-EPMC3278509 | biostudies-literature

Hybrid finite element and Brownian dynamics method for diffusion-controlled reactions.

Project description:Diffusion is often the rate determining step in many biological processes. Currently, the two main computational methods for studying diffusion are stochastic methods, such as Brownian dynamics, and continuum methods, such as the finite element method. This paper proposes a new hybrid diffusion method that couples the strengths of each of these two methods. The method is derived for a general multidimensional system, and is presented using a basic test case for 1D linear and radially symmetric diffusion systems.

| S-EPMC3350536 | biostudies-other

Finite element solution of the steady-state Smoluchowski equation for rate constant calculations.

Project description:This article describes the development and implementation of algorithms to study diffusion in biomolecular systems using continuum mechanics equations. Specifically, finite element methods have been developed to solve the steady-state Smoluchowski equation to calculate ligand binding rate constants for large biomolecules. The resulting software has been validated and applied to mouse acetylcholinesterase. Rates for inhibitor binding to mAChE were calculated at various ionic strengths with several different reaction criteria. The calculated rates were compared with experimental data and show very good agreement when the correct reaction criterion is used. Additionally, these finite element methods require significantly less computational resources than existing particle-based Brownian dynamics methods.

| S-EPMC1304055 | biostudies-literature

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

Project description: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.

| S-EPMC1853150 | biostudies-literature

Role of secondary sialic acid binding sites in influenza N1 neuraminidase.

Project description:Within influenza viral particles, the intricate balance between host cell binding and sialic acid receptor destruction is carefully maintained by the hemagglutinin (HA) and neuraminidase (NA) glycoproteins, respectively. A major outstanding question in influenza biology is the function of a secondary sialic acid binding site on the NA enzyme. Through a series of Brownian dynamics (BD) simulations of the avian N1, human pandemic N2, and currently circulating pandemic (H1)N1 enzymes, we have probed the role of this secondary sialic acid binding site in the avian N1 subtype. Our results suggest that electrostatic interactions at the secondary and primary sites in avian NA may play a key role in the recognition process of the sialic acid receptors and catalytic efficiency of NA. This secondary site appears to facilitate the formation of complexes with the NA protein and the sialic acid receptors, as well as provide HA activity to a lesser extent. Moreover, this site is able to steer inhibitor binding as well, albeit with reduced capacity in N1, and may have potential implications for drug resistance or optimal inhibitor design. Although the secondary sialic acid binding site has previously been shown to be nonconserved in swine NA strains, our investigations of the currently circulating pandemic H1N1 strain of swine origin appears to have retained some of the key features of the secondary sialic acid binding site. Our results indicate possible lowered HA activity for this secondary sialic acid site, which may be an important event in the emergence of the current pandemic strain.

| S-EPMC2835425 | biostudies-literature

Pseudovirus-based neuraminidase inhibition assays reveal potential H5N1 drug-resistant mutations.

Project description:The use of antiviral drugs such as influenza neuraminidase (NA) inhibitors is a critical strategy to prevent and control flu pandemic, but this strategy faces the challenge of emerging drug-resistant strains. For a highly pathogenic avian influenza (HPAI) H5N1 virus, biosafety restrictions have significantly limited the efforts to monitor its drug responses and mechanisms involved. In this study, a rapid and biosafe assay based on NA pseudovirus was developed to study the resistance of HPAI H5N1 virus to NA inhibitor drugs. The H5N1 NA pseudovirus was comprehensively tested using oseltamivir-sensitive strains and their resistant mutants. Results were consistent with those in previous studies, in which live H5N1 viruses were used. Several oseltamivir-resistant mutations reported in human H1N1 were also identified to cause decreased oseltamivir sensitivity in H5N1 NA by using the H5N1 NA pseudovirus. Thus, H5N1 NA pseudoviruses could be used to monitor HPAI H5N1 drug resistance rapidly and safely.

| S-EPMC4875545 | biostudies-literature

Impact of calcium on N1 influenza neuraminidase dynamics and binding free energy.

Project description:The highly pathogenic influenza strains H5N1 and H1N1 are currently treated with inhibitors of the viral surface protein neuraminidase (N1). Crystal structures of N1 indicate a conserved, high affinity calcium binding site located near the active site. The specific role of this calcium in the enzyme mechanism is unknown, though it has been shown to be important for enzymatic activity and thermostability. We report molecular dynamics (MD) simulations of calcium-bound and calcium-free N1 complexes with the inhibitor oseltamivir (marketed as the drug Tamiflu), independently using both the AMBER FF99SB and GROMOS96 force fields, to give structural insight into calcium stabilization of key framework residues. Y347, which demonstrates similar sampling patterns in the simulations of both force fields, is implicated as an important N1 residue that can "clamp" the ligand into a favorable binding pose. Free energy perturbation and thermodynamic integration calculations, using two different force fields, support the importance of Y347 and indicate a +3 to +5 kcal/mol change in the binding free energy of oseltamivir in the absence of calcium. With the important role of structure-based drug design for neuraminidase inhibitors and the growing literature on emerging strains and subtypes, inclusion of this calcium for active site stability is particularly crucial for computational efforts such as homology modeling, virtual screening, and free energy methods.

| S-EPMC2902668 | biostudies-literature

Steel benchmark frames for structural analysis and validation studies: Finite element models and numerical simulation data.

Project description:A sample of twenty-two steel benchmark frames was formed and used to investigate the accuracy of a single increment predictor-corrector (SIPC) solution scheme; an approximate method of second-order elastic analysis. Each of the frames is based on a structure from the literature and, as a collection, the set includes a diverse assortment of practical planar geometries and a wide range of sensitivities to second-order effects. This data article presents the details of these frames, including finite element models, relevant nodal coordinates and element connectivity, and detailed information regarding member sizes, support conditions, and applied loading. In addition, this article presents simulation data obtained from testing the SIPC method using the benchmark frames, and assessing its accuracy and precision. Error analysis results, based on comparisons of joint displacements and member design moments simulated using the SIPC method with those obtained using the more exact and computationally expensive work-control (WC) method, are summarized. The finite element modelling and subsequent structural analysis utilized the software package MASTAN2, which provides user-friendly features to execute both SIPC and WC methods. A detailed description of these analysis methods and the algorithms used to generate data is provided in "Efficient geometric nonlinear elastic analysis for design of steel structures: Benchmark studies" Ziemian and Ziemian [1]. The benchmark frame models are more generally useful for any researcher interested in testing and validating structural analysis and design methods, and the simulation data allow for comparisons with the results of other proposed solution schemes.

| S-EPMC8637320 | biostudies-literature

Volumetric finite-element modelling of biological growth.

Project description:Differential growth is the driver of tissue morphogenesis in plants, and also plays a fundamental role in animal development. Although the contributions of growth to shape change have been captured through modelling tissue sheets or isotropic volumes, a framework for modelling both isotropic and anisotropic volumetric growth in three dimensions over large changes in size and shape has been lacking. Here, we describe an approach based on finite-element modelling of continuous volumetric structures, and apply it to a range of forms and growth patterns, providing mathematical validation for examples that admit analytic solution. We show that a major difference between sheet and bulk tissues is that the growth of bulk tissue is more constrained, reducing the possibility of tissue conflict resolution through deformations such as buckling. Tissue sheets or cylinders may be generated from bulk shapes through anisotropic specified growth, oriented by a polarity field. A second polarity field, orthogonal to the first, allows sheets with varying lengths and widths to be generated, as illustrated by the wide range of leaf shapes observed in nature. The framework we describe thus provides a key tool for developing hypotheses for plant morphogenesis and is also applicable to other tissues that deform through differential growth or contraction.

| S-EPMC6544983 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data