Structural basis of hierarchical multiple substates of a protein. II: Monte Carlo simulation of native thermal fluctuations and energy minimization

Conformational fluctuations in a globular protein, bovine pancreatic trypsin inhibitor, in the time range between picoseconds and nanoseconds are studied by a Monte Carlo simulation method. Multiple energy minima are derived from sampled conformations by minimizing their energy. They are distributed in clusters in the conformational space. A hierarchical structure is observed in the simulated dynamics. In the time range between 10(-14) and 10(-10) seconds dynamics is well represented by a superposition of vibrational motions within an energy well with transitions among minima within each cluster. Transitions among clusters take place in the time range of nanoseconds or longer.     

The role of fluoride on the process of titanium corrosion in oral cavity

Titanium is known to possess excellent biocompatibility as a result of corrosion resistance, lack of allergenicity when compared with many other metals. Fluoride is well known as a specific and effective caries prophylactic agent and its systemic application has been recommended widely over recent decades. Nevertheless, high fluoride concentrations impair the corrosion resistance of titanium. The purpose of this article is to summarize the current data regarding the influence of fluoride on titanium corrosion process in the last 5?years. These data demonstrate noxious effects induced by high fluoride concentration as well as low pH in the oral cavity. Therefore, such conditions should be considered when prophylactic actions are administrated in patients containing titanium implants or other dental devices.     

Structural basis of hierarchical multiple substates of a protein. I: Introduction

A computer experiment of protein dynamics is carried out, which consists of two steps: (1) A Monte Carlo simulation of thermal fluctuations in the native state of a globular protein, bovine pancreatic trypsin inhibitor; and (2) a simulation of the quick freezing of fluctuating conformations into energy minima by minimization of the energy of a number of conformations sampled in the Monte Carlo simulation. From the analysis of results of the computer experiment is obtained the following picture of protein dynamics: multiple energy minima exist in the native state, and they are distributed in clusters in the conformational space. The dynamics has a hierarchical structure which has at least two levels. In the first level, dynamics is restricted within one of the clusters of minima. In the second, transitions occur among the clusters. Local parts of a protein molecule, side chains and local main chain segments, can take multiple locally stable conformations in the native state. Many minima result from combinations of these multiple local conformations. The hierarchical structure in the dynamics comes from interactions among the local parts. Protein molecules have two types of flexibility, each associated with elastic and plastic deformations, respectively.     

Hydrogen bond network of cytochrome P-450cam: a network connecting the heme group with helix K.

During investigations of the structural character of a mutant P-450cam where Glu-286 is replaced with lysine, we obtained evidence of a hydrogen bond network between helix K and the heme group via helix L of P-450cam. This mutant protein loses the ability to maintain the heme group in a proper position, possibly due to a break in the hydrogen bond network.     

Genotyping of Mycobacterium tuberculosis isolates from a low-endemic setting in northwestern state of Paraná in Southern Brazil

The purpose of this study was to provide information about the genetic diversity and prevalent genotype of Mycobacterium tuberculosis in a low-endemic setting in northwestern state of Paraná in Southern Brazil. We employed spoligotyping and mycobacterial interspersed repetitive units-variable number tandem repeat (MIRU-VNTR) techniques to genotype M. tuberculos isisolates from patients with pulmonary tuberculosis (TB). The 93 isolates analyzed by spoligotyping were divided into 36 different patterns, 30 of which were described in the SITVIT database. Latin American and Mediterranean, Haarlem and T families were responsible for 26.9%, 17.2% and 11.8% of TB cases, respectively. From the 84 isolates analyzed by MIRU-VNTR, 58 shared a unique pattern and the remaining 26 belonged to nine clusters. The MIRU loci 40, 23, 10 and 16 were the most discriminatory. A combination of MIRU-VNTR and spoligotyping resulted in 85.7% discriminatory power (Hunter-Gaston index = 0.995). Thus, combining spoligotyping and MIRU-VNTR typing proved to be most useful for epidemiological study in this low-endemic setting in Southern Brazil. The current study demonstrated that there is significant diversity in circulating strains in the city of Maringá and the surrounding regions, with no single genotype of M. tuberculosis predominating.     

