Search PubMed⌕ Search

Biomedical subjects

M Karplus

Publications and source records attributed to M Karplus.

At least 217 records · Page 12Linked to original sources

X-ray structure and refinement of carbon-monoxy (Fe II)-myoglobin at 1.5 A resolution.

The structure of carbon-monoxy (Fe II) myoglobin at 260 K has been solved at a resolution of 1.5 A by X-ray diffraction and a model refined against the X-ray data by restrained least-squares. The CO ligand is disordered and distorted from the linear conformation seen in model compounds. At least two conformations, with Fe--C--O angles of 140 degrees and 120 degrees, are required to model the system. The heme pocket is significantly larger than in deoxy-myoglobin because the distal residues have relaxed around the ligand; the largest displacement occurs for the distal histidine side-chain, which moves more than 1.4 A on ligand binding. The side-chain of Arg45 (CD3) is disordered and apparently exists in two equally populated conformations. One of these does not block the motion of the distal histidine out of the binding pocket, suggesting a mechanism for ligand entry. The heme group is planar (root-mean-square deviation from planarity is 0.08 A) with no doming of the pyrrole groups. The Fe--N epsilon 2 (His93) bond length is 2.2 A and the Fe--C bond length in the CO complex is 1.9 A. The iron is the least-squares plane of the heme, and this leads to the proximal histidine moving by 0.4 A relative to its position in deoxy-myoglobin. This shift correlates with a global structural change, with the proximal part of the molecule translated towards the heme plane.

Binding Sites↗

Application of molecular dynamics with interproton distance restraints to three-dimensional protein structure determination. A model study of crambin.

The applicability of restrained molecular dynamics for the determination of three-dimensional protein structures on the basis of short interproton distances (less than 4 A) that can be realistically determined from nuclear magnetic resonance measurements in solution is assessed. The model system used is the 1.2 A resolution crystal structure of the 46 residue protein crambin, from which a set of 240 approximate distance restraints, divided into three ranges (2.5 +/- 0.5, 3.0+0.5(-1.0) and 4 +/- 1 A), is derived. This interproton distance set comprises 159 short-range ([i-j] less than or equal to 5) and 56 ([i-j] greater than 5) long-range inter-residue distances and 25 intra-residue distances. Restrained molecular dynamics are carried out using a number of different protocols starting from two initial structures: a completely extended beta-strand; and an extended structure with two alpha-helices in the same positions as in the crystal structure (residues 7 to 19, and 23 to 30) and all other residues in the form of extended beta-strands. The root-mean-square (r.m.s.) atomic differences between these two initial structures and the crystal structure are 43 A and 23 A, respectively. It is shown that, provided protocols are used that permit the secondary structure elements to form at least partially prior to folding into a tertiary structure, convergence to the correct final structure, both globally and locally, is achieved. The r.m.s. atomic differences between the converged restrained dynamics structures and the crystal structure range from 1.5 to 2.2 A for the backbone atoms and from 2.0 to 2.8 A for all atoms. The r.m.s. atomic difference between the X-ray structure and the structure obtained by first averaging the co-ordinates of the converged restrained dynamics structures is even smaller: 1.0 A for the backbone atoms and 1.6 A for all atoms. These results provide a measure with which to judge future experimental results on proteins whose crystal structures are unknown. In addition, from an examination of the dynamics trajectories, it is shown that the convergence pathways followed by the various simulations are different.

Hydrogen Bonding↗

Molecular dynamics simulations of native and substrate-bound lysozyme. A study of the average structures and atomic fluctuations.

Molecular dynamics simulations of hen egg-white lysozyme in the free and substrate-bound states are reported and the nature of the average structures and atomic fluctuations are analyzed. Crystallographic water molecules of structural importance, as determined by hydrogen-bonding, were included in the simulations. Comparisons are made between the dynamics and the X-ray results for the atomic positions, the main-chain and side-chain dihedral angles, and the hydrogen-bonding geometry. Improvements over earlier simulations in the potential energy function and methodology resulted in stable trajectories with the C alpha co-ordinates within 1.5 A of the starting X-ray structure. Structural features analyzed in the simulations agreed well with the X-ray results except for some surface residues. The Asx chi 2 dihedral distribution and the geometry of hydrogen bonding at reverse turns show differences; possible causes are discussed. The relation between the magnitudes and time-scales of the residue fluctuations and secondary structural features, such as helices beta-sheets and coiled loops, is examined. Significant differences in the residue mobilities between the simulations of the free and substrate-bound states were found in a region of the enzyme that is in direct contact with the substrate and in a region that is distant from the active-site cleft. The dynamic behavior of the structural water molecules is analyzed by examining the correlation between the fluctuations of the water oxygens and the lysozyme heavy-atoms to which they are hydrogen-bonded.

