PubMed Health⌕ Search

Biomedical subjects

M Karplus

Publications and source records attributed to M Karplus.

At least 37 records · Page 2Linked to original sources
PubMed1991

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.

Algorithms↗
PubMed1991

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↗
PubMed1991

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↗
PubMed1991

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↗
PubMed1990

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↗
PubMed1990

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↗
PubMed1990

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↗
PubMed1990

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↗
PubMed1990

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↗
PubMed1990

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↗
PubMed1990

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↗
PubMed1990

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↗
PubMed1990

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↗
PubMed1990

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↗
PubMed1989

X-ray refinement of protein structures by simulated annealing: test of the method on myohemerythrin.

The recently developed method of structure factor refinement by molecular dynamics with simulated annealing [Brünger, Kuriyan & Karplus (1987). Science, 235, 458-460] is tested on the 118 residue protein myohemerythrin. A highly refined structure for this protein at 1.3/1.7 A resolution has recently been published [Sheriff, Hendrickson & Smith (1987). J. Mol. Biol. 197, 273-296]. This is compared with the results of simulated annealing refinement (with no manual intervention) starting from an earlier model for the protein from a stage in the refinement when conventional least-squares methods could not improve the structure. Simulated annealing reduces the R factor at 2.5 A from 39 to 31%, with uniform temperature factors and no solvent molecules and with similar stereochemistry; the comparable value for the manually refined structure is 27.9%. Errors in backbone and sidechain positions up to about 3 A are corrected by the method. The error in backbone positions for roughly 85% of the initial structure is within this range, and in these regions the r.m.s. backbone error is reduced from 1.1 to 0.4 A. For the rest of the structure, including a region which was incorrectly built due to a sequence error, the procedure does not yield any improvement and manual intervention appears to be required. Nevertheless, the overall improvement in the structure results in electron density maps that are easier to interpret and permit identification of the errors in the structure. The general utility of the simulated annealing methodology in X-ray refinement is discussed.

Evaluation Studies as Topic↗
PubMed1988

Inelastic neutron scattering analysis of picosecond internal protein dynamics. Comparison of harmonic theory with experiment.

The experimental inelastic neutron scattering spectrum of a protein, the bovine pancreatic trypsin inhibitor (BPTI), in a powder sample is presented together with the generalized density of states, G(omega), as a function of the frequency, omega, derived from the scattering data. The experimental results are compared with calculations from two different normal mode analyses of BPTI. One of these, based on an improved model, gives a calculated spectrum and density of states in general agreement with those obtained experimentally; the other, based on an earlier model, shows considerable disagreement. The important improvements in the newer normal mode analysis are the explicit treatment of all atoms (non-polar as well as polar hydrogens are included) and a modified truncation scheme for the long-range electrostatic interactions. The fact that the inelastic neutron scattering measurements can distinguish between the two theoretical models makes clear their utility for the analysis of protein dynamics.

Animals↗