Multiscale method for modeling binding phenomena involving large objects: application to kinesin motor domains motion along microtubules

Electrostatics in Computational Biophysics and Its Implications for Disease Effects

This review outlines the role of electrostatics in computational molecular biophysics and its implication in altering wild-type characteristics of biological macromolecules, and thus the contribution of electrostatics to disease mechanisms. The work is not intended to review existing computational approaches or to propose further developments. Instead, it summarizes the outcomes of relevant studies and provides a generalized classification of major mechanisms that involve electrostatic effects in both wild-type and mutant biological macromolecules. It emphasizes the complex role of electrostatics in molecular biophysics, such that the long range of electrostatic interactions causes them to dominate all other forces at distances larger than several Angstroms, while at the same time, the alteration of short-range wild-type electrostatic pairwise interactions can have pronounced effects as well. Because of this dual nature of electrostatic interactions, being dominant at long-range and being very specific at short-range, their implications for wild-type structure and function are quite pronounced. Therefore, any disruption of the complex electrostatic network of interactions may abolish wild-type functionality and could be the dominant factor contributing to pathogenicity. However, we also outline that due to the plasticity of biological macromolecules, the effect of amino acid mutation may be reduced, and thus a charge deletion or insertion may not necessarily be deleterious.

Long-range electrostatic interactions significantly modulate the affinity of dynein for microtubules

The dynein family of microtubule minus-end-directed motor proteins drives diverse functions in eukaryotic cells, including cell division, intracellular transport, and flagellar beating. Motor protein processivity, which characterizes how far a motor walks before detaching from its filament, depends on the interaction between its microtubule-binding domain (MTBD) and the microtubule. Dynein's MTBD switches between high- and low-binding affinity states as it steps. Significant structural and functional data show that specific salt bridges within the MTBD and between the MTBD and the microtubule govern these affinity state shifts. However, recent computational work suggests that nonspecific, long-range electrostatic interactions between the MTBD and the microtubule may also play an important role in the processivity of dynein. To investigate this hypothesis, we mutated negatively charged amino acids remote from the dynein MTBD-microtubule-binding interface to neutral residues and measured the binding affinity using microscale thermophoresis and optical tweezers. We found a significant increase in the binding affinity of the mutated MTBDs for microtubules. Furthermore, we found that charge screening by free ions in solution differentially affected the binding and unbinding rates of MTBDs to microtubules. Together, these results demonstrate a significant role for long-range electrostatic interactions in regulating dynein-microtubule affinity. Moreover, these results provide insight into the principles that potentially underlie the biophysical differences between molecular motors with various processivities and protein-protein interactions more generally.

Computational Study on E-Hooks of Tubulins in the Binding Process with Kinesin

Cargo transport within cells is essential to healthy cells, which requires microtubules-based motors, including kinesin. The C-terminal tails (E-hooks) of alpha and beta tubulins of microtubules have been proven to play important roles in interactions between the kinesins and tubulins. Here, we implemented multi-scale computational methods in E-hook-related analyses, including flexibility investigations of E-hooks, binding force calculations at binding interfaces between kinesin and tubulins, electrostatic potential calculations on the surface of kinesin and tubulins. Our results show that E-hooks have several functions during the binding process: E-hooks utilize their own high flexibilities to increase the chances of reaching a kinesin; E-hooks help tubulins to be more attractive to kinesin. Besides, we also observed the differences between alpha and beta tubulins: beta tubulin shows a higher flexibility than alpha tubulin; beta tubulin generates stronger attractive forces (about twice the strengths) to kinesin at different distances, no matter with E-hooks in the structure or not. Those facts may indicate that compared to alpha tubulin, beta tubulin contributes more to attracting and catching a kinesin to microtubule. Overall, this work sheds the light on microtubule studies, which will also benefit the treatments of neurodegenerative diseases, cancer treatments, and preventions in the future.

Computational Study on DNA Repair: The Roles of Electrostatic Interactions Between Uracil-DNA Glycosylase (UDG) and DNA

Uracil-DNA glycosylase (UDG) is one of the most important base excision repair (BER) enzymes involved in the repair of uracil-induced DNA lesion by removing uracil from the damaged DNA. Uracil in DNA may occur due to cytosine deamination or deoxy uridine monophosphate (dUMP) residue misincorporation during DNA synthesis. Medical evidences show that an abnormal expression of UDG is related to different types of cancer, including colorectal cancer, lung cancer, and liver cancer. Therefore, the research of UDG is crucial in cancer treatment and prevention as well as other clinical activities. Here we applied multiple computational methods to study UDG in several perspectives: Understanding the stability of the UDG enzyme in different pH conditions; studying the differences in charge distribution between the pocket side and non-pocket side of UDG; analyzing the field line distribution at the interfacial area between UDG and DNA; and performing electrostatic binding force analyses of the special region of UDG (pocket area) and the target DNA base (uracil) as well as investigating the charged residues on the UDG binding pocket and binding interface. Our results show that the whole UDG binding interface, and not the UDG binding pocket area alone, provides the binding attractive force to the damaged DNA at the uracil base.

