Search PubMed⌕ Search

Biomedical subjects

M Karplus

Publications and source records attributed to M Karplus.

At least 199 records · Page 11Linked to original sources
PubMed1987

Analysis of side-chain orientations in homologous proteins.

The side-chain conformations of topologically equivalent residues in seven pairs of proteins ranging in sequence homology from 16% to 60% are compared. Both identical and mutated residues are included. For proteins with greater than 40% homology, it is found that at least 80% of the side-chain orientations of identical residues and 75% or more of the mutated residues in each pair of proteins have matching gamma atom dihedral angles (+/- 40 degrees); the comparison is not based strictly on chi 1 angles. Further, if a match is obtained at the gamma position, there is a high probability of matching for the delta atom(s) of the side-chain. For proteins with less than 25% homology the percentages are somewhat lower. Trends observed for conservative substitutions are essentially the same as those noted for mutated residues in general. Side-chain accessibility does not affect the probability of matches of identical residues; however, less accessible pairs of mutated residues have 10 to 20% higher matching probabilities than do exposed residues. Mismatches can frequently be related to large B-factors, certain types of amino acid substitutions, or the appearance of multiple minima on the side-chain potential energy surfaces and are most likely to occur for certain small residues (Ser, Thr, Val). Analysis of all the results makes possible the formulation of a set of rules for side-chain positioning in the modeling of homologous proteins.

Amino Acid Sequence↗
PubMed1987

Refinement of the solution structure of the DNA hexamer 5'd(GCATGC)2: combined use of nuclear magnetic resonance and restrained molecular dynamics.

The solution structure of the self-complementary DNA hexamer 5'd(GCATGC)2 comprising the specific target site for the restriction endonuclease Sph 1 is investigated by using nuclear magnetic resonance spectroscopy and restrained molecular dynamics. All the nonexchangeable proton resonances are assigned sequentially, and from time-dependent nuclear Overhauser enhancement measurements a set of 158 approximate interproton distances are determined. These distances are used as the basis of a structure refinement using restrained molecular dynamics in which the interproton distances are incorporated into the total energy function of the system in the form of an effective potential term. Two restrained molecular dynamics simulations are carried out, starting from classical B- and A-DNA [atomic root mean square (rms) difference 3.3 A]. In both cases convergence is achieved to essentially identical structures satisfying the experimental restraints and having a root mean square difference of only 0.3 A between them, which is within the rms fluctuations of the atoms about their average positions. These results suggest that the restrained molecular dynamics structures represent reasonable approximations of the solution structure. The converged structures are of the B type and exhibit clear sequence-dependent variations of helical parameters, some of which follow Calladine's rules and can be attributed to the relief of interstrand purine-purine clash at adjacent base pairs. In addition, the converged restrained dynamics structures appear bent with a radius of curvature of approximately 20 A. This bending appears to be due almost entirely to the large positive base roll angles, particularly at the Pyr-Pur steps. Further, the global and local helix axes are not coincident, and the global helix axis represents a superhelical axis which the bent DNA, when extended into an "infinite" helix by repeated translation and rotation, wraps around.

Calorimetry↗
PubMed1987

Conformations of intermediates in the folding of the pancreatic trypsin inhibitor.

Intermediates in the folding pathway of the bovine pancreatic trypsin inhibitor (PTI) have been examined by 1H nuclear magnetic resonance (n.m.r.). The intermediates were trapped during the reoxidation and consequent refolding of reduced PTI by alkylating free thiols; each intermediate contained different disulphide linkages. The n.m.r. spectra reveal that conformational features of the native protein are present in the intermediate containing just one of the three normal disulphide linkages (30-51). As additional normal disulphide bonds are formed, the conformation becomes more similar to that of the native protein. Introduction of additional but incorrect disulphide bonds does not lead to an increase in observable globular structure. A description of the folding process in terms of the conformations of the different intermediates is proposed. The significance of these results for the general mechanism of protein folding is outlined.

Amino Acid Sequence↗
PubMed1987

A comparison of the restrained molecular dynamics and distance geometry methods for determining three-dimensional structures of proteins on the basis of interproton distances.

A direct comparison of the metric matrix distance geometry and restrained molecular dynamics methods for determining three-dimensional structures of proteins on the basis of interproton distances is presented using crambin as a model system. It is shown that both methods reproduce the overall features of the secondary and tertiary structure (shape and polypeptide fold). The region of conformational space sampled by the converged structures generated by the two methods is similar in size, and in both cases the converged structures are distributed about mean structures which are closer to the X-ray structure than any of the individual structures. The restrained molecular dynamics structures are superior to those obtained from distance geometry as regards local backbone conformation, side chain positions and non-bonding energies.

Chemical Phenomena↗
PubMed1987

1H NMR aromatic spectrum of the operator binding domain of the lambda repressor: resonance assignment with application to structure and dynamics.

