Search PubMed⌕ Search

Biomedical subjects

B M Pettitt

Publications and source records attributed to B M Pettitt.

At least 19 recordsLinked to original sources
PubMed2001

Simulations of the bis-penicillamine enkephalin in sodium chloride solution: a parameter study.

A simulation study of DPDPE in sodium chloride solution has been performed and compared with previous simulations using a different interaction potential for the ions. Both global thermodynamics as well as a characterization of association to DPDPE have been calculated. We show that the parameters used for the ions have a profound effect on the association to the peptide in 1M NaCl. The observed differences suggest that individual associations in these and previous simulations are sensitive to parameters.

Computer Simulation↗
PubMed2001

Structural basis for the activity of pp60(c-src) protein tyrosine kinase inhibitors.

Conformational searches on three closely related pp60(c-src) protein tyrosine kinase inhibitors of varying potencies were performed to determine a structural basis for their activity. The first was a linear peptide (PDNEYAFFQf), the second its 10-membered cyclic analogue, and the third a cyclic analogue with a para carboxyphenylalanine in place of one the F residues. A common backbone conformation with an antiparallel beta-sheet-like geometry capped by similar beta-turns was found for all three peptides, which may be a binding conformation and gives a candidate pharmacophore for further testing. The interaction between some polar side chains and between some of the aromatic rings may be important for maintaining the correct conformation. The differences in potencies of these inhibitors may be attributed to certain thermodynamic and chemical reasons.

Amino Acid Sequence↗
PubMed2001

Conformations of an adenine bulge in a DNA octamer and its influence on DNA structure from molecular dynamics simulations.

Molecular dynamics simulations have been applied to the DNA octamer d(GCGCA-GAAC). d(GTTCGCGC), which has an adenine bulge at the center to determine the pathway for interconversion between the stacked and extended forms. These forms are known to be important in the molecular recognition of bulges. From a total of ~35 ns of simulation time with the most recent CHARMM27 force field a variety of distinct conformations and subconformations are found. Stacked and fully looped-out forms are in excellent agreement with experimental data from NMR and x-ray crystallography. Furthermore, in a number of conformations the bulge base associates with the minor groove to varying degrees. Transitions between many of the conformations are observed in the simulations and used to propose a complete transition pathway between the stacked and fully extended conformations. The effect on the surrounding DNA sequence is investigated and biological implications of the accessible conformational space and the suggested transition pathway are discussed, in particular for the interaction of the MS2 replicase operator RNA with its coat protein.

Adenine↗
PubMed2000

Residence times of water molecules in the hydration sites of myoglobin.

Hydration sites are high-density regions in the three-dimensional time-averaged solvent structure in molecular dynamics simulations and diffraction experiments. In a simulation of sperm whale myoglobin, we found 294 such high-density regions. Their positions appear to agree reasonably well with the distributions of waters of hydration found in 38 x-ray and 1 neutron high-resolution structures of this protein. The hydration sites are characterized by an average occupancy and a combination of residence time parameters designed to approximate a distribution of residence times. It appears that although the occupancy and residence times of the majority of sites are rather bulk-like, the residence time distribution is shifted toward the longer components, relative to bulk. The sites with particularly long residence times are located only in the cavities and clefts of the protein. This indicates that other factors, such as hydrogen bonds and hydrophobicity of underlying protein residues, play a lesser role in determining the residence times of the longest-lived sites.

Binding Sites↗
PubMed1999

Comparison of the potentials of mean force for alanine tetrapeptide between integral equation theory and simulation.

The dielectrically consistent reference interaction site model (DRISM) integral equation theory is applied to determine the potential of mean force (PMF) for an alanine tetramer. A stochastic dynamics simulation of the alanine tetramer using this PMF is then compared with an explicit water molecular dynamics simulation. In addition, comparison is also done with simulations using other solvent models like the extended reference interaction site model (XRISM) theory, constant dielectric and linear distance-dependent dielectric models. The results show that the DRISM method offers a fairly accurate and computationally inexpensive alternative to explicit water simulations for studies on small peptides.

Algorithms↗
PubMed1999

Modeling high-resolution hydration patterns in correlation with DNA sequence and conformation.

