Search PubMed⌕ Search

Biomedical subjects

Jeremy C Smith

Publications and source records attributed to Jeremy C Smith.

At least 55 records · Page 3Linked to original sources
PubMed2003

Fluorescence quenching of dyes by tryptophan: interactions at atomic detail from combination of experiment and computer simulation.

Fluorescence spectroscopy and molecular dynamics (MD) simulation are combined to characterize the interaction of two organic fluorescent dyes, rhodamine 6G (R6G) and an oxazine derivative (MR121), with the amino acid tryptophan in aqueous solution. Steady-state and time-resolved fluorescence quenching experiments reveal the formation of essentially nonfluorescent ground-state dye/Trp complexes. The MD simulations are used to elucidate the molecular interaction geometries involved. The MD-derived probability distribution of the distance r between the centers of geometry of the dye and quencher ring systems, P(r), extends to higher distances for R6G than for MR121 due to population in the R6G/Trp system of fluorescent interaction geometries between Trp and the phenyl ring and ester group of the dye. The consequence of this is the experimental finding that under the conditions used in the simulations about 25% of the R6G dye is fluorescent in comparison with 10% of the MR121. Combining the above findings allows determination of the "quenching distance", r, above which no quenching occurs. r is found to be very similar (approximately 5.5 A) for both dye/Trp systems, corresponding to close to van der Waals contact. Both experimental dynamic Stern-Volmer analysis and the MD trajectories demonstrate that the main determinant of the fluorescence intensity is static quenching. The approach presented is likely to be useful in the structural interpretation of data obtained from fluorescent conjugates commonly used for monitoring the binding and dynamics of biomolecular systems.

Computer Simulation↗
PubMed2003

Principal components of the protein dynamical transition.

Proteins exhibit a solvent-driven dynamical transition at 180-220 K, manifested by nonlinearity in the temperature dependence of the average mean-square displacement. Here, molecular dynamics simulations of hydrated myoglobin show that the onset of the transition at approximately 180 K is characterized by the appearance of a single double-well principal component mode involving a global motion of two groups of helices. As the temperature is raised a few more quasiharmonic and multiminimum components successively appear. The results indicate an underlying simplicity in the protein dynamical transition.

Computer Simulation↗
PubMed2003

How well does charge reparametrisation account for solvent screening in molecular mechanics calculations? The example of myosin.

The rate constant of an enzyme-catalysed reaction is one of the major target properties to understand protein function. Atomic-detail computer simulations can in principle be used to estimate rate constants from the energy profile along the reaction coordinate. For such simulations, molecular mechanics is combined with a quantum description of the reaction process. In molecular mechanics calculations, the electrostatic field is represented by the Coulomb potential of partial atomic charges which have been parametrised for small building blocks in vacuum and transferred to the macromolecule. In aqueous solution, however, the electrostatic interactions are affected by the solvent polarization. While this can be described by numerically solving the Poisson-Boltzmann equation, it is computationally expensive. A simple approximation to this is to optimally reproduce the electrostatic potential in solution by reparametrising the partial atomic charges in such a way that a simple Coulomb potential can still be used. Such a procedure would allow to perform fast calculations of reaction processes in proteins while accounting for the solvent screening effect. Here, this method is tested on myosin, a motor protein that is both an enzyme and exists in very different conformations.

Adenosine Triphosphate↗
PubMed2003

Molecular mechanics force field parameterization of the fluorescent probe rhodamine 6G using automated frequency matching.

Novel single-molecule fluorescence experimental techniques have prompted a growing need to develop refined computational models of dye-tagged biomolecules. As a necessary first step towards useful molecular simulations of fluorescence-labeled biomolecules, we have derived a force field for the commonly used dye, rhodamine 6G (R6G). A novel automated method is used that includes fitting the molecular mechanics potential to both vibrational frequencies and eigenvector projections derived from quantum chemical calculations. The method is benchmarked on a series of aromatic molecules then applied to derive new parameters for R6G. The force field derived reproduces well the crystal structure of R6G.

Algorithms↗
PubMed2003

Energy resolution and dynamical heterogeneity effects on elastic incoherent neutron scattering from molecular systems.

