Search PubMed⌕ Search

Biomedical subjects

M Karplus

Publications and source records attributed to M Karplus.

At least 235 records · Page 13Linked to original sources
PubMed1985

Solution conformation of a heptadecapeptide comprising the DNA binding helix F of the cyclic AMP receptor protein of Escherichia coli. Combined use of 1H nuclear magnetic resonance and restrained molecular dynamics.

A nuclear magnetic resonance study on a heptadecamer (17-mer) peptide comprising the DNA binding helix F of the cyclic AMP receptor protein of Escherichia coli is presented under solution conditions (viz. 40% (v/v) trifluorethanol) where it adopts an ordered helical structure as judged by circular dichroism. Using a combination of two-dimensional nuclear magnetic resonance techniques, complete resonance assignments are obtained in a sequential manner. From the two-dimensional nuclear Overhauser enhancement spectra, a set of 87 approximate distance restraints is derived and used as the basis for three-dimensional structure determination with a restrained molecular dynamics algorithm in which the interproton distances are incorporated into the total energy function of the system in the form of an additional effective potential term. The convergence properties of this approach are tested by starting from three different initial structures, namely an alpha-helix, a beta-strand and a 3-10 helix. In all three cases, convergence to an alpha-helical structure is achieved with a root mean square difference of less than 3 A for all atoms and less than 2 A for the backbone atoms.

Amino Acid Sequence↗
PubMed1985

Transition from B to Z DNA: contribution of internal fluctuations to the configurational entropy difference.

The internal motions of the double-stranded DNA oligomer (dCdG)3 (dC, deoxycytidylate; dG, deoxyguanylate) in the B and Z forms have been calculated in the harmonic approximation. A complete vibrational analysis has been made, and the resulting normal mode frequencies have been used to evaluate the vibrational entropy of B and Z DNA. The greater flexibility of the B DNA hexamer leads to an entropic stabilization relative to the stiffer Z DNA hexamer of 22 calories per mole per kelvin at 300 K. The calculated value is of the same order as that (21 to 27 calories per mole per kelvin) obtained from nuclear magnetic resonance measurements on the methylated duplexes (m5dCdG)3 and (dCdGm5dCdGdCdG). This result demonstrates the importance of internal motions, which have been neglected in earlier studies of the transition from B to Z DNA, in the stability of different nucleic acid conformers.

DNA↗
PubMed1985

Vicinal coupling constants and protein dynamics.

The effects of motional averaging on the analysis of vicinal spin-spin coupling constants derived from proton NMR studies of proteins have been examined. Trajectories obtained from molecular dynamics simulations of bovine pancreatic trypsin inhibitor and of hen egg white lysozyme were used in conjunction with an expression for the dependence of the coupling constant on the intervening dihedral angle to calculate the time-dependent behavior of the coupling constants. Despite large fluctuations, the time-average values of the coupling constants are not far from those computed for the average structure in the cases where fluctuations occur about a single potential well. The calculated differences show a high correlation with the variation in the magnitude of the fluctuations of individual dihedral angles. For the cases where fluctuations involve multiple sites, large differences are found between the time-average values and the average structure values for the coupling constants. Comparison of the simulation results with the experimental trends suggests that side chains with more than one position are more common in proteins than is inferred from X-ray results. It is concluded that for the main chain, motional effects do not introduce significant errors where vicinal coupling constants are used in structure determinations; however, for side chains, the motional average can alter deductions about the structure. Accurately measured coupling constants are shown to provide information concerning the magnitude of dihedral angle fluctuations.

Amino Acids↗
PubMed1985

A de novo translocation, 14q21q, with a microchromosome-14p21p.

A familial translocation, t(14;21)(14p21p;14q21q), in a mother and her child is described. The translocation was ascertained through the birth of a Down syndrome baby with the chromosome constitution 47,XX,-14, +der 14, +der 21,t(14;21)(q11;p12) mat. A 1:3 segregation in the maternal meiosis is suggested for the evolution of the unbalanced chromosome state. The main translocated chromosome 14q21q mimics the product of a Robertsonian translocation, while the 14p21p chromosome has the morphology of a satellited microchromosome. The cytogenetic nature of this translocation is discussed.

