Stanley K. Burt

CALCULATION OF THE AQUEOUS SOLVATION FREE ENERGY OF THE PROTON

The value of the proton hydration free energy, ΔGhyd(H+), has been quoted in the literature to be from −252.6 to −262.5 kcal/mol. In this article, we present a theoretical model for calculating the hydration free energy of ions in aqueous solvent and use this model to calculate the proton hydration free energy, ΔGhyd(H+), in an effort to resolve the uncertainty concerning its exact value. In the model we define ΔGhyd(H+) as the free energy change associated with the following process: ΔG[H+(gas)+H2nOn(aq)→H+(H2nOn)(aq)], where the solvent is represented by a neutral n-water cluster embedded in a dielectric continuum and the solvated proton is represented by a protonated n-water cluster also in the continuum. All solvated species are treated as quantum mechanical solutes coupled to a dielectric continuum using a self consistent reaction field cycle. We investigated the behavior of ΔGhyd(H+) as the number of explicit waters of hydration is increased from n=1 to n=6. As n increases from 1 to 3, the hydration...

Flap opening in HIV-1 protease simulated by ‘activated’ molecular dynamics

We have used an ‘activated’ molecular dynamics approach to simulate flap opening in HIV-1 protease. An initial impulse for flap opening was provided by applying harmonic constraints to non-flap residues. After an initial ‘melting’ phase, the two β-hairpin structures that constitute the flaps opened to a 25 Å gap within 200 ps of simulation. Analysis of backbone torsion angles suggests that flap opening is related to conformational changes at Lys 45, Met 46, Gly 52 and Phe 53. In contrast, similar molecular dynamics simulations on the M461 mutant, which is associated with drug resistance, indicates that this mutation stabilizes the flaps in a closed conformation.

On the structure and thermodynamics of solvated monoatomic ions using a hybrid solvation model

The hydration free energies relative to that of the proton are calculated for a representative set of monatomic ions Z±. These include cationic forms of the alkali earth elements Li, Na, and K, and anionic forms of the halogens F, Cl, and Br. In the current model the relative ion hydration free energy is defined as Δ[ΔGhyd(Z±)]=G(Z±[H2O]n(aq))−G(H+[H2O]n(aq))−G(Z±(gas))−G(H+(gas)), where the solvated ions are represented by ion–water clusters coupled to a dielectric continuum using a self-consistent reaction field cycle. An investigation of the behavior of Δ[ΔGhyd(Z±)] as the number of explicit waters of hydration is increased reveals convergence by n=4. This convergence indicates that the free energy change for the addition of water to a solvated proton–water complex is the same as the free energy change associated with the addition of water to a solvated Z±–water complex. This is true as long as there are four explicitly solvating waters associated with the ion. This convergence is independent of the ty...

QM/MM modeling the Ras–GAP catalyzed hydrolysis of guanosine triphosphate

The mechanism of the hydrolysis reaction of guanosine triphosphate (GTP) by the protein complex Ras–GAP (p21ras – p120GAP) has been modeled by the quantum mechanical—molecular mechanical (QM/MM) and ab initio quantum calculations. Initial geometry configurations have been prompted by atomic coordinates of a structural analog (PDBID:1WQ1). It is shown that the minimum energy reaction path is consistent with an assumption of two‐step chemical transformations. At the first stage, a unified motion of Arg789 of GAP, Gln61, Thr35 of Ras, and the lytic water molecule results in a substantial spatial separation of the γ‐phosphate group of GTP from the rest of the molecule (GDP). This phase of hydrolysis process proceeds through the low‐barrier transition state TS1. At the second stage, Gln61 abstracts and releases protons within the subsystem including Gln61, the lytic water molecule and the γ‐phosphate group of GTP through the corresponding transition state TS2. Direct quantum calculations show that, in this particular environment, the reaction GTP + H2O → GDP + H2PO  4− can proceed with reasonable activation barriers of less than 15 kcal/mol at every stage. This conclusion leads to a better understanding of the anticatalytic effect of cancer‐causing mutations of Ras, which has been debated in recent years. Proteins 2005.

Mechanisms of guanosine triphosphate hydrolysis by Ras and Ras-GAP proteins as rationalized by ab initio QM/MM simulations

The hydrolysis reaction of guanosine triphosphate (GTP) by p21ras (Ras) has been modeled by using the ab initio type quantum mechanical–molecular mechanical simulations. Initial geometry configurations have been prompted by atomic coordinates of the crystal structure (PDBID: 1QRA) corresponding to the prehydrolysis state of Ras in complex with GTP. Multiple searches of minimum energy geometry configurations consistent with the hydrogen bond networks have been performed, resulting in a series of stationary points on the potential energy surface for reaction intermediates and transition states. It is shown that the minimum energy reaction path is consistent with an assumption of a two‐step mechanism of GTP hydrolysis. At the first stage, a unified action of the nearest residues of Ras and the nearest water molecules results in a substantial spatial separation of the γ‐phosphate group of GTP from the rest of the molecule (GDP). This phase of hydrolysis process proceeds through the low barrier (16.7 kcal/mol) transition state TS1. At the second stage, the inorganic phosphate is formed in consequence of proton transfers mediated by two water molecules and assisted by the Gln61 residue from Ras. The highest transition state at this segment, TS3, is estimated to have an energy 7.5 kcal/mol above the enzyme–substrate complex. The results of simulations are compared to the previous findings for the GTP hydrolysis in the Ras‐GAP (p21ras–p120GAP) protein complex. Conclusions of the modeling lead to a better understanding of the anticatalytic effect of cancer causing mutation of Gln61 from Ras, which has been debated in recent years. Proteins 2007.

