Biomedical subjects

Donald Bashford

Publications and source records attributed to Donald Bashford.

14 recordsLinked to original sources

PubMed2006

Model for proton transport coupled to protein conformational change: application to proton pumping in the bacteriorhodopsin photocycle.

A modeling method is presented for protein systems in which proton transport is coupled to conformational change, as in proton pumps and in motors driven by the proton-motive force. Previously developed methods for calculating pKa values in proteins using a macroscopic dielectric model are extended beyond the equilibrium case to a master-equation model for the time evolution of the system through states defined by ionization microstate and a discrete set of conformers. The macroscopic dielectric model supplies free energy changes for changes of protonation microstate, while the method for obtaining the energetics of conformational change and the relaxation rates, the other ingredients needed for the master equation, are system dependent. The method is applied to the photoactivated proton pump, bacteriorhodopsin, using conformational free energy differences from experiment and treating relaxation rates through three adjustable parameters. The model is found to pump protons with an efficiency relatively insensitive to parameter choice over a wide range of parameter values, and most of the main features of the known photocycle from very early M to the return to the resting state are reproduced. The boundaries of these parameter ranges are such that short-range proton transfers are faster than longer-range ones, which in turn are faster than conformational changes. No relaxation rates depend on conformation. The results suggest that an "accessibility switch", while not ruled out, is not required and that vectorial proton transport can be achieved through the coupling of the energetics of ionization and conformational states.

Bacteriorhodopsins↗

PubMed2005

Disordered p27Kip1 exhibits intrinsic structure resembling the Cdk2/cyclin A-bound conformation.

p27Kip1 (p27) influences cell division by regulating nuclear cyclin-dependent kinases. Before binding, p27 is at least partially disordered and folds upon binding its Cdk/cyclin targets. 30-40% of human proteins, including p27, are predicted to contain disordered segments, and have been termed intrinsically unstructured proteins (IUPs). Unfortunately, the inherent dynamics of IUPs hamper detailed analysis of their structure/function relationships. Here, we describe the use of molecular dynamics (MD) computations and solution NMR spectroscopy to reveal that several segments of the p27 kinase inhibitory domain (p27-KID), in addition to the previously characterized helical segment, exist as highly populated, intrinsically folded structural units (IFSUs). Several IFSUs resemble structural features of bound p27-KID, while another exhibits alternative conformations. Interestingly, the highly conserved, specificity determining segment of p27 is shown to be highly disordered. Elucidation of IFSUs within p27-KID allows consideration of their influences on the thermodynamics and kinetics of Cdk/cyclin binding. The degree to which IFSUs are populated within p27-KID is surprising and suggests that other putative IUPs contain IFSUs that may be studied using similar techniques.

Computer Simulation↗

PubMed2004

On the role of the conserved aspartate in the hydrolysis of the phosphocysteine intermediate of the low molecular weight tyrosine phosphatase.

The usual rate-determining step in the catalytic mechanism of the low molecular weight tyrosine phosphatases involves the hydrolysis of a phosphocysteine intermediate. To explain this hydrolysis, general base-catalyzed attack of water by the anion of a conserved aspartic acid has sometimes been invoked. However, experimental measurements of solvent deuterium kinetic isotope effects for this enzyme do not reveal a rate-limiting proton transfer accompanying dephosphorylation. Moreover, base activation of water is difficult to reconcile with the known gas-phase proton affinities and solution phase pK(a)'s of aspartic acid and water. Alternatively, hydrolysis could proceed by a direct nucleophilic attack by a water molecule. To understand the hydrolysis mechanism, we have used high-level density functional methods of quantum chemistry combined with continuum electrostatics models of the protein and the solvent. Our calculations do not support a catalytic activation of water by the aspartate. Instead, they indicate that the water oxygen directly attacks the phosphorus, with the aspartate residue acting as a H-bond acceptor. In the transition state, the water protons are still bound to the oxygen. Beyond the transition state, the barrier to proton transfer to the base is greatly diminished; the aspartate can abstract a proton only after the transition state, a result consistent with experimental solvent isotope effects for this enzyme and with established precedents for phosphomonoester hydrolysis.

Animals↗

PubMed2004

Exploring protein native states and large-scale conformational changes with a modified generalized born model.