Hydration around the DNA fragment d(C5T5).(A5G5) is presented from two molecular dynamics simulations of 10 and 12 ns total simulation time. The DNA has been simulated as a flexible molecule with both the CHARMM and AMBER force fields in explicit solvent including counterions and 0.8 M additional NaCl salt. From the previous analysis of the DNA structure B-DNA conformations were found with the AMBER force-field and A-DNA conformations with CHARMM parameters. High-resolution hydration patterns are compared between the two conformations and between C.G and T.A base-pairs from the homopolymeric parts of the simulated sequence. Crystallographic results from a statistical analysis of hydration sites around DNA crystal structures compare very well with the simulation results. Differences between the crystal sites and our data are explained by variations in conformation, sequence, and limitations in the resolution of water sites by crystal diffraction. Hydration layers are defined from radial distribution functions and compared with experimental results. Excellent agreement is found when the measured experimental quantities are compared with the equivalent distribution of water molecules in the first hydration shell. The number of water molecules bound to DNA was found smaller around T.A base-pairs and around A-DNA as compared to B-DNA. This is partially offset by a larger number of water molecules in hydrophobic contact with DNA around T.A base-pairs and around A-DNA. The numbers of water molecules in minor and major grooves have been correlated with helical roll, twist, and inclination angles. The data more fully explain the observed B-->A transition at low humidity.

Computer Simulation↗
PubMed1999

Hydration effects on the electrostatic potential around tuftsin.

The electrostatic potential and component dielectric constants from molecular dynamics (MD) trajectories of tuftsin, a tetrapeptide with the amino acid sequence Thr-Lys-Pro-Arg in water and in saline solution are presented. The results obtained from the analysis of the MD trajectories for the total electrostatic potential at points on a grid using the Ewald technique are compared with the solution to the Poisson-Boltzmann (PB) equation. The latter was solved using several sets of dielectric constant parameters. The effects of structural averaging on the PB results were also considered. Solute conformational mobility in simulations gives rise to an electrostatic potential map around the solute dominated by the solute monopole (or lowest order multipole). The detailed spatial variation of the electrostatic potential on the molecular surface brought about by the compounded effects of the distribution of water and ions close to the peptide, solvent mobility, and solute conformational mobility are not qualitatively reproducible from a reparametrization of the input solute and solvent dielectric constants to the PB equation for a single structure or for structurally averaged PB calculations. Nevertheless, by fitting the PB to the MD electrostatic potential surfaces with the dielectric constants as fitting parameters, we found that the values that give the best fit are the values calculated from the MD trajectories. Implications of using such field calculations on the design of tuftsin peptide analogues are discussed.

Computer Simulation↗
PubMed1999

Structure and dynamics of alpha-MSH using DRISM integral equation theory and stochastic dynamics.

The structural and dynamical features of the hormone alpha-MSH in solution have been examined over a 100 ns time scale by using free energy molecular mechanics models at room temperature. The free energy surface has been modeled using methods from integral equation theory and the dynamics by the Langevin equation. A modification of the accessible surface area friction drag model was used to calculate the atomic friction coefficients. The molecule shows a stable beta-turn conformation in the message region and a close interaction between the side chains of His6, Phe7, and Trp9. A salt bridge between Glu5 and Arg8 was found not to be a preferred interaction, whereas a Glu5 and Lys11 salt bridge was not sampled, presumably due to relatively high free energy barriers. The message region was more conformationally rigid than the N-terminal region. Several structural features observed here agree well with experimental results. The conformational features suggest a receptor-hormone interaction model where the hydrophobic side chains of Phe7 and Trp9 interact with the transmembrane portion of the MC1 receptor. Also, the positively charged side chain of Arg8 and the imidazole side chain of His6 may interact with the negatively charged portions of the receptor which may even be on the receptor's extracellular loops.

Algorithms↗
PubMed1999

Sodium and chlorine ions as part of the DNA solvation shell.

