Next Article in Journal
Visual Cortical Plasticity: Molecular Mechanisms as Revealed by Induction Paradigms in Rodents
Next Article in Special Issue
Impact of E484Q and L452R Mutations on Structure and Binding Behavior of SARS-CoV-2 B.1.617.1 Using Deep Learning AlphaFold2, Molecular Docking and Dynamics Simulation
Previous Article in Journal
β-Adrenoreceptors as Therapeutic Targets for Ocular Tumors and Other Eye Diseases—Historical Aspects and Nowadays Understanding
Previous Article in Special Issue
Atomistic Study for the Tantalum and Tantalum–Tungsten Alloy Threshold Displacement Energy under Local Strain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Intrinsic Radius as a Key Parameter in the Generalized Born Model to Adjust Protein-Protein Electrostatic Interaction

1
Research Institute for Science and Engineering, Waseda University, Okubo 3-4-1, Sinjuku-ku, Tokyo 169-8555, Japan
2
Department of Pure and Applied Physics, Waseda University, Okubo 3-4-1, Sinjuku-ku, Tokyo 169-8555, Japan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(5), 4700; https://doi.org/10.3390/ijms24054700
Submission received: 7 January 2023 / Revised: 30 January 2023 / Accepted: 19 February 2023 / Published: 28 February 2023
(This article belongs to the Special Issue Research on Molecular Dynamics)

Abstract

:
The generalized Born (GB) model is an extension of the continuum dielectric theory of Born solvation energy and is a powerful method for accelerating the molecular dynamic (MD) simulations of charged biological molecules in water. While the effective dielectric constant of water that varies as a function of the separation distance between solute molecules is incorporated into the GB model, adjustment of the parameters is indispensable for accurate calculation of the Coulomb (electrostatic) energy. One of the key parameters is the lower limit of the spatial integral of the energy density of the electric field around a charged atom, known as the intrinsic radius ρ . Although ad hoc adjustment of ρ has been conducted to improve the Coulombic (ionic) bond stability, the physical mechanism by which ρ affects the Coulomb energy remains unclear. Via energetic analysis of three differently sized systems, here, we clarify that the Coulomb bond stability increases with increasing ρ and that the increased stability is caused by the interaction energy term, not by the self-energy (desolvation energy) term, as was supposed previously. Our results suggest that the use of larger values for the intrinsic radii of hydrogen and oxygen atoms, together with the use of a relatively small value for the spatial integration cutoff in the GB model, can better reproduce the Coulombic attraction between protein molecules.

1. Introduction

The Coulomb interaction is important for the assembly of biomolecules and their fast association in aqueous solution [1,2]. The stability of a Coulomb bond (i.e., an ionic or hydrogen bond) in water, however, is marginal compared to that in a vacuum, which is usually understood by the large dielectric constant of water. For the Coulomb interaction between two oppositely charged molecules that have their own water-excluded volumes, the effective dielectric constant is lowered as the two molecules approach each other because of the presence of the water-excluded (and hence, low-dielectric) regions. The lowering of the effective dielectric constant causes augmentation of the attractive Coulomb interaction energy of the two molecules and, simultaneously, causes the increase in the Coulomb self-energy, which is always positive (in other words, the reduction in the Born solvation energy [3], which is always negative). Thus, the favorable change in the Coulomb interaction energy upon the Coulomb bond formation is counterbalanced by the unfavorable change in the Coulomb self-energy, which makes the stability of the Coulomb bond in water marginal.
The lowering of the effective dielectric constant and the counterbalance between the Coulomb interaction and the self-energies, as mentioned above, are taken into consideration in the generalized Born (GB) solvation model [4] that has been applied to the molecular dynamics (MD) simulation of biomolecular systems in water [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23]. The obvious benefit of using the GB model is the speed-up of the MD simulation due to replacing a large number of the solvent degrees of freedom with a continuum dielectric. On the other hand, the adjustment of the parameters employed in the GB model is indispensable to calculating the Coulomb energies accurately [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23]. Among the key parameters of the GB model, we previously studied the effect of the upper limit (cutoff distance), denoted by r max , for the spatial integral of the energy density of the electric field around a charged atom [24]. The spatial integral of the energy density of the electric field is used to obtain the effective Born radius of the atom (the most important quantity in the GB model that measures the effective dielectric constant around the atom). We found that r max affects the electrostatic interaction between charged molecules; the attractive Coulomb interaction between oppositely charged molecules becomes repulsive when a large r max value is used, which was shown to be caused by the overestimation of the unfavorable Coulomb self-energy (or desolvation penalty) upon binding so that the use of a relatively small r max value was recommended [24].
As another important parameter of the GB model, the lower limit of the above-mentioned integral, was recognized early [4]. The lower limit of the integral defines the dielectric boundary of the solute region and the solvent region, and the van der Waals radius of the atom (called the Bondi radius [25]) is basically used for this lower limit, which is referred to as the GB atomic radius or the intrinsic radius, denoted by ρ [11,12]. The parameter ρ has been used for the ad hoc adjustment of the stability of the Coulomb bond [11,12,16,19,20,23,26,27]. For example, Tsui and Case augmented the stability of hydrogen bonds between DNA strands by increasing the ρ of the hydrogen atom [11,12]. Im and coworkers reduced the stability of the backbone hydrogen bonds of helix-forming peptides by reducing the ρ of the nitrogen atom [19]. Although it is commonly observed that increasing the ρ leads to the enhancement of the Coulombic (i.e, ionic or hydrogen) bond stability, the underlying physical mechanism has been only speculated from the viewpoint of the reduced desolvation penalty.
In the present study, which forms a counterpart to the previous study [24], we systematically analyze the dependence of the Coulomb interaction energy and the Coulomb self-energy on ρ and show that the enhancement of the Coulomb bond stability with increasing ρ is caused by the interaction energy term, not by the self-energy (desolvation energy) term, as was previously supposed. We also show that ρ (the lower limit of the spatial integration) affects the Coulomb interaction at short distances, whereas r max (the upper limit) affects it at long distances, thereby enabling us to use these two parameters in combination. Indeed, we show that larger intrinsic radii, in combination with a relatively small cutoff distance for the spatial integral, can better reproduce the attractive interaction between oppositely charged protein molecules.

2. Results and Discussion

