Characterization of the restricted rotation of the dimethyl groups in chemically N-terminal 13C-labeled antifreeze glycoproteins: a temperature-dependent study in water to ice through the supercooled state.

Site-specific chemical modification, especially with isotopically enriched groups, allows one to study the structure and dynamics of proteins for which uniform enrichment is difficult. When the N-terminal alanine in antifreeze glycoprotein (AFGP) is replaced with an N,N-dimethyl alanine the methyl groups show signatures of slow rotation about the C-N bond. In order to separate the local dynamics of the N-terminus from the overall protein dynamics, we present a complete characterization of this dynamics. Temperature-dependent nuclear magnetic-resonance experiments from room temperature to subzero temperatures, including the supercooled state and in the presence of ice, are presented. Quantum chemical calculations are also performed on a localized N-terminus of the AFGP. Our results show that in the solution state at room temperature and in the super cooled regime, the dimethyl groups undergo a slow, restricted rotation with an unequal distribution of population between two major conformations. At lower temperatures in the presence of ice, the dynamics become much more complex due to freezing out of several conformational states. Based on these results, we conclude that the segmental dynamics of the N-terminus are local to the first residue and do not affect the overall dynamics of the protein.

Brownian dynamics simulations of sequence-dependent duplex denaturation in dynamically superhelical DNA.

The topological state of DNA in vivo is dynamically regulated by a number of processes that involve interactions with bound proteins. In one such process, the tracking of RNA polymerase along the double helix during transcription, restriction of rotational motion of the polymerase and associated structures, generates waves of overtwist downstream and undertwist upstream from the site of transcription. The resulting superhelical stress is often sufficient to drive double-stranded DNA into a denatured state at locations such as promoters and origins of replication, where sequence-specific duplex opening is a prerequisite for biological function. In this way, transcription and other events that actively supercoil the DNA provide a mechanism for dynamically coupling genetic activity with regulatory and other cellular processes. Although computer modeling has provided insight into the equilibrium dynamics of DNA supercoiling, to date no model has appeared for simulating sequence-dependent DNA strand separation under the nonequilibrium conditions imposed by the dynamic introduction of torsional stress. Here, we introduce such a model and present results from an initial set of computer simulations in which the sequences of dynamically superhelical, 147 base pair DNA circles were systematically altered in order to probe the accuracy with which the model can predict location, extent, and time of stress-induced duplex denaturation. The results agree both with well-tested statistical mechanical calculations and with available experimental information. Additionally, we find that sites susceptible to denaturation show a propensity for localizing to supercoil apices, suggesting that base sequence determines locations of strand separation not only through the energetics of interstrand interactions, but also by influencing the geometry of supercoiling.

Photodissociation of CCH: classical trajectory calculations involving seven electronic states.

The photodissociation dynamics of ethynyl radical, C(2)H, involving seven electronic states is studied by classical trajectory calculations. Initial values of the trajectories are selected based on relative absorption intensities calculated by Mebel et al. The energies and the derivatives are interpolated by three-dimensional cubic spline interpolator using an extended data pool. Mean square errors and standard deviations in interpolation of energies for 450 data points are found to be in the range 3.1 x 10(-6)-1.4 x 10(-5) and 1.7 x 10(-3)-3.8 x 10(-3) hartrees, respectively. The photofragments of C(2) and H are produced mainly in the X (1)Sigma(g) (+), a (3)Pi(u), b (3)Sigma(g) (-), c (3)Sigma(u) (+), A (1)Pi(u), B (1)Delta(g) electronic states of C(2) as product. The avoided crossings do not appear to be in the main dissociation pathways. The internal distributions are in good accord with the experimental results where comparison is possible, suggesting that the fragmentation mechanism of C(2)H(2) into C(2) and H is a two step process involving C(2)H radical as an intermediate with a life time long enough to allow complete collection of the phase space in the experiments.

Transcription-driven twin supercoiling of a DNA loop: a Brownian dynamics study.

The torque generated by RNA polymerase as it tracks along double-stranded DNA can potentially induce long-range structural deformations integral to mechanisms of biological significance in both prokaryotes and eukaryotes. In this paper, we introduce a dynamic computer model for investigating this phenomenon. Duplex DNA is represented as a chain of hydrodynamic beads interacting through potentials of linearly elastic stretching, bending, and twisting, as well as excluded volume. The chain, linear when relaxed, is looped to form two open but topologically constrained subdomains. This permits the dynamic introduction of torsional stress via a centrally applied torque. We simulate by Brownian dynamics the 100 micros response of a 477-base pair B-DNA template to the localized torque generated by the prokaryotic transcription ensemble. Following a sharp rise at early times, the distributed twist assumes a nearly constant value in both subdomains, and a succession of supercoiling deformations occurs as superhelical stress is increasingly partitioned to writhe. The magnitude of writhe surpasses that of twist before also leveling off when the structure reaches mechanical equilibrium with the torsional load. Superhelicity is simultaneously right handed in one subdomain and left handed in the other, as predicted by the "transcription-induced twin-supercoiled-domain" model [L. F. Liu and J. C. Wang, Proc. Natl. Acad. Sci. U.S.A. 84, 7024 (1987)]. The properties of the chain at the onset of writhing agree well with predictions from theory, and the generated stress is ample for driving secondary structural transitions in physiological DNA.