Implicit solvation models provide, for many applications, a reasonably accurate and computationally effective way to describe the electrostatics of aqueous solvation. Here, a popular analytical Generalized Born (GB) solvation model is modified to improve its accuracy in calculating the solvent polarization part of free energy changes in large-scale conformational transitions, such as protein folding. In contrast to an earlier GB model (implemented in the AMBER-6 program), the improved version does not overstabilize the native structures relative to the finite-difference Poisson-Boltzmann continuum treatment. In addition to improving the energy balance between folded and unfolded conformers, the algorithm (available in the AMBER-7 and NAB molecular modeling packages) is shown to perform well in more than 50 ns of native-state molecular dynamics (MD) simulations of thioredoxin, protein-A, and ubiquitin, as well as in a simulation of Barnase/Barstar complex formation. For thioredoxin, various combinations of input parameters have been explored, such as the underlying gas-phase force fields and the atomic radii. The best performance is achieved with a previously proposed modification to the torsional potential in the Amber ff99 force field, which yields stable native trajectories for all of the tested proteins, with backbone root-mean-square deviations from the native structures being approximately 1.5 A after 6 ns of simulation time. The structure of Barnase/Barstar complex is regenerated, starting from an unbound state, to within 1.9 A relative to the crystal structure of the complex.

Bacterial Proteins↗

PubMed2004

Macroscopic electrostatic models for protonation states in proteins.

The use of macroscopic electrostatic models to calculate the relative energetics of protonation states and the pH-titration properties of ionizable groups in proteins is described. These methods treat the protein as an irregularly-shaped low-dielectric object containing embedded atomic charges immersed in a high-dielectric (solvent) medium. The energetics of altering protonation states then involves the electrostatic work of altering the embedded atomic charges. The governing electrostatic equation is either the Poisson or linearized Poisson-Boltzmann equation, which generally requires numerical solution. A tutorial approach is taken, the main aim of which is a thorough understanding of the method.

Models, Chemical↗

PubMed2003

A theoretical study of the UV/visible absorption and emission solvatochromic properties of solvent-sensitive dyes.

Using the density-functional vertical self-consistent reaction field (VSCRF) solvation model, incorporated with the conductor-like screening model (COSMO) and the self-consistent reaction field (SCRF) methods, we have studied the solvatochromic shifts of both the absorption and emission bands of four solvent-sensitive dyes in different solutions. The dye molecules studied here are: S-TBA merocyanine, Abdel-Halim's merocyanine, the rigidified amino-coumarin C153, and Nile red. These dyes were selected because they exemplify different structural features likely to impact the solvent-sensitive fluorescence of "push-pull", or merocyanine, fluorophores. All trends of the blue or red shifts were correctly predicted, comparing with the experimental observations. Explict H-bonding interactions were also considered in several protic solutions like H2O, methanol and ethanol, showing that including explicit H-bonding solvent molecule(s) in the calculations is important to obtain the correct order of the excitation and emission energies. The geometries, electronic structures, dipole moments, and intra- and intermolecular charge transfers of the dyes in different solvents are also discussed.

Absorption↗

PubMed2003

Proton affinity changes driving unidirectional proton transport in the bacteriorhodopsin photocycle.

Bacteriorhodopsin is the smallest autonomous light-driven proton pump. Proposals as to how it achieves the directionality of its trans-membrane proton transport fall into two categories: accessibility-switch models in which proton transfer pathways in different parts of the molecule are opened and closed during the photocycle, and affinity-switch models, which focus on changes in proton affinity of groups along the transport chain during the photocycle. Using newly available structural data, and adapting current methods of protein protonation-state prediction to the non-equilibrium case, we have calculated the relative free energies of protonation microstates of groups on the transport chain during key conformational states of the photocycle. Proton flow is modeled using accessibility limitations that do not change during the photocycle. The results show that changes in affinity (microstate energy) calculable from the structural models are sufficient to drive unidirectional proton transport without invoking an accessibility switch. Modeling studies for the N state relative to late M suggest that small structural re-arrangements in the cytoplasmic side may be enough to produce the crucial affinity change of Asp96 during N that allows it to participate in the reprotonation of the Schiff base from the cytoplasmic side. Methodologically, the work represents a conceptual advance compared to the usual calculations of pK(a) using macroscopic electrostatic models. We operate with collective states of protonation involving all key groups, rather than the individual-group pK(a) values traditionally used. When combined with state-to-state transition rules based on accessibility considerations, a model for non-equilibrium proton flow is obtained. Such methods should also be applicable to other active proton-transport systems.

Aspartic Acid↗

PubMed2003

Photoisomerization and proton transfer in photoactive yellow protein.