The aromatic 1H NMR resonances of the operator binding domain of lambda repressor are completely assigned. Since the resonances of this 23-kilodalton domain are too broad for the application of two-dimensional strategies for sequence-specific assignment, an alternative approach has been used. Assignments are obtained by a combination of one- and two-dimensional NMR methods, by the study of genetically altered domains, and by the biosynthetic incorporation of deuterium labels. The resulting assignments provide sensitive markers for tertiary and quaternary structure. Nuclear Overhauser enhancements demonstrate that the major features of the crystal structure, including the dimer contacts, are retained in solution. The rates of aromatic ring rotation indicate that the globular domain is not rigid; significant barriers to ring rotation are observed only in the dimer contact.

Bacteriophage lambda↗
PubMed1987

Dimerization of the operator binding domain of phage lambda repressor.

Dimerization of lambda repressor is required for its binding to operator DNA. As part of a continuing study of the structural basis of the coupling between dimer formation and operator binding, we have undertaken 1H NMR and gel filtration studies of the dimerization of the N-terminal domain of lambda repressor. Five protein fragments have been studied: three are wild-type fragments of different length (1-102, 1-92, and 1-90), and two are fragments bearing single amino acid substitutions in residues involved in the dimer interface (1-102, Tyr-88----Cys; 1-92, Ile-84----Ser). The tertiary structure of each species is essentially the same, as monitored by the 1H NMR resonances of internal aromatic groups. However, significant differences are observed in their dimerization properties. 1H NMR resonances of aromatic residues that are involved in the dimer contact allow the monomer-dimer equilibrium to be monitored in solution. The structure of the wild-type dimer contact appears to be similar to that deduced from X-ray crystallography and involves the hydrophobic packing of symmetry-related helices (helix 5) from each monomer. Removal of two contact residues, Val-91 and Ser-92, by limited proteolysis disrupts this interaction and also prevents crystallization. The Ile-84----Ser substitution also disrupts this interaction, which accounts for the severely reduced operator affinity of this mutant protein.

Amino Acids↗
PubMed1987

Multiple conformational states of proteins: a molecular dynamics analysis of myoglobin.

A molecular dynamics simulation of myoglobin provides the first direct demonstration that the potential energy surface of a protein is characterized by a large number of thermally accessible minima in the neighborhood of the native structure (for example, approximately 2000 minima were sampled in a 300-picosecond trajectory). This is expected to have important consequences for the interpretation of the activity of transport proteins and enzymes. Different minima correspond to changes in the relative orientation of the helices coupled with side-chain rearrangements that preserve the close packing of the protein interior. The conformational space sampled by the simulation is similar to that found in the evolutionary development of the globins. Glasslike behavior is expected at low temperatures. The minima obtained from the trajectory do not satisfy certain criteria for ultrametricity.

Models, Structural↗
PubMed1987

Thermal expansion of a protein.

The thermal expansion of a protein, metmyoglobin, was investigated by analysis of the refined X-ray crystal structures at 80 and 255-300 K. On heating from 80 to 300 K, the volume occupied by myoglobin increases by approximately 3%. The linear thermal expansion coefficient is estimated to be 115 X 10(-6) K-1. This value is more than twice as large as that of liquid water but less than that of benzene. As the temperature is raised, the internal volume change does not come from the large, atom-sized internal cavities in the structure but from an increase in the small, subatomic free volumes between atoms. The largest expansion occurs in the region of the CD and GH corners; both these regions move away from the center of the protein. The remainder of the expansion results from the lengthening of contacts between segments of secondary structure.

Animals↗
PubMed1987

Estimation of uncertainties in X-ray refinement results by use of perturbed structures.

The uncertainties in the refined parameters for a 1.5-A X-ray structure of carbon-monoxy (FeII) myoglobin are estimated by combining energy minimization with least-squares refinement against the X-ray data. The energy minimizations, done without reference to the X-ray data, provide perturbed structures which are used to restart conventional X-ray refinement. The resulting refined structures have the same, or better, R-factor and stereochemical parameters as the original X-ray structure, but deviate from it by 0.13 A rms for the backbone atoms and 0.31 A rms for the sidechain atoms. Atoms interacting with a disordered sidechain, Arg 45 CD3, are observed to have larger positional uncertainties. The uncertainty in the B-factors, within the isotropic harmonic motion approximation, is estimated to be 15%. The resulting X-ray structures are more consistent with the energy parameters used in simulations.

Arginine↗
PubMed1987

Anisotropy and anharmonicity of atomic fluctuations in proteins: analysis of a molecular dynamics simulation.

