Targeted molecular dynamics simulation studies of calcium binding and conformational change in the C-terminal half of gelsolin

Gelsolin consists of six related domains (G1-G6) and the C-terminal half (G4-G6) acts as a calcium sensor during the activation of the whole molecule, a process that involves large domain movements. In this study, we used targeted molecular dynamics simulations to elucidate the conformational transitions of G4-G6 at an atomic level. Domains G4 and G6 are initially ruptured, followed by a rotation of G6 by approximately 90 degrees , which is the dominant conformational change. During this period, local conformational changes occur at the G4 and G5 calcium-binding sites, facilitating large changes in interdomain distances. Alterations in the binding affinities of the calcium ions in these three domains appear to be related to local conformational changes at their binding sites. Analysis of the relative stabilities of the G4-G6-bound calcium ions suggests that they bind first to G6, then to G4, and finally to G5.     

Low temperature optical absorption spectroscopy: an approach to the study of stereodynamic properties of hemeproteins

In this short review we show how suitable analysis of the temperature dependence of the optical absorption spectra of metalloproteins can give insight into their stereodynamic properties in the region of the chromophore. To this end, the theory of coupling between an intense allowed electronic transition of a chromophore and Franck-Condon active vibrations of the nearby atoms is applied to the Soret band of hemeproteins to obtain an analytical expression suitable for fitting the spectral profile at various temperatures. The reported approach enables one to separate the various contributions to the overall bandwidth together with the parameters that characterize the vibrational coupling. The thermal behavior of these quantities gives information on the dynamic properties of the active site and on their dependence upon protein structure and ligation state. The Soret band of hemeproteins appears to be coupled to high frequency vibrational modes of the heme group (as already shown by resonance Raman spectroscopy) and to a bath of low frequency modes most likely deriving from the bulk of the protein. For the deoxy derivatives inhomogeneous broadening arising from conformational heterogeneity appears to contribute substantially to the linewidth. The data indicate the onset; at temperatures near 180 K, of large scale anharmonic motions that can be attributed to jumping among different conformational substates of the protein.Abbreviations MbCO Carbonmonoxy-myoglobin - Mb Deoxymyoglobin - Mb³⁺ Aquomet-myoglobin - SWMbCO Spermwhale carbonmonoxy-myoglobin - SWMb Spermwhale deoxy-myoglobin Correspondence to: A. Cupane     

Crystallization and preliminary x-ray studies of I-CreI: a group I intron-encoded endonuclease from C. reinhardtii

Group I intron endonuclease I-CreI is encoded by an open reading frame contained within a self-splicing intron in the Chlamydomonas reinhardtii chloroplast 23S rRNA gene. I-CreI initiates the lateral transfer or homing of this intron by specifically recognizing and cleaving a pseudopalindromic 19–24 bp homing site in chloroplast 23S rRNA genes that lack the intron. The gene encoding this enzyme has been subcloned, and the protein product has been purified and crystallized. The crystals belong to space group P321, with unit cell dimensions a = b = 78.2 Å, c = 67.4 Å. The crystal unit cell is consistent with an asymmetric unit consisting of the enzyme monomer. The specific volume of this unit cell is 3.3 ų/Da. The crystals diffract to at least 3.0 Å resolution after flash-cooling, when using a rotating anode x-ray source and an RAXIS image plate detector. © 1997 Wiley-Liss Inc.     

Picosecond timescale rigid-helix and side-chain motions in deoxymyoglobin

The contribution of rigidbody motions to the atomic trajectories in a 100 ps molecular dynamics simulation of deoxymyoglobin is examined. Two typesof rigid-body motions are considered: one in which the helices are rigid units and one in which the side-chains are rigid units. Using a quaternionbased algorithm, fits of the rigid reference structures are made to each time frame of the simulation to derive trajectories of the rigid-body motions. The fitted trajectories are analysed in terms of atomic position fluctuations, mean-square displacements as a function of time, velocity autocorrelation functions and densities of states. The results are compared with the corresponding quantities calculated from the full trajectory. The relative contribution of the rigid helix motions to the helix atom dynamics depends on which quantity is examined and on which subset of atoms is chosen: rigid-helix motions contribute 86% of the rms helix backbone atomic position fluctuations, but 30% of the helix,: atom (backbone and side-chain) mean square displacements and only 1.1% of total kinetic energy. Only very low-frequency motions contribute to the rigid-helix dynamics; the rigid-body analysis allows characteristic rigid-helix vibrations to be identified and described. Treating the side-chains as rigid bodies is foundto be an excellent approximation to both their diffusive and vibrationalmean-square displacements: 96% of side-chain atom mean-square displacements originate from rigid side-Chain motions. However, the errors in theside-chain atomic positional fits are not always small. An analysis is madeof factors contributing to the positional error for different types of side-chain. © Wiley-Liss, Inc.     

