How well does Poisson-Boltzmann implicit solvent agree with explicit solvent? A quantitative analysis.

Chunhu Tan; Lijiang Yang; Ray Luo

doi:10.1021/jp063479b

Home/Archive/PubMed

PubMed · 16970499

How well does Poisson-Boltzmann implicit solvent agree with explicit solvent? A quantitative analysis.

Chunhu Tan·Lijiang Yang·Ray Luo

View original source ↗Export metadata ↓

Abstract

We have quantitatively studied the performance of a finite-difference Poisson-Boltzmann implicit solvent with respect to the TIP3P explicit solvent in a range of systems of biochemical interest. An overall agreement was found between the tested implicit and explicit solvents for hydrogen-bonding/salt-bridging dimers and peptide monomers and dimers of different conformations and different lengths. These comparative analyses also indicate a good transferability of empirically optimized parameters for the implicit solvent from small training molecules to large testing peptides. However, deviations between the two tested solvents are also apparent. Specifically, a consistent deviation was observed when hydrogen-bonding or salt-bridging dimers are within 4-6 A. The deviation reaches a maximum at about 5.5 A, the so-called water-bridging distance. The tested implicit solvent, even with optimized parameters, cannot capture the subtle fluctuation in the distance-dependent reaction field energy profiles, although smoothed profiles can still be obtained and are in overall agreement with those in the explicit solvent. Interestingly, the same mechanism underlining the above discrepancy is also responsible for the larger deviations of certain peptide conformations, such as parallel beta-strand dimers. It is likely that the observed discrepancy may cause improper conformational distributions in simulations with the implicit solvent when hydrogen-bonding or salt-bridging interactions are crucial, such as secondary structure populations in proteins. Validation of the implicit solvent with optimized parameters in dynamics simulations will be the next step to study the influences of the observed discrepancy at biological conditions.

Explore related subjects

Alanine Computer Simulation Dimerization Hydrogen Bonding Models, Chemical Models, Molecular Models, Statistical Models, Theoretical Molecular Conformation Peptides Protein Structure, Secondary Salts Solvents Thermodynamics Water Journal Article Research Support, N.I.H., Extramural

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Chunhu Tan, Lijiang Yang, Ray Luo. 2006-09-21. How well does Poisson-Boltzmann implicit solvent agree with explicit solvent? A quantitative analysis.. https://doi.org/10.1021/jp063479b

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

PubMed2007

An improved algorithm for analytical gradient evaluation in resolution-of-the-identity second-order Møller-Plesset perturbation theory: application to alanine tetrapeptide conformational analysis.

We present a new algorithm for analytical gradient evaluation in resolution-of-the-identity second-order Møller-Plesset perturbation theory (RI-MP2) and thoroughly assess its computational performance and chemical accuracy. This algorithm addresses the potential I/O bottlenecks associated with disk-based storage and access of the RI-MP2 t-amplitudes by utilizing a semi-direct batching approach and yields computational speed-ups of approximately 2-3 over the best conventional MP2 analytical gradient algorithms. In addition, we attempt to provide a straightforward guide to performing reliable and cost-efficient geometry optimizations at the RI-MP2 level of theory. By computing relative atomization energies for the G3/99 set and optimizing a test set of 136 equilibrium molecular structures, we demonstrate that satisfactory relative accuracy and significant computational savings can be obtained using Pople-style atomic orbital basis sets with the existing auxiliary basis expansions for RI-MP2 computations. We also show that RI-MP2 geometry optimizations reproduce molecular equilibrium structures with no significant deviations (>0.1 pm) from the predictions of conventional MP2 theory. As a chemical application, we computed the extended-globular conformational energy gap in alanine tetrapeptide at the extrapolated RI-MP2/cc-pV(TQ)Z level as 2.884, 4.414, and 4.994 kcal/mol for structures optimized using the HF, DFT (B3LYP), and RI-MP2 methodologies and the cc-pVTZ basis set, respectively. These marked energetic discrepancies originate from differential intramolecular hydrogen bonding present in the globular conformation optimized at these levels of theory and clearly demonstrate the importance of long-range correlation effects in polypeptide conformational analysis.

Alanine↗

PubMed2007

Structural and functional effects of disease-causing amino acid substitutions affecting residues Ala72 and Glu76 of the protein tyrosine phosphatase SHP-2.

Mutations of the protein tyrosine phosphatase SHP-2 are implicated in human diseases, causing Noonan syndrome (NS) and related developmental disorders or contributing to leukemogenesis depending on the specific amino acid substitution involved. SHP-2 is composed by a catalytic (PTP) and two regulatory (N-SH2 and C-SH2) domains that bind to signaling partners and control the enzymatic activity by limiting the accessibility of the catalytic site. Wild type SHP-2 and four disease-associated mutants recurring in hematologic malignancies (Glu76Lys and Ala72Val) or causing NS (Glu76Asp and Ala72Ser), with affected residues located in the PTP-interacting region of the N-SH2 domain, were analyzed by molecular dynamics simulations and in vitro biochemical assays. Simulations demonstrate that mutations do not affect significantly the conformation of the N-SH2 domain. Rather they destabilize the interaction of this domain with the catalytic site, with more evident effects in the two leukemia associated mutants. Consistent with this structural evidence, mutants exhibit an increased level of basal phosphatase activity in the order Glu76Lys > Ala72Val > Glu76Asp > Ala72Ser > WT. The experimental data also show that the mutants with higher basal activity are more responsive to an activating phosphopeptide. A thermodynamic analysis demonstrates that an increase in the overall phosphopeptide affinity of mutants can be explained by a shift in the equilibrium between the inactive and active SHP-2 structure. These data support the view that an increase in the affinity of SHP-2 for its binding partners, caused by destabilization of the closed, inactive conformation, rather than protein basal activation per se, would represent the molecular mechanism, leading to pathogenesis in these mutants.

Alanine↗

PubMed2007

How a small change in retinal leads to G-protein activation: initial events suggested by molecular dynamics calculations.

Rhodopsin is the prototypical G-protein coupled receptor, coupling light activation with high efficiency to signaling molecules. The dark-state X-ray structures of the protein provide a starting point for consideration of the relaxation from initial light activation to conformational changes that may lead to signaling. In this study we create an energetically unstable retinal in the light activated state and then use molecular dynamics simulations to examine the types of compensation, relaxation, and conformational changes that occur following the cis-trans light activation. The results suggest that changes occur throughout the protein, with changes in the orientation of Helices 5 and 6, a closer interaction between Ala 169 on Helix 4 and retinal, and a shift in the Schiff base counterion that also reflects changes in sidechain interactions with the retinal. Taken together, the simulation is suggestive of the types of changes that lead from local conformational change to light-activated signaling in this prototypical system.

Alanine↗