Electrostatic features for nucleocapsid proteins of SARS-CoV and SARS-CoV-2

COVID-19 is increasingly affecting human health and global economy. Understanding the fundamental mechanisms of Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) is highly demanded to develop treatments for COVID-19. SARS-CoV and SARS-CoV-2 share 92.06% identity in their N protein RBDs' sequences, which results in very similar structures. However, the SARS-CoV-2 is more easily to spread. Utilizing multi-scale computational approaches, this work studied the fundamental mechanisms of the nucleocapsid (N) proteins of SARS-CoV and SARS-CoV-2, including their stabilities and binding strengths with RNAs at different pH values. Electrostatic potential on the surfaces of N proteins show that both the N proteins of SARS-CoV and SARS-CoV-2 have dominantly positive potential to attract RNAs. The binding forces between SARS-CoV N protein and RNAs at different distances are similar to that of SARS-CoV-2, both in directions and magnitudes. The electric filed lines between N proteins and RNAs are also similar for both SARS-CoV and SARS-CoV-2. The folding energy and binding energy dependence on pH revealed that the best environment for N proteins to perform their functions with RNAs is the weak acidic environment.

Hybrid method for representing ions in implicit solvation calculations

Fast and accurate calculations of the electrostatic features of highly charged biomolecules such as DNA, RNA, and highly charged proteins are crucial and challenging tasks. Traditional implicit solvent methods calculate the electrostatic features quickly, but these methods are not able to balance the high net biomolecular charges effectively. Explicit solvent methods add unbalanced ions to neutralize the highly charged biomolecules in molecular dynamic simulations, which require more expensive computing resources. Here we report developing a novel method, Hybridizing Ions Treatment (HIT), which hybridizes the implicit solvent method with an explicit method to realistically calculate the electrostatic potential for highly charged biomolecules. HIT utilizes the ionic distribution from an explicit method to predict the bound ions. The bound ions are then added in the implicit solvent method to perform the electrostatic potential calculations. In this study, two training sets were developed to optimize parameters for HIT. The performance on the testing set demonstrates that HIT significantly improves the electrostatic calculations. Results on molecular motors myosin and kinesin reveal some mechanisms and explain some previous experimental findings. HIT can be widely used to study highly charged biomolecules, including DNA, RNA, molecular motors, and other highly charged biomolecules. The HIT package is available at http://compbio.utep.edu/static/downloads/download_hit.zip.

StructureMan: A Structure Manipulation Tool to Study Large Scale Biomolecular Interactions

Studying biomolecular interactions is a crucial but challenging task. Due to their large scales, many biomolecular interactions are difficult to be simulated via all atom models. An effective approach to investigate the biomolecular interactions is highly demanded in many areas. Here we introduce a Structure Manipulation (StructureMan) program to operate the structures when studying the large-scale biomolecular interactions. This novel StructureMan tool provides comprehensive operations which can be utilized to study the interactions in various large biological systems. Combining with electrostatic calculation programs such as DelPhi and DelPhiForce, StructureMan was implemented to reveal the detailed electrostatic features in two large biological examples, the viral capsid and molecular motor-microtubule complexes. Applications on these two examples revealed interesting binding mechanisms in the viral capsid and molecular motor. Such applications demonstrated that the StructureMan can be widely used when studying the biomolecular interactions in large scale biological problems. This novel tool provides an alternative approach to efficiently study the biomolecular interactions, especially for large scale biology systems. The StructureMan tool is available at our website: http://compbio.utep.edu/static/downloads/script-for-munipulation2.zip.

Spike Proteins of SARS-CoV and SARS-CoV-2 Utilize Different Mechanisms to Bind With Human ACE2

