Search PubMed⌕ Search

Biomedical subjects

M Karplus

Publications and source records attributed to M Karplus.

At least 73 records · Page 4Linked to original sources
PubMed1997

Dynamics and unfolding pathways of a hyperthermophilic and a mesophilic rubredoxin.

Molecular dynamics simulations in solution are performed for a rubredoxin from the hyperthermophilic archaeon Pyrococcus furiosus (RdPf) and one from the mesophilic organism Desulfovibrio vulgaris (RdDv). The two proteins are simulated at four temperatures: 300 K, 373 K, 473 K (two sets), and 500 K; the various simulations extended from 200 ps to 1,020 ps. At room temperature, the two proteins are stable, remain close to the crystal structure, and exhibit similar dynamic behavior; the RMS residue fluctuations are slightly smaller in the hyperthermophilic protein. An analysis of the average energy contributions in the two proteins is made; the results suggest that the intraprotein energy stabilizes RdPf relative to RdDv. At 373 K, the mesophilic protein unfolds rapidly (it begins to unfold at 300 ps), whereas the hyperthermophilic does not unfold over the simulation of 600 ps. This is in accord with the expected stability of the two proteins. At 473 K, where both proteins are expected to be unstable, unfolding behavior is observed within 200 ps and the mesophilic protein unfolds faster than the hyperthermophilic one. At 500 K, both proteins unfold; the hyperthermophilic protein does so faster than the mesophilic protein. The unfolding behavior for the two proteins is found to be very similar. Although the exact order of events differs from one trajectory to another, both proteins unfold first by opening of the loop region to expose the hydrophobic core. This is followed by unzipping of the beta-sheet. The results obtained in the simulation are discussed in terms of the factors involved in flexibility and thermostability.

Amino Acid Sequence↗
PubMed1997

Computer simulations of the OmpF porin from the outer membrane of Escherichia coli.

Molecular dynamics simulations were used to study the structure and dynamics of the Escherichia coli OmpF porin, which is composed of three identical 16-stranded beta-barrels. Simulations of the full trimer in the absence of water and the membrane led to significant contraction of the channel in the interior of each beta-barrel. With very weak harmonic constraints (0.005 kcal/mol A2/atom) applied to the main-chain C alpha atoms of the beta-barrel, the structure was stabilized without alteration of the average fluctuations. The resulting distribution of the fluctuations (small for beta-strands, large for loops and turns) is in good agreement with the x-ray B factors. Dynamic cross-correlation functions showed the importance of coupling between the loop motions and barrel flexibility. This was confirmed by the application of constraints corresponding to the observed temperature factors to the barrel C alpha atoms. With these constraints, the beta-barrel fluctuations were much smaller than the experimental values because of the intrinsic restrictions on the atomic motions, and the loop motions were reduced significantly. This result indicates that considerable care is required in introducing constraints to keep proteins close to the experimental structure during simulations, as has been done in several recent studies. Loop 3, which is thought to be important in gating the pore, undergoes a displacement that shifts it away from the x-ray structure. Analysis shows that this arises from the breakdown of a hydrogen bond network, which appears to result more from the absence of solvent that from the use of standard ionization states for the side chains of certain beta-barrel residues.

Amino Acid Sequence↗
PubMed1997

The Levinthal paradox: yesterday and today.

A change in the perception of the protein folding problem has taken place recently. The nature of the change is outlined and the reasons for it are presented. An essential element is the recognition that a bias toward the native state over much of the effective energy surface may govern the folding process. This has replaced the random search paradigm of Levinthal and suggests that there are many ways of reaching the native state in a reasonable time so that a specific pathway does not have to be postulated. The change in perception is due primarily to the application of statistical mechanical models and lattice simulations to protein folding. Examples of lattice model results on protein folding are presented. It is pointed out that the new optimism about the protein folding problem must be complemented by more detailed studies to determine the structural and energetic factors that introduce the biases which make possible the folding of real proteins.

Models, Chemical↗
PubMed1997

A comparison between molecular dynamics and X-ray results for dissociated CO in myoglobin.

The distribution of carbon monoxide after photodissociation in the myoglobin haem pocket has been investigated using molecular dynamics simulations at 300 K. The results show that both intermediates (one close to the haem iron and one further away) observed in recent low temperature X-ray studies of photodissociated CO have a high probability of occurrence, even at ambient temperatures. The fact that the O of CO is oriented toward the haem iron in the closer intermediate provides an explanation for the slow rate of CO geminate rebinding. A refinement against X-ray data generated from the molecular dynamics simulations indicates that the CO has a broader distribution in the haem pocket than is apparent from the experimental electron density. This effect is likely to be general for systems containing highly mobile groups.

Algorithms↗
PubMed1996

Genetic neural networks for quantitative structure-activity relationships: improvements and application of benzodiazepine affinity for benzodiazepine/GABAA receptors.

