Protein structure prediction based on statistical potential.
Explore the source record for details and available documents.
Biomedical subjects
Publications and source records attributed to R Rein.
Explore the source record for details and available documents.
This paper suggests that for medical knowledge to be effective, it must be reliably translated into a series of synchronized activities, both clinical and operational. The goal of medical information science must be not only to identify the best course of care, but also through the use of systems for the ordering, scheduling, and deployment of staff and resources to assure that the planned course of care is reliably and efficiently delivered in actual practice settings. We present a brief history of several innovations in patient care coordination and logistics at Mayo Clinic to provide background for the main concepts presented. The principle of patient care coordination underlying this history of innovation is shown to have evolved into a new need for information systems capable of proactively effecting the delivery of care through patient, staff and resource deployment and monitoring.
The effect of G.T mispair incorporation into a double-helical environment was examined by molecular dynamics simulation. The 60-ps simulations performed on the two hexanucleotide duplexes d (G3C3)2 and d(G3TC2)2 included 10 Na+ counterions and first hydration shell waters. The resulting backbone torsional angle trajectories were analyzed to select time spans representative of conformational domains. The average backbone angles and helical parameters of the last time span for both duplexes are reported. During the simulation the hexamers retained B-type DNA structures that differed from typical A- or B-DNA forms. The overall helical structures for the two duplexes are vary similar. The presence of G.T mispairs did not alter the overall helical structure of the oligonucleotide duplex. Large propeller twist and buckle angles were obtained for both duplexes. The purine/pyrimidine crossover step showed a large decrease in propeller twist in the normal duplex but not in the mismatch duplex. Upon the formation of wobble mispairs in the mismatched duplex, the guanines moved into the minor groove and the thymines moved into the major groove. This helped prevent purine/purine clash and created a deformation in the relative orientation of the glycosidic bonds. It also exposed the free O4 of the thymines in the major groove and N2 of the guanines in the minor groove to interactions with solvent and counterions. These factors seemed to contribute to the apparently higher rigidity of the mismatched duplex during the simulation.
The origin of torsional potentials in H3CSSCH3, H3CSSH, and HOOH and the anisotropy of the local charge distribution has been analyzed in terms of atomic multipoles calculated from the ab initio LCAO-MO-SCF wave function in the 6-31G* basis set. The results indicate that for longer -S-S-bonds the major contribution to these torsional barriers are electrostatic interactions of the atomic multipoles located on two atoms forming the rotated bond. This finding demonstrates the important role of electrostatic 1-2 interatomic interactions, usually neglected in conformational studies. It also opens the possibility to derive directly from accurate ab initio wave functions a simple nonempirical torsional potential involving atomic multipoles of two bonded atoms defining the torsional angle. For shorter -O-O- bonds, use of more precise models and inclusion of 1-3 interactions seems to be necessary.
A systematic analysis is presented of the algorithm for converting a virtual-bond chain, defined by the coordinates of the alpha-carbons of a given protein, into a complete polypeptide backbone. An alternative algorithm, based upon the same set of geometric parameters used in the Purisima-Scheraga algorithm but with a different "linkage map" of the algorithmic procedures, is proposed. The global virtual-bond chain geometric constraints are more easily separable from the loal peptide geometric and energetic constraints derived from, for example, the Ramachandran criterion, within the framework of this approach.
The changes in the catalytic activity resulting from amino acid substitutions in the active site region have been theoretically modeled for tyrosyl tRNA synthetase (Tyr-RS). The catalytic activity was calculated as the differential stabilization of the transition state using electrostatic approximation. The results indicate that charged residues His45, His48, Asp78, Asp176, Asp194, Lys225, Lys230, Lys233, Arg265, and Lys268 play essential roles in catalysis of aminoacyl adenylate formation in Tyr-RS, which is in general agreement with previously known experimental data for residues 45, 48, 194, 230, and 233. These catalytic residues have also been used to search for sequence homology patterns among class I aminoacyl RSs of which HIGH and KMSKS conserved sequence motifs are well known. His45 and His48 belong to the HIGH signature sequence of class I aminoacyl tRNA synthetases (aRSs), whereas Arg265 and Lys268 can constitute a part of the KMSKS charge pattern. Lys225, Lys230, and Lys233 may be part of the conservative substitution pattern [HKR]-X(4)-[HKR]-X(2)-[HKR], and Asp194 is part of the new GSDQ motif. This demonstrates that the three dimensional charge distribution near the active site is an essential feature of the catalytic activity of aRS and that the theoretical technique used in this work can be utilized in searches for the catalytically important residues that may provide a clue for a charge residue pattern conserved in evolution. The appearance of patterns I-IV in Arg-, Gln-, Met-, Ile-, Leu-, Trp-, Val-, Glu-, Cys-, and Tyr-RS indicates that all these enzymes could have the same ancestor.
The phospholipid and Ca2+ dependent protein kinase (PKC) plays an essential role in a variety of cellular events. Inhibition of PKC was shown to arrest growth in tumor cell cultures making it a target for possible antitumor therapy. Calphostins are potent inhibitors of PKC with high affinity for the enzyme regulatory site. Structural characteristics of calphostins, which confer the inhibitory activity, are investigated by comparing their optimized structures with the existing models for PKC activation. The resulting model of inhibitory activity assumes interaction with two out of the three electrostatic interaction sites postulated for activators. The model shows two sites of hydrophobic interaction and enables the inhibitory activity of gossypol to be accounted for.
In an effort to understand the conformational and structural changes in DNA brought about by thymine photodimer, computer modeling and molecular mechanics energy calculations were performed on DNA hexamer and dodecamer duplexes containing a cis-syn photodimer. The conformation of the crystal structure of the cyanoethyl phosphate ester of the thymine dimer (Hruska et al., Biopolymers 25, 1399-1417 (1986)) was used in modeling the photodimer portion. Various starting conformations were used in the modeling procedure and the structures were minimized both retaining and later relaxing the crystallographic geometry of the cyclobutane ring. The results indicate that most of the deformation is restricted to the thymine dimer region, and that the conformational changes decrease rapidly on either side of the region containing the photodimer. The structural changes brought about by the introduction of the photodimer can be accommodated within six base paired duplex without significant bend in the DNA. More conformational changes are observed on the 5'-side of the photodimer than on the 3'-side. The conformational features, such as backbone torsion angles and sugar puckers, of the energy minimized structures are discussed in the context of the solution structures determined by NMR on a series of oligomers containing photodimers (Rycyna et al., Biochemistry 27, 3152-3163 (1988)).
Two methods of side-chain placement on a modeled protein have been examined. Two molecular models of calmodulin were constructed that differ in the treatment of side chains prior to optimization of the molecule. A virtual bond analysis program developed by Purisima and Scheraga was used to determine the backbone conformation based on 2.2 angstroms resolution C alpha coordinates for the molecules. In the first model, side chains were initially constructed in an extended conformation. In the second model, a conformational grid search technique was employed. Calcium ions were treated explicitly during energy optimization using CHARMM. The models are compared to a recently published refined crystal structure of calmodulin. The results indicate that the initial choices for side-chains, but also significant effects on the main-chain conformation and supersecondary structure. The conformational differences are discussed. Analysis of these and other methods makes possible the formulation of a methodology for more appropriate side-chain placement in modeled proteins.
Explore the source record for details and available documents.
About 20 years subsequent to a thyroid operation, syncopes occurred in a 42-year-old male patient. A prolonged QT interval was seen in the ECG. Both calcium and magnesium concentrations were decreased (calcium 1.1 mmol/l, magnesium 0.6 mmol/l). There was hyperphosphataemia of 2.8 mmol/l with a parathormone level at the lower limit of normal. After substituting vitamin D and calcium the patient was discharged symptom-free.
A molecular dynamics (MD) study of (dG)6.(dC)6 including counter ions and 292 water molecules was made. The hydrogen bonding pattern and propeller twist angles for the mini-helix are reported as averages for times spanning 21-30, 31-40, 41-50, and 51-60 ps. The propeller twist angles range from 18 degrees to 38 degrees. Bifurcated and interstrand neighboring base (twisted) hydrogen bonding patterns were found.
Recent discoveries have established the fact that RNA is capable of acting as an enzyme. In this study two different types of molecular orbital calculations, INDO and ab initio, were used in an attempt to assess the structural/functional role of the Mg2+ hydrated complex in ribozyme reactions. Preliminary studies indicate that the reaction is multistep and that the Mg2+ complex exerts a stabilizing effect on the intermediate or midpoint of the reaction.
The human genomic clone pb2m13 contains a functional beta 2-microglobulin (B2m) gene, which upon transfection is readily expressed in murine fibroblasts. Here we report the nucleotide sequence of the human beta 2m gene and of a nearly full length cDNA clone. A comparison with the murine beta 2m gene reveals that exon/intron boundaries are absolutely conserved. In the protein-coding regions the similarity is 70%. As far as intron sequences of the murine beta 2m gene are available, no significant similarity between human and murine genes is observed. The transcriptional start site of the human beta 2m gene was determined by S1 mapping, and comparison with the nearly full length cDNA clone now defines the transcriptional unit of the beta 2m gene. In the 5' region of the gene strong clustering of the usually underrepresented CpG dinucleotide is found resembling a similar overrepresentation in the 5' regions of the major histocompatibility complex class I genes.
The scope of multistep modeling (MSM) is expanding by adding a least-squares minimization step in the procedure to fit backbone reconstruction consistent with a set of C-alpha coordinates. The analytical solution of Phi and Psi angles, that fits a C-alpha x-ray coordinate is used for tyr-tRNA synthetase. Phi and Psi angles for the region where the above mentioned method fails, are obtained by minimizing the difference in C-alpha distances between the computed model and the crystal structure in a least-squares sense. We present a stepwise application of this part of MSM to the determination of the complete backbone geometry of the 321 N terminal residues of tyrosine tRNA synthetase to a root mean square deviation of 0.47 angstroms from the crystallographic C-alpha coordinates.
The molecular mechanisms underlying the primitive translational apparatus have been studied in light of present day protein biosynthesis. Using the structural information available from the contemporary system as a key to its function, both the structural necessities for an early adaptor and the multipoint recognition properties of such adaptors have been investigated. This was done by first critically examining the potential feasibility of right- and left-handed hairpin adaptor models. Second, a molecular model of the contemporary transpeptidation complex has been constructed in order to ascertain the structural requirements of the adaptor molecule needed for peptidyl transfer. Third, a model of the tRNATyr-tyrosyl tRNA synthetase complex including the positioning of the disordered region is proposed. This model is used to illustrate those required recognition properties of aminoacyl synthetase which lead to a perspective on the structure of the ancestor synthetase.
Molecular models are proposed to describe the relative arrangement of aminoacyl and peptidyl tRNAs when bound to their respective A and P sites on the ribosome. The crystallographically determined structures of tRNAasp and tRNAphe have served as the models for these bound structures, while the imposed steric constraints for the model complexes were based on the results of published experimental data. The constructed models satisfy the stereochemical requirements needed for codon-anticodon interaction and for peptide bond formation. In this paper, the results of the complex containing tRNAphe as the A and P site bound transfer RNAs, is compared to a similarly constructed model which uses tRNAasp as the ribosome-bound transfer RNAs. The models have the following three major features: 1) the aminoacyl and peptidyl transfer RNAs assume an angle of approximately 45 degrees relative to each other; 2) in providing the proper stereochemistry for peptide bond condensation, a significant kink must be present in the messenger RNA between the A site and P site codons; and 3) a comparison of the two model complexes indicates that structural variations between the tRNAs or any allosteric transitions of the transfer RNAs associated with codon-anticodon recognition may be accommodated in the model by way of freedom of rotation about the phosphate backbone bonds in the mRNA between consecutive codons.
Adaptor properties of linear hairpin helices have been examined. The analysis suggests that neither right nor left handed hairpin helices can simultaneously read a comma free messenger and align aminoacyl residues for peptide condensation. Comparison of these studies with the model of the present day peptidyl transfer intermediate suggests that the "L" shaped folding of the present day tRNAs may be a prerequisite for adaptor function. Therefore, the three-dimensional organization of the ancestral adaptor molecule must have had structural features similar to its present day counterpart.