Incoherent neutron scattering is widely used to probe picosecond-nanosecond time scale dynamics of molecular systems. In systems of spatially confined atoms the relatively high intensity of elastic incoherent neutron scattering is often used to obtain a first estimate of the dynamics present. For many complex systems, however, experimental elastic scattering is difficult to interpret unambiguously using analytical dynamical models that go beyond the determination of an average mean-square displacement. To circumvent this problem a description of the scattering is derived here that encompasses a variety of analytical models in a common framework. The framework describes the time-converged part of the dynamic structure factor [the elastic incoherent scattering function (EISF)] and lends itself to practical use by explicitly incorporating effects due to the finite energy resolution of the instrument used. The dependence of the elastic scattering on wave vector is examined, and it is shown how heterogeneity in the distribution of mean-square displacements can be related to deviations of the scattering from Gaussian behavior. In this case, a correction to fourth order in the scattering vector can be used to extract the variance of the distribution of mean-square displacements. The formalism is used in a discussion of measurements on dynamics accompanying the glass transition in molecular systems. By fitting to experimental data obtained on a protein solution the present methodology is used to show how the existence of a temperature-dependent relaxation frequency can lead to a transition in the measured mean-square displacement in the absence of an EISF change.

Journal Article↗
PubMed2003

Molecular dynamics simulation reveals a surface salt bridge forming a kinetic trap in unfolding of truncated Staphylococcal nuclease.

Surface salt bridges are ubiquitous in globular proteins. Their contribution to protein stability has been extensively debated in the past decade. Here, molecular dynamics simulations are performed starting from a non-equilibrium state of Staphylococcal nuclease (SNase) with C-terminal truncation (SNaseDelta). The results indicate a key role in the unfolding of the surface salt bridge between arginine 105 and glutamate 135. Experimentally, SNaseDelta is known to be partially unfolded. However, in simulations over 1 ns at 300 K and over 500 ps at 400 K, SNaseDelta remains stable in the native-like folded conformation, the salt bridge hindering unfolding. When the potential function is altered so as to selectively weaken the salt bridge, which then breaks rapidly at 430 K, the protein starts to unfold. The results suggest that breaking of this salt bridge presents a significant barrier to the unfolding transition of SNaseDelta from a native-like state to the unfolded state. Potential of mean force calculations indicate that the barrier height for this transition is approximately 7 kcal/mol.

Arginine↗
PubMed2003

Convergence in peptide folding simulation: multiple trajectories of a potential AIDS pharmacophore.

To examine the conformational properties in aqueous solution of a 15-residue peptide that is a potential pharmacophore for AIDS vaccine development, molecular dynamics simulations were performed in water starting from structures determined experimentally in three different organic solvents. Convergence characteristics of the simulation are examined in Cartesian and conformational spaces. In addition, novel analysis tools are employed including a multidimensional scaling method to represent the distance between trajectory frames. As these methods are based on a variety of physical parameters, they provide a useful cross-check on the structural convergence. Theoretical two-dimensional (2D) 1H-NMR spectra are also generated. These are superficially quite different in appearance, demonstrating that backbone similarities difficult to identify by visual inspection of 2D NMR data can be revealed using the methods described here.

AIDS Vaccines↗
PubMed2003

Molecular dynamics decomposition of temperature-dependent elastic neutron scattering by a protein solution.

Molecular dynamics simulations are performed of bovine pancreatic trypsin inhibitor in a cryosolution over a range of temperatures from 80 to 300 K and the origins identified of elastic dynamic neutron scattering from the solution. The elastic scattering and mean-square displacement calculated from the molecular dynamics trajectories are in reasonable agreement with experiments on a larger protein in the same solvent. The solvent and protein contributions to the scattering from the simulation model are determined. At lower temperatures (< approximately 200 K) or on shorter timescales ( approximately 10 ps) the scattering contributions are proportional to the isotopic nuclear scattering cross-sections of each component. However, for T > 200 K marked deviations from these cross-sections are seen due to differences in the dynamics of the components of the solution. Rapid activation of solvent diffusion leads to the variation with temperature of the total elastic intensity being determined largely by that of the solvent. At higher temperatures (>240 K) and longer times ( approximately 100 ps) the protein makes the only significant contribution to the scattering, the solvent scattering having moved out of the accessible time-space window. Decomposition of the protein mean-square displacement shows that the observed dynamical transition in the solution at 200-220 K involves activation of both internal motions and external whole-molecule rotational and translational diffusion. The proportion that the external dynamics contributes to the protein mean-square displacement increases to approximately 30 and 60% at 300 K on the 10- and 100-ps timescales, respectively.

Aprotinin↗
PubMed2003

Translational hydration water dynamics drives the protein glass transition.

