Physical Sciences and Mathematics Commons^™

On The Modeling Of Polar Component Of Solvation Energy Using Smooth Gaussian-Based Dielectric Function, Lin Li, Chuan Li, Emil Alexov

Publications

Traditional implicit methods for modeling electrostatics in biomolecules use a two-dielectric approach: a biomolecule is assigned low dielectric constant while the water phase is considered as a high dielectric constant medium. However, such an approach treats the biomolecule-water interface as a sharp dielectric border between two homogeneous dielectric media and does not account for inhomogeneous dielectric properties of the macromolecule as well. Recently we reported a new development, a smooth Gaussian-based dielectric function which treats the entire system, the solute and the water phase, as inhomogeneous dielectric medium (J Chem Theory Comput. 2013 Apr 9; 9(4): 2126-2136.). Here we examine …

On The Electrostatic Properties Of Homodimeric Proteins, Brandon Campbell, Marharyta Petukh, Emil Alexov, Chuan Li

Publications

A large fraction of proteins function as homodimers, but it is not always clear why the dimerization is important for functionality since frequently each monomer possesses a distinctive active site. Recent work (PLoS Computational Biology, 9(2), e1002924) indicates that homodimerization may be important for forming an electrostatic funnel in the spermine synthase homodimer which guides changed substrates toward the active centers. This prompted us to investigate the electrostatic properties of a large set of homodimeric proteins and resulted in an observation that in a vast majority of the cases the dimerization indeed results in specific electrostatic features, although not necessarily …

Continuous Development Of Schemes For Parallel Computing Of The Electrostatics In Biological Systems: Implementation In Delphi, Chuan Li, Marharyta Petukh, Lin Li, Emil Alexov

Publications

Due to the enormous importance of electrostatics in molecular biology, calculating the electrostatic potential and corresponding energies has become a standard computational approach for the study of biomolecules and nano-objects immersed in water and salt phase or other media. However, the electrostatics of large macromolecules and macromolecular complexes, including nano-objects, may not be obtainable via explicit methods and even the standard continuum electrostatics methods may not be applicable due to high computational time and memory requirements. Here, we report further development of the parallelization scheme reported in our previous work (J Comput Chem. 2012 Sep 15; 33(24):1960–6.) to include parallelization …

The Role Of Protonation States In Ligand-Receptor Recognition And Binding, Marharyta Petukh, Shannon Stefl, Emil Alexov

Publications

In this review we discuss the role of protonation states in receptor-ligand interactions, providing experimental evidences and computational predictions that complex formation may involve titratable groups with unusual pKa’s and that protonation states frequently change from unbound to bound states. These protonation changes result in proton uptake/release, which in turn causes the pHdependence of the binding. Indeed, experimental data strongly suggest that almost any binding is pH-dependent and to be correctly modeled, the protonation states must be properly assigned prior to and after the binding. One may accurately predict the protonation states when provided with the structures of the unbound …

Delphi Web Server: A Comprehensive Online Suite For Electrostatic Calculations Of Biological Macromolecules And Their Complexes., Subhra Sarkar, Shawn Witham, Jie Zhang, Maxim Zhenirovskyy, Walter Rocchia

Publications

Here we report a web server, the DelPhi web server, which utilizes DelPhi program to calculate electrostatic energies and the corresponding electrostatic potential and ionic distributions, and dielectric map. The server provides extra services to fix structural defects, as missing atoms in the structural file and allows for generation of missing hydrogen atoms. The hydrogen placement and the corresponding DelPhi calculations can be done with user selected force field parameters being either Charmm22, Amber98 or OPLS. Upon completion of the calculations, the user is given option to download fixed and protonated structural file, together with the parameter and Delphi output …

In Silico Investigation Of Ph-Dependence Of Prolactin And Human Growth Hormone Binding To Human Prolactin Receptor, Lin Wang, Shawn Witham, Zhe Zhang, Michael E. Hodsdon, Emil Alexov

Publications