Positional probability density functions (pdf) for the atomic fluctuations are determined from a molecular dynamics simulation for hen egg-white lysozyme. Most atoms are found to have motions that are highly anisotropic but only slightly anharmonic. The largest deviations from harmonic motion are in the direction of the largest rms fluctuations in the local principal axis frame. Backbone atoms tend to be more nearly harmonic than sidechain atoms. The atoms with the largest anharmonicities tend to have pdfs with multiple peaks, each of which is close to harmonic. Several model pdfs are evaluated on the basis of how well they fit probability densities from the dynamics simulations when parameterized in terms of the moments of the distribution. Gram-Charlier and Edgeworth perturbation expansions, which have been successful in describing the motions of small molecules in crystals, are shown to be inadequate for the distributions found in the dynamics of proteins. Multipeaked distribution functions are found to be more appropriate.

Animals↗
PubMed1987

Configurational entropy of native proteins.

Simulations of the residual configurational entropy of a protein in the native state suggest that it is nearly an order of magnitude larger than the entropy of denaturation. The implications of this result are discussed.

Mathematics↗
PubMed1987

Molecular and stochastic dynamics of proteins.

The rapid fluctuations of protein atoms derived from molecular dynamics simulations can be extrapolated to longer-time motions by effective single-particle stochastic models. This is demonstrated by an analysis of velocity autocorrelation functions for the atoms of lysine side chains in the active site of RNase A. The atomic motions are described by a bounded stochastic model with the friction and noise parameters determined from a molecular dynamics simulation. The low-frequency relaxation behavior is shown to result from collisional damping rather than dephasing. Extrapolation of these results to the quasistochastic motion of the heme group in myoglobin provides an explanation of 57Fe Mössbauer spectroscopic data.

Binding Sites↗
PubMed1987

Quaternary structure and function in phage lambda repressor: 1H-NMR studies of genetically altered proteins.

The quaternary structure and dynamics of phage lambda repressor are investigated in solution by 1H-NMR methods. lambda repressor contains two domains separable by proteolysis: an N-terminal domain that mediates sequence-specific DNA-A binding, and a C-terminal domain that contains strong dimer and higher-order contacts. The active species in operator recognition is a dimer. Although the crystal structure of an N-terminal fragment has been determined, the intact protein has not been crystallized, and there is little evidence concerning its structure. 1H-NMR data indicate that the N-terminal domain is only loosely tethered to the C-terminal domain, and that its tertiary structure is unperturbed by proteolysis of the "linker" polypeptide. It is further shown that in the intact repressor structure a quaternary interaction occurs between N-terminal domains. This domain-domain interaction is similar to the dimer contact observed in the crystal structure of the N-terminal fragment and involves the hydrophobic packing of symmetry-related helices (helix 5). In the intact structure this interaction is disrupted by the single amino-acid substitution, Ile84----Ser, which reduces operator affinity at least 100-fold. We conclude that quaternary interactions between N-terminal domains function to appropriately orient the DNA-binding surface with respect to successive major grooves of B-DNA.

Bacteriophage lambda↗
PubMed1987

Models for the three-dimensional structure of renin inhibitors bound in the active site of human renin: an analysis of the properties that produce tight binding.

Three-dimensional models of the octapeptide segment of renin substrate and of several inhibitors characterized by a modified scissile bond were constructed at the active site of human renin. The substrate and inhibitor models were based on the structure of pepstatin solved by x-ray crystallography to a resolution of 1.8 A. The renin structure, previously reported by us, was based on the structures of homologous aspartyl proteases solved by x-ray diffraction techniques to a resolution of 2.1 A or higher. An energy minimization program, CHARMM, was used to refine these structures with respect to the optimal interaction of the enzyme and its inhibitor, and to calculate the energy of each inhibitor-enzyme interaction. The results indicate that the most significant difference between the binding of inhibitors and of the natural substrate is caused by a reduction in electrostatic rather than conformational strain.

Binding Sites↗
PubMed1986

Interdomain motion in liver alcohol dehydrogenase. Structural and energetic analysis of the hinge bending mode.

A study of the hinge bending mode in the enzyme liver alcohol dehydrogenase is made by use of empirical energy functions. The enzyme is a dimer, with each monomer composed of a coenzyme binding domain and a catalytic domain with a large cleft between the two. Superposition of the apoenzyme and holoenzyme crystal structures is used to determine a rigid rotation axis for closing of the cleft. It is shown that a rigid body transformation of the apoenzyme to the holoenzyme structure corresponds to a 10 degrees rotation of the catalytic domain about this axis. The rotation is not along the least-motion path for closing of the cleft but instead corresponds to the catalytic domain coming closer to the coenzyme binding domain by a sliding motion. Estimation of the energy associated with the interdomain motion of the apoenzyme over a range of 90 degrees (-40 to 50 degrees, where 0 degrees corresponds to the minimized crystal structure) demonstrates that local structural relaxation makes possible large-scale rotations with relatively small energy increments. A variety of structural rearrangements associated with the domain motion are characterized. They involve the hinge region residues that provide the covalent connections between the two domains and certain loop regions that are brought into contact by the rotation. Differences between the energy minimized and the holoenzyme structures point to the existence of alternative conformations for loops and to the importance of the ligands in the structural rearrangements.

Alcohol Dehydrogenase↗