Search PubMed⌕ Search

PubMed · 12594794

Transferable atom equivalent multicentered multipole expansion method.

Abstract

The transferability of atomic and functional group properties is an implicit concept in chemistry. The work presented here describes the use of Transferable Atom Equivalents (TAE) to represent molecular electrostatic potential fields through the use of integrated atomic multipole moments that are associated with each TAE atom type used in the reconstruction. TAE molecular surface distributions of electrostatic potentials are compared with analytical ab initio and empirical (Gasteiger) partial charge reference models for several conformations of test peptides. Surface electrostatic potential distributions computed using TAE multipole representations were found to converge at the octopole level, with incremental improvement observed when hexadecapoles were included. Molecular electrostatic potential fields that were produced using the TAE method were observed to be responsive to conformational changes and to compare well with ab initio reference distributions. Generation of TAE atom types and their associated multipoles does not involve fitting to sample electrostatic potential fields, but rather utilizes integrated AIM atomic electron density distributions within representative chemical environments. The RECON program was used for TAE reconstruction. RECON is capable of processing 5,000 drug-sized molecules or 25 proteins per minute per 1.7 GHz P4 Linux processor.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

C E Whitehead, C M Breneman, N Sukumar, M D Ryan. 2003. Transferable atom equivalent multicentered multipole expansion method.. https://doi.org/10.1002/jcc.10240

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

KEEP EXPLORING

Related citations

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↗

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↗

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↗