Biomedical subjects

Rob D Coalson

Publications and source records attributed to Rob D Coalson.

12 recordsLinked to original sources

PubMed2006

Condensed-phase relaxation of multilevel quantum systems. I. An exactly solvable model.

An analytically solvable model of multilevel condensed-phase quantum dynamics relevant to vibrational relaxation and electron transfer is presented. Exact solutions are derived for the reduced system density matrix dynamics of a degenerate N-level quantum system characterized by nearest-neighbor hopping and off-diagonal coupling (which is linear in the bath coordinates) to a harmonic oscillator bath. We demonstrate that for N> 2 the long-time steady-state system site occupation probabilities are not the same for all sites; that is, they are distributed in a non-Boltzmann manner, which depends on the initial conditions and the number of levels in the system. Although the system-bath Hamiltonian considered here is restricted in form, the availability of an exact solution enables us to study the model in all regions of an extensive parameter space.

Journal Article↗

PubMed2006

Condensed-phase relaxation of multilevel quantum systems. II. Comparison of path integral calculations and second-order relaxation theory for a nondegenerate three-level system.

An exactly solvable model of multisite condensed-phase vibrational relaxation was studied in Paper I (Peter, S.; Evans, D. G.; Coalson, R. D. J. Phys. Chem. B 2006, 110, 18758.), where it was shown that long-time steady-state site populations of a degenerate N-level system are not equal (hence, they are non-Boltzmann) and depend on the initial preparation of the system and the number of sites that it comprises. Here we consider a generalization of the model to the case of a nondegenerate three-level system coupled to a high-dimensional bath: such a model system has direct relevance to a large class of donor-bridge-acceptor electron transfer processes. Because the quantum dynamics of this system cannot be computed analytically, we compare numerically exact path integral calculations to the predictions of second-order time-local relaxation theory. For modest system-bath coupling strengths, the two sets of results are in excellent agreement. They show that non-Boltzmann long-time steady-state site populations are obtained when the level splitting is small but nonzero, whereas at larger values of the system bias (asymmetry) these populations become Boltzmann distributed.

Journal Article↗

PubMed2006

A rigid-body Newtonian propagation scheme based on instantaneous decomposition into rotation and translation blocks.

The rotation and translation block (RTB) method of Durand et al. [Biopolymers 34, 759 (1994)] and Tama et al. [Proteins 41, 1 (2000)] provides an appealing way to calculate low-frequency normal modes of large biomolecules by restricting the space of motions to exclude internal motions of preselected rigid fragments within the molecule. These fragments are modeled essentially as rigid bodies and the need to calculate high-frequency relative motions of the atoms that form them is obviated in a natural way. Here we extend the RTB approach into a method for computing the classical (Newtonian) dynamics of a biomolecule, or any large molecule, with effective rigid-body constraints applied to a prechosen set of internal molecular fragments. This method, to be termed RTB dynamics, is easy to implement, conserves the total energy of the system, does not require the construction of the matrix of second spatial derivatives of the potential-energy function (Hessian matrix), and can be used to compute the classical dynamics of a system moving in an arbitrary anharmonic force field. An elementary numerical application to signal propagation in the small membrane-bound polypeptide gramicidin-A is presented for illustration purposes.

Journal Article↗

PubMed2006

Diffusion constant of K+ inside Gramicidin A: a comparative study of four computational methods.

The local diffusion constant of K(+) inside the Gramicidin A (GA) channel has been calculated using four computational methods based on molecular dynamics (MD) simulations, specifically: Mean Square Displacement (MSD), Velocity Autocorrelation Function (VACF), Second Fluctuation Dissipation Theorem (SFDT) and analysis of the Generalized Langevin Equation for a Harmonic Oscillator (GLE-HO). All methods were first tested and compared for K(+) in bulk water--all predicted the correct diffusion constant. Inside GA, MSD and VACF methods were found to be unreliable because they are biased by the systematic force exerted by the membrane-channel system on the ion. SFDT and GLE-HO techniques properly unbias the influence of the systematic force on the diffusion properties and predicted a similar diffusion constant of K(+) inside GA, namely, ca. 10 times smaller than in the bulk. It was found that both SFDT and GLE-HO methods require extensive MD sampling on the order of tens of nanoseconds to predict a reliable diffusion constant of K(+) inside GA.

Cations, Monovalent↗

PubMed2005

Theoretical studies of the M2 transmembrane segment of the glycine receptor: models of the open pore structure and current-voltage characteristics.

