PubMed Health⌕ Search

Biomedical subjects

M Karplus

Publications and source records attributed to M Karplus.

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

Three key residues form a critical contact network in a protein folding transition state.

Determining how a protein folds is a central problem in structural biology. The rate of folding of many proteins is determined by the transition state, so that a knowledge of its structure is essential for understanding the protein folding reaction. Here we use mutation measurements--which determine the role of individual residues in stabilizing the transition state--as restraints in a Monte Carlo sampling procedure to determine the ensemble of structures that make up the transition state. We apply this approach to the experimental data for the 98-residue protein acylphosphatase, and obtain a transition-state ensemble with the native-state topology and an average root-mean-square deviation of 6 A from the native structure. Although about 20 residues with small positional fluctuations form the structural core of this transition state, the native-like contact network of only three of these residues is sufficient to determine the overall fold of the protein. This result reveals how a nucleation mechanism involving a small number of key residues can lead to folding of a polypeptide chain to its unique native-state structure.

Acid Anhydride Hydrolases↗
PubMed2001

Evaluation of designed ligands by a multiple screening method: application to glycogen phosphorylase inhibitors constructed with a variety of approaches.

Glycogen phosphorylase (GP) is an important enzyme that regulates blood glucose level and a key therapeutic target for the treatment of type II diabetes. In this study, a number of potential GP inhibitors are designed with a variety of computational approaches. They include the applications of MCSS, LUDI and CoMFA to identify additional fragments that can be attached to existing lead molecules; the use of 2D and 3D similarity-based QSAR models (HQSAR and SMGNN) and of the LUDI program to identify novel molecules that may bind to the glucose binding site. The designed ligands are evaluated by a multiple screening method, which is a combination of commercial and in-house ligand-receptor binding affinity prediction programs used in a previous study (So and Karplus, J. Comp.-Aid. Mol. Des., 13 (1999), 243-258). Each method is used at an appropriate point in the screening, as determined by both the accuracy of the calculations and the computational cost. A comparison of the strengths and weaknesses of the ligand design approaches is made.

Binding Sites↗
PubMed2001

Functional group placement in protein binding sites: a comparison of GRID and MCSS.

One approach to combinatorial ligand design begins by determining optimal locations (i.e., local potential energy minima) for functional groups in the binding site of a target macromolecule. MCSS and GRID are two methods, based on significantly different algorithms, which are used for this purpose. A comparison of the two methods for the same functional groups is reported. Calculations were performed for nonpolar and polar functional groups in the internal hydrophobic pocket of the poliovirus capsid protein, and on the binding surface of the src SH3 domain. The two approaches are shown to agree qualitatively; i.e., the global characteristics of the functional group maps generated by MCSS and GRID are similar. However, there are significant differences in the relative interaction energies of the two sets of minima, a consequence of the different functional form used to evaluate polar interactions (electrostatics and hydrogen bonding) in the two methods. The single sphere representation used by GRID affords only positional information, supplemented by the identification of hydrogen bonding interactions. By contrast, the multi-atom representation of most MCSS groups yields in both positional and orientational information. The two methods are most similar for small functional groups, while for larger functional groups MCSS yields results consistent with GRID but superior in detail. These results are in accord with the somewhat different purposes for which the two methods were developed. GRID has been used mainly to introduce functionalities at specific positions in lead compounds, in which case the orientation is predetermined by the structure of the latter. The orientational information provided by MCSS is important for its use in the de novo design of large, multi-functional ligands, as well as for improving lead compounds.

Algorithms↗
PubMed2000

A dynamic model for the allosteric mechanism of GroEL.

GroEL-assisted protein folding is regulated by a cycle of large coordinated domain movements in the 14-subunit double-ring assembly. The transition path between the closed (unliganded) and the open (liganded) states, calculated with a targeted molecular dynamics simulation, shows the highly complex subunit displacements required for the allosteric transition. The early downward motion of the small intermediate domain induced by nucleotide binding emerges as the trigger for the larger movements of the apical and equatorial domains. The combined twisting and upward displacement of the apical domain determined for a single subunit is accommodated easily in the heptamer ring only if its opening is concerted. This is a major source of cooperative ligand binding within a ring. It suggests also that GroEL has evolved so that the motion required for heptamer cooperativity is encoded in the individual subunits. A calculated model for a di-cis 14-subunit assembly is found to be destabilized by strong steric repulsion between the equatorial domains of the two rings, the source of negative cooperativity. The simulation results, which indicate that transient interactions along the transition path are essential for GroEL function, provide a detailed structural description of the motions that are involved in the GroEL allosteric cycle.