A novel tool, called a genetic neural network (GNN), has been developed for obtaining quantitative structure-activity relationships (QSAR) for high-dimensional data sets (J. Med. Chem. 1996, 39, 1521-1530). The GNN method uses a neural network to correlate activity with descriptors that are preselected by a genetic algorithm. To provide an extended test of the GNN method, the data on 57 benzodiazepines given by Maddalena and Johnston (MJ; J. Med. Chem. 1995, 38, 715-724) have been examined with an enhanced version of GNN, and the results are compared with the excellent QSAR of MJ. The problematic steepest descent training has been replaced by the scaled conjugate gradient algorithm. This leads to a substantial gain in performance in both robustness of prediction and speed of computation. The cross-validation GNN simulation and the subsequent run based on an unbiased and more efficient protocol led to the discovery of other 10-descriptor QSARs that are superior to the best model of MJ based on backward elimination selection and neural network training. Results from a series of GNNs with a different number of inputs showed that a neural network with fewer inputs can produce QSARs as good as or even better than those with higher dimensions. The top-ranking models from a GNN simulation using only six input descriptors are presented, and the chemical significance of the chosen descriptors is discussed. The statistical significance of these GNN QSARs is validated. The best QSARs are used to provide a graphical tool that aids the design of new drug analogues. By replacing functional groups at the 7- and 2'-positions with ones that have optimal substituent parameters, a number of new benzodiazepines with high potency are predicted.

Algorithms↗
PubMed1996

Comment on a "fluctuation and cross correlation analysis of protein motions observed in nanosecond molecular dynamics simulations".

Long range positive correlations in the internal atomic motions have been observed in molecular dynamics simulations of the bovine pancreatic trypsin inhibitor. Since the results depend on the removal of overall translation and rotation (which is not unique), some care is required. It is pointed out that original use of all of the C atoms as a frame of reference is more appropriate than the selection of a particular subset. The former comes closest to reflecting the zero rotation condition in a normal mode analysis, which also shows the existence of long range correlations.

Animals↗
PubMed1996

The folding mechanism of larger model proteins: role of native structure.

The folding mechanism of a 125-bead heteropolymer model for proteins is investigated with Monte Carlo simulations on a cubic lattice. Sequences that do and do not fold in a reasonable time are compared. The overall folding behavior is found to be more complex than that of models for smaller proteins. Folding begins with a rapid collapse followed by a slow search through the semi-compact globule for a sequence-dependent stable core with about 30 out of 176 native contacts which serves as the transition state for folding to a near-native structure. Efficient search for the core is dependent on structural features of the native state. Sequences that fold have large amounts of stable, cooperative structure that is accessible through short-range initiation sites, such as those in anti-parallel sheets connected by turns. Before folding is completed, the system can encounter a second bottleneck, involving the condensation and rearrangement of surface residues. Overly stable local structure of the surface residues slows this stage of the folding process. The relation of the results from the 125-mer model studies to the folding of real proteins is discussed.

Models, Theoretical↗
PubMed1996

Evolutionary optimization in quantitative structure-activity relationship: an application of genetic neural networks.

A new hybrid method (GNN) combining a genetic algorithm and an artificial neural network has been developed for quantitative structure-activity relationship (QSAR) studies. A suitable set of molecular descriptors are selected by a genetic algorithm. This set serves as input to a neural network, in which model-free mapping of multivariate data is performed. Multiple predictors are generated that are superior to results obtained from previous studies of the Selwood data set, which is used to test the method. The neural network technique provides a graphical description of the functional form of the descriptors that play an important role in determining drug activity. This can serve as an aid in future design of drug analogues. The effectiveness of GNN is tested by comparing its results with a benchmark obtained by exhaustive enumeration. Different fitness strategies that tune the evolution of genetic models are examined, and QSARs with higher predictiveness are found. From these results, a composite model is constructed by averaging predictions from several high-ranking models. The predictions of the resulting QSAR should be more reliable than those derived from a single predictor because it makes greater use of information and also permits error estimation. An analysis of the sets of descriptors selected by GNN shows that it is essential to have one each for the steric, electrostatic, and hydrophobic attributes of a drug candidate to obtain a satisfactory QSAR for this data set. This type of result is expected to be of general utility in designing and understanding QSAR.

Algorithms↗
PubMed1996

The importance of larger data sets for protein secondary structure prediction with neural networks.