Reaction path and free energy calculations of the transition between alternate conformations of HIV-1 protease†

Two different structures of ligand‐free HIV protease have been determined by X‐ray crystallography. These structures differ in the position of two 12 residue, β‐hairpin regions (or “flaps”) which cap the active site. The movements of the flaps must be involved in the binding of substrates since, in either conformation, the flaps block the binding site. One of these structures is similar to structures of the ligand‐bound enzyme; however, the importance of both structures to enzyme function is unclear. This transformation takes place on a time scale too long for conventional molecular dynamics simulations, so the process was studied by first identifying a reaction path between the two structures and then calculating the free energy along this path using umbrella sampling. For the ligand‐free enzyme, it is found that the two structures are nearly equally stable, with the ligand‐bound‐type structure being less stable, consistent with X‐ray crystallography data. The more stable open structure does not have a lower potential energy, but is stabilized by entropy. The transition occurs through a collapse and reformation of the β‐sheet structure of the conformationally flexible, glycine‐rich flap ends. Additionally, some problems in studying conformational changes in proteins through the use of a single reaction path are addressed. Proteins 32:7–16, 1998.

Mechanism of the myosin catalyzed hydrolysis of ATP as rationalized by molecular modeling

The intrinsic chemical reaction of adenosine triphosphate (ATP) hydrolysis catalyzed by myosin is modeled by using a combined quantum mechanics and molecular mechanics (QM/MM) methodology that achieves a near ab initio representation of the entire model. Starting with coordinates derived from the heavy atoms of the crystal structure (Protein Data Bank ID code 1VOM) in which myosin is bound to the ATP analog ADP·VO4−, a minimum-energy path is found for the transformation ATP + H2O → ADP + Pi that is characterized by two distinct events: (i) a low activation-energy cleavage of the PγOβγ bond and separation of the γ-phosphate from ADP and (ii) the formation of the inorganic phosphate as a consequence of proton transfers mediated by two water molecules and assisted by the Glu-459–Arg-238 salt bridge of the protein. The minimum-energy model of the enzyme–substrate complex features a stable hydrogen-bonding network in which the lytic water is positioned favorably for a nucleophilic attack of the ATP γ-phosphate and for the transfer of a proton to stably bound second water. In addition, the PγOβγ bond has become significantly longer than in the unbound state of the ATP and thus is predisposed to cleavage. The modeled transformation is viewed as the part of the overall hydrolysis reaction occurring in the closed enzyme pocket after ATP is bound tightly to myosin and before conformational changes preceding release of inorganic phosphate.

Flexible effective fragment QM/MM method: Validation through the challenging tests

A new version of the QM/MM method, which is based on the effective fragment potential (EFP) methodology [Gordon, M. et al., J Phys Chem A 2001, 105, 293] but allows flexible fragments, is verified through calculations of model molecular systems suggested by different authors as challenging tests for QM/MM approaches. For each example, the results of QM/MM calculations for a partitioned system are compared to the results of an all‐electron ab initio quantum chemical study of the entire system. In each case we were able to achieve approximately similar or better accuracy of the QM/MM results compared to those described in original publications. In all calculations we kept the same set of parameters of our QM/MM scheme. A new test example is considered when calculating the potential of internal rotation in the histidine dipeptide around the CαCβ side chain bond.

An exhaustive DNA micro-satellite map of the human genome using high performance computing.

The current pace of the generation of sequence data requires the development of software tools that can rapidly provide full annotation of the data. We have developed a new method for rapid sequence comparison using the exact match algorithm without repeat masking. As a demonstration, we have identified all perfect simple tandem repeats (STR) within the draft sequence of the human genome. The STR elements (chromosome, position, length and repeat subunit) have been placed into a relational database. Repeat flanking sequence is also publicly accessible at http://grid.abcc.ncifcrf.gov. To illustrate the utility of this complete set of STR elements, we documented the increased density of potentially polymorphic markers throughout the genome. The new STR markers may be useful in disease association studies because so many STR elements manifest multiallelic polymorphism. Also, because triplet repeat expansions are important for human disease etiology, we identified trinucleotide repeats that exist within exons of known genes. This resulted in a list that includes all 14 genes known to undergo polynucleotide expansion, and 48 additional candidates. Several of these are non-polyglutamine triplet repeats. Other examinations of the STR database demonstrated repeats spanning splice junctions and identified SNPs within repeat elements.

The use of urine proteomic and metabonomic patterns for the diagnosis of interstitial cystitis and bacterial cystitis

The advent of systems biology approaches that have stemmed from the sequencing of the human genome has led to the search for new methods to diagnose diseases. While much effort has been focused on the identification of disease-specific biomarkers, recent efforts are underway toward the use of proteomic and metabonomic patterns to indicate disease. We have developed and contrasted the use of both proteomic and metabonomic patterns in urine for the detection of interstitial cystitis (IC). The methodology relies on advanced bioinformatics to scrutinize information contained within mass spectrometry (MS) and high-resolution proton nuclear magnetic resonance (1H-NMR) spectral patterns to distinguish IC-affected from non-affected individuals as well as those suffering from bacterial cystitis (BC). We have applied a novel pattern recognition tool that employs an unsupervised system (self-organizing-type cluster mapping) as a fitness test for a supervised system (a genetic algorithm). With this approach, a training set comprised of mass spectra and 1H-NMR spectra from urine derived from either unaffected individuals or patients with IC is employed so that the most fit combination of relative, normalized intensity features defined at precise m/z or chemical shift values plotted in n-space can reliably distinguish the cohorts used in training. Using this bioinformatic approach, we were able to discriminate spectral patterns associated with IC-affected, BC-affected, and unaffected patients with a success rate of approximately 84%.