The photoactive yellow protein (PYP) is a bacterial photosensor containing a para-coumaryl thioester chromophore that absorbs blue light, initiating a photocycle involving a series of conformational changes. Here, we present computational studies to resolve uncertainties and controversies concerning the correspondence between atomic structures and spectroscopic measurements on early photocycle intermediates. The initial nanoseconds of the PYP photocycle are examined using time-dependent density functional theory (TDDFT) to calculate the energy profiles for chromophore photoisomerization and proton transfer, and to calculate excitation energies to identify photocycle intermediates. The calculated potential energy surface for photoisomerization matches key, experimentally determined, spectral parameters. The calculated excitation energy of the photocycle intermediate cryogenically trapped in a crystal structure by Genick et al. [Genick, U. K.; Soltis, S. M.; Kuhn, P.; Canestrelli, I. L.; Getzoff, E. D. Nature 1998, 392, 206-209] supports its assignment to the PYP(B) (I(0)) intermediate. Differences between the time-resolved room temperature (298 K) spectrum of the PYP(B) intermediate and its low temperature (77 K) absorbance are attributed to a predominantly deprotonated chromophore in the former and protonated chromophore in the latter. This contrasts with the widely held belief that chromophore protonation does not occur until after the PYP(L) (I(1) or pR) intermediate. The structure of the chromophore in the PYP(L) intermediate is determined computationally and shown to be deprotonated, in agreement with experiment. Calculations based on our PYP(B) and PYP(L) models lead to insights concerning the PYP(BL) intermediate, observed only at low temperature. The results suggest that the proton is more mobile between Glu46 and the chromophore than previously realized. The findings presented here provide an example of the insights that theoretical studies can contribute to a unified analysis of experimental structures and spectra.

Bacterial Proteins↗

PubMed2003

Structural details, pathways, and energetics of unfolding apomyoglobin.

Protein folding is often difficult to characterize experimentally because of the transience of intermediate states, and the complexity of the protein-solvent system. Atomistic simulations, which could provide more detailed information, have had to employ highly simplified models or high temperatures, to cope with the long time scales of unfolding; direct simulation of folding is even more problematic. We report a fully atomistic simulation of the acid-induced unfolding of apomyoglobin in which the protonation of acidic side-chains to simulate low pH is sufficient to induce unfolding at room temperature with no added biasing forces or other unusual conditions; and the trajectory is validated by comparison to experimental characterization of intermediate states. Novel insights provided by their analysis include: characterization of a dry swollen globule state forming a barrier to initial unfolding or final folding; observation of cooperativity in secondary and tertiary structure formation and its explanation in terms of dielectric environments; and structural details of the intermediate and the completely unfolded states. These insights involve time scales and levels of structural detail that are presently beyond the range of experiment, but come within reach through the simulation methods described here. An implicit solvation model is used to analyze the energetics of protein folding at various pH and ionic strength values, and a reasonable estimate of folding free energy is obtained. Electrostatic interactions are found to disfavor folding.

Apoproteins↗

PubMed2002

Effective Born radii in the generalized Born approximation: the importance of being perfect.

Generalized Born (GB) models provide, for many applications, an accurate and computationally facile estimate of the electrostatic contribution to aqueous solvation. The GB models involve two main types of approximations relative to the Poisson equation (PE) theory on which they are based. First, the self-energy contributions of individual atoms are estimated and expressed as "effective Born radii." Next, the atom-pair contributions are estimated by an analytical function f(GB) that depends upon the effective Born radii and interatomic distance of the atom pairs. Here, the relative impacts of these approximations are investigated by calculating "perfect" effective Born radii from PE theory, and enquiring as to how well the atom-pairwise energy terms from a GB model using these perfect radii in the standard f(GB) function duplicate the equivalent terms from PE theory. In tests on several biological macromolecules, the use of these perfect radii greatly increases the accuracy of the atom-pair terms; that is, the standard form of f(GB) performs quite well. The remaining small error has a systematic and a random component. The latter cannot be removed without significantly increasing the complexity of the GB model, but an alternative choice of f(GB) can reduce the systematic part. A molecular dynamics simulation using a perfect-radii GB model compares favorably with simulations using conventional GB, even though the radii remain fixed in the former. These results quantify, for the GB field, the importance of getting the effective Born radii right; indeed, with perfect radii, the GB model gives a very good approximation to the underlying PE theory for a variety of biomacromolecular types and conformations.

Macromolecular Substances↗

PubMed2002

Density functional study of the mechanism of a tyrosine phosphatase: I. Intermediate formation.

The first step in the catalytic mechanism of a protein tyrosine phosphatase, the transfer of a phosphate group from the phosphotyrosine substrate to a cysteine side chain of the protein to form a phosphoenzyme intermediate, has been studied by combining density functional calculations of an active-site cluster with continuum electrostatic descriptions of the solvent and the remainder of the protein. This approach provides the high level of quantum chemical methodology needed to adequately model phosphotransfer reactions with a reasonable description of the environment around the active site. Energy barriers and geometries along a reaction pathway are calculated. In the literature, mechanisms assuming both a monoanionic and a dianionic substrate have been proposed; this disagreement is addressed by performing calculations for both possibilities. For the dianionic substrate, a dissociative reaction pathway with early proton transfer to the leaving group and a 9 kcal/mol energy barrier is predicted (the experimental estimate is ca. 14 kcal/mol), while for the monoanionic substrate, an associative pathway with late proton transfer and a 22 kcal/mol energy barrier is predicted. These results, together with a review of experimental evidence, support the dianionic-substrate/dissociative-pathway alternative. The relationship between a dianionic or monoanionic substrate and a dissociative or associative pathway, respectively, can be understood in terms of classical organic chemical reaction pathways.