The ongoing outbreak of COVID-19 has been a serious threat to human health worldwide. The virus SARS-CoV-2 initiates its infection to the human body via the interaction of its spike (S) protein with the human Angiotensin-Converting Enzyme 2 (ACE2) of the host cells. Therefore, understanding the fundamental mechanisms of how SARS-CoV-2 S protein receptor binding domain (RBD) binds to ACE2 is highly demanded for developing treatments for COVID-19. Here we implemented multi-scale computational approaches to study the binding mechanisms of human ACE2 and S proteins of both SARS-CoV and SARS-CoV-2. Electrostatic features, including electrostatic potential, electric field lines, and electrostatic forces of SARS-CoV and SARS-CoV-2 were calculated and compared in detail. The results demonstrate that SARS-CoV and SARS-CoV-2 S proteins are both attractive to ACE2 by electrostatic forces even at different distances. However, the residues contributing to the electrostatic features are quite different due to the mutations between SARS-CoV S protein and SARS-CoV-2 S protein. Such differences are analyzed comprehensively. Compared to SARS-CoV, the SARS-CoV-2 binds with ACE2 using a more robust strategy: The electric field line related residues are distributed quite differently, which results in a more robust binding strategy of SARS-CoV-2. Also, SARS-CoV-2 has a higher electric field line density than that of SARS-CoV, which indicates stronger interaction between SARS-CoV-2 and ACE2, compared to that of SARS-CoV. Key residues involved in salt bridges and hydrogen bonds are identified in this study, which may help the future drug design against COVID-19.

Revealing the mechanism of SARS-CoV-2 spike protein binding with ACE2

A large population in the world has been infected by COVID-19. Understanding the mechanisms of Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) is important for the management and treatment of COVID-19. When it comes to the infection process, one of the most important proteins in SARS-CoV-2 is the spike (S) protein, which is able to bind to human Angiotensin-Converting Enzyme 2 (ACE2) and initializes the entry of the host cell. In this study, we implemented multiscale computational approaches to study the electrostatic features of the interfaces of the SARS-CoV-2 S protein receptor binding domain and ACE2. The simulations and analyses were performed on high-performance computing resources in the Texas Advanced Computing Center. Our study identified key residues on SARS-CoV-2, which can be used as targets for future drug design. The results shed light on future drug design and therapeutic targets for COVID-19.

Ab-initio binding of barnase–barstar with DelPhiForce steered Molecular Dynamics (DFMD) approach

Receptor–ligand interactions are involved in various biological processes, therefore understanding the binding mechanism and ability to predict the binding mode are essential for many biological investigations. While many computational methods exist to predict the 3D structure of the corresponding complex provided the knowledge of the monomers, here we use the newly developed DelPhiForce steered Molecular Dynamics (DFMD) approach to model the binding of barstar to barnase to demonstrate that first-principles methods are also capable of modeling the binding. Essential component of DFMD approach is enhancing the role of long-range electrostatic interactions to provide guiding force of the monomers toward their correct binding orientation and position. Thus, it is demonstrated that the DFMD can successfully dock barstar to barnase even if the initial positions and orientations of both are completely different from the correct ones. Thus, the electrostatics provides orientational guidance along with pulling force to deliver the ligand in close proximity to the receptor.

A computational model of ESAT-6 complex in membrane

One quarter of the world's population are infected by Mycobacterium tuberculosis (Mtb), which is a leading death-causing bacterial pathogen. Recent evidence has demonstrated that two virulence factors, ESAT-6 and CFP-10, play crucial roles in Mtb's cytosolic translocation. Many efforts have been made to study the ESAT-6 and CFP-10 proteins. Some studies have shown that ESAT-6 has an essential role in rupturing phagosome. However, the mechanisms of how ESAT-6 interacts with the membrane have not yet been fully understood. Recent studies indicate that the ESAT-6 disassociates with CFP-10 upon their interaction with phagosome membrane, forming a membrane-spanning pore. Based on these observations, as well as the available structure of ESAT-6, ESAT-6 is hypothesized to form an oligomer for membrane insertion as well as rupturing. Such an ESAT-6 oligomer may play a significant role in the tuberculosis infection. Therefore, deeper understanding of the oligomerization of ESAT-6 will establish new directions for tuberculosis treatment. However, the structure of the oligomer of ESAT-6 is not known. Here, we proposed a comprehensive approach to model the complex structures of ESAT-6 oligomer inside a membrane. Several computational tools, including MD simulation, symmetrical docking, MM/PBSA, are used to obtain and characterize such a complex structure. Results from our studies lead to a well-supported hypothesis of the ESAT-6 oligomerization as well as the identification of essential residues in stabilizing the ESAT-6 oligomer which provide useful insights for future drug design targeting tuberculosis. The approach in this research can also be used to model and study other cross-membrane complex structures.

DelPhi Suite: New Developments and Review of Functionalities