Intermolecular interaction studies of winter flounder antifreeze protein reveal the existence of thermally accessible binding state.

The physical nature underlying intermolecular interactions between two rod-like winter flounder antifreeze protein (AFP) molecules and their implication for the mechanism of antifreeze function are examined in this work using molecular dynamics simulations, augmented with free energy calculations employing a continuum solvation model. The energetics for different modes of interactions of two AFP molecules is examined in both vacuum and aqueous phases along with the water distribution in the region encapsulated by two antiparallel AFP backbones. The results show that in a vacuum two AFP molecules intrinsically attract each other in the antiparallel fashion, where their complementary charge side chains face each other directly. In the aqueous environment, this attraction is counteracted by both screening and entropic effects. Therefore, two nearly energetically degenerate states, an aggregated state and a dissociated state, result as a new aspect of intermolecular interaction in the paradigm for the mechanism of action of AFP. The relevance of these findings to the mechanism of function of freezing inhibition in the context of our work on Antarctic cod antifreeze glycoprotein (Nguyen et al., Biophysical Journal, 2002, Vol. 82, pp. 2892-2905) is discussed.

Translational dynamics of antifreeze glycoprotein in supercooled water.

Structure and dynamics of biomolecules in supercooled water assume a particular and distinct importance in the case of antifreeze glycoproteins (AFGPs), which function at sub-zero temperatures. To investigate whether any large-scale structural digressions in the supercooled state are correlated to the function of AFGPs, self-diffusion behavior of the AFGP8, the smallest AFGP is monitored as a function of temperature from 243 to 303 K using nuclear magnetic resonance (NMR) spectroscopy. The experimental results are compared with the hydrodynamic calculations using the viscosity of water at the same temperature range. In order to evaluate results on AFGP8, the smallest AFGP, constituting approximately two-thirds of the total AFGP fraction in fish blood serum, similar experimental and computational calculations were also performed on a set of globular proteins. These results show that even though the general trend of translational dynamics of AFGP is similar to that of the other globular proteins, AFGP8 appears to be more hydrated (approximately 30% increase in the bead radius) than the others over the temperature range studied. These results also suggest that local conformational changes such as segmental librations or hydrogen bond dynamics that are closer to the protein surface are more likely the determining dynamic factors for the function of AFGPs rather than any large-scale structural rearrangements.

Synthesis and characterization of the monomeric imides Ar'MNAr' ' (M = Ga or In; Ar' or Ar' ' = terphenyl ligands) with two-coordinate gallium and indium.

