A super-Gaussian Poisson-Boltzmann model for electrostatic solvation energy calculation: smooth dielectric distributions for protein cavities and in both water and vacuum states

dc.contributorZhao, Shan
dc.contributorHalpern, David
dc.contributorMoen, Kabe
dc.contributorAmes, Brendan
dc.contributorSu, Weihua
dc.contributor.advisorZhao, Shan
dc.contributor.authorHazra, Tania
dc.contributor.otherUniversity of Alabama Tuscaloosa
dc.descriptionElectronic Thesis or Dissertationen_US
dc.description.abstractCalculations of electrostatic potential and solvation energy of macromolecules are essential for understanding the mechanism of many biological processes. In the classical implicit solvent Poisson-Boltzmann (PB) model, the macromolecule and water are modeled as two-dielectric media with a sharp border. However, the dielectric property of interior cavities and ion-channels is difficult to model in a two-dielectric setting. In fact, whether there are water molecules or cavity-fluid inside a protein cavity remains to be an experimental challenge. Physically, this uncertainty affects the subsequent solvation free energy calculation. In order to compensate this uncertainty, a novel super-Gaussian dielectric PB model is introduced in this work, which devices an inhomogeneous dielectric distribution to represent the compactness of atoms and characterize empty cavities via a gap dielectric value. Moreover, the minimal molecular surface level set function is adopted so that the dielectric profile remains to be smooth when the protein is transfer from water phase to vacuum. A nice feature of this new model is that as the order of super-Gaussian function approaches the infinity, the dielectric distribution reduces a piecewise constant of the two-dielectric model. Mathematically, a simple effective dielectric constant analysis is introduced in this work to benchmark the dielectric model and select optimal parameter values. Computationally, a pseudo-time alternative direction implicit (ADI) algorithm is utilized for solving the super-Gaussian PB equation, which is found to be unconditionally stable in a smooth dielectric setting. Solvation free energy calculation of a Kirkwood sphere and various proteins is carried out to validate the super-Gaussian model and ADI algorithm. One macromolecule with both cavity-fluids and empty cavities is employed to demonstrate how the cavity uncertainty in protein structure can be bypassed through dielectric modeling in the biomolecular electrostatic analysis.en_US
dc.format.extent75 p.
dc.publisherUniversity of Alabama Libraries
dc.relation.hasversionborn digital
dc.relation.ispartofThe University of Alabama Electronic Theses and Dissertations
dc.relation.ispartofThe University of Alabama Libraries Digital Collections
dc.rightsAll rights reserved by the author unless otherwise indicated.en_US
dc.subjectApplied mathematics
dc.titleA super-Gaussian Poisson-Boltzmann model for electrostatic solvation energy calculation: smooth dielectric distributions for protein cavities and in both water and vacuum statesen_US
etdms.degree.departmentUniversity of Alabama. Department of Mathematics
etdms.degree.grantorThe University of Alabama
Original bundle
Now showing 1 - 1 of 1
Thumbnail Image
4.85 MB
Adobe Portable Document Format