The pentameric glycine receptor (GlyR), a member of the nicotinicoid superfamily of ligand-gated ion channels, is an inhibitory Cl(-) channel that is gated by glycine. Using recently published NMR data of the second transmembrane segment (M2) of the human alpha1 GlyR, structural models of pentameric assemblies embedded in a lipid bilayer were constructed using a combination of experimentally determined constraints coupled with all-atom energy minimization. Based on this structure of the pentameric M2 "pore", Brownian dynamics simulations of ion permeation through this putative conducting open state of the channel were carried out. Simulated I-V curves were in good agreement with published experimental current-voltage curves and the anion/cation permeability ratio, suggesting that our open-state model may be representative of the conducting channel of the full-length receptor. These studies also predicted regions of chloride occupancy and suggested residues critical to anion permeation. Calculations of the conductance of the cation-selective mutant A251E channel are also consistent with experimental data. In addition, both rotation and untilting of the pore helices of our model were found to be broadly consistent with closing of the channel, albeit at distinct regions that may reflect alternate gates of the receptor.

Algorithms↗

PubMed2005

An accurate and efficient empirical approach for calculating the dielectric self-energy and ion-ion pair potential in continuum models of biological ion channels.

This paper presents empirical formulas for calculating the dielectric self-energy and ion-ion pair interactions in cylindrical ion channels. The proposed approach can be extended to more complex channel structures, for example, (i) a "straight" channel with variable radius and (ii) a "curved" channel with constant radius. For calibration purposes, we compare results obtained based on the approximate effective potentials developed herein to exact electrostatic calculations obtained via the algorithm of Graf et al.: the agreement is satisfactory. A dynamic lattice Monte Carlo (DLMC) technique is used to further assess the accuracy and efficiency of the proposed empirical potentials. The concentration profiles and current-voltage curves produced with our simple empirical energy formulas are in excellent agreement with numerical results obtained using the algorithm of Graf et al., which calculates all relevant electrostatic forces exactly. The use of effective ion-ion potentials greatly reduces the computer memory required to perform DLMC ion permeation simulations in dielectrically inhomogeneous environments, thus enabling treatment of larger systems than can be handled by numerically exact techniques.

Algorithms↗

PubMed2005

Poisson-Nernst-Planck theory approach to the calculation of current through biological ion channels.

The Poisson-Nernst-Planck (PNP) theory of electro-diffusion is reviewed. Techniques for numerical solution of the three-dimensional PNP equations are summarized, and several illustrative applications to ion transport through protein channels are presented. Strengths and weaknesses of the theory are discussed, as well as attempts to improve it via increasingly realistic evaluation of the force acting on each ion due to the protein/membrane environment.

Animals↗

PubMed2004

Diffraction in crystalline colloidal-array photonic crystals.

We characterized the diffraction and crystal structure of a crystalline colloidal array (CCA) photonic crystal composed of 270 nm diameter polystyrene spheres which have a nearest neighbor spacing of approximately 540 nm. This CCA diffracts light in first order at approximately 1200 nm and shows strong diffraction in the visible spectral region from higher order planes. We quantitatively examined the relative diffraction intensities of the putative fcc (111), (200), (220), and (311) planes. Comparing these intensities to those calculated theoretically we find that the crystal structure is fcc with significant stacking faults. Essentially, no light transmits at the Bragg angle for the fcc (111) planes even through thin approximately 40 microm thick CCA. However, much of this light is diffusely scattered about the Bragg angle due to crystal imperfections. Significant transmission occurs from thin samples oriented at the Bragg condition for the fcc (200), (220), and (311) planes. We also observe moderately intense two-dimensional diffraction from the first few layers at the crystal surfaces. We also examined the sample thickness dependence of diffraction from CCA photonic crystals prepared from approximately 120 nm polystyrene spheres whose fcc (111) planes diffract in the visible spectral region. These experimental observations, aided by calculations based upon a simple but flexible model of light scattering from an arbitrary collection of colloidal spheres, make clear that fabrication of three-dimensional photonic band gap crystals will be challenged by crystal imperfections.

Biophysics↗

PubMed2004

Flexible polyelectrolyte simulations at the Poisson-Boltzmann level: a comparison of the kink-jump and multigrid configurational-bias Monte Carlo methods.

We present a new approach for simulating the motions of flexible polyelectrolyte chains based on the continuous kink-jump Monte Carlo technique coupled to a lattice field theory based calculation of the Poisson-Boltzmann (PB) electrostatic free energy "on the fly." This approach is compared to the configurational-bias Monte Carlo technique, in which the chains are grown on a lattice and the PB equation is solved for each configuration with a linear scaling multigrid method to obtain the many-body free energy. The two approaches are used to calculate end-to-end distances of charged polymer chains in solutions with varying ionic strengths and give similar numerical results. The configurational-bias Monte Carlo/multigrid PB method is found to be more efficient, while the kink-jump Monte Carlo method shows potential utility for simulating nonequilibrium polyelectrolyte dynamics.