Localization of hydrogen-bonds within modules in barnase

Proteins in eukaryotes are composed of structural units, each encoded by discrete exons. The protein module is one such structural unit; it has been defined as the least extended or the most compact contiguous segment in a globular domain. To elucidate roles of modules in protein evolution and folding, we examined roles of hydrogen bonds and hydrophobic cores, as related to the stability of these modules. For this purpose we studied barnase, a bacterial Rnase from Bacillus amylolique-faciens. Barnase is decomposed into at least six modules, M1–M6; the module boundaries are identified at amino acid residues 24, 52, 73, 88, and 98. Hydrogen bonds are localized mainly within each of the modules, with only a few between them, thereby indicating that their locations are designed to primarily stabilize each individual module. To obtain support for this notion, an analysis was made of hypothetical modules defined as segments starting at a center of one module and ending at the center of the following one. We found that the hydrogen bonds did not localize in each hypothetical module and that many formed between the hypothetical modules. The native conformations of modules of barnase may be specified predominantly by interactions within the modules. © 1993 Wiley-Liss, Inc.     

Revisiting gap locations in amino acid sequence alignments and a proposal for a method to improve them by introducing solvent accessibility

In comparative modeling, the quality of amino acid sequence alignment still constitutes a major bottleneck in the generation of high quality models of protein three-dimensional (3D) structures. Substantial efforts have been made to improve alignment quality by revising the substitution matrix, introducing multiple sequences, replacing dynamic programming with hidden Markov models, and incorporating 3D structure information. Improvements in the gap penalty have not been a major focus, however, following the development of the affine gap penalty and of the secondary structure dependent gap penalty. We revisited the correlation between protein 3D structure and gap location in a large protein 3D structure data set, and found that the frequency of gap locations approximated to an exponential function of the solvent accessibility of the inserted residues. The nonlinearity of the gap frequency as a function of accessibility corresponded well to the relationship between residue mutation pattern and residue accessibility. By introducing this relationship into the gap penalty calculation for pairwise alignment between template and target amino acid sequences, we were able to obtain a sequence alignment much closer to the structural alignment. The quality of the alignments was substantially improved on a pair of sequences with identity in the "twilight zone" between 20 and 40%. The relocation of gaps by our new method made a significant improvement in comparative modeling, exemplified here by the Bacillus subtilis yitF protein. The method was implemented in a computer program, ALAdeGAP (ALignment with Accessibility dependent GAp Penalty), which is available at http://cib.cf.ocha.ac.jp/target_protein/.     

Efficient Monte Carlo method for simulation of fluctuating conformations of native proteins

A powerful Monte Carlo method is described to simulate thermal conformational fluctuations in native proteins by using an empirical conformational energy function in which bond lengths and bond angles are kept fixed and only dihedral angles are independent variables. In this method, collective variables corresponding to eigenvectors of the second-derivative matrix of the energy function at its minimum point are scaled according to corresponding eigenvalues in such a way that the energy function in terms of the scaled collective variables is isotropic at the minimum point. Simulation is carried out with an isotropic step size in the space of these scaled collective variables. This simulation method is applied to a small protein, bovine pancreatic trypsin inhibitor (BPTI), and its model harmonic system defined by a quadratic energy function with the same second-derivative matrix as that of BPTI at its minimum point. Efficiency of the simulation method with an isotropic step size in the space of the scaled collective variables is found to be about 500–50 times greater than the conventional method with with an isotropic step in the space of the usual nonscaled variables. One step of this new method generates conformational changes that occur in the real-time range of 0.05 ps. In a record of 5 × 10⁵ step simulation, the BPTI molecule is observed to migrate beyond a single minimum-energy region.