Adult↗
PubMed1985

Molecular dynamics of myoglobin at 298 degrees K. Results from a 300-ps computer simulation.

We have carried out a very long (300 ps) molecular dynamics simulation of the protein myoglobin. This trajectory is approximately three times longer than the longest previous molecular dynamics simulation of a protein, and ten times longer than protein simulations of comparable size (1,423 atoms in our model). Here we report results from this long simulation concerning the average structure, the mean square fluctuations of atoms about the average structure, and the nuclear magnetic resonance order parameters for various groups in myoglobin. The results demonstrate that the average coordinates change very slowly during the simulation. The relative atomic mobilities are well described by the simulation. For both the mean square atomic fluctuations and the order parameters, however, there are significant quantitative differences when values calculated using shorter portions of the trajectory are compared with results obtained for the entire 300-ps simulation. The implications of this result for obtaining converged properties from protein molecular dynamics simulations for comparison with experiment are discussed.

Leucine↗
PubMed1985

Normal modes for specific motions of macromolecules: application to the hinge-bending mode of lysozyme.

A method is presented for finding particular normal modes for large molecules such as proteins and nucleic acids. The method is based on an iterative approach that extracts eigenvectors of interest from the full second-derivative matrix. Application of the method to the interdomain (hinge-bending) motion of lysozyme yields a frequency of 3.6 cm-1. This is similar to the value obtained from earlier adiabatic-energy-minimization studies. Analysis of the mode shows that the relaxation associated with the hinge bending is highly delocalized; that is, the dihedral angle and energy changes are distributed over many residues, including some (e.g., Trp-28) that are distant from the cleft and hinge region.

Animals↗
PubMed1985

Active site dynamics of ribonuclease.

The stochastic boundary molecular dynamics method is used to study the structure, dynamics, and energetics of the solvated active site of bovine pancreatic ribonuclease A. Simulations of the native enzyme and of the enzyme complexed with the dinucleotide substrate CpA and the transition-state analog uridine vanadate are compared. Structural features and dynamical couplings for ribonuclease residues found in the simulation are consistent with experimental data. Water molecules, most of which are not observed in crystallographic studies, are shown to play an important role in the active site. Hydrogen bonding of residues with water molecules in the free enzyme is found to mimic the substrate-enzyme interactions of residues involved in binding. Networks of water stabilize the cluster of positively charged active site residues. Correlated fluctuations between the uridine vanadate complex and the distant lysine residues are mediated through water and may indicate a possible role for these residues in stabilizing the transition state.

Animals↗
PubMed1985

Structure of proteins with single-site mutations: a minimum perturbation approach.

A large number of mutant proteins with single amino acid substitutions are now being produced. The ability to predict the structural changes expected from such mutations would aid greatly in the efficient utilization of the mutagenic techniques and in the interpretation of the changes in stability and function that result. A minimum perturbation approach is suggested as a first step in such structural predictions and is tested by application to a recently isolated variant of the hemagglutinin glycoprotein. The agreement between the predicted structure and that inferred from the x-ray refinement is encouraging and provides support for the proposed modeling procedure.

Hemagglutinin Glycoproteins, Influenza Virus↗
PubMed1984

Amino-terminal arm of the lambda repressor: a 1H NMR study.

The N-terminal arm of the lambda repressor is shown to be flexible in solution by one- and two-dimensional 1H NMR methods. In particular, the relaxation of Thr-2 is largely independent of macromolecular tumbling. The conformation of the operator-binding domain is not affected by the removal of the first three residues nor by a point mutation, Lys-4----Gln. These results support a proposed model of the lambda repressor-operator complex in which the N-terminal arm of the repressor is assumed to be flexible and to wrap around the operator double helix.

DNA-Binding Proteins↗
PubMed1984

An analysis of incorrectly folded protein models. Implications for structure predictions.

Proteins with homologous amino acid sequences have similar folds and it has been assumed that an unknown three-dimensional structure can be obtained from a known homologous structure by substituting new side-chains into the polypeptide chain backbone, followed by relatively small adjustment of the model. To examine this approach of structure prediction and, more generally, to isolate the characteristics of native proteins, we constructed two incorrectly folded protein models. Sea-worm hemerythrin and the variable domain of mouse immunoglobulin K-chain, two proteins with no sequence homology, were chosen for study; the former is composed of a bundle of four alpha-helices and the latter consists of two 4-stranded beta-sheets. Using an automatic computer procedure, hemerythrin side-chains were substituted into the immunoglobulin domain and vice versa. The structures were energy-minimized with the program CHARMM and the resulting structures compared with the correctly folded forms. It was found that the incorrect side-chains can be incorporated readily into both types of structures (alpha-helices, beta-sheets) with only small structural adjustments. After constrained energy-minimization, which led to an average atomic co-ordinate shift of no more than 0.7 to 0.9 A, the incorrectly folded models arrived at potential energy values comparable to those of the correct structures. Detailed analysis of the energy results shows that the incorrect structures have less stabilizing electrostatic, van der Waals' and hydrogen-bonding interactions. The difference is particularly pronounced when the electrostatic and van der Waals' energy terms are calculated by modified equations that include an approximate representation of solvent effects. The incorrectly folded structures also have a significantly larger solvent-accessible surface and a greater fraction of non-polar side-chain atoms exposed to solvent. Examination of their interior shows that the packing of side-chains at the secondary structure interfaces, although corresponding to sterically allowed conformations, deviates from the characteristics found in normal proteins. The analysis of incorrectly folded structures has made it clear that the absence of bad non-bonded contacts, though necessary, is not sufficient to demonstrate the validity of model-built structures and that modeling of homologous structures has to be accompanied by a thorough quantitative evaluation of the results. Further, certain features that characterize native proteins are made evident by their absence in misfolded models.

Amino Acid Sequence↗
PubMed1984

1H-NMR study of the lambda operator site OL1: assignment of the imino and adenine H2 resonances.

One- and two-dimensional proton NMR methods are being used to study the synthetic lambda operator site O-L1, a 17 base-pair DNA duplex recognized by lambda repressor and Cro protein. The complete assignment of the 17 imino protons, which participate in Watson-Crick hydrogen bonding, and of the eight adenine H2 protons, which lie in the minor groove of the double helix, is presented.

Adenine↗
PubMed1984

A new two-disulphide intermediate in the refolding of reduced bovine pancreatic trypsin inhibitor.

The apparently complete refolding of reduced bovine pancreatic trypsin inhibitor (BPTI) is shown to produce a mixture of two species. One of these is native BPTI, but the other lacks the disulphide bond between cysteines 30 and 51. The latter species has a folded conformation very like that of native BPTI, and is oxidized by air to native BPTI on warming in aqueous solution. The two unreactive cysteine thiol groups appear to be buried in the interior of the molecule, which restricts access by reagents that can alkylate them or oxidize them to form the disulphide bond. The implications of this intermediate and its conformation for the understanding of protein folding are discussed.

Amino Acid Sequence↗
PubMed1984

Dynamics of proteins.

Theoretical methods for studying protein motions have been outlined and some examples have been given. Simulation methods have shown that structural fluctuations in proteins cover a wide range in their magnitudes and time scales. Some indication of their biological role has been given, but the most important and interesting applications to protein function have to be made.

Cytochrome c Group↗
PubMed1984

Diffusion-collision model for the folding kinetics of the lambda-repressor operator-binding domain.

The operator-binding domain of the lambda-repressor contains five alpha-helices and an extended N-terminal arm in the crystal structure determined by Pabo and Lewis reported in Nature 298, 443, 1982 (1). The four helices form a "box" enclosing a hydrophobic core with the fifth helix interacting with the equivalent helix in a dimer. With a small number of well-defined secondary structure elements (microdomains), the repressor is well suited for an analysis of its folding pathways and kinetics by use of the diffusion-collision model. In this paper, the basic elements of the model appropriate to a several microdomain protein are formulated and applied to a set of folding pathways consistent with the crystal structure of the operator-binding domain. The overall kinetics, as well as the time-dependence of intermediate states are determined as a function of the microdomain stability parameter.

Binding Sites↗