Myoglobin-CO substate structures and dynamics: multidimensional vibrational echoes and molecular dynamics simulations

Spectrally resolved infrared stimulated vibrational echo data were obtained for sperm whale carbonmonoxymyoglobin (MbCO) at 300 K. The measured dephasing dynamics of the CO ligand are in agreement with dephasing dynamics calculated with molecular dynamics (MD) simulations for MbCO with the residue histidine-64 (His64) having its imidazole epsilon nitrogen protonated (N(epsilon)-H). The two conformational substate structures B(epsilon) and R(epsilon) observed in the MD simulations are assigned to the spectroscopic A(1) and A(3) conformational substates of MbCO, respectively, based on the agreement between the experimentally measured and calculated dephasing dynamics for these substates. In the A(1) substate, the N(epsilon)-H proton and N(delta) of His64 are approximately equidistant from the CO ligand, while in the A(3) substate, the N(epsilon)-H of His64 is oriented toward the CO, and the N(delta) is on the surface of the protein. The MD simulations show that dynamics of His64 represent the major source of vibrational dephasing of the CO ligand in the A(3) state on both femtosecond and picosecond time scales. Dephasing in the A(1) state is controlled by His64 on femtosecond time scales, and by the rest of the protein and the water solvent on longer time scales.     

Prediction of Complexation Properties of Crown Ethers Using Computational Neural Networks

A computational neural network method was used for the prediction of stability constants of simple crown ether complexes. The essence of the method lies in the ability of a computer neural network to recognize the structure-property relationships in these host-guest systems. Testing of the computational method has demonstrated that stability constants of alkali metal cation (Na+, K+, Cs+)-crown ether complexes in methanol at 25 °C can be predicted with an average error of ±0.3 log K units based on the chemical structure of the crown ethers alone. The computer model was then used for the preliminary analysis of trends in the stabilities of the above complexes.     

Classical and quantum mechanical infrared echoes from resonantly coupled molecular vibrations

The nonlinear response function associated with the infrared vibrational echo is calculated for a quantum mechanical model of resonantly coupled, anharmonic oscillators at zero temperature. The classical mechanical response function is determined from the quantum response function by setting variant Planck's over 2pi-->0, permitting the comparison of the effects of resonant vibrational coupling among an arbitrary number of anharmonic oscillators on quantum and classical vibrational echoes. The quantum response function displays a time dependence that reflects both anharmonicity and resonant coupling, while the classical response function depends on anharmonicity only through a time-independent amplitude, and shows a time dependence controlled only by the resonant coupling. In addition, the classical response function grows without bound in time, a phenomenon associated with the nonlinearity of classical mechanics, and absent in quantum mechanics. This unbounded growth was previously identified in the response function for a system without resonant vibrational energy transfer, and is observed to persist in the presence of resonant coupling among vibrations. Quantitative agreement between classical and quantum response functions is limited to a time scale of duration inversely proportional to the anharmonicity.     

Development of Nanosensors and Bioprobes

We describe the development and application of nanosensors having bioreceptor probes for bioanalysis. The nanoprobes were fabricated with optical fibers pulled down to tips having distal end sizes of approximately 30–60nm. The use of two different types of receptors was investigated. Fiberoptic nanoprobes were covalently bound either with bioreceptors, such as antibodies, or with other receptors, such as cyclodextrins that are selective for the size and chemical structure of the analyte molecules. Theoretical calculations were performed to model the binding of beta-cyclodextrin with pyrene and 5,6-benzoquinoline, and to illustrate the possibility of comparing experimental data with theoretical data. The antibody-based nanoprobe was used for in situ measurements of benzopyrene tetrol in single cells. The performance of the nanosensor is illustrated by intracellular measurements performed on a rat liver epithelial cell line (Clone 9) used as the model cell system. The usefulness and potential of these nanotechnology-based biosensors in biological research and applications are discussed.     

Spectra of atomic Rydberg states in strong magnetic fields

Transition frequencies and intensities for exited states of atomic hydrogen in a strong magnetic field are calculated by Fourier transformation of dipole moments computed using classical trajectories. We consider states with n = 30 and L_z = 1 for fields of up to 7 T. Spectra for states labeled as librators and rotators are found to be qualitatively different, especially for the spectral component perpendicular to the field. In addition to the zero-field Kepler line, a new intra-manifold transition is located at low frequency. The frequencies and intensities are found to be sensitive functions of the field strength, and of the particular state of the Coulomb manifold involved.     

Scattering in the quasiperiodic and chaotic regimes

The results of a quasiclassical investigation of the scattering process for a collinear collision between a particle and two coupled anharmonic oscillators is reported. The energy transfer was found to fall at the critical energy where the oscillators made the transition from quasiperiodic to chaotic motion.     

Thermal analysis via molecular dynamics simulation

Thermal analysis by classical molecular dynamics simulations is discussed on hand of heat capacity of crystals of 9600 atoms. The differences between quantum mechanical and classical mechanical calculations are shown. Anharmonicity is proven to be an important factor. Finally, it is found that defects contribute to an increase in heat capacity before melting. The energy of conformational gauche defects within the crystal is only about 10% due to internal rotation. The other energy must be generated by cooperative strain. The conclusion is that the next generation of faster computers may permit wider use of molecular dynamics simulations in support of the interpretation of thermal analysis.Dedicated to Professor Bernhard Wunderlich on the occasion of his 65th birthday     

Multiscale coarse-graining and structural correlations: connections to liquid-state theory

A statistical mechanical framework elucidates the significance of structural correlations between coarse-grained (CG) sites in the multiscale coarse-graining (MS-CG) method (Izvekov, S.; Voth, G. A. J. Phys. Chem. B 2005, 109, 2469; J. Chem. Phys. 2005, 123, 134105). If no approximations are made, then the MS-CG method yields a many-body multidimensional potential of mean force describing the interactions between CG sites. However, numerical applications of the MS-CG method typically employ a set of pair potentials to describe nonbonded interactions. The analogy between coarse-graining and the inverse problem of liquid-state theory clarifies the general significance of three-particle correlations for the development of such CG pair potentials. It is demonstrated that the MS-CG methodology incorporates critical three-body correlation effects and that, for isotropic homogeneous systems evolving under a central pair potential, the MS-CG equations are a discretized representation of the well-known Yvon-Born-Green equation. Numerical calculations validate the theory and illustrate the role of these structural correlations in the MS-CG method.