We analyzed the effect of the intrinsic radii, ρ , on the Coulomb self-energy (i.e., the desolvation energy defined by Equation (13)) and the Coulomb interaction energy (Equation (15)) upon the Coulomb bond (ionic bond) formation. We employed three systems for the energy analysis: (i) a two-atom system (H–O system), (ii) a two amino acid system (Arg–Glu system), and (iii) a two-protein system (barstar–barnase system).

2.1. Two-Atom System (H–O System)

We first show the results of the energetic analysis of the Coulomb bond formation in the two-atom (H–O) system in water, as shown in Figure 1a. The total energy (the sum of the Coulomb and van der Waals energies) is shown in Figure 1b as a function of d (the separation distance between the hydrogen and oxygen atoms). We can see that the energy minimum appears at d 1.6 Å and is lowered as the intrinsic radius of the hydrogen atom, ρ H , is increased from 1.1 Å (red) to 1.5 Å (blue), indicating that the Coulomb bond is more stabilized with increasing ρ H . This increase in the Coulomb bond stability with increasing ρ H is in accord with what has been observed in the previous studies [11,12,16,23,27]. The effect of ρ H on the total energy is limited to the short distance range ( d 3 Å).
In Figure 1c,d, we show how the Coulomb self-energy (i.e., the desolvation energy) and the Coulomb interaction energy are affected by ρ H . The energy profile of Δ G O self [Figure 1c] indicates that the self-energy of the oxygen is increased as d is decreased, indicating that the unfavorable desolvation energy, called the “desolvation penalty”, arises upon the Coulomb bond formation. The desolvation penalty is seen in the total energy [Figure 1b] as an energy barrier at d 3 Å, and the energy barrier is increased as ρ H is increased. The increase in Δ G O self with increasing ρ H is understood by inspecting the spatial integral in Equation (11), that is, the I-value of the oxygen atom increased with increasing ρ H ; see the dark blue region in Figure 1a that indicates the increased volume to be integrated in Equation (11). The increase in the I-value reflects the lowering of the effective dielectric constant around the oxygen (red). Figure 1d shows that the desolvation penalty upon the Coulomb bond formation arises in Δ G H self as well. However, the ρ H dependence of Δ G H self is opposite to that of Δ G O self ; it decreases with increasing ρ H . This opposite behavior can also be understood by inspecting the spatial integral in Equation (11), that is, the I-value of the hydrogen atom decreases with increasing ρ H ; see the dark red region in the right part of Figure 1a that indicates the reduced volume to be integrated. As displayed in Figure 1f, the increase in Δ G O self with increasing ρ H is largely compensated by the decrease in Δ G H self , so the desolvation energy in total remains almost unchanged (with a slight increase at a larger ρ H ).
On the other hand, the Coulomb interaction energy, G int , as a function of d [Figure 1e] decreases as d is decreased, stabilizing the Coulomb bond, and the amount of the decrease (stabilization) is enhanced with increasing ρ H , which is clearly seen in Figure 1g, in which G int at the total energy minimum ( d = 1.57 Å) is shown as a function of ρ H . Therefore, we can tell that it is the Coulomb interaction energy G int , not the self (desolvation) energy, as was previously supposed, that contributes to the enhancement of the Coulomb bond stability at a larger ρ H . The enhancement of the attractive interaction at a larger ρ H is due to the increased water-excluded (and hence, low-dielectric) region encompassing the two interacting charges.

2.2. Two Amino Acid System (Arg–Glu System)

The Coulomb interaction between basic amino acid residues and acidic amino acid residues is of vital importance in the binding of the protein molecules and their biological function. As one example among many, we recently found the Coulomb interaction between the rigorously conserved arginine and glutamic acid in the membrane-embedded portion of ATP synthase plays a vital role in driving the proton-motive rotary motion [28,29]. Here, we study the Coulomb interaction between arginine (Arg) and glutamic acid (Glu) [Figure 2a]. We first calculated the potential of mean force (PMF) for the Coulomb bond formation between the Arg and Glu [Figure 2b]. We obtained the PMF as a function of the distance between the guanidyl carbon (C ζ atom) of the Arg and the carboxyl carbon (C δ atom), which was calculated by R T ln [ n ( d ) / N ref ] , where R is the gas constant, T is the temperature (300 K), n ( d ) is the number of the snapshots that fall in the range of d ± 0.05 Å, and N ref is the number of the snapshots in the reference state ( d = 10 Å). In Figure 2b, we can see the minimum PMF at d 4 Å both for the GB model and the explicit water model, indicating the thermodynamic stability of the Coulomb bond. It is also clearly seen that the Coulomb bond stability in the GB model increases with increasing ρ H . We then examined the Coulomb energies in detail. In Figure 2c, we show the total Coulomb energy as a function of d, which was calculated by using the snapshot structures obtained from the MD simulations with the explicit water model (a total of 40,000 snapshot structures sampled at a 0.1 ns interval). The stability of the Coulomb bond is increased with increasing ρ H , which is essentially the same as what was observed in the two-atom (H–O) system described above. In addition, by comparing the PMF [Figure 2b] with the total Coulomb energy [Figure 2c], we can find that the stability enhancement in the PMF is almost the same amount as that in the total Coulomb energy; the amount of the stability enhancement at d 4 Å from ρ H = 1.1 Å (red) to 1.5 Å (blue) is about 4 kcal/mol for both the PMF and the total Coulomb energy. We then scrutinized the self- (desolvation) energy term (Equation (16)) and the interaction energy term (Equation (17)). As seen in Figure 2d,e, the desolvation energy increases as d is decreased, reflecting the desolvation penalty upon Coulomb bond formation, and the amount of the increase in Δ G self is enhanced with increasing ρ H . On the other hand, the interaction energy decreases as d is decreased, bringing about the stability of the Coulomb bond, and the amount of the decrease in Δ G int is enhanced with increasing ρ H . The ρ H -enhanced desolvation penalty is counterbalanced by the ρ H -enhanced Coulomb bond stability, as is clearly seen in Figure 2f,g, where the average values averaged in the Coulomb bond formation state ( 3.9 d 4.0 Å) are shown as a function of ρ H . Here, the ρ H -enhancement in G int is larger than that in Δ G self , so the former surpasses the latter, and, consequently, the Coulomb bond stability is enhanced as ρ H is increased. The physical reason for the weaker ρ H -enhancement in Δ G self is essentially the same as explained in the two-atom system (the counterbalance between the increase in Δ G O self and the decrease in Δ G H self with increasing ρ H ), and, hence, the increased Coulomb bond stability with increasing ρ H is also observed in other basic residue–acidic residue pairs, such as the Lys–Asp pair (see Supplementary Figure S1).