Allosteric Regulation↗
PubMed2000

Structures of scrambled disulfide forms of the potato carboxypeptidase inhibitor predicted by molecular dynamics simulations with constraints.

The structures of two species of potato carboxypeptidase inhibitor with nonnative disulfide bonds were determined by molecular dynamics simulations in explicit solvent using disulfide bond constraints that have been shown to work for the native species. Ten structures were determined; five for scrambled A (disulfide bonds between Cys8-Cys27, Cys12-Cys18, and Cys24-Cys34) and five for the scrambled C (disulfide bonds Cys8-Cys24, Cys12-Cys18, and Cys27-Cys34). The two scrambled species were both more solvent exposed than the native structure; the scrambled C species was more solvent exposed and less compact than the scrambled A species. Analysis of the loop regions indicates that certain loops in scrambled C are more nativelike than in scrambled A. These factors, combined with the fact that scrambled C has one native disulfide bond, may contribute to the observed faster conversion to the native structure from scrambled C than from scrambled A. Results from the PROCHECK program using the standard parameter database and a database specially constructed for small, disulfide-rich proteins indicate that the 10 scrambled structures have correct stereochemistry. Further, the results show that a characteristic feature of small, disulfide-rich proteins is that they score poorly using the standard PROCHECK parameter database. Proteins 2000;40:482-493.

Carboxypeptidases↗
PubMed2000

Dynamic ligand design and combinatorial optimization: designing inhibitors to endothiapepsin.

The dynamic ligand design (DLD) algorithm, an automated method for the creation of novel ligands, links up small functional groups that have been placed in energetically favorable positions in the binding site of a target molecule. The positions and orientations of the small functional groups can be determined using the multi-copy simultaneous search approach (MCSS) or experimental data. In this work the original DLD methodology is extended by using a modified version of the pseudo-potential energy function. A novel simulated annealing protocol is presented for optimizing the pseudo-potential energy of ligands in the binding site; the protocol is expected to be applicable to other optimization problems. The utility of the method is illustrated by designing an inhibitor for endothiapepsin. The binding affinity of the inhibitor is assessed using a thermodynamic cycle that decomposes the binding free energy into a sum of translational, rotational, configurational, hydrophobic, and electrostatic contributions. The calculations suggest that the designed molecule will bind endothiapepsin with high affinity.

Algorithms↗
PubMed2000

Modeling of the TCR-MHC-peptide complex.

The methodology for generating a homology model of the T1 TCR-PbCS-K(d) class I major histocompatibility complex (MHC) class I complex is presented. The resulting model provides a qualitative explanation of the effect of over 50 different mutations in the region of the complementarity determining region (CDR) loops of the T cell receptor (TCR), the peptide and the MHC's alpha(1)/alpha(2) helices. The peptide is modified by an azido benzoic acid photoreactive group, which is part of the epitope recognized by the TCR. The construction of the model makes use of closely related homologs (the A6 TCR-Tax-HLA A2 complex, the 2C TCR, the 14.3.d TCR Vbeta chain, the 1934.4 TCR Valpha chain, and the H-2 K(b)-ovalbumine peptide), ab initio sampling of CDR loops conformations and experimental data to select from the set of possibilities. The model shows a complex arrangement of the CDR3alpha, CDR1beta, CDR2beta and CDR3beta loops that leads to the highly specific recognition of the photoreactive group. The protocol can be applied systematically to a series of related sequences, permitting the analysis at the structural level of the large TCR repertoire specific for a given peptide-MHC complex.

Algorithms↗
PubMed2000

Unfolding proteins by external forces and temperature: the importance of topology and energetics.

Unfolding of proteins by forced stretching with atomic force microscopy or laser tweezer experiments complements more classical techniques using chemical denaturants or temperature. Forced unfolding is of particular interest for proteins that are under mechanical stress in their biological function. For beta-sandwich proteins (a fibronectin type III and an immunoglobulin domain), both of which appear in the muscle protein titin, the results of stretching simulations show important differences from temperature-induced unfolding, but there are common features that point to the existence of folding cores. Intermediates detected by comparing unfolding with a biasing perturbation and a constant pulling force are not evident in temperature-induced unfolding. For an alpha-helical domain (alpha-spectrin), which forms part of the cytoskeleton, there is little commonality in the pathways from unfolding induced by stretching and temperature. Comparison of the forced unfolding of the two beta-sandwich proteins and two alpha-helical proteins (the alpha-spectrin domain and an acyl-coenzyme A-binding protein) highlights important differences within and between protein classes that are related to the folding topologies and the relative stability of the various structural elements.

