Biomedical subjects

Mark A Olson

Publications and source records attributed to Mark A Olson.

8 recordsLinked to original sources

PubMed2005

Calculation of absolute protein-ligand binding affinity using path and endpoint approaches.

A comparative analysis is provided of rigorous and approximate methods for calculating absolute binding affinities of two protein-ligand complexes: the FKBP protein bound with small molecules 4-hydroxy-2-butanone and FK506. Our rigorous approach is an umbrella sampling technique where a potential of mean force is determined by pulling the ligand out of the protein active site over several simulation windows. The results of this approach agree well with experimentally observed binding affinities. Also assessed is a commonly used approximate endpoint approach, which separately estimates enthalpy, solvation free energy, and entropy. We show that this endpoint approach has numerous variations, all of which are prone to critical shortcomings. For example, conventional harmonic and quasiharmonic entropy estimation procedures produce disparate results for the relatively simple protein-ligand systems studied in this work.

Binding Sites↗

PubMed2005

Evaluation of Poisson solvation models using a hybrid explicit/implicit solvent method.

Implicit solvent methods have become popular tools in the field of protein dynamics simulations, yet evaluation of their validity has been primarily limited to comparisons with experimental and theoretical data for small molecules. In this paper, we use a recently developed hybrid explicit/implicit solvent methodology to evaluate the accuracy of several Poisson-based implicit solvent models. Specifically, we focus on the calculation of electrostatic solvation free energies of various fixed conformations for two proteins. We show that, among various dielectric boundary definitions, the Lee-Richards molecular surface has the best agreement with hybrid solvent results. Furthermore, certain modifications of the molecular surface Poisson protocol provide varied results. For instance, simple modifications of atomic radii on charged residues generally improve absolute errors but do not significantly reduce relative errors among conformations. On the other hand, using a water-probe radius of 1.0 A, as opposed to the standard value of 1.4 A, to generate the molecular surface, moderately improves both absolute and relative results.

Calibration↗

PubMed2004

Modeling loop reorganization free energies of acetylcholinesterase: a comparison of explicit and implicit solvent models.

The treatment of hydration effects in protein dynamics simulations varies in model complexity and spans the range from the computationally intensive microscopic evaluation to simple dielectric screening of charge-charge interactions. This paper compares different solvent models applied to the problem of estimating the free-energy difference between two loop conformations in acetylcholinesterase. Molecular dynamics (MD) simulations were used to sample potential energy surfaces of the two basins with solvent treated by means of explicit and implicit methods. Implicit solvent methods studied include the generalized Born (GB) model, atomic solvation potential (ASP), and the distance-dependent dieletric constant. By using the linear response approximation (LRA), the explicit solvent calculations determined a free-energy difference that is in excellent agreement with the experimental estimate, while rescoring the protein conformations with GB or the Poisson equation showed inconsistent and inferior results. While the approach of rescoring conformations from explicit water simulations with implicit solvent models is popular among many applications, it perturbs the energy landscape by changing the solvent contribution to microstates without conformational relaxation, thus leading to non-optimal solvation free energies. Calculations applying MD with a GB solvent model produced results of comparable accuracy as observed with LRA, yet the electrostatic free-energy terms were significantly different due to optimization on a potential energy surface favored by an implicit solvent reaction field. The simpler methods of ASP and the distance-dependent scaling of the dielectric constant both produced considerable distortions in the protein internal free-energy terms and are consequently unreliable.

Acetylcholinesterase↗

PubMed2004

Conformational energy landscape of the acyl pocket loop in acetylcholinesterase: a Monte Carlo-generalized Born model study.

The X-ray crystal structure of the reaction product of acetylcholinesterase (AChE) with the inhibitor diisopropylphosphorofluoridate (DFP) showed significant structural displacement in a loop segment of residues 287-290. To understand this conformational selection, a Monte Carlo (MC) simulation study was performed of the energy landscape for the loop segment. A computational strategy was applied by using a combined simulated annealing and room temperature Metropolis sampling approach with solvent polarization modeled by a generalized Born (GB) approximation. Results from thermal annealing reveal a landscape topology of broader basin opening and greater distribution of energies for the displaced loop conformation, while the ensemble average of conformations at 298 K favored a shift in populations toward the native by a free-energy difference in good agreement with the estimated experimental value. Residue motions along a reaction profile of loop conformational reorganization are proposed where Arg-289 is critical in determining electrostatic effects of solvent interaction versus Coulombic charging.

Acetylcholinesterase↗

PubMed2004

Finding a new vaccine in the ricin protein fold.

