Journal article
Numerical framework models of single proton conduction through gramicidin
Frontiers in bioscience, Vol.8(6), pp.s982-s991
09/01/2003
Handle:
https://hdl.handle.net/2376/101203
PMID: 12957852
Abstract
A framework model of single-proton conduction through gramicidin was previously designed to incorporate potentials of mean force and diffusion coefficients computed by the molecular dynamics simulations of Pom s and Roux (1). The resulting diffusion model was solved analytically using the lumped state approximation (LSA), allowing a detailed comparison to be made with conductance data from gramicidin A and two Trp--> Phe analogs (2). The comparison included a sensitivity analysis which required over 1 million current evaluations. A numerical method for constructing framework models is now introduced which involves finding the steady states of random walks using a trapezoid rule closely related to the rule for numerical integration. The method is described and then applied directly to the LSA. Convergence of the results to the analytical solution is seen as the number of random walk sites increase. The numerical method is then used to construct a more elaborate framework model which avoids the LSA. This is also in very good agreement with the analytical solution under the experimental conditions, confirming the accuracy of the LSA. The numerical method remains fast enough to allow an extensive comparison with conductance data.
Metrics
6 Record Views
Details
- Title
- Numerical framework models of single proton conduction through gramicidin
- Creators
- Mark F Schumaker - Department of Mathematics and Center for Integrated Biotechnology, Washington State University, Pullman WA 99164-3113, USA. schumaker@wsu.edu
- Publication Details
- Frontiers in bioscience, Vol.8(6), pp.s982-s991
- Academic Unit
- Mathematics and Statistics, Department of
- Publisher
- United States
- Identifiers
- 99900546532801842
- Language
- English
- Resource Type
- Journal article