Collective motions in proteins: a covariance analysis of atomic fluctuations in molecular dynamics and normal mode simulations.

A method is described for identifying collective motions in proteins from molecular dynamics trajectories or normal mode simulations. The method makes use of the covariances of atomic positional fluctuations. It is illustrated by an analysis of the bovine pancreatic trypsin inhibitor. Comparison of the covariance and cross-correlation matrices shows that the relative motions have many similar features in the different simulations. Many regions of the protein, especially regions of secondary structure, move in a correlated manner. Anharmonic effects, which are included in the molecular dynamics simulations but not in the normal analysis, are of some importance in determining the larger scale collective motions, but not the more local fluctuations. Comparisons of molecular dynamics simulations in the present and absence of solvent indicate that the environment is of significance for the long-range motions.     

Versatile method employing basic techniques of genetic engineering to study the ability of low-molecular-weight compounds to bind covalently with DNA in cell-free systems

Numerous antitumor and carcinogenic compounds and free radicals are able to modify DNA by forming covalent bonds, mainly with nucleophilic centers in nucleobases. Such a binding is usually of utmost importance for the biological outcome. The level of DNA adducts formed by a given agent is in most cases extremely low; hence their detection is very difficult. Here we propose a simple approach, exploiting techniques widely used in genetic engineering, to demonstrate and characterize the covalent modification of a DNA fragment by any low-molecular-weight compound of interest in a cell-free system. The specifically designed, several-hundred-base-pairs-long double-stranded deoxyoligonucleotide (PCR amplified)--subject to modification--includes two restriction sites: one containing only GC base pairs recognized by restriction endonuclease MspI and the other including only AT base pairs recognized by restriction endonuclease Tru1I. The covalent modification of the restriction sites abolishes their recognition and thus cleavage by the endonucleases applied. The formation of DNA adducts is induced by incubating the oligonucleotide with increasing concentrations of a studied compound, in the appropriate activating system if required. Then, the modified oligonucleotide is submitted to digestion by the above-mentioned restriction endonucleases and the DNA fragments are separated by polyacrylamide gel electrophoresis. The inhibition of cleavage indicates the occurrence of covalent modification of the restriction site(s) while simultaneously pointing at the kind of base pairs involved in DNA adduct formation. The validation of the method was performed for two DNA binding antitumor compounds, cisplatin and CC-1065, which form adducts preferentially with guanine and adenine, respectively.     

Capecitabine as a minor groove binder of DNA: molecular docking,molecular dynamics,and multi-spectroscopic studies

The interaction mechanism and binding mode of capecitabine with ctDNA was extensively investigated using docking and molecular dynamics simulations, fluorescence and circular dichroism (CD) spectroscopy, DNA thermal denaturation studies, and viscosity measurements. The possible binding mode and acting forces on the combination between capecitabine and DNA had been predicted through molecular simulation. Results indicated that capecitabine could relatively locate stably in the G-C base-pairs-rich DNA minor groove by hydrogen bond and several weaker nonbonding forces. Fluorescence spectroscopy and fluorescence lifetime measurements confirmed that the quenching was static caused by ground state complex formation. This phenomenon indicated the formation of a complex between capecitabine and ctDNA. Fluorescence data showed that the binding constants of the complex were approximately 2 × 10⁴ M^?1. Calculated thermodynamic parameters suggested that hydrogen bond was the main force during binding, which were consistent with theoretical results. Moreover, CD spectroscopy, DNA melting studies, and viscosity measurements corroborated a groove binding mode of capecitabine with ctDNA. This binding had no effect on B-DNA conformation.