Electrostatic potential, energies, and forces affect virtually any process in molecular biology, however, computing these quantities is a difficult task due to irregularly shaped macromolecules and the presence of water. Here, we report a new edition of the popular software package DelPhi along with describing its functionalities. The new DelPhi is a C++ object-oriented package supporting various levels of multiprocessing and memory distribution. It is demonstrated that multiprocessing results in significant improvement of computational time. Furthermore, for computations requiring large grid size (large macromolecular assemblages), the approach of memory distribution is shown to reduce the requirement of RAM and thus permitting large-scale modeling to be done on Linux clusters with moderate architecture. The new release comes with new features, whose functionalities and applications are described as well. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Modeling Electrostatic Force in Protein-Protein Recognition

Electrostatic interactions are important for understanding molecular interactions, since they are long-range interactions and can guide binding partners to their correct binding positions. To investigate the role of electrostatic forces in molecular recognition, we calculated electrostatic forces between binding partners separated at various distances. The investigation was done on a large set of 275 protein complexes using recently developed DelPhiForce tool and in parallel, evaluating the total electrostatic force via electrostatic association energy. To accomplish the goal, we developed a method to find an appropriate direction to move one chain of protein complex away from its bound position and then calculate the corresponding electrostatic force as a function of separation distance. It is demonstrated that at large distances between the partners, the electrostatic force (magnitude and direction) is consistent among the protocols used and the main factors contributing to it are the net charge of the partners and their interfaces. However, at short distances, where partners form specific pair-wise interactions or de-solvation penalty becomes significant, the outcome depends on the precise balance of these factors. Based on the electrostatic force profile (force as a function of distance), we group the cases into four distinctive categories, among which the most intriguing is the case termed "soft landing." In this case, the electrostatic force at large distances is favorable assisting the partners to come together, while at short distance it opposes binding, and thus slows down the approach of the partners toward their physical binding.

Microtubules: From understanding their dynamics to using them as potential therapeutic targets

Microtubules (MT) and actin microfilaments are dynamic cytoskeleton components involved in a range of intracellular processes. MTs play a role in cell division, beating of cilia and flagella, and intracellular transport. Over the past decades, much knowledge has been gained regarding MT function and structure, and its role in underlying disease progression. This makes MT potential therapeutic targets for various disorders. Disturbances in MT and their associated proteins are the underlying cause of diseases such as Alzheimer's disease, cancer, and several genetic diseases. Some of the advances in the field of MT research, as well as the potenti G beta gamma, is needed al uses of MT-targeting agents in various conditions have been reviewed here.

Prediction of pKa Values for Neutral and Basic Drugs based on Hybrid Artificial Intelligence Methods

The pKa value of drugs is an important parameter in drug design and pharmacology. In this paper, an improved particle swarm optimization (PSO) algorithm was proposed based on the population entropy diversity. In the improved algorithm, when the population entropy was higher than the set maximum threshold, the convergence strategy was adopted; when the population entropy was lower than the set minimum threshold the divergence strategy was adopted; when the population entropy was between the maximum and minimum threshold, the self-adaptive adjustment strategy was maintained. The improved PSO algorithm was applied in the training of radial basis function artificial neural network (RBF ANN) model and the selection of molecular descriptors. A quantitative structure-activity relationship model based on RBF ANN trained by the improved PSO algorithm was proposed to predict the pKa values of 74 kinds of neutral and basic drugs and then validated by another database containing 20 molecules. The validation results showed that the model had a good prediction performance. The absolute average relative error, root mean square error, and squared correlation coefficient were 0.3105, 0.0411, and 0.9685, respectively. The model can be used as a reference for exploring other quantitative structure-activity relationships.

Divide-and-conquer strategy for large-scale Eulerian solvent excluded surface

MOTIVATION

Surface generation and visualization are some of the most important tasks in biomolecular modeling and computation. Eulerian solvent excluded surface (ESES) software provides analytical solvent excluded surface (SES) in the Cartesian grid, which is necessary for simulating many biomolecular electrostatic and ion channel models. However, large biomolecules and/or fine grid resolutions give rise to excessively large memory requirements in ESES construction. We introduce an out-of-core and parallel algorithm to improve the ESES software.

RESULTS

The present approach drastically improves the spatial and temporal efficiency of ESES. The memory footprint and time complexity are analyzed and empirically verified through extensive tests with a large collection of biomolecule examples. Our results show that our algorithm can successfully reduce memory footprint through a straightforward divide-and-conquer strategy to perform the calculation of arbitrarily large proteins on a typical commodity personal computer. On multi-core computers or clusters, our algorithm can reduce the execution time by parallelizing most of the calculation as disjoint subproblems. Various comparisons with the state-of-the-art Cartesian grid based SES calculation were done to validate the present method and show the improved efficiency. This approach makes ESES a robust software for the construction of analytical solvent excluded surfaces.

