Showing papers by "David A. Case published in 1998"

Continuum Solvent Studies of the Stability of DNA, RNA, and Phosphoramidate−DNA Helices

Abstract: We apply continuum solvent models to investigate the relative stability of A- and B-form helices for three DNA sequences, d(CCAACGTTGG)2, d(ACCCGCGGGT)2, and d(CGCGAATTCGCG)2, a phosphoramidate-mod...

Continuum solvent studies of the stability of RNA hairpin loops and helices.

Abstract: We apply continuum solvent models to investigate the relative stability of various conformational forms for two RNA sequences, GGAC(UUCG)GUCC and GGUG(UGAA)CACC. In the first part, we compare alternate hairpin conformations to explore the reliability of these models to discriminate between different local conformations. A second part looks at the hairpin-duplex conversion for the UUCG sequence, identifying major contributors to the thermodynamics of a much large scale transition. Structures were taken as snapshots from multi-nanosecond molecular dynamics simulations computed in a consistent fashion using explicit solvent and with long-range electrostatics accounted for using the Particle-Mesh Ewald procedure. The electrostatic contribution to solvation energies were computed using both a finite-difference Poisson-Boltzmann (PB) model and a pairwise Generalized Born model; non-electrostatic contributions were estimated with a surface-area dependent term. To these solvation free energies were added the mean solute internal energies (determined from a molecular mechanics potential) and estimates of the solute entropy (from a harmonic analysis). Consistent with experiment and with earlier solvated molecular dynamics simulations, the UUCG hairpin was found to prefer conformers close to a recent NMR structure determination in preference to those from an earlier NMR study. Similarly, results for the UGAA hairpin favored an NMR-derived structure over that to be expected for a generic GNRA hairpin loop. Experimental free energies are not known for the hairpin/duplex conversion, but must be close to zero since hairpins are seen in solution and duplexes in crystals; out calculations find a value near zero and illustrate the expected interplay of solvation, salt effects and entropy in affecting this equilibrium.

Molecular Dynamics and Continuum Solvent Studies of the Stability of PolyG-PolyC and PolyA-PolyT DNA Duplexes in Solution

Abstract: Molecular dynamics simulation in explicit solvent and continuum solvent models are applied to investigate the relative stability of A- and B-form helices for two DNA sequences, dA10-dT10 and dG10-dC10 in three structural forms. One structural form is based on an unrestrained molecular dynamics (MD) trajectory starting from a canonical B-DNA structure, the second is based on a MD trajectory starting in a canonical B-DNA structure with the sugars constrained to be C2'-endo and the third simulation started from a canonical A-DNA structure with the sugars constrained to C3'-endo puckers. For the energetic analysis, structures were taken as snapshots from nanosecond length molecular dynamics simulations computed in a consistent fashion in explicit solvent, applying the particle mesh Ewald method and the Cornell et al. force field. The electrostatic contributions to solvation free energies are computed using both a finite-difference Poisson-Boltzmann model and a pairwise Generalized Born model. The non-electrostatic contributions to the solvation free energies are estimated with a solvent accessible surface area dependent term. To estimate the gas phase component of the relative free energy between the various structures, the mean solute internal energies (determined with the Cornell et al. molecular mechanics potential including all pairwise interactions within the solute) and estimates of the solute entropy (using a harmonic approximation) were used. Consistent with experiment, the polyG-polyC (GC) structures are found to be much more A-phillic than the polyA-polyT (AT) structures, the latter being quite A-phobic. The dominant energy components responsible for this difference comes from the internal and van der Waal energies. A perhaps less appreciated difference between the GC and AT rich sequences is suggested by the calculated salt dependence which demonstrates a significantly enhanced ability to drive GC rich sequences towards an A-form structure compared to AT rich sequences. In addition to being A-phobic, the AT structure also has a noticably larger helical repeat than GC and other mixed sequence duplexes, consistent with experiment. Analysis of the average solvent density from the trajectories shows hydration patterns in qualitative agreement with experiment and previous theoretical treatments.