Catalysis↗

PubMed2002

Efficient aminoacylation of the tRNA(Ala) acceptor stem: dependence on the 2:71 base pair.

Specific aminoacylation by aminoacyl-tRNA synthetases requires accurate recognition of cognate tRNA substrates. In the case of alanyl-tRNA synthetase (AlaRS), RNA duplexes that mimic the acceptor stem of the tRNA are efficient substrates for aminoacylation in vitro. It was previously shown that recognition by AlaRS is severely affected by a simple base pair transversion of the G2:C71 pair at the second position in the RNA helix. In this study, we determined the aminoacylation efficiencies of 50 variants of the tRNA(Ala) acceptor stem containing substitutions at the 2:71 position. We find that there is not a single functional group of the wild-type G2:C71 base pair that is critical for positive recognition. Rather, we observed that base-pair orientation plays an important role in recognition. In particular, pyrimidine2:purine71 combinations generally resulted in decreased aminoacylation efficiency compared to the corresponding purine:pyrimidine pair. Moreover, the activity of a pyrimidine:purine variant could be partially restored by the presence of a major groove amino group at position 71. In an attempt to understand this result further, dielectric continuum electrostatic calculations were carried out, in some cases with additional inclusion of van der Waals interaction energies, to determine interaction potentials of the wild-type duplexAla and seven 2:71 variants. This analysis revealed a positive correlation between major groove negative electrostatic potential in the vicinity of the 3:70 base pair and measured aminoacylation efficiency.

Acylation↗

PubMed1999

CuZn Superoxide Dismutase Geometry Optimization, Energetics, and Redox Potential Calculations by Density Functional and Electrostatic Methods.

The structures, energetics, and orbital- and charge-dependent properties of copper zinc superoxide dismutase (CuZnSOD) have been studied using density functional and electrostatic methods. The CuZnSOD was represented with a model consisting of copper and zinc sites connected by a bridging histidine ligand. In addition to the bridge, three histidine ligands and one water molecule were bonded to the Cu ion in the copper site as first-shell ligands. Two histidine ligands and an aspartate were coordinated to the zinc ion in the zinc site. Full optimization of the model was performed using different functionals, both local and nonlocal. Geometrical parameters calculated with the nonlocal functionals agree well with the experimental X-ray data. In our calculated results, the His61 Nepsilon-Cu bond in the active site breaks during the reduction and protonation, consistent with a number of X-ray structures and with EXAFS and NMR evidence. The reduction potential and pK(a) of the coupled electron/proton reaction catalyzed by CuZnSOD were determined using different models for the extended environment-from an electrostatic representation of continuum solvent, to the full protein/solvent environment using a Poisson-Boltzmann method. The predicted redox potential and pK(a) values determined using the model with the full protein/solvent environment are in excellent agreement with experiment. Inclusion of the full protein environment is essential for an accurate description of the redox process. Although the zinc ion does not play a direct redox role in the dismutation, its electronic contribution is very important for the catalytic mechanism.

Journal Article↗

PubMed1999

Density Functional and Electrostatic Calculations of Manganese Superoxide Dismutase Active Site Complexes in Protein Environments.

Density functional and electrostatic methods have been applied to calculate active site geometries and the redox potential of manganese superoxide dismutase (MnSOD). The initial active site clusters were built up by including only first-shell side chain ligands and then augmented by second-shell ligands. The density functional optimized Mn-ligand bond lengths for the reduced complexes in general compared fairly well with protein crystallography data; however, large deviations for calculated Mn-OH distances were found for the oxidized active site clusters. Our calculations suggest that this deviation can be attributed to the redox heterogeneity of the oxidized protein in X-ray crystallography studies. The redox potential was calculated by treating the protein environment and the solvent bulk by a semimacroscopic electrostatic model. The protein structures were taken from the Thermus thermophilus enzyme. The calculated coupled redox potentials converge toward experimental values with increasing size of the active site cluster models, and the final calculated value was +0.06 V, compared to experimental values of +0.26 V determined for Bacillus stearothermophilus and +0.31 V in Escherichia coli enzymes. Using an energy decomposition scheme, the effects of the second-shell ligands and the protein and reaction fields have been analyzed.

Journal Article↗