AVAILABILITY AND IMPLEMENTATION

http://weilab.math.msu.edu/ESES.

Forces and Disease: Electrostatic force differences caused by mutations in kinesin motor domains can distinguish between disease-causing and non-disease-causing mutations

The ability to predict if a given mutation is disease-causing or not has enormous potential to impact human health. Typically, these predictions are made by assessing the effects of mutation on macromolecular stability and amino acid conservation. Here we report a novel feature: the electrostatic component of the force acting between a kinesin motor domain and tubulin. We demonstrate that changes in the electrostatic component of the binding force are able to discriminate between disease-causing and non-disease-causing mutations found in human kinesin motor domains using the receiver operating characteristic (ROC). Because diseases may originate from multiple effects not related to kinesin-microtubule binding, the prediction rate of 0.843 area under the ROC plot due to the change in magnitude of the electrostatic force alone is remarkable. These results reflect the dependence of kinesin’s function on motility along the microtubule, which suggests a precise balance of microtubule binding forces is required.

Accurate, robust, and reliable calculations of Poisson-Boltzmann binding energies

Poisson-Boltzmann (PB) model is one of the most popular implicit solvent models in biophysical modeling and computation. The ability of providing accurate and reliable PB estimation of electrostatic solvation free energy, ΔGel, and binding free energy, ΔΔGel, is important to computational biophysics and biochemistry. In this work, we investigate the grid dependence of our PB solver (MIBPB) with solvent excluded surfaces for estimating both electrostatic solvation free energies and electrostatic binding free energies. It is found that the relative absolute error of ΔGel obtained at the grid spacing of 1.0 Å compared to ΔGel at 0.2 Å averaged over 153 molecules is less than 0.2%. Our results indicate that the use of grid spacing 0.6 Å ensures accuracy and reliability in ΔΔGel calculation. In fact, the grid spacing of 1.1 Å appears to deliver adequate accuracy for high throughput screening. © 2017 Wiley Periodicals, Inc.

A New DelPhi Feature for Modeling Electrostatic Potential around Proteins: Role of Bound Ions and Implications for Zeta-Potential

A new feature of the popular software DelPhi is developed and reported, allowing for computing the surface averaged electrostatic potential (SAEP) of macromolecules. The user is given the option to specify the distance from the van der Waals surface where the electrostatic potential will be outputted. In conjunction with DelPhiPKa and the BION server, the user can adjust the charges of titratable groups according to specific pH values, and add explicit ions bound to the macromolecular surface. This approach is applied to a set of four proteins with "experimentally" delivered zeta (ζ)-potentials at different pH values and salt concentrations. It has been demonstrated that the protocol is capable of predicting ζ-potentials in the case of proteins with relatively large net charges. This protocol has been less successful for proteins with low net charges. The work demonstrates that in the case of proteins with large net charges, the electrostatic potential should be collected at distances about 4 Å away from the vdW surface and explicit ions should be added at a binding energy cutoff larger than 1-2kT, in order to accurately predict ζ-potentials. The low salt conditions substantiate this effect of ions on SAEP.

Cytoplasmic dynein binding, run length, and velocity are guided by long-range electrostatic interactions

Dyneins are important molecular motors involved in many essential biological processes, including cargo transport along microtubules, mitosis, and in cilia. Dynein motility involves the coupling of microtubule binding and unbinding to a change in the configuration of the linker domain induced by ATP hydrolysis, which occur some 25 nm apart. This leaves the accuracy of dynein stepping relatively inaccurate and susceptible to thermal noise. Using multi-scale modeling with a computational focusing technique, we demonstrate that the microtubule forms an electrostatic funnel that guides the dynein's microtubule binding domain (MTBD) as it finally docks to the precise, keyed binding location on the microtubule. Furthermore, we demonstrate that electrostatic component of the MTBD's binding free energy is linearly correlated with the velocity and run length of dynein, and we use this linearity to predict the effect of mutating each glutamic and aspartic acid located in MTBD domain to alanine. Lastly, we show that the binding of dynein to the microtubule is associated with conformational changes involving several helices, and we localize flexible hinge points within the stalk helices. Taken all together, we demonstrate that long range electrostatic interactions bring a level of precision to an otherwise noisy dynein stepping process.