Density Functional Study of Ribose and Deoxyribose Chemical Shifts

Abstract: Density functional chemical shielding calculations are reported for methyl β-d-2-deoxyribofuranoside and for methyl β-d-ribofuranoside, models for the deoxyribose and ribose sugars in nucleic acids. The variation of the chemical shielding as a function of the sugar ring conformation is reported, as well as the influence of the ring conformation on the chemical shift anisotropy. The calculated chemical shieldings are sensitive to the puckering of the sugar ring. The value of the exocyclic torsion angles, particularly γ(O5‘−C5‘−C4‘−C3‘), are also found to influence the chemical shielding of the ring atoms. The chemical shielding of the C3‘ carbon is the most sensitive to the sugar ring pucker, with a variation of 10 ppm between the C3‘ endo and C2‘ endo conformations. H3‘ and H4‘ hydrogen shieldings vary by 0.4−0.6 ppm between the C3‘ endo and C2‘ endo conformations. Chemical shift anisotropies at C1‘ and C3‘ are strongly influenced by sugar pucker. Our results agree well with experimentally reported values...

Calculations of proton-binding thermodynamics in proteins.

Abstract: Computational models of proton binding can range from the chemically complex and statistically simple (as in the quantum calculations) to the chemically simple and statistically complex. Much progress has been made in the multiple-site titration problem. Calculations have improved with the inclusion of more flexibility in regard to both the geometry of the proton binding and the larger scale protein motions associated with titration. This article concentrated on the principles of current calculations, but did not attempt to survey their quantitative performance. This is (1) because such comparisons are given in the cited papers and (2) because continued developments in understanding conformational flexibility and interaction energies will be needed to develop robust methods with strong predictive power. Nevertheless, the advances achieved over the past few years should not be underestimated: serious calculations of protonation behavior and its coupling to conformational change can now be confidently pursued against a backdrop of increasing understanding of the strengths and limitations of such models. It is hoped that such theoretical advances will also spur renewed experimental interest in measuring both overall titration curves and individual pKa values or pKa shifts. Exploration of the shapes of individual titration curves (as measured by Hill coefficients and other parameters) would also be useful in assessing the accuracy of computations and in drawing connections to functional behavior.

Distributed torsion angle grid search in high dimensions: A systematic approach to NMR structure determination

Abstract: Two complementary approaches for systematic search in torsion angle space are described for the generation of all conformations of polypeptides which satisfy experimental NMR restraints, hard-sphere van der Waals radii, and rigid covalent geometry. The first procedure is based on a recursive, tree search algorithm for the examination of linear chains of torsion angles, and uses a novel treatment to propagate the search results to neighboring regions so that the structural consequences of the restraints are fully realized. The second procedure is based on a binary combination of torsion vector spaces for connected submolecules, and produces intermediate results in Cartesian space for a more robust restraint analysis. Restraints for NMR applications include bounds on torsion angles and internuclear distances, including relational and degenerate restraints involving equivalent and nonstereoassigned protons. To illustrate these methods, conformation search results are given for the tetrapeptide APGA restrained to an idealized β-turn conformation, an alanine octapeptide restrained to a right-handed helical conformation, and the structured region of the peptide SYPFDV.

High-resolution solution structure of Bacillus subtilis IIAglc.

Abstract: The high-resolution solution structure of the phosphocarrier protein IIAglc from Bacillus subtilis is determined using 3D and 4D heteronuclear NMR methods. B. subtilis IIAglc contains 162 amino acid residues and is one of the larger proteins for which high-resolution solution structure has been determined by NMR methods. The structures have been calculated from a total of 2,232 conformational constraints. Comparison with the X-ray crystal structure indicates that the overall fold is the same in solution and in crystalline environments, although some local structural differences are observed. These occur largely in turns and loops, and mostly correspond to regions with high-temperature factors in the crystal structure. The N-terminus of IIAglc is disordered in solution. The active site is located in a concave region of the protein surface. The histidine, which accepts the phosphoryl group (His 83), interacts with a neighboring histidine (His 68) and is surrounded by hydrophobic residues.