The distribution of sodium and chlorine ions around DNA is presented from two molecular dynamics simulations of the DNA fragment d(C(5)T(5)). (A(5)G(5)) in explicit solvent with 0.8 M additional NaCl salt. One simulation was carried out for 10 ns with the CHARMM force field that keeps the DNA structure close to A-DNA, the other for 12 ns with the AMBER force field that preferentially stabilizes B-DNA conformations (, Biophys. J. 75:134-149). From radial distributions of sodium and chlorine ions a primary ion shell is defined. The ion counts and residence times of ions within this shell are compared between conformations and with experiment. Ordered sodium ion sites were found in minor and major grooves around both A and B-DNA conformations. Changes in the surrounding hydration structure are analyzed and implications for the stabilization of A-DNA and B-DNA conformations are discussed.

Base Pairing↗
PubMed1999

Modeling of alpha-MSH conformations with implicit solvent.

A conformational search for the most probable structures of the hormone alpha-MSH in aqueous solution was performed in order to help determine the structural features necessary for biological activity. The free-energy surface was modeled using methods from integral equation theory, and high-temperature molecular dynamics was used to enhance conformational sampling. Families of low free-energy structures have been found. The minimum energy structure shows a stable beta-turn conformation in the putative message region that is stabilized by a salt bridge between Glu5 and Lys11. The orientation of the side chains reflects the amphiphilic nature of the peptide, and a close interaction between the side chains of the His6, Phe7 and Trp9 was observed. Several structural features observed in the minimum energy structure agree well with experimental results. The conformational features led to a hypothesis of a receptor-hormone interaction model in which the hydrophobic side chains of Phe7 and Trp9 interact with the transmembrane portion of the human melanocortin (MC1) receptor. Also, the positively charged side chain of Arg8 and the imidazole side chain of His6 may interact with the negatively charged portions of the receptor which may even be on the receptor's extracellular loops.

Computer Simulation↗
PubMed1998

Characterizing global substates of myoglobin.

BACKGROUND: The massive amount of information generated from current molecular dynamics simulations makes the data difficult to analyze efficiently. Principal component analysis has been used for almost a century to detect and characterize data relationships and to reduce the dimensionality for problems in many fields. Here, we present an adaptation of principal component analysis using a partial singular value decomposition (SVD) for investigating both the localized and global motions of macromolecules. RESULTS: Configuration space projections from the SVD analysis of a variety of myoglobin simulations are used to characterize the dynamics of the protein. This technique reveals new dynamical motifs, which quantify proposed hierarchical structures of conformational substates for proteins and provide a means by which configuration space sampling efficiency may be probed. The SVD clearly shows that solvent effects facilitate transitions between global conformational substates for myoglobin molecular dynamics simulations. Lyapunov exponents calculated from the configuration space divergence of 15 trajectories agree with previous predictions for the chaotic behavior of complex protein systems. CONCLUSIONS: Configuration space projections provide invaluable information about protein motions that would be extremely difficult to obtain otherwise. While the configuration space for myoglobin is quite large, it does have structure. Our analysis of this structure shows that the protein hops between a number of distinct global conformational states, much like the local behavior observed for an individual residue.

Computer Simulation↗
PubMed1998

A molecular simulation picture of DNA hydration around A- and B-DNA.

Recent results from molecular dynamics (MD) simulations on hydration of DNA with respect to conformation are reviewed and compared with experimental data. MD simulations of explicit solvent around DNA can now give a detailed model of DNA that not only matches well with the experimental data but provides additional insight beyond current experimental limitations. Such simulation results are analyzed with a focus on differential hydration properties between A- and B-DNA and between C/G and A/T base pairs. The extent of hydration is determined from the number of waters in the primary shell and compared to experimental numbers from different measurements. High-resolution hydration patterns around the whole DNA are shown and correlated with the conformations. The role of ions associating with DNA is discussed with respect to changes in the hydration structure correlating with DNA conformation.

Animals↗
PubMed1998

Reconstructing the protein-water interface.

Using molecular dynamics simulations of fully hydrated proteins and analysis of crystal structures contained in the Protein Data Bank, we develop a transferable set of perpendicular radial distribution functions for water molecules around globular proteins. These universal functions may be used to reconstruct the unique three-dimensional solvent density distribution around every individual protein with a modest error. We discuss potential applications of this solvent treatment in protein x-ray crystallographic refinements and in theoretical modeling. We also present a fast, grid-based algorithm for construction of the perpendicular solvent density distributions.

