Biomedical subjects

B M Pettitt

Publications and source records attributed to B M Pettitt.

At least 19 recordsLinked to original sources

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

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↗

PubMed1998

NMR and quenched molecular dynamics studies of superpotent linear and cyclic alpha-melanotropins.

Conformational searching, computer simulations, synthesis and NMR are used on a variety of alpha melanocyte-stimulating hormone (alpha-MSH) analogues to understand the physical characteristics required for biological potency. Peptides I (Ac-[Nle4,Asp5,D-Phe7,Lys10]alpha-MSH(4-10)-NH2), II (Ac-c[Nle4,Asp5,D-Phe7,Lys10]alpha-MSH(4-10)-NH2) and III (Ac-[Nle4,Asp5,D-Phe7,Dap10]alpha-MSH(4-10)-NH2 all show very similar conformational properties (backbone and side-chain torsional angles), and all display high biological potencies. The modeling results for these compounds are supported by the NMR data. Peptide IV (Ac-c[Nle4,Asp5,D-Phe7,Dap10]alpha-MSH(4-10)-NH2) appears to have a markedly different conformation and has decreased biological potency.

Animals↗

PubMed1997

Modeling the DNA-solvent interface.

We extend the technique of using perpendicular distribution functions to salt solutions around nucleic acids. Both solute density averaged and nonaveraged reference frames are considered and compared. Using a previous simulation of DNA in salt water of over a nanosecond in duration, the aqueous distribution functions were found to be well coveraged, whereas the salt perpendicular distribution functions were less well determined. Three-dimensional density reconstructions reliably showed the prominent solvation features with transferable functions. The number of solute atom types needed for reconstructions of a given precision was determined in the context of the reference simulation data set with the goal of achieving a required level of reconstruction quality.

Computer Simulation↗

PubMed1997

A simple two-dimensional representation for the common secondary structural elements of polypeptides and proteins.

A simple method is presented for projecting the conformation of extended secondary structure elements of peptides and proteins that extend over four C alpha atoms onto a simple two-dimensional surface. A new set of two degrees of freedom is defined, a pseudodihedral involving four sequential C alpha atoms, as well as the triple scalar product for the vectors describing the orientation of the three intervening peptide groups. The method provides a reduction in dimensionality, from the usual combination of multiple phi,psi pairs to a single pair, yielding valuable information concerning the structure and dynamics of these important elements. The new two-dimensional surface is explored by reference to 63 selected protein crystal structures together with a comparison of model built peptides representing the common secondary structural elements. Dynamical aspects on this new surface are examined using a molecular dynamics trajectory of Basic Pancreatic Trypsin Inhibitor.

Aprotinin↗

PubMed1996

Salt effects on peptide conformers: a dielectric study of tuftsin.

Four 1-ns molecular dynamics computer simulations of tuftsin, Thr-Lys-Pro-Arg, are analyzed: (1) cis tuftsin in water, (2) trans tuftsin in water, (3) cis tuftsin in 1 M NaCl, and (4) trans tuftsin in 1 M NaCl. Independently of the salt concentration, the trans conformer has a higher dielectric constant than the cis conformer because the former exhibits a more widely distributed charge distribution in space. Independently of the peptide conformation, the presence of salt reduces the dielectric constants of both the peptide and the solvating water molecules because ions, on binding, restrict the motion of other atoms. In contrast to the dielectric constants, neither the peptide conformation nor the salt concentration shows a significant influence on the dielectric relaxation time of water molecules.

Electrochemistry↗

PubMed1995

Structure and stability of a model pyrimidine-purine-purine DNA triple helix with a GC.T mismatch by simulation.

A 1.5 ns long molecular dynamics simulation was conducted to compare the structure and stability of a model DNA triplex in saline solution with that found from experiments. The model DNA was an antiparallel py.pu.pu (CG.G) 7-mer structure which contained a GC.T mismatch triplet at the middle of the sequence. The local conformation of the mismatch triplet and the effects of this triplet on the global helical structure suggest that the GC.T triplet forms stable hydrogen bonds and shows distortions from an in-plane alignment. The overall rms deviation of the triplex is similar to one without a mismatch, although the thymine base in the mismatch triplet shows significantly higher mobility. A high coordination probability for water between the G and T bases in the mismatch triplet was observed to have an effect on the stability of non-hydrogen-bonded base pairs. Average helical parameters, sugar pucker, and backbone dihedral angles indicate that the CG.G triplets on the 3' side of the mismatch triplet possess different structural and dynamical properties than that of the 5' side. These observations are consistent with recently available experimental results and provide an interpretation of the observed experimental structure. They also suggest that inclusion of explicit water molecules is necessary in order to understand and predict the interaction between the third strand and duplex DNA.

Computer Simulation↗

PubMed1995

A sampling problem in molecular dynamics simulations of macromolecules.

Correlations in low-frequency atomic displacements predicted by molecular dynamics simulations on the order of 1 ns are undersampled for the time scales currently accessible by the technique. This is shown with three different representations of the fluctuations in a macromolecule: the reciprocal space of crystallography using diffuse x-ray scattering data, real three-dimensional Cartesian space using covariance matrices of the atomic displacements, and the 3N-dimensional configuration space of the protein using dimensionally reduced projections to visualize the extent to which phase space is sampled.

Computer Simulation↗

PubMed1995

Solvent effects on model d(CG.G)7 and d(TA.T)7 DNA triple helices.

Using free energy molecular mechanics, we find that the molecular effects of solvent are critical in determining relative stabilities in DNA triple helices or triplexes. The continuum solvent model is unable to differentiate the thermodynamics reflecting the basic solvation differences around the occupied major groove in triplexes. In order to avoid the local minimum problem, which is a major limitation of any modeling study, we started our computations with multiple structures rather than relying on the optimization of a single reference structure. By constructing triplex models with different initial helical twists, helical rises, and sugar-pucker permutations, we explore the potential surface and the structural preference with respect to these variations. We find that in order to accommodate a third strand in triplex formation, the backbone geometry of the B-DNA duplex target has to be adjusted into A-DNA-like form with a deep major groove. This is achieved by concerted adjustment in torsions beta, epsilon, and zeta around the phosphate groups. However, the sugar pucker displays a more rich variation, resulting in conformations not usually associated with the canonical duplex structures.

DNA↗

PubMed1995

Conformational states governing the rates of spontaneous transition mutations.

We report the results of a theoretical study, combining the results of sequence analysis and integral equation structural methods for nucleic acids in aqueous solutions, on the effects of nearest neighbors on the (T.G) mispair in solution, for 12 nearest neighbor contexts. Attempts have been made to classify the structural and energetic effects of the 5' and 3' neighbors with respect to the observed spontaneous mutation rates in vertebrates. It is found that 5' nearest neighbor is probably the most critical structural factor in facilitating or discouraging mutations. Local conformational states correlate with discrimination of bases to be excised in mispairs. Our study confirms the role of the flexibility of the DNA molecule in governing the rates of spontaneous mutations.

Base Composition↗

PubMed1995

Structure and dynamics of the water around myoglobin.

The interplay between simulations at various levels of hydration and experimental observables has led to a picture of the role of solvent in thermodynamics and dynamics of protein systems. One of the most studied protein-solvent systems is myoglobin, which serves as a paradigm for the development of structure-function relationships in many biophysical studies. We review here some aspects of the solvation of myoglobin and the resulting implications. In particular, recent theoretical and simulation studies unify much of the diverse set of experimental results on water near proteins.

Crystallography, X-Ray↗

PubMed1995

Dielectric response of triplex DNA in ionic solution from simulations.

We have analyzed a 1.2-ns molecular dynamics simulation of 51 mM d(CG.G)7 with 21 Na+ counter-ions and 1 M NaCl in water. Via the dipole fluctuations, the dielectric constant for the DNA is found to be around 16, whereas that for the bases and sugars combined is only 3. The dielectric constant for water in this system is 41, which is much smaller than 71 for pure SPC/E water, because of the strong restriction imposed on the motion of water molecules by the DNA and the ions. Also addressed in the present work are several technical issues related to the calculation of the dipole moment of an ionic solution from molecular dynamics simulations using periodic boundary conditions.

DNA↗

PubMed1995

Simulations of conformers of tuftsin and a cyclic tuftsin analog.

The conformational properties of the configurational isomers of tuftsin, a linear tetrapeptide with the sequence Thr-Lys-Pro-Arg, were investigated with six 1 ns molecular dynamics simulations in explicit water and in a 1.0 m NaCl solution. The average conformation of the cis isomer is a type VI beta-turn. Our results indicate that water-peptide hydrogen bonding, in addition to intramolecular hydrogen bonds, stabilizes the cis conformer. The trans isomer is neither a beta- nor a gamma-turn. Results are compared with parallel studies on a cyclic analog of tuftsin, cyclo(Thr-Lys-Pro-Arg-Gly). The addition of salt does, not influence the backbone conformation of the peptide. Differences between the structures are confined to the side-chain orientations of the Lys and Arg residues.

Amino Acid Sequence↗

PubMed1994

A connected-cluster of hydration around myoglobin: correlation between molecular dynamics simulations and experiment.

An analysis of a molecular dynamics simulation of metmyoglobin in an explicit solvent environment of 3,128 water molecules has been performed. Both statics and dynamics of the protein-solvent interface are addressed in a comparison with experiment. Three-dimensional density distributions, temperature factors, and occupancy weights are computed for the solvent by using the trajectory coordinates. Analysis of the hydration leads to the localization of more than 500 hydration sites distributed into multiple layers of solvation located between 2.6 and 6.8 A from the atomic protein surface. After locating the local solvent density maxima or hydration sites we conclude that water molecules of hydration positions and hydration sites are distinct concepts. Both global and detailed properties of the hydration cluster around myoglobin are compared with recent neutron and X-ray data on myoglobin. Questions arising from differences between X-ray and neutron data concerning the locations of the protein-bound water are investigated. Analysis of water site differences found from X-ray and neutron experiments compared with our simulation shows that the simulation gives a way to unify the hydration picture given by the two experiments.

Amino Acid Sequence↗