Protein Folding↗
PubMed2000

Protein motions at zero-total angular momentum: the importance of long-range correlations.

A constant-energy molecular dynamics simulation is used to monitor protein motion at zero-total angular momentum. With a simple protein model, it is shown that overall rotation is possible at zero-total angular momentum as a result of flexibility. Since the rotational motion is negligible on a time scale of 1000 reduced time units, the essentially rotation-free portion of the trajectory provides an unbiased test of the common approximate methods for separating overall rotation from internal motions by optimal superposition. Removing rotation by minimizing the root-mean-square deviation (RMSD) for the entire system is found to be more appropriate than using the RMSD for only the more rigid part of the system. The results verify the existence of positive cross-correlation in the motions of atoms separated by large distances.

Models, Molecular↗
PubMed2000

Congenital dyserythropoietic anemia type I presenting as persistent pulmonary hypertension of the newborn.

Congenital dyserythropoietic anemia (CDA) is a rare group of inherited bone marrow disorders characterized by anemia with ineffective erythropoiesis. We report 3 siblings from a family known to have CDA type I who presented with persistent pulmonary hypertension of the newborn (PPHN). We suggest that the diagnosis of CDA type I should be considered in any neonate with PPHN and anemia.

Anemia, Dyserythropoietic, Congenital↗
PubMed2000

Effective energy functions for protein structure prediction.

Protein structure prediction, fold recognition, homology modeling and design rely mainly on statistical effective energy functions. Although the theoretical foundation of such functions is not clear, their usefulness has been demonstrated in many applications. Molecular mechanics force fields, particularly when augmented by implicit solvation models, provide physical effective energy functions that are beginning to play a role in this area.

Animals↗
PubMed2000

Understanding protein folding via free-energy surfaces from theory and experiment.

The ability of protein molecules to fold into their highly structured functional states is one of the most remarkable evolutionary achievements of biology. In recent years, our understanding of the way in which this complex self-assembly process takes place has increased dramatically. Much of the reason for this advance has been the development of energy surfaces (landscapes), which allow the folding reaction to be described and visualized in a meaningful manner. Analysis of these surfaces, derived from the constructive interplay between theory and experiment, has led to the development of a unified mechanism for folding and a recognition of the underlying factors that control the rates and products of the folding process.

Animals↗
PubMed2000

Design of dimerization inhibitors of HIV-1 aspartic proteinase: a computer-based combinatorial approach.

Inhibition of dimerization to the active form of the HIV-1 aspartic proteinase (HIV-1 PR) may be a way to decrease the probability of escape mutations for this viral protein. The Multiple Copy Simultaneous Search (MCSS) methodology was used to generate functionality maps for the dimerization interface of HIV-1 PR. The positions of the MCSS minima of 19 organic fragments, once postprocessed to take into account solvation effects, are in good agreement with experimental data on peptides that bind to the interface. The MCSS minima combined with an approach for computational combinatorial ligand design yielded a set of modified HIV-1 PR C-terminal peptides that are similar to known nanomolar inhibitors of HIV-1 PR dimerization. A number of N-substituted 2,5-diketopiperazines are predicted to be potential dimerization inhibitors of HIV-1 PR.

Aspartic Acid Endopeptidases↗
PubMed2000

Solvent mobility and the protein 'glass' transition.

Proteins and other biomolecules undergo a dynamic transition near 200 K to a glass-like solid state with small atomic fluctuations. This dynamic transition can inhibit biological function. To provide a deeper understanding of the relative importance of solvent mobility and the intrinsic protein energy surface in the transition, a novel molecular dynamics simulation procedure with the protein and solvent at different temperatures has been used. Solvent mobility is shown to be the dominant factor in determining the atomic fluctuations above 180 K, although intrinsic protein effects become important at lower temperatures. The simulations thus complement experimental studies by demonstrating the essential role of solvent in controlling functionally important protein fluctuations.

Computer Simulation↗
PubMed2000

Early postnatal dexamethasone treatment and increased incidence of cerebral palsy.