Experimental and computer simulation studies have revealed the presence of a glass-like transition in the internal dynamics of hydrated proteins at approximately 200 K involving an increase of the amplitude of anharmonic dynamics. This increase in flexibility has been correlated with the onset of protein activity. Here, we determine the driving force behind the protein transition by performing molecular dynamics simulations of myoglobin surrounded by a shell of water. A dual heat bath method is used with which, in any given simulation, the protein and solvent are held at different temperatures, and sets of simulations are performed varying the temperature of the components. The results show that the protein transition is driven by a dynamical transition in the hydration water that induces increased fluctuations primarily in side chains in the external regions of the protein. The water transition involves activation of translational diffusion and occurs even in simulations where the protein atoms are held fixed.

Animals↗
PubMed2003

The position of QB in the photosynthetic reaction center depends on pH: a theoretical analysis of the proton uptake upon QB reduction.

Electrostatics-based calculations have been performed to examine the proton uptake upon reduction of the terminal electron acceptor Q(B) in the photosynthetic reaction center of Rhodobacter sphaeroides as a function of pH and the associated conformational equilibrium. Two crystal structures of the reaction center were considered: one structure was determined in the dark and the other under illumination. In the two structures, the Q(B) was found in two different positions, proximal or distal to the nonheme iron. Because Q(B) was found mainly in the distal position in the dark and only in the proximal position under illumination, the two positions have been attributed mostly to the oxidized and the reduced forms of Q(B), respectively. We calculated the proton uptake upon Q(B) reduction by four different models. In the first model, Q(B) is allowed to equilibrate between the two positions with either oxidation state. This equilibrium was allowed to vary with pH. In the other three models the distribution of Q(B) between the proximal position and the distal position was pH-independent, with Q(B) occupying only the distal position or only the proximal position or populating the two positions with a fixed ratio. Only the first model, which includes the pH-dependent conformational equilibrium, reproduces both the experimentally measured pH dependence of the proton uptake and the crystallographically observed conformational equilibrium at pH 8. From this model, we find that Q(B) occupies only the distal position below pH 6.5 and only the proximal position above pH 9.0 in both oxidation states. Between these pH values both positions are partially occupied. The reduced Q(B) has a higher occupancy in the proximal position than the oxidized Q(B). In summary, the present results indicate that the conformational equilibrium of Q(B) depends not only on the redox state of Q(B), but also on the pH value of the solution.

Computer Simulation↗
PubMed2003

The dynamic transition in proteins may have a simple explanation.

The transition that has been observed in the dynamics of hydrated proteins at low temperatures (180-230 K) is normally interpreted as a change from vibrational, harmonic motion at low temperatures to anharmonic motions as the temperature is raised. It is taken to be an intrinsic property of proteins and has been associated with the onset of protein functions. Examination of the dynamic behaviour of proteins in solution within a defined timescale window suggests that certain observations can be explained without the need to invoke a discontinuity in the dynamics of proteins with temperature, i.e. the existence of a dynamical transition is not required. This is discussed in the context of recent evidence that enzyme activity is independent of the activation of anharmonic picosecond dynamics and declines steadily with temperature through the apparent dynamic transition, in accordance with the Arrhenius relationship. That similar timescale dependent dynamical behaviour has been observed experimentally in chain polymers, and seen also in computer simulations of silica glasses, suggests that the phenomenon may be of wide general relevance in both simple glassy and more complex polymeric systems.

Enzymes↗
PubMed2003

Time-resolved computational protein biochemistry: solvent effects on interactions, conformational transitions and equilibrium fluctuations.

Solvent plays an important role in modulating internal motions of proteins. Here we present a computational method for including solvent effects on charge-charge interactions and on pathways between functional protein conformations, and examine solvent effects on equilibrium internal fluctuations in proteins. A computationally efficient charge reparametrisation method is presented that satisfactorily reproduces the electrostatic interactions present in a full continuum Poisson-Boltzmann representation. The application of charge reparametrisation in the calculation of a large-scale conformational transition pathway in a protein, annexin V, is illustrated. We also examine solvent effects on fast (picosecond timescale) internal protein dynamics. Nosé-Hoover dual heatbath molecular dynamics simulations are performed. These simulations allow the solvent region to be fixed at one temperature and the protein at another. The results of the Nosé-Hoover simulations on hydrated myoglobin confirm that the solvent temperature strongly influences the protein fluctuations. We consider to what extent the solvent can be considered to determine the high temperature protein dynamics.

Animals↗
PubMed2003

Computational tools for analysing structural changes in proteins in solution.