Reaction of Ar'MMAr' (M = Ga or In) with N3Ar' ' (Ar' = C6H3-2,6-Dipp2, Dipp = C6H3-2,6-Pri2, Ar' ' = C6H3-2,6(Xyl-4-But)2) afforded the first monomeric imides of heavier group 13 elements with two-coordinate metals. Planar, trans-bent structures with short M-N bond distances were observed, which are consistent with lone pair character at both M and N and a bond order less than the formally expected triple one.

Binding to an RNA aptamer changes the charge distribution and conformation of malachite green.

RNA plays a central role in many biological processes and is therefore an important target for drug development. In recent years an increasing wealth of structural and functional information about RNA-ligand complexes has been obtained using in vitro selected RNAs (aptamers). However, all those studies focused on structure and changes of the nucleic acid and mostly considered the ligand as a rigid target. To develop a detailed picture of ligand structure and dynamics in RNA-small molecule complexes, the malachite green binding aptamer was studied. Isotopically labeled ligand in complex with RNA was analyzed by NMR spectroscopy in solution. The surprisingly asymmetric changes in the (13)C chemical shift of the ligand methyl groups indicate that the dye undergoes changes in its conformation and charge distribution upon binding. The role of the RNA electrostatic field in this interaction was explored using ab initio calculations of the ligand structure and charge distribution. The results indicate that the uneven charge distribution in the RNA binding pocket provides a major contribution to the driving force of the ligand structural changes. The observation that not only the RNA adapts to the ligand, in what is called adaptive binding, but that the ligand itself also undergoes conformational changes ("induced fit") is crucial for the rational design of RNA ligands and for understanding the properties of RNA-ligand complexes.

The dynamics, structure, and conformational free energy of proline-containing antifreeze glycoprotein.

Recent NMR studies of the solution structure of the 14-amino acid antifreeze glycoprotein AFGP-8 have concluded that the molecule lacks long-range order. The implication that an apparently unstructured molecule can still have a very precise function as a freezing inhibitor seems startling at first consideration. To gain insight into the nature of conformations and motions in AFGP-8, we have undertaken molecular dynamics simulations augmented with free energy calculations using a continuum solvation model. Starting from 10 different NMR structures, 20 ns of dynamics of AFGP were explored. The dynamics show that AFGP structure is composed of four segments, joined by very flexible pivots positioned at alanine 5, 8, and 11. The dynamics also show that the presence of prolines in this small AFGP structure facilitates the adoption of the poly-proline II structure as its overall conformation, although AFGP does adopt other conformations during the course of dynamics as well. The free energies calculated using a continuum solvation model show that the lowest free energy conformations, while being energetically equal, are drastically different in conformations. In other words, this AFGP molecule has many structurally distinct and energetically equal minima in its energy landscape. In addition, conformational, energetic, and hydrogen bond analyses suggest that the intramolecular hydrogen bonds between the N-acetyl group and the protein backbone are an important integral part of the overall stability of the AFGP molecule. The relevance of these findings to the mechanism of freezing inhibition is discussed.

Rapidly converging lattice sums for nonelectrostatic interactions.

The relative energies of one-, two-, and three-dimensional Bravais lattice Lennard-Jones particles can be calculated by lattice sums. The expression of lattice sums over a Lennard-Jones potential can be manipulated into a form that converges rapidly. A formalism capable of calculating the lattice potential at arbitrary points of a completely general lattice has been developed. This method provides an alternative way to calculate the relative energies from the surface and the interior bulk sites of many chemical systems. The method is illustrated with application to hcp and fcc Lennard-Jonesium, both for the relative binding energy and for calculating the potential along the geometric diffusion pathway between tetrahedral and octahedral interstitial sites. Diffusion from the tetrahedral site to the octahedral site experiences a barrier of 752.600 in units of 4 epsilon. The reverse pathway experiences a barrier of 1035.614 in units of 4 epsilon.

Dynamics of antifreeze glycoproteins in the presence of ice.

Antifreeze glycoproteins from the Greenland cod Boreogadus saida were dimethylated at the N-terminus (m*AFGP) and their dynamics and conformational properties were studied in the presence of ice using (13)C-NMR and FTIR spectroscopy. (13)C-NMR experiments of m*AFGP in D(2)O, in H(2)O, and of freeze-dried m*AFGP were performed as a function of temperature. Dynamic parameters ((1)H T(1 rho) and T(CH)) obtained by varying the contact time revealed notable differences in the motional properties of AFGP between the different states. AFGP/ice dynamics was dominated by fast-scale motions (nanosecond to picosecond time scale), suggesting that the relaxation is markedly affected by the protein hydration. The data suggest that AFGP adopts a similar type of three-dimensional fold both in the presence of ice and in the freeze-dried state. FTIR studies of the amide I band did not show a single prevailing secondary structure in the frozen state. The high number of conformers suggests a high flexibility, and possibly reflects the necessity to expose more ice-binding groups. The data suggest that the effect of hydration on the local mobility of AFGP and the lack of significant change in the backbone conformation in the frozen state may play a role in inhibiting the ice crystal growth.

Stable, Monomeric Imides of Aluminum and Gallium: Synthesis and Characterization of [{HC(MeCDippN)2 }MN-2,6-Trip2 C6 H3 ] (M=Al or Ga; Dipp=2,6-iPr2 C6 H3 ; Trip=2,4,6-iPr3 C6 H2 ).

A short Ga-N bond with double-bond character is displayed by the first monomeric imide of gallium, which was obtained by the reaction of [{HC(MeCDippN)2 }M:] (Dipp=2,6-iPr2 C6 H3 , M=Ga; see picture) with N3 -2,6-Trip2 C6 H3 (Trip=2,4,6-iPr3 C6 H2 ). The analogous aluminum (M=Al) compound is also readily available.

Comparison of pi-Bond Strengths in M-E (M = B, Al, Ga; E = O, N, S) Compounds. Ab Initio Calculation of Rotational Barriers.

Results of ab initio calculations of the electronic structure of compounds of the type R(2)MER'(x) and R(2)MEMR'(2) with R = H, Me, M = Al, Ga, and E = O, N, S are reported at the Hartree-Fock level with split-valence, polarization basis sets for all atoms except hydrogen where a split-valence basis set is used. Full optimizations for the equilibrium geometry and partial optimization at constrained rotational transition states have been performed to evaluate the barriers to rotation as a measure of the pi interactions in these compounds. We conclude that, although important for determining the final conformational equilibrium geometries, pi interactions are weak in these compounds as the rotational barriers are smaller than that for ethylene by 1 or 2 orders of magnitude.