Animals↗

Effect of anisotropy and anharmonicity on protein crystallographic refinement. An evaluation by molecular dynamics.

Molecular dynamics simulations are employed to determine the errors introduced by anharmonicity and anisotropy in the structure and temperature factors obtained for proteins by refinement of X-ray diffraction data. Simulations (25 ps and 300 ps) of metmyoglobin are used to generate time-averaged diffraction data at 1.5 A resolution. The crystallographic restrained-parameter least-squares refinement program PROLSQ is used to refine models against these simulated data. The resulting atomic positions and isotropic temperature factors are compared with the average structure and fluctuations calculated directly from the simulations. It is found that significant errors in the atomic positions and fluctuations are introduced by the refinement, and that the errors increase with the magnitude of the atomic fluctuations. Of particular interest is the fact that the refinement generally underestimates the atomic motions. Moreover, while the actual fluctuations go up to a mean-square value of about 5 A2, the X-ray results never go above approximately 2 A2. This systematic deviation in the motional parameters appears to be due to the use of a single-site isotropic model for the atomic fluctuations. Many atoms have multiple peaks in their probability distribution functions. For some atoms, the multiple peaks are seen in difference electron density maps and it is possible to include these in the refinement as disordered residues. However, for most atoms the refinement fits only one peak and neglects the rest, leading to the observed errors in position and temperature factor. The use of strict stereochemical restraints is inconsistent with the average dynamical structure; nevertheless, refinement with tight restraints results in structures that are comparable to those obtained with loose restraints and better than those obtained with no restraints. The results support the use of tight stereochemical restraints, but indicate that restraints on the variation of temperature factors are too restrictive.

Amino Acid Sequence↗

Structure refinement of oligonucleotides by molecular dynamics with nuclear Overhauser effect interproton distance restraints: application to 5' d(C-G-T-A-C-G)2.

The solution structure of the self-complementary DNA hexamer 5' d(C-G-T-A-C-G)2 is refined by restrained molecular dynamics in which 192 interproton distances, determined from pre-steady-state nuclear Overhauser enhancement measurements, are incorporated into the total energy of the system in the form of effective potentials. First the method is tested by applying an idealized set of distance restraints taken from classical B-DNA to a simulation starting off from A-DNA and vice versa. It is shown that in both cases the expected transition between A- and B-DNA occurs. Second, a set of restrained molecular dynamics calculations is carried out starting from both A- and B-DNA with the experimental interproton distances for 5' d(C-G-T-A-C-G)2 as restraints. Convergence to the same B-type structure is achieved with the interproton distances equal to the measured values within experimental error. The root-mean-square atomic difference between the two average restrained dynamics structures (less than 1 A) is approximately the same as the root-mean-square fluctuations of the atoms.

DNA↗

The three-dimensional structure of alpha1-purothionin in solution: combined use of nuclear magnetic resonance, distance geometry and restrained molecular dynamics.

The determination of the three-dimensional solution structure of alpha1-purothionin using a combination of metric matrix distance geometry and restrained molecular dynamics calculations based on n.m.r. data is presented. The experimental data comprise complete sequence-specific proton resonance assignments, a set of 310 approximate interproton distance restraints derived from nuclear Overhauser effects, 27 Ø backbone torsion angle restraints derived from vicinal coupling constants, 4 distance restraints from hydrogen bonds and 12 distance restraints from disulphide bridges. The average atomic rms difference between the final nine converged structures and the mean structure obtained by averaging their coordinates is 1.5 +/- 0.1 å for the backbone atoms and 2.0 +/- 0.1 å for all atoms. The overall shape of alpha1-purothionin is that of the capital letter L, similar to that of crambin, with the longer arm comprising two approximately parallel alpha-helices and the shorter arm a strand and a mini anti-parallel beta sheet.

Journal Article↗

The impact of prenatal care on the outcome of pregnancy.