2.3. Two-Protein System (Bn-Bs System)

The binding between the bacterial ribonuclease barnase (Bn) and its inhibitor barstar (Bs) has been intensively studied both experimentally [30] and theoretically [1,2,31,32,33] because of their strong binding affinities and fast binding rates, which are largely due to the Coulomb interaction between Bn and Bs. In fact, a clear charge complementarity is seen at the binding interface [Figure 3a], where the surface of Bn is positively charged and the surface of Bs is negatively charged. We first calculated the potential of mean force (PMF) for the binding between Bn and Bs as a function of the distance between the center of mass of Bn and that of Bs, denoted by d [Figure 3b]. The minimum PMF is found at d = 23.7 Å, which corresponds to the bound state and should be stabilized by the Coulomb interaction. Note that the non-Coulomb interaction should also contribute to the stability of the bound state in the case of the Bn–Bs binding, in which the buried surface area upon binding is much larger than that in the case of the Arg–Glu system. Accurate evaluation of the contribution of the non-Coulomb energy to the binding stability, which is needed to obtain the total binding free energy in addition to the Coulomb energy, is another important issue [23]. Here, we focus on the contribution of the Coulomb energy (we will address the contribution of the non-Coulomb energy in the next study). In contrast to the case of the Arg–Glu system, in which the water-excluded volume for an amino acid monomer is small, and, hence, the effect of the upper limit (cutoff) of the spatial integral in Equation (11), r max , is negligible, the water-excluded volume of a protein molecule is much larger so that the effect of r max on the Coulomb energy should become much greater in the Bn–Bs system. In Figure 3c, we show the total Coulomb energy as a function of d at three different values of r max (10, 15, and 25 Å). We can see that a large energy barrier arises when r max is large, which is caused by the overestimation of the desolvation penalty in the GB model, in which the so-called Coulomb field approximation is used [24]. The overestimation of the desolvation penalty can be corrected by more elaborate, yet computationally demanding, GB models (e.g., the one developed by Onufriev and coworkers [22]); alternatively, it can also be remedied to a large extent by using a smaller r max [24] in the plain and the computationally less demanding GB model [16] that we used here. However, even with r max = 10 Å, the energy barrier remains, and the stabilization of the bound state is only a little. Therefore, keeping r max at the small value of 10 Å, we studied the effect of the lower limit of the spatial integral, i.e., the intrinsic radius ρ , on the total Coulomb energy [Figure 3d]. As was the case in the H–O system and the Arg–Glu system, increasing ρ H (the intrinsic radius of the hydrogen atoms covalently bonded to a nitrogen atom in the side chain of Arg or Lys) has the effect of stabilizing the binding between Bn and Bs. The stabilizing effect is short ranged and localized in the close contact region with d 25 Å. Furthermore, Figure 3e,f clearly show that the physical mechanism of the enhancement of the stabilization with increasing ρ H is the same as that in the H–O system and in the Arg–Glu system, i.e., the decrease in Δ G int (the enhancement of the Coulomb attraction) surpasses the increase in Δ G self (Coulomb desolvation penalty) with increasing ρ H . The same physical mechanism of the stability enhancement of the Coulomb bonds holds true for the intrinsic radius of the carboxyl oxygen atoms in the side chain of Glu or Asp, so that Bn–Bs binding was further stabilized by increasing ρ O from 1.5 Å (the Bondi radius of the oxygen atom [25]) to 1.6 Å [Figure 3d–f].

3. Methods

3.1. Intrinsic Radius in GB Model

We first explain where the intrinsic radius (also referred to as the GB atomic radius) appears in the basic framework of the GB model [4].
Suppose a system composed of two spherical particles solvated in water with the effective dielectric constant ε ; the two particles possess the charges q 1 and q 2 at their centers, the sphere radii are ρ 1 and ρ 2 , respectively, and the distance between the centers of the particles is d. When we take the infinite separation of the two particles in a vacuum as the reference state, the total Coulomb energy of the above system in water at the separation distance of d, G C ( d ) is represented by the sum of the Coulomb (Born) solvation energy G sol self at the infinite separation and the Coulomb interaction energy in water G wat int , i.e.,
G C ( d ) = G sol self ( ) + G wat int ( d ) ,
where G sol self and G wat int are given by
G sol self = 1 ε w 1 q 1 2 2 ρ 1 + 1 ε w 1 q 2 2 2 ρ 2 ,
G wat int ( d ) = q 1 q 2 ε d .
The effective dielectric constant, ε , in Equation (3) should take the dielectric constant of bulk water ( ε w ) when d is large enough, and it should be lowered as d becomes smaller because the water-excluded regions come closer. Here, if we consider solvating the system at the separation distance d (instead of solvating the system at infinite separation, as assumed in Equation (1)), then Equation (1) can be replaced by
G C ( d ) = G sol self ( d ) + G sol int ( d ) + E vac int ( d ) ,
where G sol self , G sol int , and E vac int are given by
G sol self ( d ) = 1 ε ( d ) 1 q 1 2 2 ρ 1 + 1 ε ( d ) 1 q 2 2 2 ρ 2 ,
G sol int ( d ) = 1 ε ( d ) 1 q 1 q 2 d ,
E vac int ( d ) = q 1 q 2 d .
The sum of G sol self and G sol int is called the generalized Born (GB) solvation energy [4,35], which we designate as G gb . The separation-distance dependent change in the effective dielectric constant of water, as mentioned above, is taken into consideration in the GB model proposed by Still et al. [4], where G sol self and G sol int are approximated by
G sol self ( d ) G gb self ( d ) = 1 ε w 1 q 1 2 2 R 1 + 1 ε w 1 q 2 2 2 R 2 ,
G sol int ( d ) G gb int ( d ) = 1 ε w 1 q 1 q 2 d 2 + R 1 R 2 exp d 2 / ( 4 R 1 R 2 ) 1 2 ,
where R i ( i = 1 , 2 ) is called the effective Born radius and represents the d-dependent behavior of the effective dielectric constant ε around the particle i in such a manner that
1 R i = 1 ρ i I ,
I = 1 4 π ρ i | r i | r max , r i V excl 1 | r i | 4 d r i .
The right-hand side of Equation (10) originates from the spatial integral of the energy density of the electric field due to the particle i [10]. The variable of the integration in Equation (11), r i , represents the position vector, with the center of the particle i set at the origin; r max is the cutoff distance (the upper limit) of the integral; and V excl denotes the water-excluded volume due to the other particle over which the spatial integral runs (in fact, the integration is replaced by the summation of analytical pairwise functions [5,6]). The intrinsic radius, ρ i , on which we focus in this study, appears here as the lower limit for the spatial integral of Equation (11). The quantity I acts as the measure of the surrounding dielectric environment of the particle i, changing from 0 for the completely solvated state to 1 / ρ i for the completely desolvated state. Accordingly, R i increases from ρ i to infinity as the desolvation proceeds (i.e., as the effective dielectric constant ε around the particle i decreases from ε w to one). Using Equations (8) and (9), G gb can be written in a unified form as
G gb ( d ) = G gb self ( d ) + G gb int ( d ) = 1 2 1 ε w 1 i , j 2 q i q j d 2 + R i R j exp d 2 / ( 4 R i R j ) 1 2 .