Journal Article↗

PubMed2003

The wide-angle equation and its solution through the short-time iterative Lanczos method.

Properties of the wide-angle equation (WAEQ), a nonparaxial scalar wave equation used to propagate light through media characterized by inhomogeneous refractive-index profiles, are studied. In particular, it is shown that the WAEQ is not equivalent to the more complicated but more fundamental Helmholtz equation (HEQ) when the index of refraction profile depends on the position along the propagation axis. This includes all nonstraight waveguides. To study the quality of the WAEQ approximation, we present a novel method for computing solutions to the WAEQ. This method, based on a short-time iterative Lanczos (SIL) algorithm, can be applied directly to the full three-dimensional case, i.e., systems consisting of the propagation axis coordinate and two transverse coordinates. Furthermore, the SIL method avoids series-expansion procedures (e.g., Padé approximants) and thus convergence problems associated with such procedures. Detailed comparisons of solutions to the HEQ, WAEQ, and the paraxial equation (PEQ) are presented for two cases in which numerically exact solutions to the HEQ can be obtained by independent analysis, namely, (i) propagation in a uniform dielectric medium and (ii) propagation along a straight waveguide that has been tilted at an angle to the propagation axis. The quality of WAEQ and PEQ, compared with exact HEQ results, is investigated. Cases are found for which the WAEQ actually performs worse than the PEQ.

Journal Article↗

PubMed2003

The role of the dielectric barrier in narrow biological channels: a novel composite approach to modeling single-channel currents.

A composite continuum theory for calculating ion current through a protein channel of known structure is proposed, which incorporates information about the channel dynamics. The approach is utilized to predict current through the Gramicidin A ion channel, a narrow pore in which the applicability of conventional continuum theories is questionable. The proposed approach utilizes a modified version of Poisson-Nernst-Planck (PNP) theory, termed Potential-of-Mean-Force-Poisson-Nernst-Planck theory (PMFPNP), to compute ion currents. As in standard PNP, ion permeation is modeled as a continuum drift-diffusion process in a self-consistent electrostatic potential. In PMFPNP, however, information about the dynamic relaxation of the protein and the surrounding medium is incorporated into the model of ion permeation by including the free energy of inserting a single ion into the channel, i.e., the potential of mean force along the permeation pathway. In this way the dynamic flexibility of the channel environment is approximately accounted for. The PMF profile of the ion along the Gramicidin A channel is obtained by combining an equilibrium molecular dynamics (MD) simulation that samples dynamic protein configurations when an ion resides at a particular location in the channel with a continuum electrostatics calculation of the free energy. The diffusion coefficient of a potassium ion within the channel is also calculated using the MD trajectory. Therefore, except for a reasonable choice of dielectric constants, no direct fitting parameters enter into this model. The results of our study reveal that the channel response to the permeating ion produces significant electrostatic stabilization of the ion inside the channel. The dielectric self-energy of the ion remains essentially unchanged in the course of the MD simulation, indicating that no substantial changes in the protein geometry occur as the ion passes through it. Also, the model accounts for the experimentally observed saturation of ion current with increase of the electrolyte concentration, in contrast to the predictions of standard PNP theory.

Cell Membrane Permeability↗

PubMedDate not supplied

An introduction to simulation and visualization of biological systems at multiple scales: a summer training program for interdisciplinary research.

Advances in biomedical research require a new generation of researchers having a strong background in both the life and physical sciences and a knowledge of computational, mathematical, and engineering tools for tackling biological problems. The NIH-NSF Bioengineering and Bioinformatics Summer Institute at the University of Pittsburgh (BBSI @ Pitt; www.ccbb.pitt.edu/bbsi) is a multi-institutional 10-week summer program hosted by the University of Pittsburgh, Duquesne University, the Pittsburgh Supercomputing Center, and Carnegie Mellon University, and is one of nine Institutes throughout the nation currently participating in the NIH-NSF program. Each BBSI focuses on a different area; the BBSI @ Pitt, entitled "Simulation and Computer Visualization of Biological Systems at Multiple Scales", focuses on computational and mathematical approaches to understanding the complex machinery of molecular-to-cellular systems at three levels, namely, molecular, subcellular (microphysiological), and cellular. We present here an overview of the BBSI @ Pitt, the objectives and focus of the program, and a description of the didactic training activities that distinguish it from other traditional summer research programs. Furthermore, we also report several challenges that have been identified in implementing such an interdisciplinary program that brings together students from diverse academic programs for a limited period of time. These challenges notwithstanding, presenting an integrative view of molecular-to-system analytical models has introduced these students to the field of computational biology and has allowed them to make an informed decision regarding their future career prospects.

Computational Biology↗