Weihua Geng

Weihua Geng

In this paper, we present a parallel higher-order boundary integral method to solve the linear Poisson-Boltzmann (PB) equation. In our method, a well-posed boundary integral formulation is used to ensure the fast convergence of Krylov subspace linear solver such as GMRES. The molecular surfaces are first discretized with flat triangles and then converted to curved triangles with the assistance of normal information at vertices. To maintain the desired accuracy, four-point Gauss-Radau quadratures are used on regular triangles and sixteen-point Gauss-Legendre quadratures together with regularization transformations are applied on singular triangles. To speed up our method, we take advantage of the embarrassingly parallel feature of boundary integral formulation, and parallelize the schemes with the message passing interface (MPI) implementation. Numerical tests show significantly improved accuracy and convergence of the proposed higher-order boundary integral Poisson-Boltzmann (HOBI-PB) solver compared with boundary integral PB solver using often-seen centroid collocation on flat triangles. The higher-order accuracy results achieved by present method are important to sensitive solvation analysis of biomolecules, particularly when accurate electrostatic surface potentials are critical in the molecular simulation. In addition, the higher-order boundary integral schemes presented here and their associated parallelization potentially can be applied to solving boundary integral equations in a general sense.

Elyssa Sliheet, Md Abu Talha, Weihua Geng

In this project, we provide a deep-learning neural network (DNN) based biophysics model to predict protein properties. The model uses multi-scale and uniform topological and electrostatic features generated with protein structural information and force field, which governs the molecular mechanics. The topological features are generated using the element specified persistent homology (ESPH) while the electrostatic features are fast computed using a Cartesian treecode. These features are uniform in number for proteins with various sizes thus the broadly available protein structure database can be used in training the network. These features are also multi-scale thus the resolution and computational cost can be balanced by the users. The machine learning simulation on over 4000 protein structures shows the efficiency and fidelity of these features in representing the protein structure and force field for the predication of their biophysical properties such as electrostatic solvation energy. Tests on topological or electrostatic features alone and the combination of both showed the optimal performance when both features are used. This model shows its potential as a general tool in assisting biophysical properties and function prediction for the broad biomolecules using data from both theoretical computing and experiments.

Jiahui Chen, Weihua Geng, Guo-Wei Wei

Monte Carlo (MC) methods are important computational tools for molecular structure optimizations and predictions. When solvent effects are explicitly considered, MC methods become very expensive due to the large degree of freedom associated with the water molecules and mobile ions. Alternatively implicit-solvent MC can largely reduce the computational cost by applying a mean field approximation to solvent effects and meanwhile maintains the atomic detail of the target molecule. The two most popular implicit-solvent models are the Poisson-Boltzmann (PB) model and the Generalized Born (GB) model in a way such that the GB model is an approximation to the PB model but is much faster in simulation time. In this work, we develop a machine learning-based implicit-solvent Monte Carlo (MLIMC) method by combining the advantages of both implicit solvent models in accuracy and efficiency. Specifically, the MLIMC method uses a fast and accurate PB-based machine learning (PBML) scheme to compute the electrostatic solvation free energy at each step. We validate our MLIMC method by using a benzene-water system and a protein-water system. We show that the proposed MLIMC method has great advantages in speed and accuracy for molecular structure optimization and prediction.