3.2. Energy Analysis

We employed the following three differently sized systems in the energy analysis: (i) a two-atom system (H–O system), (ii) a two amino acid system (Arg–Glu system), and (iii) a two protein system (barstar–barnase system).
(i) Two-atom system (H–O system): The oppositely charged two-atom system [see Figure 1a], to which Equations (7)–(12) can be directly applied, was studied as the simplest system. Here, by using the charge and the intrinsic radius parameters in AMBER ff99SB [36,37] for the O ϵ atom of glutamic acid [ q O = 0.82 e (e is the elementary charge) and ρ O = 1.5 Å] and those for the H η atom of arginine ( q H = + 0.45 e and ρ H = 1.2 Å), we calculated the energy components in Equations (7)–(11) to obtain the Coulomb self-energy and interaction energy. To study the energy change accompanying the Coulomb bond (ionic bond) formation, we first calculated the Coulomb self-energy at a separation distance of d relative to that at infinite separation in water, i.e., the “desolvation energy”, for each atom (H or O), defined by
Δ G X self ( d ) = G gb , X self ( d ) G gb , X self ( ) ( X = O , H ) ,
G gb , X self ( d ) = 1 ε w 1 q X 2 2 R X ,
where R X depends on d, and R X = ρ X at d = . We also studied the change in the interaction energy accompanying the Coulomb bond formation (i.e., the sum of the second and third terms of Equation (4)),
G int ( d ) = E vac int ( d ) + G sol int ( d ) E vac int ( d ) + G gb int ( d ) = q H q O d + 1 ε w 1 q H q O d 2 + R H R O exp d 2 / ( 4 R H R O ) 1 2 .
Then, we calculated Δ G X self ( d ) and G int ( d ) by systematically changing ρ H . The r max in Equation (11) was set to 999 Å (or effectively no cutoff for the spatial integral). The bulk dielectric constant of water, ε w , was set to 78.5 throughout this study.
(ii) Two amino acid system (Arg–Glu system): As an example of the Coulomb bond (salt bridge) formation between amino acid residues, we studied the interaction between the side chain of arginine (Arg) and that of glutamic acid (Glu) [see Figure 2a]. The main-chain amino terminal and the carboxyl terminal of the Arg and Glu were neutralized by acetylation and methylamidation, respectively, and the Arg and Glu were aligned so that the main-chain orientation (the vector connecting the methyl carbons in the terminal acetyl and methylamide groups) of the Arg becomes parallel to that of the Glu, and the side-chain orientation of the Arg (the vector connecting the C α and C β atoms) faces that of the Glu. The distance between the C α atoms of the Arg and Glu was set to 14 Å. After energy minimization and subsequent equilibration at 300 K, in which the main-chain heavy atoms were positionally restrained with a spring constant of 1000 kcal/mol/Å 2 to relax the side-chain conformation, four independent 1- μ s MD runs were conducted at 300 K, with the main-chain heavy atoms positionally restrained by a spring constant of 20 kcal/mol/Å 2 . The temperature was controlled using the Langevin thermostat with a collision frequency of 1.0 ps 1 . Throughout the present study, we employed the GB model refined by Onufriev et al. [16], using the GB OBC (II) parameter set. In this study, the intrinsic radius of the hydrogen atom covalently bonded to a nitrogen atom, denoted by ρ H , was systematically changed in the range from 1.1 to 1.5 Å, in which the default value (the mBondi2 value) of 1.3 Å [16] is included.
For the purpose of comparison, we used the explicit water model for the MD simulation, in addition to the GB solvation model. We immersed the Arg and Glu in a water box containing 2337 water molecules (the TIP3P water model [38] was used), to which the periodic boundary condition and the particle mesh Ewald method [39] (with a direct space cutoff of 8 Å) were applied. After energy minimization and subsequent equilibration at 300 K and 0.1 MPa (using the Berendsen barostat [40]), in which the main-chain heavy atoms were positionally restrained with a spring constant of 1000 kcal/mol/Å 2 to relax the side-chain conformation and the system volume, four independent 1- μ s MD runs were conducted at 300 K under a constant volume condition, with the main-chain heavy atoms positionally restrained by a spring constant of 20 kcal/mol/Å 2 . We employed AMBER [41] (version 12 or later) with the ff99SB force field [36,37] for both MD simulations with the GB solvation model and the explicit solvent model. For each solvation model, we obtained the potential of mean force (PMF) as a function of the distance between the C ζ atom (guanidyl carbon) of the Arg and the C δ atom (carboxyl carbon) of the Glu to measure the stability of the Coulomb bond between the Arg and Glu.
For each of the snapshot structures obtained from the MD runs conducted in explicit water (a total of 40,000 snapshot structures sampled at a 0.1 ns interval), we calculated the total self- (or desolvation) energy and the total interaction energy, defined by
Δ G self ( d ) = i N [ G gb , i self ( d ) G gb , i self ( ) ] ,
G int ( d ) = k < l E vac , k l int ( d ) + G gb , k l int ( d ) ,
where G gb , i self conforms to Equation (14), N is the total number of the atoms in the Arg and Glu ( N = 63 ), and E vac , k l int and G gb , k l int conform to Equation (15) (k and l are the atom identifiers that belong to the Arg and Glu, respectively). Then, we calculated Δ G self ( d ) and G int ( d ) by systematically changing ρ H . The r max was set to 999 Å (or effectively no cutoff for the spatial integral).
(iii) Two-protein system (barstar–barnase system): As an example of the Coulomb bond (salt bridge) formation between protein molecules, we studied the binding between barnase (Bn) and barstar (Bs) [see Figure 3a]. We solvated the Bn–Bs complex (PDB ID: 1BRS) [30] in a water box containing 15,611 water molecules (the TIP3P [38] model was used) and four sodium ions [42] to neutralize the system, and we employed the AMBER ff03 force field [43]. We set all of the basic and acidic residues in the charged state and all histidines in the single-protonation state ( ϵ -protonated state). The missing residues (N-terminal two residues in Bn) were complemented by Modeller [44]. We then conducted the umbrella sampling in explicit water to obtain the potential of mean force of the Bn–Bs binding (as a function of the separation distance between the center of mass of the Bn and that of the Bs, denoted by d). We conducted the MD simulation using AMBER [41] (version 12 or later) at a constant temperature (set to 300 K using the Langevin thermostat with a collision frequency of 1.0 ps 1 ) and at a constant volume (the system volume was relaxed in the preparation run conducted at 0.1 MPa using the Berendsen barostat [40]).
For the umbrella sampling, we applied the umbrella potential, 1 2 a i ( d d i ) 2 ( i = 1 , 2 , . . . , 19 ) for the i-th window and changed d i from 22 to 32 Å at an interval of 0.5 Å for d i 28 with a i = 9.6 kcal/mol/Å 2 (two windows were added at d i = 25.75 and 26.75 Å to improve the overlap between the distributions of the neighboring windows). For d i > 28 Å, we changed d i at an interval of 1 Å with a i = 2.4 kcal/mol/Å 2 . During the umbrella sampling, to maintain the relative orientation between the Bn and Bs and the side-chain conformations of the residue involved in the binding interface, we employed the same restraints as those in the study by Gumbart et al [33]. For each of the 19 windows, we conducted five independent 5 ns MD runs and used the weighted histogram analysis method (WHAM) [45] to obtain the PMF as a function of d.
Using the snapshot structures obtained from the umbrella sampling [snapshots sampled at a 0.1 ns interval from the last 1 ns of the five 5 ns MD runs (i.e., total 50 snapshots) in each window, we calculated the total desolvation energy, Δ G self ( d ) , and the total interaction energy, G int ( d ) , in the same way as in Equations (16) and (17). Then, we studied the effect of ρ H (the intrinsic radius of the hydrogen atom covalently bonded to a nitrogen atom in the side chain of Arg or Lys) on Δ G self and G int . For the Bn–Bs system, we also studied the effect of the upper limit of the spatial integral, r max , on Δ G self and G int by using different values for r max (10, 15, and 25 Å).

4. Concluding Remarks

We clarified the physical mechanism of how the Coulomb bond stability is increased as the intrinsic radius, ρ , is increased in the generalized Born (GB) model. As ρ is increased, the desolvation penalty upon Coulomb bond formation is increased. However, due both to the counterbalance between the desolvation energies of the two interacting atoms and to the enhancement of the Coulomb interaction as a result of the enlarged region with a low effective dielectric constant, the desolvation penalty is more than compensated for by the enhanced stabilization of the Coulomb interaction energy. We note that increasing the intrinsic radius, the effect of which is the alleviation of the desolvation penalty and the enhancement of the attractive interaction, can be regarded as an effective modeling of the intercalating water molecules at the binding interface, which have a crucial role in strengthening the binding affinity [46]. Practically, to study the binding between proteins using the GB model (with the Coulomb field approximation), which is fast yet tends to overestimate the Coulomb desolvation energy, it is worthwhile to consider adjusting the two parameters in combination, i.e., the lower and the upper limits for the spatial integral of the energy density of the electric field around a charged atom, ρ and r max , respectively. The use of a larger ρ for the Coulomb-bond-forming hydrogen and oxygen atoms strengthens the Coulomb attraction in the short range. The use of a relatively small r max alleviates the overestimated desolvation penalty in the long range. In the case of the binding between barnase and barstar, the above-mentioned two-way adjustment was found to work well for the intermolecular binding driven by the Coulomb attraction. Investigating whether the same two-way strategy can be applied to other systems, as well as the investigation of the contribution of the non-Coulomb energy to the binding, remains for a future work.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24054700/s1.

Author Contributions

Design of this research, D.P. and M.T.; development of computational method, D.P. and M.T.; validation, D.P. and M.T.; conducting calculations, D.P.; analysis and data interpretation, D.P. and M.T.; writing the manuscript, D.P. and M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by Grant-in-Aids from JSPS (JP18H05264 and JP21H02451), by the Top Global University Project from MEXT, and by the Waseda University Grant for Special Research Project.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in Supplementary Material.

Acknowledgments

We thank Y. Mizuhara for the valuable discussion and M. Yaguchi for the assistance with figure drawing. D.P. was supported by the Leading Graduate Program in Science and Engineering at Waseda University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sheinerman, F.B.; Norel, R.; Honig, B. Electrostatic Aspects of Protein-Protein Interactions. Curr. Opin. Struct. Biol. 2000, 10, 153–159. [Google Scholar] [CrossRef] [PubMed]
  2. Schreiber, G.; Haran, G.; Zhou, H.-X. Fundamental Aspects of Protein-Protein Association Kinetics. Chem. Rev. 2009, 109, 839–860. [Google Scholar] [CrossRef] [Green Version]
  3. Born, M. Volumen und hydratationswärme der Ionen. Z. Phys. 1920, 1, 45–48. [Google Scholar] [CrossRef] [Green Version]
  4. Still, W.C.; Tempczyk, A.; Hawley, R.C.; Hendrickson, T. Semianalytical Treatment of Solvation for Molecular Mechanics and Dynamics. J. Am. Chem. Soc. 1990, 112, 6127–6129. [Google Scholar] [CrossRef]
  5. Hawkins, G.D.; Cramer, C.J.; Truhlar, D.G. Pairwise Solute Descreening of Solute Charges from a Dielectric Medium. Chem. Phys. Lett. 1995, 246, 122–129. [Google Scholar] [CrossRef]
  6. Hawkins, G.D.; Cramer, C.J.; Truhlar, D.G. Parametrized Models of Aqueous Free Energies of Solvation Based on Pairwise Descreening of Solute Atomic Charges from a Dielectric Medium. J. Phys. Chem. 1996, 100, 19824–19839. [Google Scholar] [CrossRef]
  7. Qiu, D.; Shenkin, P.S.; Hollinger, F.P.; Still, W.C. The GB/SA Continuum Model for Solvation. A Fast Analytical Method for the Calculation of Approximate Born Radii. J. Phys. Chem. A 1997, 101, 3005–3014. [Google Scholar] [CrossRef]
  8. Schaefer, M.; Karplus, M. A Comprehensive Analytical Treatment of Continuum Electrostatics. J. Phys. Chem. 1996, 100, 1578–1599. [Google Scholar] [CrossRef]
  9. Dominy, B.N.; Brooks, C.L., III. Development of a Generalized Born Model Parametrization for Proteins and Nucleic Acids. J. Phys. Chem. B 1999, 103, 3765–3773. [Google Scholar] [CrossRef]
  10. Bashford, D.; Case, D.A. Generalized Born Models of Macromolecular Solvation Effects. Annu. Rev. Phys. Chem. 2000, 51, 129–152. [Google Scholar] [CrossRef]
  11. Tsui, V.; Case, D.A. Molecular Dynamics Simulations of Nucleic Acids with a Generalized Born Solvation Model. J. Am. Chem. Soc. 2000, 122, 2489–2498. [Google Scholar] [CrossRef]
  12. Tsui, V.; Case, D.A. Theory and Applications of The Generalized Born Solvation Model in Macromolecular Simulations. Biopolymers 2001, 56, 275–291. [Google Scholar] [CrossRef]
  13. Lee, M.S.; Salsbury, F.R.H.; Brooks, C.L., III. Novel Generalized Born Methods. J. Chem. Phys. 2002, 116, 10606–10614. [Google Scholar] [CrossRef]
  14. Lee, M.S.; Feig, M.; Salsbury, F.R., Jr.; Brooks, C.L., III. New Analytic Approximation to the Standard Molecular Volume Definition and Its Application to Generalized Born Calculations. J. Comput. Chem. 2003, 24, 1348–1356. [Google Scholar] [CrossRef]
  15. Im, W.; Lee, M.S.; Brooks, C.L., III. Generalized Born Model with a Simple Smoothing Function. J. Comput. Chem. 2003, 24, 1691–1702. [Google Scholar] [CrossRef]
  16. Onufriev, A.; Bashford, D.; Case, D.A. Exploring Protein Native States and Large-scale Conformational Changes with a Modified Generalized Born Model. Proteins Struct. Funct. Bioinf. 2004, 55, 383–394. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Feig, M.; Onufriev, A.; Lee, M.S.; Im, W.; Case, D.A.; Brooks, C.L., III. Performance Comparison of Generalized Born and Poisson Methods in the Calculation of Electrostatic Solvation Energies for Protein Structures. J. Comput. Chem. 2004, 25, 265–284. [Google Scholar] [CrossRef]
  18. Gallicchio, E.; Levy, R.M. AGBNP: An Analytic Implicit Solvent Model Suitable for Molecular Dynamics Simulations and High-Resolution Modeling. J. Comput. Chem. 2004, 25, 479–499. [Google Scholar] [CrossRef]
  19. Im, W.; Chen, J.; Brooks, C.L. Peptide and Protein Folding and Conformational Equilibria: Theoretical Treatment of Electrostatics and Hydrogen Bonding with Implicit Solvent Models. Adv. Protein Chem. 2005, 72, 173–198. [Google Scholar]
  20. Chen, J.; Im, W.; Brooks, C.L. Balancing Solvation and Intramolecular Interactions: Toward a Consistent Generalized Born Force Field. J. Am. Chem. Soc. 2006, 128, 3728–3736. [Google Scholar] [CrossRef] [Green Version]
  21. Mongan, J.; Simmerling, C.; McCammon, J.A.; Case, D.A.; Onufriev, A. Generalized Born Model with a Simple, Robust Molecular Volume Correction. J. Chem. Theory Comput. 2007, 3, 156–169. [Google Scholar] [CrossRef]
  22. Aguilar, B.; Shadrach, R.; Onufriev, A.V. Reducing the Secondary Structure Bias in the Generalized Born Model via R6 Effective Radii. J. Chem. Theory Comput. 2010, 6, 3613–3630. [Google Scholar] [CrossRef]
  23. Nguyen, H.; Roe, D.R.; Simmerling, C. Improved Generalized Born Solvent Model Parameters for Protein Simulations. J. Chem. Theory Comput. 2013, 9, 2020–2034. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Mizuhara, Y.; Parkin, D.; Umezawa, K.; Ohnuki, J.; Takano, M. Over-Destabilization of Protein-Protein Interaction in Generalized Born Model and Utility of Energy Density Integration Cutoff. J. Phys. Chem. B 2017, 121, 4669–4677. [Google Scholar] [CrossRef] [PubMed]
  25. Bondi, A. van der Waals volumes and radii. J. Phys. Chem. 1964, 68, 441–451. [Google Scholar] [CrossRef]
  26. Felts, A.K.; Gallicchio, E.; Chekmarev, D.; Paris, K.A.; Friesner, R.A.; Levy, R.M. Prediction of protein loop conformations using the AGBNP implicit solvent model and torsion angle sampling. J. Chem. Theory Comput. 2008, 4, 855–868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Geney, R.; Layten, M.; Gomperts, R.; Hornak, V.; Simmerling, C. Investigation of Salt Bridge Stability in a Generalized Born Solvent Model. J. Chem. Theory Comput. 2006, 2, 115–127. [Google Scholar] [CrossRef] [PubMed]
  28. Parkin, D.; Takano, M. Coulombic Organization in Membrane-Embedded Rotary Motor of ATP synthase. J. Phys. Chem. B 2023, 127, 1552–1562. [Google Scholar] [CrossRef]
  29. Kamiyama, Y.; Parkin, D.; Takano, M. Torque Generation Mechanism in Fo Motor of ATP Synthase Elucidated by Free-energy and Coulomb-energy Landscapes along the c-ring Rotation. Biochem. Biophys. Res. Commun. 2023, 651, 56–61. [Google Scholar] [CrossRef]
  30. Buckle, A.M.; Schreiber, G.; Fersht, A.R. Protein-Protein Recognition: Crystal Structural Analysis of Barnase-Barstar Complex at 2.0-Å Resolution. Biochemistry 1994, 33, 8878–8889. [Google Scholar] [CrossRef]
  31. Dong, F.; Vijayakumar, M.; Zhou, H.-X. Comparison of Calculation and Experiment Implicates Significant Electrostatic Contributions to the Binding Stability of Barnase and Barstar. Biophys. J. 2003, 85, 49–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Wang, T.; Tomic, S.; Gabdoulline, R.R.; Wade, R.C. How Optimal are the Binding Energetics of Barnase and Barstar? Biophys. J. 2004, 87, 1618–1630. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Gumbart, J.C.; Roux, B.; Chipot, C. Efficient Determination of Protein-Protein Standard Binding Free Energies from First Principles. J. Chem. Theory Comput. 2013, 9, 3789–3798. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Baker, N.A.; Sept, D.; Joseph, S.; Holst, M.J.; McCammon, J.A. Electrostatics of Nanosystems: Application to Microtubules and the Ribosome. Proc. Natl. Acad. Sci. USA 2001, 98, 10037–10041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Constanciel, R.; Contreras, R. Self Consistent Field Theory of Solvent Effects Representation by Continuum Models: Introduction of Desolvation Contribution. Theor. Chim. Acta 1984, 65, 1–11. [Google Scholar] [CrossRef]
  36. Cornell, W.D.; Cieplak, P.; Bayly, C.I.; Gould, I.R.; Merz, K.M.; Ferguson, D.M.; Spellmeyer, D.C.; Fox, T.; Caldwell, J.W.; Kollman, P.A. A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules. J. Am. Chem. Soc. 1995, 117, 5179–5197. [Google Scholar] [CrossRef] [Green Version]
  37. Hornak, V.; Abel, R.; Okur, A.; Strockbine, B.; Roitberg, A.; Simmerling, C. Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins: Struct. Funct. Bioinform. 2006, 65, 712–725. [Google Scholar] [CrossRef] [Green Version]
  38. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  39. Essmann, U.; Perera, L.; Berkowitz, M.L.; Darden, T.; Lee, H.; Pedersen, L.G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 8577–8593. [Google Scholar] [CrossRef] [Green Version]
  40. Berendsen, H.J.C.; Postma, J.P.M.; van Gunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular Dynamics with Coupling to an External Bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef] [Green Version]
  41. Case, D.A.; Darden, T.A.; Cheatham, T.E., III; Simmerling, C.L.; Wang, J.; Duke, R.E.; Luo, R.; Walker, R.C.; Zhang, W.; Merz, K.M.; et al. AMBER12; University of California: San Francisco, CA, USA, 2012. [Google Scholar]
  42. Joung, I.S.; Cheatham, T.E., III. Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations. J. Phys. Chem. B 2008, 112, 9020–9041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Duan, Y.; Wu, C.; Chowdhury, S.; Lee, M.C.; Xiong, G.; Zhang, W.; Yang, R.; Cieplak, P.; Luo, R.; Lee, T.; et al. A Point-Charge Force Field for Molecular Mechanics Simulations of Proteins Based on Condensed-Phase Quantum Mechanical Calculations. J. Comput. Chem. 2003, 24, 1999–2012. [Google Scholar] [CrossRef] [PubMed]
  44. Eswar, N.; Webb, B.; Marti-Renom, M.A.; Madhusudhan, M.S.; Eramian, D.; Shen, M.Y.; Pieper, U.; Šali, A. Comparative Protein Structure Modeling Using MODELLER. Curr. Protoc. Protein Sci. 2007. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kumar, S.; Bouzida, D.; Swendsen, R.H.; Kollman, P.A.; Rosenberg, J.M. The Weighted Histogram Analysis Method for Free-energy Calculations on Biomolecules. I. The Method. J. Comput. Chem. 1992, 13, 1011–1021. [Google Scholar] [CrossRef]
  46. Wong, S.; Amaro, R.E.; McCammon, J.A. MM-PBSA Captures Key Role of Intercalating Water Molecules at a Protein-Protein Interface. J. Chem. Theory Comput. 2009, 5, 422–429. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Effect Author request: Please enlarge the figure size for visibility (we changed the ration from 0.8 to 0.95 of ρ H on Coulomb energies in H–O system. (a) Illustration of the Coulomb bond formation. Dark-colored regions indicate the regions that are subjected to the integral in Equation (11) and are affected by the increase in ρ H . (b) Total energy [i.e., the sum of the Coulomb and van der Waals (Lennard–Jones) energies] is shown as a function of the separation distance d between the two atoms at different ρ H values (1.1 to 1.5 Å). The dotted line represents the total energy where Equation (3) is used for the Coulomb energy ( ε w = 78.5 ). (c) Desolvation Δ G O self , (d) desolvation energy Δ G H self , and (e) interaction energy G int are shown as a function of d. Line colors correspond to those in (b). (f) Desolvation energies [ Δ G O self (cyan), Δ G H self (magenta), and Δ G O self + Δ G H self (open square)] at d = 1.57 Å are shown as a function of ρ H . (g) Interaction energy G int at d = 1.57 Å is shown as a function of ρ H . Energy unit is kcal/mol.
Figure 1. Effect Author request: Please enlarge the figure size for visibility (we changed the ration from 0.8 to 0.95 of ρ H on Coulomb energies in H–O system. (a) Illustration of the Coulomb bond formation. Dark-colored regions indicate the regions that are subjected to the integral in Equation (11) and are affected by the increase in ρ H . (b) Total energy [i.e., the sum of the Coulomb and van der Waals (Lennard–Jones) energies] is shown as a function of the separation distance d between the two atoms at different ρ H values (1.1 to 1.5 Å). The dotted line represents the total energy where Equation (3) is used for the Coulomb energy ( ε w = 78.5 ). (c) Desolvation Δ G O self , (d) desolvation energy Δ G H self , and (e) interaction energy G int are shown as a function of d. Line colors correspond to those in (b). (f) Desolvation energies [ Δ G O self (cyan), Δ G H self (magenta), and Δ G O self + Δ G H self (open square)] at d = 1.57 Å are shown as a function of ρ H . (g) Interaction energy G int at d = 1.57 Å is shown as a function of ρ H . Energy unit is kcal/mol.
Ijms 24 04700 g001
Figure 2. Author request: Please enlarge the figure size for visibility (we changed the ration from 0.8 to 0.95 Effect of ρ H on Coulomb energies in Arg–Glu system. (a) Illustration of the system. Encircled hydrogen atoms in the guanine group of Arg are the atoms, the intrinsic radii of which ( ρ H ) were uniformly and systematically changed from 1.1 Å to 1.5 Å. (b) Potential of mean force (PMF) is shown as a function of d [see (a)] for the explicit water model (black) and for the GB solvation model at different ρ H values (1.1 to 1.5 Å). (c) Total Coulomb energy ( Δ G self + G int ), (d) total desolvation energy ( Δ G self ), and (e) total Coulomb interaction energy ( G int ) are shown as a function of d at different ρ H values [line colors correspond to those in (c)]. (f) Total desolvation energy ( Δ G self ) and (g) total interaction energy ( G int ) averaged in the range of 3.9 d 4.0 Å are shown as a function of ρ H . Energy unit is kcal/mol.
Figure 2. Author request: Please enlarge the figure size for visibility (we changed the ration from 0.8 to 0.95 Effect of ρ H on Coulomb energies in Arg–Glu system. (a) Illustration of the system. Encircled hydrogen atoms in the guanine group of Arg are the atoms, the intrinsic radii of which ( ρ H ) were uniformly and systematically changed from 1.1 Å to 1.5 Å. (b) Potential of mean force (PMF) is shown as a function of d [see (a)] for the explicit water model (black) and for the GB solvation model at different ρ H values (1.1 to 1.5 Å). (c) Total Coulomb energy ( Δ G self + G int ), (d) total desolvation energy ( Δ G self ), and (e) total Coulomb interaction energy ( G int ) are shown as a function of d at different ρ H values [line colors correspond to those in (c)]. (f) Total desolvation energy ( Δ G self ) and (g) total interaction energy ( G int ) averaged in the range of 3.9 d 4.0 Å are shown as a function of ρ H . Energy unit is kcal/mol.
Ijms 24 04700 g002
Figure 3. Effect of r max and ρ H on Coulomb energies in Bn–Bs system. (a) Illustration of the system. Coulomb potential at the surface of Bn and Bs is shown by the color map [red: 10 k T / e ; blue + 10 k T / e , where e is the elementary charge, k is the Boltzmann constant, and T is temperature (300 K)]. Binding interfaces of Bn and Bs (rotated by 90 ) are displayed beneath. Coulomb potential was calculated by APBS [34]. (b) Potential of mean force (PMF) is shown as a function of the intermolecular distance d [see (a)] obtained by conducting an umbrella sampling with the explicit water. (c) Total Coulomb energy is shown as a function of d at different r max values ( r max = 25, 15, and 10 Å). (d) Total Coulomb energy is shown as a function of d at different ρ H values ( ρ H = 1.3 and 1.4 Å). The asterisk (*) indicates that the intrinsic radius of the carboxyl oxygen ( ρ O ) in the side chain of acidic residues is also changed (from 1.5 to 1.6 Å). (e) Total desolvation energy ( Δ G self ) and (f) total interaction energy ( G int ), averaged in the range of 23.6 d 23.8 Å, are shown as a function of ρ H . Energy unit is kcal/mol.
Figure 3. Effect of r max and ρ H on Coulomb energies in Bn–Bs system. (a) Illustration of the system. Coulomb potential at the surface of Bn and Bs is shown by the color map [red: 10 k T / e ; blue + 10 k T / e , where e is the elementary charge, k is the Boltzmann constant, and T is temperature (300 K)]. Binding interfaces of Bn and Bs (rotated by 90 ) are displayed beneath. Coulomb potential was calculated by APBS [34]. (b) Potential of mean force (PMF) is shown as a function of the intermolecular distance d [see (a)] obtained by conducting an umbrella sampling with the explicit water. (c) Total Coulomb energy is shown as a function of d at different r max values ( r max = 25, 15, and 10 Å). (d) Total Coulomb energy is shown as a function of d at different ρ H values ( ρ H = 1.3 and 1.4 Å). The asterisk (*) indicates that the intrinsic radius of the carboxyl oxygen ( ρ O ) in the side chain of acidic residues is also changed (from 1.5 to 1.6 Å). (e) Total desolvation energy ( Δ G self ) and (f) total interaction energy ( G int ), averaged in the range of 23.6 d 23.8 Å, are shown as a function of ρ H . Energy unit is kcal/mol.
Ijms 24 04700 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Parkin, D.; Takano, M. The Intrinsic Radius as a Key Parameter in the Generalized Born Model to Adjust Protein-Protein Electrostatic Interaction. Int. J. Mol. Sci. 2023, 24, 4700. https://doi.org/10.3390/ijms24054700

AMA Style

Parkin D, Takano M. The Intrinsic Radius as a Key Parameter in the Generalized Born Model to Adjust Protein-Protein Electrostatic Interaction. International Journal of Molecular Sciences. 2023; 24(5):4700. https://doi.org/10.3390/ijms24054700

Chicago/Turabian Style

Parkin, Dan, and Mitsunori Takano. 2023. "The Intrinsic Radius as a Key Parameter in the Generalized Born Model to Adjust Protein-Protein Electrostatic Interaction" International Journal of Molecular Sciences 24, no. 5: 4700. https://doi.org/10.3390/ijms24054700

APA Style

Parkin, D., & Takano, M. (2023). The Intrinsic Radius as a Key Parameter in the Generalized Born Model to Adjust Protein-Protein Electrostatic Interaction. International Journal of Molecular Sciences, 24(5), 4700. https://doi.org/10.3390/ijms24054700

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop