Search PubMed⌕ Search

Biomedical subjects

Brajesh K Rai

Publications and source records attributed to Brajesh K Rai.

6 recordsLinked to original sources

MMM: a sequence-to-structure alignment protocol.

MOTIVATION: Accurate alignment of a target sequence to a template structure continues to be a bottleneck in producing good quality comparative protein structure models. RESULTS: Multiple Mapping Method (MMM) is a comparative protein structure modeling server with an emphasis on a novel alignment optimization protocol. MMM takes inputs from five profile-to-profile based alignment methods. The alternatively aligned regions from the input alignment set are combined according to their fit in the structural environment of the template structure. The resulting, optimally spliced MMM alignment is used as input to an automated comparative modeling module to produce a full atom model. AVAILABILITY: The MMM server is freely accessible at http://www.fiserlab.org/servers/mmm

Algorithms↗

Multiple mapping method: a novel approach to the sequence-to-structure alignment problem in comparative protein structure modeling.

A major bottleneck in comparative protein structure modeling is the quality of input alignment between the target sequence and the template structure. A number of alignment methods are available, but none of these techniques produce consistently good solutions for all cases. Alignments produced by alternative methods may be superior in certain segments but inferior in others when compared to each other; therefore, an accurate solution often requires an optimal combination of them. To address this problem, we have developed a new approach, Multiple Mapping Method (MMM). The algorithm first identifies the alternatively aligned regions from a set of input alignments. These alternatively aligned segments are scored using a composite scoring function, which determines their fitness within the structural environment of the template. The best scoring regions from a set of alternative segments are combined with the core part of the alignments to produce the final MMM alignment. The algorithm was tested on a dataset of 1400 protein pairs using 11 combinations of two to four alignment methods. In all cases MMM showed statistically significant improvement by reducing alignment errors in the range of 3 to 17%. MMM also compared favorably over two alignment meta-servers. The algorithm is computationally efficient; therefore, it is a suitable tool for genome scale modeling studies.

Algorithms↗

Single-atom test of all-atom empirical potentials: Fe in myoglobin.

The measured Fe vibrational density of states in deoxy-myoglobin, obtained from nuclear resonance vibrational spectroscopy, is compared to results from a normal-mode analysis using an all-atom empirical potential. Substantial disagreement reveals that for this one atom, the empirical potential does not accurately describe the actual forces. A Green function technique is developed to calculate the iron vibrational spectrum of deoxy-myoglobin by coupling the independently calculated heme and globin normal modes; nonbonded interactions between the heme molecule and the protein are essential for a good fit to the measurements. A projection of the eigenvectors from this potential onto the displacements induced by binding of CO demonstrates that normal modes over a broad range centered around 50-150 cm(-1) may drive the ligand-induced structural changes.

Binding Sites↗

Direct determination of the complete set of iron normal modes in a porphyrin-imidazole model for carbonmonoxy-heme proteins: [Fe(TPP)(CO)(1-MeIm)].

Detailed Fe vibrational spectra have been obtained for the heme model complex [Fe(TPP)(CO)(1-MeIm)] using a new, highly selective and quantitative technique, Nuclear Resonance Vibrational Spectroscopy (NRVS). This spectroscopy measures the complete vibrational density of states for iron atoms, from which normal modes can be calculated via refinement of the force constants. These data and mode assignments can reveal previously undetected vibrations and are useful for validating predictions based on optical spectroscopies and density functional theory, for example. Vibrational modes of the iron porphyrin-imidazole compound [Fe(TPP)(CO)(1-MeIm)] have been determined by refining normal mode calculations to NRVS data obtained at an X-ray synchrotron source. Iron dynamics of this compound, which serves as a useful model for the active site in the six-coordinate heme protein, carbonmonoxy-myoglobin, are discussed in relation to recently determined dynamics of a five-coordinate deoxy-myoglobin model, [Fe(TPP)(2-MeHIm)]. For the first time in a six-coordinate heme system, the iron-imidazole stretch mode has been observed, at 226 cm(-)(1). The heme in-plane modes with large contributions from the nu(42), nu(49), nu(50), and nu(53) modes of the core porphyrin are identified. In general, the iron modes can be attributed to coupling with the porphyrin core, the CO ligand, the imidazole ring, and/or the phenyl rings. Other significant findings are the observation that the porphyrin ring peripheral substituents are strongly coupled to the iron doming mode and that the Fe-C-O tilting and bending modes are related by a negative interaction force constant.

Biomimetic Materials↗

Iron normal mode dynamics in a porphyrin-imidazole model for deoxyheme proteins.

Iron vibrational modes of a deoxyheme protein model (2-methylimidazole)(tetraphenylporphinato)iron(II), [Fe(TPP)(2-MeImH)], have been studied by refining normal mode calculations to nuclear resonance vibrational spectroscopy (NRVS) data. The NRVS measurements give quantitative frequencies and iron amplitudes of all modes with significant Fe vibrational motion. Modes with in-plane displacement of iron are distinguished from those involving out-of-plane motion by measurements on oriented single-crystal samples. Normal modes having large overlaps with in-plane nu(42), nu(50), and nu(53) modes of the porphyrin core are identified, as well as several modes with large iron-imidazole stretch components. An out-of-plane mode at 78 cm(-1) shows significant doming of the porphyrin core, but the largest Fe doming motion arises from the coupling of phenyls and imidazole at 25 cm(-1).

Biophysical Phenomena↗

Iron normal mode dynamics in (nitrosyl)iron(II)tetraphenylporphyrin from X-ray nuclear resonance data.

The complete iron atom vibrational spectrum has been obtained by refinement of normal mode calculations to nuclear inelastic x-ray absorption data from (nitrosyl)iron(II)tetraphenylporphyrin, FeTPP(NO), a useful model for heme dynamics in myoglobin and other heme proteins. Nuclear resonance vibrational spectroscopy (NRVS) provides a direct measurement of the frequency and iron amplitude for all normal modes involving significant displacement of (57)Fe. The NRVS measurements on isotopically enriched single crystals permit determination of heme in-plane and out-of-plane modes. Excellent agreement between the calculated and experimental values of frequency and iron amplitude for each mode is achieved by a force-field refinement. Significantly, we find that the presence of the phenyl groups and the NO ligand leads to substantial mixing of the porphyrin core modes. This first picture of the entire iron vibrational density of states for a porphyrin compound provides an improved model for the role of iron atom dynamics in the biological functioning of heme proteins.

Biophysical Phenomena↗