Many important structural changes in proteins involve long-time dynamics, which are outside the timescale presently accessible by a straightforward integration of Newton's equations of motion. This problem is addressed with minimisation-based algorithms, which are applied on possible reaction pathways using atomic-detail models. For reasons of efficiency, an implicit treatment of solvent is imperative. We present the charge reparameterisation protocol, which is a method that approximates the interaction energies obtained by a numerical solution of the Poisson-Boltzmann equation. Furthermore, we present a number of methods that can be used to compute possible reaction pathways associated with a particular conformational change. Two of them, the self-penalty walk and the nudged elastic band method, define an objective function, which is minimised to find optimal paths. A third method, conjugate peak refinement, is a heuristic method, which finds minimum energy paths without the use of an explicit objective function. Finally, we discuss problems and limitations with these methods and give a perspective on future research.

Algorithms↗
PubMed2002

Is the first hydration shell of lysozyme of higher density than bulk water?

Characterization of the physical properties of protein surface hydration water is critical for understanding protein structure and folding. Here, using molecular dynamics simulation, we provide an explanation of recent x-ray and neutron solution scattering data that indicate that the density of water on the surface of lysozyme is significantly higher than that of bulk water. The simulation-derived scattering profiles are in excellent agreement with the experiment. In the simulation, the 3-A-thick first hydration layer is 15% denser than bulk water. About two-thirds of this increase is the result of a geometric contribution that would also be present if the water was unperturbed from the bulk. The remaining third arises from modification of the water structure and dynamics, involving approximately equal contributions from shortening of the average water-water O-O distance and an increase in the coordination number. Variation in the first hydration shell density is shown to be determined by topographical and electrostatic properties of the protein surface. On average, denser water is found in depressions on the surface in which the water dipoles tend to be aligned parallel to each other by the electrostatic field generated by the protein atoms.

Biophysical Phenomena↗
PubMed2002

Solution structure of a putative HIV1 immunogenic peptide: computer simulation of the principal CD4 binding domain of gp120.

A 15-residue segment within the principal CD4 binding domain of gp120 from HIV1 is a natural subject for vaccine development. Knowledge of a three-dimensional structure of this peptide in aqueous solution can in principle be used for pharmacophoric footprint identification. However, the peptide is resistant to structural characterization using nuclear magnetic resonance (NMR) due to aggregation properties. Here we examine the conformational properties of the peptide using molecular dynamics simulation. Three simulations were performed starting from widely different structures of the peptide determined by NMR in organic solvents. The three structures converge during the simulations to a common structure in both Cartesian and in backbone dihedral conformational space. Preliminary examination of the average three-dimensional structure identifies, after sequence threading, characteristics of potential use as a pharmacophoric footprint.

Binding Sites↗
PubMed2002

SASSIM: a method for calculating small-angle X-ray and neutron scattering and the associated molecular envelope from explicit-atom models of solvated proteins.

A method is presented to calculate efficiently small-angle neutron and X-ray solution scattering intensities from explicit-atom models of macromolecules and the surrounding solvent. The method is based on a multipole expansion of the scattering amplitude. It is particularly appropriate for extensive configurational averaging, as is required for calculations based on computer-simulation results. In test calculations, excellent agreement with experiment is found between neutron and X-ray scattering profiles calculated from a molecular-dynamics simulation of lysozyme in water. The question of definition of the protein surface is also addressed. For comparison with the continuum model, an analytical envelope around the protein is defined in terms of spherical harmonics and is calculated using a Lebedev grid. The analytical surface thus defined is shown to reproduce well the scattering profile calculated from the explicit-atom model of the protein.

Animals↗
PubMed2002

Can the calculation of ligand binding free energies be improved with continuum solvent electrostatics and an ideal-gas entropy correction?

The prediction of a ligand binding constant requires generating three-dimensional structures of the complex concerned and reliably scoring these structures. Here, the scoring problem is investigated by examining benzamidine-like inhibitors of trypsin, a system for which errors in the structures are small. Precise and consistent binding free energies for the inhibitors are determined experimentally for this test system. To examine possible improvement of scoring methods, we test the suitability of continuum electrostatics to account for solvation effects and use an ideal-gas entropy correction to account for the changes in the degrees of freedom of the ligand. The small observed root-mean-square deviation of 0.55 kcal/mol of the calculated relative to the experimental values indicates that the essentials of the binding process have been captured. Even though all six ligands make the same salt bridge and H-bonds to the protein, the electrostatic contribution varies among the ligands by as much as 2 kcal/mol. Moreover, although the ligands are rigid and similar in size, the entropic terms also significantly affect the relative binding affinities (by up to 2.7 kcal/mol). The present approach to solvation and entropy may allow the ranking of the ligands to be considerably improved at a cost that makes the method applicable to the optimization of lead compounds or to the screening of small collections of ligands.

Journal Article↗