In a 12-months' prospective project all women delivering at the Soroka Medical Center, Beer-Sheba, Israel, were studied with regard to the type of prenatal care (P.C.) and the outcome of pregnancy. The Soroka Medical Center is the sole medical facility providing obstetric and neonatal services to the 300,000 inhabitants of the region. Prenatal care is delivered by a uniform network of 70 community-based stations. The availability, type and quality of services were uniform for all women and were not altered throughout the study period. The project encompassed 7308 deliveries. 2154 Bedouin women were excluded in order to avoid possible bias due to cultural and genetic characteristics. Data regarding 5154 Jewish women were analysed. Perinatal mortality was inversely proportional to the number of prenatal contacts. The uncorrected mortality rates were 12.7% in women entirely lacking P.C. and 6.2%, 1.9% and 1.0% in patients who had 1-6, 7-10 and 11 or more prenatal contacts respectively. The low-birthweight rate was significantly increased in women lacking prenatal care (22.8%) or having rudimentary care (17.4%) as compared to those who had 11 or more prenatal contacts with medical personnel (5.9%). Prenatal care reduced neonatal morbidity as expressed by the length of hospitalization, the frequency of infants with multiple diagnoses and the incidence of specific pathologies such as respiratory distress syndrome, light for dates and asphyxia. The incidence of some complications of pregnancy (abruptio placentae, premature labor, PROM) was significantly increased in women lacking P.C. or having inadequate P.C. (1-6 prenatal contacts). Moreover, lack of prenatal care was clearly connected with high-risk delivery, thus increasing the danger to the baby. It is suggested that many adverse effects of various socio-economic, genetic and general health factors may be diminished by proper prenatal care coupled with adequate obstetric and neonatal services.

Adult↗

Three-dimensional structure of proteins determined by molecular dynamics with interproton distance restraints: application to crambin.

Model calculations are performed to evaluate the utility of molecular dynamics with NMR interproton distance restraints for determining the three-dimensional structure of proteins. The system used for testing the method is the 1.5-A resolution crystal structure of crambin (a protein of 46 residues) from which a set of 240 approximate interproton distances of less than 4 A are derived. The convergence properties of the method are examined by using different dynamics protocols and starting from two initial structures; one is a completely extended beta-strand, and the other has residues 7-19 and 23-30 in the form of alpha-helices (as in the crystal structure) with the remaining residues in the form of extended beta-strands. In both cases global and local convergence to the correct final structure is achieved with rms atomic differences between the restrained dynamics structures and the crystal structure of 1.5-2.1 A for the backbone atoms and 2.1-2.8 A for all atoms; the averaged structure has backbone and all atom rms deviations of 1.3 and 1.9 A, respectively. Further, it is shown that a restrained dynamics structure with significantly larger deviations (i.e., 5.7 A for the backbone atoms) can be characterized as incorrect, independent of a knowledge of the crystal structure.

Magnetic Resonance Spectroscopy↗

Aplasia cutis congenita in one of monozygotic twins.

Aplasia cutis congenita is an uncommon disorder characterized by a patchy defect of the scalp and other areas of the skin of the newborn. Several familial cases have been reported, including a set of identical twin girls. The disease was unique in our patient, one monozygotic twin.

Diseases in Twins↗

Domain association in immunoglobulin molecules. The packing of variable domains.

We have analyzed the structure of the interface between VL and VH domains in three immunoglobulin fragments: Fab KOL, Fab NEW and Fab MCPC 603. About 1800 A2 of protein surface is buried between the domains. Approximately three quarters of this interface is formed by the packing of the VL and VH beta-sheets in the conserved "framework" and one quarter from contacts between the hypervariable regions. The beta-sheets that form the interface have edge strands that are strongly twisted (coiled) by beta-bulges. As a result, the edge strands fold back over their own beta-sheet at two diagonally opposite corners. When the VL and VH domains pack together, residues from these edge strands form the central part of the interface and give what we call a three-layer packing; i.e. there is a third layer composed of side-chains inserted between the two backbone side-chain layers that are usually in contact. This three-layer packing is different from previously described beta-sheet packings. The 12 residues that form the central part of the three observed VL-VH packings are absolutely or very strongly conserved in all immunoglobulin sequences. This strongly suggests that the structure described here is a general model for the association of VL and VH domains and that the three-layer packing plays a central role in forming the antibody combining site.

Amino Acid Sequence↗