OBJECTIVE: To study the long term neurodevelopmental outcome of children who participated in a randomised, double blind, placebo controlled study of early postnatal dexamethasone treatment for prevention of chronic lung disease. METHODS: The original study compared a three day course of dexamethasone (n = 132) with a saline placebo (n = 116) administered from before 12 hours of age in preterm infants, who were ventilated for respiratory distress syndrome and had received surfactant treatment. Dexamethasone treatment was associated with an increased incidence of hypertension, hyperglycaemia, and gastrointestinal haemorrhage and no reduction in either the incidence or severity of chronic lung disease or mortality. A total of 195 infants survived to discharge and five died later. Follow up data were obtained on 159 of 190 survivors at a mean (SD) age of 53 (18) months. RESULTS: No differences were found between the groups in terms of perinatal or neonatal course, antenatal steroid administration, severity of initial disease, or major neonatal morbidity. Dexamethasone treated children had a significantly higher incidence of cerebral palsy than those receiving placebo (39/80 (49%) v. 12/79 (15%) respectively; odds ratio (OR) 4.62, 95% confidence interval (95% CI) 2.38 to 8.98). The most common form of cerebral palsy was spastic diplegia (incidence 22/80 (28%) v. 5/79 (6%) in dexamethasone and placebo treated infants respectively; OR 4.45, 95% CI 1.95 to 10.15). Developmental delay was significantly more common in the dexamethasone treated group (44/80 (55%)) than in the placebo treated group (23/79 (29%); OR 2. 87, 95% CI 1.53 to 5.38). Dexamethasone treated infants had more periventricular leucomalacia and less intraventricular haemorrhage in the neonatal period than those in the placebo group, although these differences were not statistically significant. Eleven children with cerebral palsy had normal ultrasound scans in the neonatal period; all 11 had received dexamethasone. Logistic regression analysis showed both periventricular leucomalacia and drug assignment to dexamethasone to be highly significant predictors of abnormal neurological outcome. CONCLUSIONS: A three day course of dexamethasone administered shortly after birth in preterm infants with respiratory distress syndrome is associated with a significantly increased incidence of cerebral palsy and developmental delay.

Anti-Inflammatory Agents↗
PubMed1999

MCSS functionality maps for a flexible protein.

The Multiple Copy Simultaneous Search (MCSS) methodology for finding energetically favorable positions and orientations of small functional groups in a binding site is extended to include flexibility of the target. This makes possible the finding of novel minima not present in a fixed structure and so extends the diversity of inhibitors that can be constructed starting with the MCSS procedure. Quenched molecular dynamics is used to generate energetically favorable positions and orientations of the functional groups in the field of a flexible protein. The method is applied to the viral protein HIV-1 protease with methanol and methyl ammonium as a test case. If the protein is quenched with many copies of functional groups randomly distributed in the binding site, the resulting minima have ligand-protein interaction energies that are, on average, less favorable than those obtained with standard MCSS. This is a consequence of the renormalized potential function employed in the Locally Enhanced Sampling (LES) approximation. However, local optimizations of existing MCSS minima with a flexible protein results in lower energy minima in regions of the protein that are of particular interest. Their use in constructing a consensus protein model for ligand design is discussed.

Algorithms↗
PubMed1999

Folding of a model three-helix bundle protein: a thermodynamic and kinetic analysis.

The kinetics and thermodynamics of an off-lattice model for a three-helix bundle protein are investigated as a function of a bias gap parameter that determines the energy difference between native and non-native contacts. A simple dihedral potential is used to introduce the tendency to form right-handed helices. For each value of the bias parameter, 100 trajectories of up to one microsecond are performed. Such statistically valid sampling of the kinetics is made possible by the use of the discrete molecular dynamics method with square-well interactions. This permits much faster simulations for off-lattice models than do continuous potentials. It is found that major folding pathways can be defined, although ensembles with considerable structural variation are involved. The large gap models generally fold faster than those with a smaller gap. For the large gap models, the kinetic intermediates are non-obligatory, while both obligatory and non-obligatory intermediates are present for small gap models. Certain large gap intermediates have a two-helix microdomain with one helix extended outward (as in domain-swapped dimers); the small gap intermediates have more diverse structures. The importance of studying the kinetic, as well as the thermodynamics, of folding for an understanding of the mechanism is discussed and the relation between kinetic and equilibrium intermediates is examined. It is found that the behavior of this model system has aspects that encompass both the "new" view and the "old" view of protein folding.

Algorithms↗