Experimental data shows that the binding of human prolactin (hPRL) to human prolactin receptor (hPRLr-ECD) is strongly pH-dependent, while the binding of the same receptor to human growth hormone (hGH) is pH-independent. Here we carry in silico analysis of the molecular effects causing such a difference and reveal the role of individual amino acids. It is shown that the computational modeling correctly predicts experimentally determined pKa’s of histidine residues in an unbound state in the majority of the cases and the pH-dependence of the binding free energy. Structural analysis carried in conjunction with calculated pH-dependence of the binding revealed that …

Highly Efficient And Exact Method For Parallelization Of Grid-Based Algorithms And Its Implementation In Delphi, Chuan Li, Lin Li, Jie Zhang, Emil Alexov

Publications

The Gauss–Seidel (GS) method is a standard iterative numerical method widely used to solve a system of equations and, in general, is more efficient comparing to other iterative methods, such as the Jacobi method. However, standard implementation of the GS method restricts its utilization in parallel computing due to its requirement of using updated neighboring values (i.e., in current iteration) as soon as they are available. Here, we report an efficient and exact (not requiring assumptions) method to parallelize iterations and to reduce the computational time as a linear/nearly linear function of the number of processes or computing units. In …

Delphi: A Comprehensive Suite For Delphi Software And Associated Resources, Lin Li, Chuan Li, Subhra Sarkar, Jie Zhang, Shawn Witham, Zhe Zhang, Lin Wang, Nicholas Smith, Marharyta Petukh, Emil Alexov

Publications

Background

Accurate modeling of electrostatic potential and corresponding energies becomes increasingly important for understanding properties of biological macromolecules and their complexes. However, this is not an easy task due to the irregular shape of biological entities and the presence of water and mobile ions.

Results

Here we report a comprehensive suite for the well-known Poisson-Boltzmann solver, DelPhi, enriched with additional features to facilitate DelPhi usage. The suite allows for easy download of both DelPhi executable files and source code along with a makefile for local installations. The users can obtain the DelPhi manual and parameter files required for the corresponding …

Developing Hybrid Approaches To Predict Pka Values Of Ionizable Groups, Shawn Witham, Kemper Talley, Lin Wang, Zhe Zhang, Daquan Gao, Wei Yang, Emil Alexov

Publications

Accurate predictions of pKa values of titratable groups require taking into account all relevant processes associated with the ionization/deionization. Frequently, however, the ionization does not involve significant structural changes and the dominating effects are purely electrostatic in origin allowing accurate predictions to be made based on the electrostatic energy difference between ionized and neutral forms alone using a static structure. On another hand, if the change of the charge state is accompanied by a structural reorganization of the target protein, then the relevant conformational changes have to be taken into account in the pKa calculations. Here we report a hybrid …

A Missense Mutation In Clic2 Associated With Intellectual Disability Is Predicted By In Silico Modeling To Affect Protein Stability And Dynamics, Shawn Witham, Kyoko Takano, Charles Schwartz, Emil Alexov

Publications

Large-scale next generation resequencing of X chromosome genes identified a missense mutation in the CLIC2 gene on Xq28 in a male with X-linked intellectual disability (XLID) and not found in healthy individuals. At the same time, numerous nsSNPs (nonsynonomous SNP) have been reported in the CLIC2 gene in healthy individuals indicating that the CLIC2 protein can tolerate amino acid substitutions and be fully functional. To test the possibility that p.H101Q is a disease-causing mutation, we performed in silico simulations to calculate the effects of the p.H101Q mutation on CLIC2 stability, dynamics, and ionization states while comparing the effects obtained for …

In Silico Modeling Of Ph-Optimum Of Protein-Protein Binding, Rooplekha C. Mitra, Zhe Zhang, Emil Alexov

Publications

Protein-protein association is a pH-dependent process and thus the binding affinity depends on the local pH. In vivo the association occurs in a particular cellular compartment, where the individual monomers are supposed to meet and form a complex. Since the monomers and the complex exist in the same micro environment, it is plausible that they coevolved toward its properties, in particular, toward the characteristic subcellular pH. Here we show that the pH at which the monomers are most stable (pH-optimum) or the pH at which stability is almost pH-independent (pH-flat) of monomers are correlated with the pH-optimum of maximal affinity …