A neural network algorithm is applied to secondary structure and structural class prediction for a database of 318 nonhomologous protein chains. Significant improvement in accuracy is obtained as compared with performance on smaller databases. A systematic study of the effects of network topology shows that, for the larger database, better results are obtained with more units in the hidden layer. In a 32-fold cross validated test, secondary structure prediction accuracy is 67.0%, relative to 62.6% obtained previously, without any evolutionary information on the sequence. Introduction of sequence profiles increases this value to 72.9%, suggesting that the two types of information are essentially independent. Tertiary structural class is predicted with 80.2% accuracy, relative to 73.9% obtained previously. The use of a larger database is facilitated by the introduction of a scaled conjugate gradient algorithm for optimizing the neural network. This algorithm is about 10-20 times as fast as the standard steepest descent algorithm.

Algorithms↗
PubMed1996

Simulations of lysozyme: internal motions and the reaction mechanism.

Simulations of lysozyme by molecular dynamics have greatly increased our understanding of this enzyme. It has been shown how the internal motions are related to the structural elements (helices, sheets and loops) of the molecule. Comparisons of the motions in the free and substrate-bound form reveal that most are similar but that there are significant differences. Comparisons of the theoretical results with X-ray and nuclear magnetic resource data show good agreement. The hinge-bending motion, which opens and closes the binding site cleft between the two domains, is correlated with substrate binding. From analysis of simulations with bound substrate, an alternative mechanism for oligoglycoside hydrolysis was proposed. It involves cleavage of an endocyclic C-O bond, instead of the exocyclic cleavage proposed in the standard mechanism. Both mechanisms have been demonstrated in solution, but it is still unclear which is prevalent in lysozyme.

Binding Sites↗
PubMed1996

Functionality map analysis of the active site cleft of human thrombin.

The Multiple Copy Simultaneous Search methodology has been used to construct functionality maps for an extended region of human thrombin, including the active site. This method allows the determination of energetically favorable positions and orientations for functional groups defined by the user on the three-dimensional surface of a protein. The positions of 10 functional group sites are compared with those of corresponding groups of four thrombin-inhibitor complexes. Many, but not all features, of known thrombin inhibitors are reproduced by the method. The results indicate that certain aspects of the binding modes of these inhibitors are not optimal. In addition, suggestions are made for improving binding by interaction with functional group sites on the thrombin surface that are not used by the thrombin inhibitors.

Amino Acid Chloromethyl Ketones↗
PubMed1996

Failure of early postnatal dexamethasone to prevent chronic lung disease in infants with respiratory distress syndrome.

OBJECTIVE: To study the effect of early postnatal dexamethasone (days 1-3) on the incidence and severity of chronic lung disease in preterm infants with respiratory distress syndrome. METHODS: A multicentre, randomised, placebo controlled, blinded study was carried out in 18 neonatal intensive care units in Israel. The primary outcome measure was survival to discharge without requirement for supplemental oxygen therapy beyond 28 days of life. The secondary outcome measures were requirement for mechanical ventilation at 3 and 7 days, duration of ventilation or oxygen therapy, need for subsequent steroids for established chronic lung disease and incidence of major morbidities. RESULTS: The study consisted of 248 infants (dexamethasone n = 132; placebo n = 116). No differences were found in the outcome variables except for a reduction in requirement for mechanical ventilation at age 3 days in treated infants (dexamethasone 44%, placebo 67%; P = 0.001). Gastrointestinal haemorrhage, hypertension, and hyperglycaemia were more common in treated infants, but no life threatening complications, such as gastrointestinal perforation, were encountered. CONCLUSIONS: These data do no support the routine use of early postnatal steroids, but may justify further study in a selected, high risk group of infants.

Birth Weight↗
PubMed1996

Fundus pigmentation in retinopathy of prematurity.

OBJECTIVE: A prospective study design was used to investigate the association between different degrees of fundus pigmentation and the indicence of retinopathy of prematurity (ROP) among very low birth weight infants in a large neonatal intensive care unit. METHODS: The study group consisted of 161 infants weighing 1500 g or less at birth and included all infants born from 1988 to 1990 who survived at least 10 weeks. Presence or lack of any acute stage ROP was determined by weekly ophthalmological examination from the age of 4 weeks. The degree of fundus pigmentation was recorded for each infant during the first examination. RESULTS: Infants with dark fundus pigmentation were found to be at half the risk of developing ROP as compared with infants having light/medium fundus pigmentation (relative risk = 0.5; 95% confidence interval = 0.2-1.1). When controlled for birth weight, gestational age, length of oxygen therapy, and ethnic group in multivariate analysis, dark pigmentation was an independent and statistically significant protective factor (odds ration = 0.09, 95% confidence interval = 0.02-0.06). None of the infants with Stage III ROP or higher had dark pigmentation. CONCLUSION: We speculate that large amounts of melanin in the retinal pigment epithelium or choroid may protect the dark-pigmented very low birth weight infant from developing ROP.

Female↗
PubMed1995

The meaning of component analysis: decomposition of the free energy in terms of specific interactions.