Previous attempts to produce a vaccine for ricin toxin have been hampered by safety concerns arising from residual toxicity and the undesirable aggregation or precipitation caused by exposure of hydrophobic surfaces on the ricin A-chain (RTA) in the absence of its natural B-chain partner. We undertook a structure-based solution to this problem by reversing evolutionary selection on the 'ribosome inactivating protein' fold of RTA to arrive at a non-functional, compacted single-domain scaffold (sequence RTA1-198) for presentation of a specific protective epitope (RTA loop 95-110). An optimized protein based upon our modeling design (RTA1-33/44-198) showed greater resistance to thermal denaturation, less precipitation under physiological conditions and a reduction in toxic activity of at least three orders of magnitude compared with RTA. Most importantly, RTA1-198 or RTA1-33/44-198 protected 100% of vaccinated animals against supra-lethal challenge with aerosolized ricin. We conclude that comparative protein analysis and engineering yielded a superior vaccine by exploiting a component of the toxin that is inherently more stable than is the parent RTA molecule.

Circular Dichroism↗

PubMed2004

An efficient hybrid explicit/implicit solvent method for biomolecular simulations.

We present a new hybrid explicit/implicit solvent method for dynamics simulations of macromolecular systems. The method models explicitly the hydration of the solute by either a layer or sphere of water molecules, and the generalized Born (GB) theory is used to treat the bulk continuum solvent outside the explicit simulation volume. To reduce the computational cost, we implemented a multigrid method for evaluating the pairwise electrostatic and GB terms. It is shown that for typical ion and protein simulations our method achieves similar equilibrium and dynamical observables as the conventional particle mesh Ewald (PME) method. Simulation timings are reported, which indicate that the hybrid method is much faster than PME, primarily due to a significant reduction in the number of explicit water molecules required to model hydration effects.

Journal Article↗

PubMed2003

Prediction of the multimeric assembly of staphylococcal enterotoxin A with cell-surface protein receptors.

Staphylococcal enterotoxin A (SEA) cross-links two class II major histocompatibility complex (MHC) molecules and forms a multimeric assembly with T-cell receptors (TcRs). The X-ray crystal structure of SEA has been solved, yet details describing molecular recognition and association remain unclear. We present a structural model for the interactions of SEA with cell-surface proteins. Molecular docking calculations predicting SEA association with the class II MHC molecule HLA-DR1 were performed by using a rigid-body docking method. Docked orientations were evaluated by a Poisson-Boltzmann model for the electrostatic free energy of binding and the hydrophobic effect calculated from molecular surface areas. We found that the best-scoring SEA conformers for the DR1alpha interface display a binding mode similar to that determined crystallographically for staphylococcal enterotoxin B bound to HLA-DR1. For the zinc-binding site of SEA, docking DR1beta yielded several orientations exhibiting tetrahedral-like coordination geometries. Combining the two interfaces, tetramers were modeled by docking an alphabeta TcR with trimolecular complexes DR1beta-SEA-DR1alpha and SEA-betaDR1alpha-SEA. Our results indicate that the complex DR1beta-SEA-DR1alpha provides a more favorable assembly for the engagement of TcRs, forming SEA molecular contacts that are in accord with reported mutagenesis studies. In contrast, the cooperative association of two SEA molecules on a single DR1 molecule sterically inhibits interactions with TcRs. We suggest that signal transduction stimulated by SEA through large-scale assembly is limited to four or five TcR-(DR1beta-SEA-DR1alpha) tetramers and requires the dimerization of class II MHC molecules, while TcR dimerization is unlikely.

Enterotoxins↗

PubMed2002

Free-energy contributions to complex formation between botulinum neurotoxin type B and synaptobrevin fragment.

Free-energy terms that contribute to complex formation between the catalytic domain of botulinum neurotoxin type B (BoNT/B-L(C)) and a 36-residue synaptobrevin fragment were estimated by using a combination of microscopic simulations and continuum methods. The complex for a non-hydrolyzed substrate was calculated by optimizing an energy function applied to the X-ray co-crystal structure of BoNT/B-L(C) bound with reaction products from a cleaved synaptobrevin peptide, refined to high crystallographic thermal factors. The estimated absolute binding affinity of the simulation structure is in good qualitative agreement with the experimental free energy of Michaelis complex formation, given the approximations of the model calculations. The simulation structure revealed significant complex stabilization from the hydrophobic effect, while the electrostatic cost of releasing water molecules from the interface determined to be highly unfavorable. By partitioning the total electrostatic and hydrophobic terms into residue free-energy contributions, a binding-affinity 'signature' for synaptobrevin was developed from the optimized conformation. The results demonstrate the effect of substrate length on complex formation and identify a peripheral high-affinity binding site near the N-terminal region that might initiate cooperative activation responsible for the large minimal substrate length requirement. The so-called SNARE motif is observed to contribute negligible free energy of binding.

Binding Sites↗