Neural networks for protein structure prediction.
Explore the source record for details and available documents.
Databases, Factual↗
Biomedical subjects
M Karplus
Publications and source records attributed to M Karplus.
Dynamical search for bis-penicillamine enkephalin conformations.
Quenched molecular dynamics is used as a conformational search technique for the constrained cyclic analog [D-Pen2,D-Pen5]enkephalin (DPDPE) in a continuum solvent. The results show a Gaussianlike distribution of conformations as a function of energy, unlike the distributions found for simple liquids which have sharp bands for different crystal forms and broad glasslike states are found. The lowest energy conformers have structural features in common with those obtained from constrained searches based on energy minimization. (Hruby, V. J., L-.F. Kao, B. M. Pettitt, and M. Karplus. 1988. J. Am. Chem. Soc. 110:3351-3359). Many of the low energy configurations are amphiphilic with the carbonyl groups on one surface and the hydrophobic groups on the other. This supports the conclusions from the previous modeling study, which yielded amphiphilic structures as the most probable conformations of DPDPE when NOE data were included.
Calorimetry↗
Ion transport in a model gramicidin channel. Structure and thermodynamics.
The potential of mean force for Na+ and K+ ions as a function of position in the interior of a periodic poly(L,D)-alanine model for the gramicidin beta-helix is calculated with a detailed atomic model and realistic interactions. The calculated free energy barriers are 4.5 kcal/mol for Na+ and 1.0 kcal/mol for K+. A decomposition of the free energy demonstrates that the water molecules make a significant contribution to the free energy of activation. There is an increase in entropy at the transition state associated with greater fluctuations. Analysis reveals that the free energy profile of ions in the periodic channel is controlled not by the large interaction energy involving the ion but rather by the weaker water-water, water-peptide and peptide-peptide hydrogen bond interactions. The interior of the channel retains much of the solvation properties of a liquid in its interactions with the cations. Of particular importance is the flexibility of the helix, which permits it to respond to the presence of an ion in a fluidlike manner. The distortion of the helix is local (limited to a few carbonyls) because the structure is too flexible to transmit a perturbation to large distances. The plasticity of the structure (i.e., the property to deform without generating a large energy stress) appears to be an essential factor in the transport of ions, suggesting that a rigid helix model would be inappropriate.
Amino Acid Sequence↗
Pseudooutbreak of Candida guilliermondii fungemia in a neonatal intensive care unit.
During a 3-week period multiple blood cultures obtained from 14 Neonatal Intensive Care Unit infants and 3 Newborn Unit babies grew Candida guilliermondii, a yeast rarely associated with infections in humans. At the time of detection of positive cultures, most infants had been hospitalized for days or weeks for serious perinatal conditions and treated with antibiotics and intravenous hyperalimentation. Two critically ill premature infants from whom the yeast was isolated were given amphotericin B. In 7 other infants, however, yeasts were recovered on the day of birth, raising the question of pseudofungemia. Exhaustive interrogation on the blood culture practices revealed that when drawing blood for a culture from small infants, "butterfly" needles were often flushed with a diluted heparin solution to prevent blood clotting. Culture of a single lot of diluted heparin vials, prepared at the hospital pharmacy and distributed to the Neonatal Intensive Care Unit and Newborn Unit shortly before the onset of the epidemic, grew between 10,000 and 15,000 colony-forming units of Candida guilliermondii/ml. Removal of contaminated heparin vials and discontinuation of heparinization of needles used for blood cultures resulted in cessation of the epidemic. The present outbreak illustrates the difficulties in recognizing pseudoinfections in sick premature infants and the importance of intensive investigation and intervention during such an outbreak.
Blood↗
Modeling of globular proteins. A distance-based data search procedure for the construction of insertion/deletion regions and Pro----non-Pro mutations.
A distance-based database search scheme is proposed for modeling Pro----in non-Pro and insertion/deletion regions of homologous globular proteins up to six residues in length. In the first step, geometric descriptors, the number of residues involved and target distances corresponding to the separation of C alpha atom positions adjacent to the "missing" segment, are chosen. In the second step, a database of high-resolution X-ray structures is scanned for segments with similar descriptors and selected segments are binned according to conformational type. In the third and fourth steps, the selected conformations are docked into the protein, and geometric and energetic criteria are used to determine their viability as segment models. The fifth step consists of an interaction scheme in which the geometric descriptors are redefined. This compensates for the use of a limited database and/or for the use of a poor original protein model adjacent to the missing segment. The procedure has been tested on Pro----non-Pro mutations in the homologous proteins penicillopepsin and endothiapepsin, and on the insertion/deletion regions of the homologs penicillopepsin and endothiapepsin, trypsin and gamma-chymotrypsin and hen and human lysozyme. The test cases represent a wide variety of secondary structural elements (helix, sheet, turn and coil) and insertion/deletion lengths (0 to 4 residues). It is shown that 79% of the test cases are accurately modeled (within 0.54 A root-mean-square (r.m.s.) deviation for main-chain atoms) using the proposed scheme. Failure of the scheme (main-chain atom r.m.s. deviations greater than 1.29 A) in 21% of the cases appears to be related to the presence of infrequently observed conformations or locally unique folds of the target proteins with respect to the database (18% of the test cases); the remaining 3% are unexplained. Geometric and energetic criteria are able to discriminate between trial conformations that correspond to the X-ray structures and those that are different in 97% of the conformations generated by the distance-weighted database search scheme. The scheme is shown to be relatively insensitive to uncertainty in the template co-ordinates, since the geometric descriptors were taken from the homologous protein (r.m.s. deviations in the position of descriptors range from 0.18 to 1.35 A for the accurately modeled test cases). It is demonstrated that the scheme can be used to correct local sequence misalignments.
Algorithms↗
pKa's of ionizable groups in proteins: atomic detail from a continuum electrostatic model.
A macroscopic electrostatic model is used to calculate the pKa values of the titratable groups in lysozyme. The model makes use of detailed structural information and treats solvation self-energies and interactions arising from permanent partial charges and titratable charges. Both the tetragonal and triclinic crystal structures are analyzed. Half of the experimentally observed pKa shifts (11 out of 21) are well reproduced by calculations for both structures; this includes the unusually high pKa of Glu 35 in the active site. For more than half the pKa's (13 out of 21), there is a large difference (1-3.3 pK units) between the results from the two structures. Many of these correspond to the titrating groups for which the calculations are in error. Since for an ionic strength of 0.1 M the Debye screening between titratable groups leads to a very high effective dielectric constant (the average value for all pairs of titrating groups is approximately 900), near-neighbor interactions dominate the pKa perturbations. Thus, the pKa values are very sensitive to the details of the local protein conformation, and it is likely that side-chain mobility has an important role in determining the observed pKa shifts.
Computer Simulation↗
Molecular dynamics simulations in biology.
Molecular dynamics--the science of simulating the motions of a system of particles--applied to biological macromolecules gives the fluctuations in the relative positions of the atoms in a protein or in DNA as a function of time. Knowledge of these motions provides insights into biological phenomena such as the role of flexibility in ligand binding and the rapid solvation of the electron transfer state in photosynthesis. Molecular dynamics is also being used to determine protein structures from NMR, to refine protein X-ray crystal structures faster from poorer starting models, and to calculate the free energy changes resulting from mutations in proteins.
Chemical Phenomena↗
Anatomy of a conformational change: hinged "lid" motion of the triosephosphate isomerase loop.
Triosephosphate isomerase (TIM) is used as a model system for the study of how a localized conformational change in a protein structure is produced and related to enzyme reactivity. An 11-residue loop region moves more than 7 angstroms and closes over the active site when substrate binds. The loop acts like a "lid" in that it moves rigidly and is attached by two hinges to the remainder of the protein. The nature of the motion appears to be built into the loop by conserved residues; the hinge regions, in contrast, are not conserved. Results of molecular dynamics calculations confirm the structural analysis and suggest a possible ligand-induced mechanism for loop closure.
Amino Acid Sequence↗
Temperature dependence of the structure and dynamics of myoglobin. A simulation approach.
The results of simulations of the structure and internal motions of carbonomonoxymyoglobin (MbCO) at two different temperatures (325 and 80 K) are presented and compared with experimental data. Properties calculated from the 120 ps trajectory at 325 K are used as a reference in the analysis of the motion of the protein at 80 K. Three separate 80 K molecular dynamics trajectories were calculated; they were started with different coordinate sets from the 325 K simulation and the lower temperature was achieved by scaling the velocities. The simulations yield results for the structural changes between 325 and 80 K that are in general accord with those from X-ray data. Both the experimental and calculated radii of gyration, distances from the center of mass and main-chain difference distance matrices show that there is a significant but inhomogeneous shrinkage with decreasing temperature. For the atomic fluctuations, by contrast, the calculated temperature dependence is very different from the X-ray results; i.e. the calculated root-mean-square backbone fluctuations decrease to 0.11 A at 80 K from 0.51 A at 325 K, while the fluctuations obtained from the X-ray B factors go from 0.56 A at 260 K to 0.47 A at 80 K. The smaller temperature dependence of the B factors suggests that there is significant conformational disorder in MbCO crystals at lower temperatures. This is in accord with the simulation results, which show that the protein is trapped in restricted regions of conformational space at 80 K, while at 325 K a much larger region is accessible to the protein. Analysis of the fluctuations at 325 K and 80 K shows that the room temperature flexibility of the protein is determined by the mobility of the loop regions and by side-chain torsional motions (in accord with earlier simulation results), while the low temperature fluctuations involve motion within a single well. Examination of the calculated iron atom fluctuations and comparison with Mossbauer data show good agreement. It is found that the dominant contribution to the iron motion arises from heme sliding; motion of the iron relative to the heme are much smaller.
Chemical Phenomena↗
Conformational sampling using high-temperature molecular dynamics.
High-temperature molecular dynamics as a method for conformational search was explored on the antigen combining site of McPC 603, a phosphorylcholine binding immunoglobulin. Simulations at temperatures of 500, 800, and 1500 K were run for 111.5, 101.7, and 76.3 ps, respectively. The effectiveness of the search was assessed using a variety of methods. For the shorter hypervariable loops, molecular dynamics explored an appreciable fraction of the conformational space as evidenced by a comparison to a simple theoretical model of the size of the conformational space. However, for the longer loops and the antigen combining site as a whole, the simulation times were too short for a complete search. The simulations at 500 and 800 K both generated conformations that minimized to energies 200 kcal/mole lower than the crystal structure. However, the 1500 K simulation produced higher energy structures, even after minimization; in addition, this highest temperature run had many cis-trans peptide isomerizations. This suggests that 1500 K is too high a temperature for unconstrained conformational sampling. Comparison of the results of high temperature molecular dynamics with a direct conformational search method, [R. E. Bruccoleri & M. Karplus (1987) Biopolymers 26, 137-168]. showed that the two methods did not overlap much in conformational space. Simple geometric measures of the conformational space indicated that the direct method covered more space than molecular dynamics at the lower temperature, but not at 1500 K. The results suggest that high-temperature molecular dynamics can aid in conformational searches.
Amino Acid Sequence↗
Postural development in very low birth weight and normal birth weight infants.
In this study the seven postural responses selected by Vojta to evaluate neuromotor development were applied to 68 very low birth weight (VLBW) (greater than 1500 g) infants and to 28 healthy infants of normal birth weight (less than 2500 g). Of the 68 VLBW infants, 41 were small for gestational age and 27 appropriate for gestational age. All infants were examined between 37 and 40 weeks postmenstrual age. They were all later assessed on the Griffiths Mental Developmental Scale at 12 and 18 months. There were significant differences in postural reactions between the two groups which confirmed the lower tone and greater extension previously described in very low birth weight infants. An important finding in the study was that poor head and trunk righting noted at four months corrected age in very low birth weight infants, was associated with less developed locomotion at 12 and 18 months as assessed by the Griffiths Mental Developmental Scale. Thus, a delay in maturation in very low birth weight infants which was apparent from the assessment of postural responses in early infancy was still identifiable on the locomotor subscales at 12 and 18 months. Five of Vojta's responses were shown to be useful as part of the neurological assessment of high risk infants.
Apgar Score↗
Free energy of sickling: A simulation analysis.
Molecular dynamics simulations were performed to calculate the difference between the dimerization free energies of normal human deoxyhemoglobin (HbA) and the mutant sickle-cell deoxyhemoglobin HbS (Glu-beta 6----Val) for one of the lateral contacts in the HbS x-ray structure. The simulations yield a value of--15 kcal/mol. Although there is no quantitative experimental value for comparison, this is in qualitative agreement with the experimental result that HbS self-assembles into multistranded fibers that are responsible for erythrocyte sickling, while HbA does not. The free-energy difference was decomposed into enthalpic and entropic terms, both of which are significant, and the contributions of individual protein residues and of the solvent were examined. Electrostatic effects play the dominant role in favoring dimerization of HbS compared with HbA; van der Waals interactions make a negligible contribution to the difference. Both differential solvation and protein-protein interactions are important. Interactions within the donor tetramer (i.e., that containing the Glu-beta 6 mutation site), as well as those with the acceptor tetramer, contribute to the preferential free energy of dimerization of HbS.
Calorimetry↗
Dynamics of myoglobin: comparison of simulation results with neutron scattering spectra.
Molecular dynamics simulations are used to calculate the incoherent neutron scattering spectra of myoglobin between 80 K and 325 K and compared with experimental data. There is good agreement over the entire temperature range for the elastic, quasi-elastic, and inelastic components of the scattering. This provides support for the accuracy of the simulations of the internal motions that make the dominant contributions to the atomic displacements on a time scale of 0.3-100 ps (100-0.3 cm-1). Analysis of the simulations shows that at low temperatures a harmonic description of the molecule is appropriate and that the molecule is trapped in localized regions of conformational space. At higher temperatures the scattering arises from a combination of vibrations within wells (substates) and transitions between them; the latter contribute to the quasi-elastic scattering.
Computer Simulation↗
The interpretation of site-directed mutagenesis experiments by linear free energy relations.
Fersht and co-workers have applied a linear free energy relation (Brønsted equation) to analyze site-directed mutagenesis experiments involving the enzyme tyrosyl-tRNA synthetase and have suggested that the Brønsted exponent is linearly correlated with the value of the reaction coordinate at the transition state. We point out that when the mutants differ solely through the formation or deletion of a hydrogen bond away from the reaction center, a linear free energy relation is expected only in limiting cases for which the Brønsted relation exponent is 0, 1 or infinity. The results may be correlated with a conformational coordinate but not with the development of the reaction coordinate per se.
Chemical Phenomena↗
Construction of side-chains in homology modelling. Application to the C-terminal lobe of rhizopuspepsin.
A detailed and rule-based side-chain modelling procedure for globular proteins is presented. It uses the conformational information contained in a homologous (template) structure as a starting point and includes recipes for atom placement and for checking and improving the atomic positions. The scheme does not rely on intuitive judgements or visual examination of the model during construction or refinement. It comprises four stages; the first three are relatively simple and the fourth is more complex. In the first stage, initial conformations for as many atoms as possible are transferred from the template structure based on the application of trends reported previously. Second, these trends are used to correct poor van der Waals overlaps. Third, the remaining side-chains atoms (those for which no information is contained in the template) are placed by evaluating their rigid rotation, van der Waals surfaces. The fourth stage consists of a hierarchial series of conformational checks. They involve the evaluation of individual residue energies in the absence and presence of the rest of the protein relative to statistical trends observed in the template structure, the comparison of hydrogen-bonding patterns and side-chain accessibilities in the model and template and brief energy minimization followed by an evaluation of the rigid rotation potential energy surfaces of each side-chain. The checks pinpoint "incorrectly" modelled side-chains, suggest conformational changes and provide a means for determining the portions of the model that are likely to be correct and those likely to be in error. The procedure developed in the paper is tested by modelling the side-chains of the C-terminal lobe of the aspartyl proteinase rhizopuspepsin, using the rhizopuspepsin backbone and the homologous protein, penicillopepsin, as a template for the side-chains. The resultant model was compared to the high-resolution X-ray structure of rhizopuspepsin. Using penicillopepsin data only (stage I), 58% of the chi 1 dihedrals and 44% of the chi 2 dihedrals were modelled correctly. Once poor van der Waals overlaps had been corrected and all of the atoms had been placed (stages II and III), 86% of the chi 1 dihedrals and 75% of the chi 2 dihedrals were correct. After the refinement had been completed (stage IV), 92% of the chi 1 dihedrals and 81% of the chi 2 dihedrals were correctly positioned.(ABSTRACT TRUNCATED AT 400 WORDS)
Amino Acid Sequence↗
Two-dimensional NMR and photo-CIDNP studies of the insulin monomer: assignment of aromatic resonances with application to protein folding, structure, and dynamics.
The aromatic 1H NMR resonances of the insulin monomer are assigned at 500 MHz by comparative studies of chemically modified and genetically altered variants, including a mutant insulin (PheB25----Leu) associated with diabetes mellitus. The two histidines, three phenylalanines, and four tyrosines are observed to be in distinct local environments; their assignment provides sensitive markers for studies of tertiary structure, protein dynamics, and protein folding. The environments of the tyrosine residues have also been investigated by photochemically induced dynamic nuclear polarization (photo-CIDNP) and analyzed in relation to packing constraints in the crystal structures of insulin. Dimerization involving specific B-chain interactions is observed with increasing protein concentration and is shown to depend on temperature, pH, and solvent composition. In the monomer large variations are observed in the line widths of amide resonances, suggesting intermediate exchange among conformational substates; such substates may relate to conformational changes observed in different crystal states and proposed to occur in the hormone-receptor complex. Additional evidence for multiple conformations in solution is provided by comparative studies of an insulin analogue containing a peptide bond between residues B29 and A1 (mini-proinsulin). This analogue forms dimers and higher-order oligomers under conditions in which native insulin is monomeric, suggesting that the B29-A1 peptide bond stabilizes a conformational substate favorable for dimerization. Such stabilization is not observed in corresponding studies of native proinsulin, in which a 35-residue connecting peptide joins residues B30 and A1; this extended tether is presumably too flexible to constrain the conformation of the B-chain. The differences between proinsulin and mini-proinsulin suggest a structural mechanism for the observation that the fully reduced B29-A1 analogue folds more efficiently than proinsulin to form the correct pattern of disulfide bonds. These results are discussed in relation to molecular mechanics calculations of insulin based on the available crystal structures.
Circular Dichroism↗
Solvent effects on protein motion and protein effects on solvent motion. Dynamics of the active site region of lysozyme.
The stochastic boundary molecular dynamics methodology is applied to the active site of the enzyme lysozyme. A comparison is made of in vacuo dynamics results from the stochastic boundary method and a full conventional molecular dynamics simulation of lysozyme. Excellent agreement between the two approaches is obtained. The influence of solvent on the residues in the active site region is explored and it is shown that both the structure and dynamics are affected. Of particular importance for the structure of the protein is the solvation of polar residues and the stabilization of like-charged ion pairs. The magnitude of the fluctuations is only slightly altered by the solvent; the overall increase in the root-mean-square fluctuations, relative to the vacuum run, is 11%. The solvent effect on dynamical properties is found not to be simply related to the solvent viscosity. Both the solvent exposure and dynamic aspects of protein-solvent interactions, including the relative time scales of the motions, are shown to play a role. The effects of the protein on solvent dynamics and structure are also observed to be significant. The solvent molecules around atoms in charged, polar and apolar side-chains show markedly different diffusion coefficients as well as exhibiting different solvation structures. One key example is the water around apolar groups, which is much less mobile than bulk water, or water solvating polar groups.
Animals↗
Hidden thermodynamics of mutant proteins: a molecular dynamics analysis.
A molecular dynamics simulation method is used to determine the contributions of individual amino acid residues and solvent molecules to free energy changes in proteins. Its application to the hemoglobin interface mutant Asp G1(99) beta----Ala shows that some of the contributions to the difference in the free energy of cooperativity are as large as 60 kilocalories (kcal) per mole. Since the overall free energy change is only -5.5 kcal/mole (versus the experimental value of -3.4 kcal/mole), essential elements of the thermodynamics are hidden in the measured results. By exposing the individual contributions, the free energy simulation provides new insights into the origin of thermodynamic changes in mutant proteins and demonstrates the role of effects beyond those usually considered in structural analyses.
Alanine↗