Free energy simulations are of particular interest for the interpretation of macroscopic data in terms of microscopic interactions. This can be done by expressing calculated free energies as a sum of components that correspond to the contributions of different energy terms or different parts of the system. Since the resulting components depend on the integration path, care is required for their use. We show that a linear coupling scheme for the alchemical creation of a chemical identity corresponds to a particularly useful path because it leads to a symmetric decoupling of the free energy components. The path dependence also provides an additional degree of freedom that can be used to study different processes. This latter point is illustrated by a reinterpretation of a recent simulation on wild-type and mutant azurin by Mark and van Gunsteren.

Alchemy↗
PubMed1995

Acid and thermal denaturation of barnase investigated by molecular dynamics simulations.

The transition in barnase from the native state to a partially unfolded conformation has been studied by molecular dynamics simulations with explicit water molecules at 360 K and low pH(450 ps), and at 600 K and neutral pH (three simulations of 120, 250 and 200 ps each). The use of several simulations provides evidence that the results are not sensitive to initial conditions. To mimic low pH conditions, the acidic sidechains in barnase were neutralized and the two histidine residues were doubly protonated. Runs at 300 K showed that the solvated structures at low pH (300 ps) and neutral pH (310 ps) are very similar. The main structural differences involved the acidic residues, histidine residues, and the beta-turn connecting strands 4 and 5. When the temperature is raised to 360 K at low pH and to 600 K at neutral pH the barnase molecule begins to unfold. The molecule rapidly expands (Rg changes from 13.9 A to 15.3 A in 450 ps at 360 K and from 13.7 A to between 15.1 and 15.5 A in 120 ps at 600 K). However, the expansion is not uniform. In all the simulations, the chain termini, loops and the N-terminal parts of the main alpha-helix (helix 1) show a continuous and progressive unfolding. An essential step in the denaturation process is that the major alpha-helix (helix 1) separates from the beta-sheet; this is coupled to the exposure of the principal hydrophobic core, many of whose non-polar side chains become solvated by hydrogen-bonded water molecules. The barnase-water interaction energy improves during unfolding at the expense of the barnase self-energy. The deterioration of the intramolecular van der Waals energy suggests that the rupture of the tight packing during the initial unfolding phase contributes to the energy barrier of the denaturation process. The mutationally well-analyzed Asp8-Arg110-Asp12 double salt-bridge on the barnase surface is found to be marginally stable in the folded form in the simulations. A Poisson-Boltzmann calculation indicates that the salt-bridge is unstable; this is probably due to an overestimate of the solvation energy. A detailed analysis of the main hydrophobic core reveals that increase in solvent-accessible surface area and penetration of water molecules are simultaneous in the high-temperature simulation; at lower temperatures there is significant cavity formation and the entrance of the water molecules is somewhat delayed. The cavities occur in the neighborhood of the hydrophobic sidechains; the region formed by the sidechains of Val10, Leu14, Leu20, Tyr24, Ala74, Ile76 and Tyr90 is involved. The loosening of the core packing is coupled to an increase in the number of dihedral transitions.(ABSTRACT TRUNCATED AT 400 WORDS)

Bacterial Proteins↗
PubMed1995

Packaging of proteases and proteoglycans in the granules of mast cells and other hematopoietic cells. A cluster of histidines on mouse mast cell protease 7 regulates its binding to heparin serglycin proteoglycans.

Mouse mast cell protease 7 (mMCP-7) is a tryptase stored in the secretory granules of mast cells. At the granule pH of 5.5, mMCP-7 is fully active and is bound to heparin-containing serglycin proteoglycans. to understand the interaction of mMCP-7 with heparin inside and outside the mast cell, this trytase was first studied by comparative protein modeling. The "pro" form of mMCP-7 was then expressed in insect cells and studied by site-directed mutagenesis. Although mMCP-7 lacks known linear sequences of amino acis that interact with heparin, the three-dimensional model of mMCP-7 revealed an area on the surface of the folded protein away from the substrate-binding site that exhibits a strong positive electrostatic potential at the acidic pH of the granule. In agreement with this calculation, recombinant pro-mMCP-7 bound to a heparin-affinity column at pH 5.5 and readily dissociated from the column at pH > 6.5. Site-directed mutagenesis confirmed the prediction that the conversion of His residues 8,68, and 70 in the positively charged region into Glu prevents the binding of pro-mMCP-7 to heparin. Because the binding requires positively charged His residues, native mMCP-7 is able to dissociate from the protease/proteoglycan macromolecular complex when the complex is exocytosed from bone marrow-derived mast cells into a neutral pH environment. Many hematopoietic effector cells store positively charged proteins in granules that contain serglycin proteoglycans. The heparin/mMCP-7 interaction, which depends on the tertiary structure of the tryptase, may be representative of a general control mechanism by which hematopoietic cells maximize storage of properly folded, enzymatically active proteins in their granules.

Amino Acid Sequence↗