Algorithms↗
PubMed1998

Structural equilibrium of DNA represented with different force fields.

We have recently indicated preliminary evidence of different equilibrium average structures with the CHARMM and AMBER force fields in explicit solvent molecular dynamics simulations on the DNA duplex d(C5T5) . d(A5G5) (Feig, M. and B.M. Pettitt, 1997, Experiment vs. Force Fields: DNA conformation from molecular dynamics simulations. J. Phys. Chem. B. (101:7361-7363). This paper presents a detailed comparison of DNA structure and dynamics for both force fields from extended simulation times of 10 ns each. Average structures display an A-DNA base geometry with the CHARMM force field and a base geometry that is intermediate between A- and B-DNA with the AMBER force field. The backbone assumes B form on both strands with the AMBER force field, while the CHARMM force field produces heterogeneous structures with the purine strand in A form and the pyrimidine strand in dynamical equilibrium between A and B conformations. The results compare well with experimental data for the cytosine/guanine part but fail to fully reproduce an overall B conformation in the thymine/adenine tract expected from crystallographic data, particularly with the CHARMM force field. Fluctuations between A and B conformations are observed on the nanosecond time scale in both simulations, particularly with the AMBER force field. Different dynamical behavior during the first 4 ns indicates that convergence times of several nanoseconds are necessary to fully establish a dynamical equilibrium in all structural quantities on the time scale of the simulations presented here.

Base Sequence↗
PubMed1998

Diffusion of solvent around biomolecular solutes: a molecular dynamics simulation study.

Effects of the macromolecular solute on the translational mobility of surrounding solvent water, and Na+ and Cl- ions are investigated by molecular dynamics (MD) simulation. Using MD trajectories of myoglobin and d(C5T5) . d(G5A5) DNA decamer of high quality and length, we determine the average diffusion coefficients for all solvent species as a function of distance from the closest solute atom. We examine solvent mobility in the directions parallel and perpendicular to the solute surface and in proximity to three different classes of solute atoms (oxygens, nitrogens, and carbons). The nature and the magnitude of the solute effects on water diffusion appear to be very similar for protein and DNA decamer. The overall diffusion rate at the interface is lower than in the bulk. The rate is higher than the average in the direction parallel to the solute surface, and lower in the direction normal to the surface, up to 15 A away from the solute. The rate is also lower in the solvation shells of the macromolecules, producing characteristic depressions in the radial profiles of the diffusion coefficient that can be correlated with peaks in the corresponding radial distribution functions. The magnitude of these depressions is small compared to the overall change in solvent mobility at the interface. Similar features are observed in the radial profiles of the diffusion coefficient of sodium and chlorine ions as well.

Biophysical Phenomena↗
PubMed1998

Comparison of simulated and experimentally determined dynamics for a variant of the Lacl DNA-binding domain, Nlac-P.

Recent advances in the experimentally determined structures and dynamics of the domains within LacI provide a rare context for evaluating dynamics calculations. A 1500-ps trajectory was simulated for a variant of the LacI DNA-binding domain, which consists of the first three helices in LacI and the hinge helix of the homologous PurR. Order parameters derived from dynamics simulations are compared to those obtained for the LacI DNA-binding domain with 15N relaxation NMR spectroscopy (Slijper et al., 1997. Biochemistry. 36:249-254). The MD simulations suggest that the unstructured loop between helices II and III does not exist in a discrete state under the conditions of no salt and neutral pH, but occupies a continuum of states between the DNA-bound and free structures. Simulations also indicate that the unstructured region between helix III and the hinge helix is very mobile, rendering motions of the hinge helix essentially independent of the rest of the protein. Finally, the alpha-helical hydrogen bonds in the hinge helix are broken after 1250 ps, perhaps as a prelude to helix unfolding.

Amino Acid Sequence↗
PubMed1998

Protein hydration density: theory, simulations and crystallography.

Models of protein hydration are becoming increasingly more accurate in comparison with experimental data. The recent success of these models implies that the major features of the solvation layers are dominated by local correlations and that such correlations are universal. The excellent agreement between theoretical and experimental solvent electron density radial distributions marks a significant success in our ability to accurately model macromolecular hydration.

Crystallography↗