Search PubMed⌕ Search

Biomedical subjects

M Karplus

Publications and source records attributed to M Karplus.

At least 127 records · Page 7Linked to original sources
PubMed1994

Free energy simulations: the meaning of the individual contributions from a component analysis.

A theoretical analysis is made of the decomposition into contributions from individual interactions of the free energy calculated by thermodynamic integration. It is demonstrated that such a decomposition, often referred to as "component analysis," is meaningful, even though it is a function of the integration path. Moreover, it is shown that the path dependence can be used to determine the relation of the contribution of a given interaction to the state of the system. To illustrate these conclusions, a simple transformation (Cl- to Br- in aqueous solution) is analyzed by use of the Reference Interaction Site Model-Hypernetted Chain Closure integral equation approach; it avoids the calculational difficulties of macromolecular simulation while retaining their conceptual complexity. The difference in the solvation free energy between chloride and bromide is calculated, and the contributions of the Lennard-Jones and electrostatic terms in the potential function are analyzed by the use of suitably chosen integration paths. The model is also used to examine the path dependence of individual contributions to the double free energy differences (delta delta G or delta delta A) that are often employed in free energy simulations of biological systems. The alchemical path, as contrasted with the experimental path, is shown to be appropriate for interpreting the effects of mutations on ligand binding and protein stability. The formulation is used to obtain a better understanding of the success of the Poisson-Boltzmann continuum approach for determining the solvation properties of polar and ionic systems.

Bromides↗
PubMed1994

Gestational diabetes among Bedouins in southern Israel: comparison of prevalence and neonatal outcomes with the Jewish population.

Differences in the prevalence of gestational diabetes mellitus (GDM) have recently been reported between various ethnic populations. In the Negev region of Israel, a universal free screening programme for GDM was implemented in 1985. Between 1 March 1987 and 31 July 1988 11,003 deliveries occurred at the Soroka Medical Center, which provides free delivery and postnatal care to the whole Jewish and Bedouin population of the region. GDM was found in 5.7% of Jewish and in 2.4% of Bedouin women (odds ratio, 2.3, 95% confidence interval (CI) 1.8-2.9; P < 0.0001). Ethnicity was unrelated to maternal outcome, perinatal mortality or to any of the examined morbidity conditions of the newborn. The incidence of major congenital malformations was significantly higher in Jewish than in Bedouin infants of GDM women (Fisher's Exact Test, P < 0.03). Conversely, Jewish infants had fewer minor congenital anomalies (odds ratio 0.26, 95% CI 0.09-0.73). In a multivariate logistic regression model, gestational age, mode of delivery and insulin requirement during pregnancy were the only factors independently associated with neonatal morbidity. The results of this study suggest that in our health care system, among women with GDM ethnicity is not associated with an excess of unfavourable maternal or infant outcomes.

Adult↗
PubMed1994

Conformational analysis of the backbone-dependent rotamer preferences of protein sidechains.

Amino acids have sidechain rotamer preferences dependent on the backbone dihedral angles phi and psi. These preferences provide a method for rapid structure prediction which is a significant improvement over backbone-independent rotamer libraries. We demonstrate here that simple arguments based on conformational analysis can account for many of the features of the observed backbone dependence of the sidechain rotamers. Steric repulsions corresponding to the 'butane' and 'syn-pentane' effects make certain conformers rare, as has been observed experimentally.

Amino Acids↗
PubMed1993

A mechanism for rotamase catalysis by the FK506 binding protein (FKBP).

A detailed mechanism for the catalysis of prolyl isomerization by the rotamase enzyme FKBP is proposed on the basis of a model constructed from the known structure of the FK506/FKBP complex. The model substrate is bound as a type VIa proline turn with the ends exposed to permit longer polypeptide chains (e.g., protein loops) to act as substrates. An ab initio potential for the isomerized imide bond is combined with a molecular mechanics representation of the rest of the system to calculate the reaction path. The resulting activation energy for the enzymatic cis-->trans isomerization is equal to about 6 kcal/mol, in good agreement with experiment. The lowering of the barrier relative to the solution value of 19 kcal/mol is found to arise from a combination of desolvation of the imide carbonyl, ground-state destabilization, substrate autocatalysis, and preferential transition-state binding. Minimal rearrangements are required in the enzyme and the substrate along the reaction path. The enzyme residues that participate in catalysis agree with the available mutation data. The type VIa turn model corresponds to a sequence-specific structural motif commonly found on the surface of proteins. It is likely to have a role in the formation of protein complexes with FKBP-like domains that function as foldases or chaperones.

Carrier Proteins↗
PubMed1993

Effective coupling in biological electron transfer: exponential or complex distance dependence?

Calculations for a simple model electron transfer system and tuna cytochrome c demonstrate a dichotomy in the distance dependence of the effective coupling. In one regime, the effective coupling varies exponentially with distance and depends primarily on the average properties of the bridging material; in the other regime, the effective coupling has a complex distance dependence and is more sensitive to the details of the bridging material. Experiments and theory indicate that both regimes may occur in biological systems, providing a perspective on a recent controversy over the nature of the distance dependence.

Animals↗
PubMed1993

Multiple copy simultaneous search and construction of ligands in binding sites: application to inhibitors of HIV-1 aspartic proteinase.

Rational ligand design is a complex problem that can be divided into three parts: the search for optimal positions and orientations of functional groups in the binding site, the connection of such positions to form candidate ligands, and the estimation of their binding constants. Approaches for addressing the first two parts of the problem are described in the present work. They are applied to the construction of peptide ligands in the binding site of the human immunodeficiency virus 1 (HIV-1) proteinase. The primary objective is to test the method by comparison of the results with the MVT-101 complex structure for which coordinates are available; the results obtained with the liganded and unliganded proteinase structure are used to examine the utility of the latter for binding studies. A secondary objective is to show how to find new inhibitor candidates. The multiple copy simultaneous search (MCSS) method is utilized to search for optimal positions and orientations of a set of functional groups. For peptide ligands, functional groups corresponding to the protein main chain (N-methylacetamide) and to protein side chains (e.g., methanol, ethyl guanidinium) are used. The resulting N-methylacetamide minima are connected to form hexapeptide main chains with a simple pseudoenergy function that permits a complete search of all possible ways of connecting the minima. Side chains are added to the main-chain candidates by application of the same pseudoenergy function to the appropriate functional group minima. A set of 15 hexapeptides with the sequence of MVT-101 is then minimized by a Monte Carlo scheme, which allows for escape from local minima. Comparison of the MCSS results with the structure of MVT-101 in the HIV-1 binding site showed that all of its functional group positions correspond (within 2.4 A) to some (usually more than one) MCSS minima. There were also many other low-energy MCSS minima which do not appear in any known inhibitors, e.g., methyl ammonium minima in the neighborhood of the catalytic aspartates. Among the 15 lowest minima are seven hexapeptides with the same main-chain orientation as the one found by X-ray crystallography for the inhibitor MVT-101 in the binding site and eight with the main chain oriented in the opposite direction; the latter tend to be more stable. [Addendum: These results are in agreement with recent high-resolution crystallographic data provided after the study was completed.(ABSTRACT TRUNCATED AT 400 WORDS)

Aspartic Acid Endopeptidases↗
PubMed1993

Nonexponential relaxation after ligand dissociation from myoglobin: a molecular dynamics simulation.

Molecular dynamics simulations of myoglobin after ligand photodissociation show that the out-of-plane motion of the heme iron has a rapid subpicosecond phase followed by a slower nonexponential process involving more global protein relaxation. Individual trajectories show rather different behavior, suggesting there is an inhomogeneous component to the relaxation. The calculated time dependence of the iron motion over 100 ps is in excellent agreement with the frequency shift of band III of the heme group [see Lim, M., Jackson, T. A. & Anfinrud, P. A. (1993) Proc. Natl. Acad. Sci. USA 90, 5801-5804]. If that the barrier to rebinding depends on the out-of-plane iron position, the time dependence obtained from the simulation can explain the nonexponential room-temperature geminate recombination of NO.

Kinetics↗
PubMed1993

Three-dimensional models of four mouse mast cell chymases. Identification of proteoglycan binding regions and protease-specific antigenic epitopes.

Mouse mast cell protease (mMCP) 1, mMCP-2, mMCP-4, and mMCP-5 are serine proteases which are predicted to have chymotryptic specificity (chymases). They are bound to negatively charged heparin or chondroitin sulfate proteoglycans and are stored in secretory granules. Three-dimensional (3D) models of these four proteases were constructed with a comparative molecular modeling technique based on satisfaction of spatial constraints. The models were used to predict immunogenic epitopes and surface regions that are likely to interact with proteoglycans. Nine potential antigenic segments in the four chymases were identified on the basis of solvent accessibility, protrusion, flexibility, and sequence variability. These segments are suitable epitopes for preparation of protease-specific antipeptide immunoglobulin. Two regions with net charges ranging from +6 to +10 at neutral pH were found on the surfaces of mMCP-4 and mMCP-5. The two regions are located far from the substrate binding cleft at diametrically opposite ends of the folded proteases. A strong positive electrostatic potential surrounds the two regions. Thus, they are good candidates for binding sites that interact with heparin proteoglycan in the granules of serosal mast cells. In contrast, mMCP-1 and mMCP-2, which are present in granules of mucosal mast cells that contain chondroitin sulfate, lack one of these regions and have a lower charge density in the other. The differences between the 3D models provide a structural basis for the selective localization of specific chymases within mouse mast cells that contain different proteoglycans.

Amino Acid Sequence↗
PubMed1993

Backbone-dependent rotamer library for proteins. Application to side-chain prediction.

A backbone-dependent rotamer library for amino acid side-chains is developed and used for constructing protein side-chain conformations from the main-chain co-ordinates. The rotamer library is obtained from 132 protein chains in the Brookhaven Protein Database. A grid of 20 degrees by 20 degrees blocks for the main-chain angles phi, psi is used in the rotamer library. Significant correlations are found between side-chain dihedral angle probabilities and backbone phi, psi values. These probabilities are used to place the side-chains on the known backbone in test applications for six proteins for which high-resolution crystal structures are available. A minimization scheme is used to reorient side-chains that conflict with the backbone or other side-chains after the initial placement. The initial placement yields 59% of both chi 1 and chi 2 values in the correct position (to within 40 degrees) for thermolysin to 81% for crambin. After refinement the values range from 61% (lysozyme) to 89% (crambin). It is evident from the results that a single protein does not adequately test a prediction scheme. The computation time required by the method scales linearly with the number of side-chains. An initial prediction from the library takes only a few seconds of computer time, while the iterative refinement takes on the order of hours. The method is automated and can easily be applied to aid experimental side-chain determinations and homology modeling. The high degree of correlation between backbone and side-chain conformations may introduce a simplification in the protein folding process by reducing the available conformational space.

Amino Acid Sequence↗
PubMed1993

The contribution of cross-links to protein stability: a normal mode analysis of the configurational entropy of the native state.

The vibrational entropy of native BPTI, with three disulfide bonds, was determined by use of normal mode calculations and compared with that of folded variants having either one less disulfide bond or lacking a peptide bond at the trypsin-reactive site. Favorable contributions to the free energy of 2.5-5.1 kcal/mol at 300 K were calculated for the reduction of disulfide bonds in the folded state, whereas no favorable contribution was found for the hydrolysis of the peptide bond cleaved by trypsin. This is on the order of the effect of disulfides in the unfolded state. The implications of these results for the stabilization of a folded protein by the introduction of cross-links are discussed.

Aprotinin↗
PubMed1993

Active site dynamics of acyl-chymotrypsin.

The motions of water molecules, the acyl moiety, the catalytic triad, and the oxyanion binding site of acyl-chymotrypsin were studied by means of a stochastic boundary molecular dynamics simulation. A water molecule that could provide the nucleophilic OH- for the deacylation stage of the catalysis was found to be trapped between the imidazole ring of His-57 and the carbonyl carbon of the acyl group. It makes a hydrogen bond with the N epsilon 2 of His-57 and is held in place through a network of hydrogen-bonded water molecules in the active site. The water molecule was found as close as 2.8 A to the carbonyl carbon. This appears to be due to the constraints imposed by nonbonded interaction in the active site. Configurations were found in which one hydrogen of the trapped water shared a bifurcated hydrogen bond with His-57-N epsilon 2 and Ser-195-O gamma, with the water oxygen very close to the carbonyl carbon. The existence of such a water molecule suggests that large movement of the His-57 imidazole ring between positions suitable for providing general-base catalyzed assistance and for providing general-acid catalyzed assistance may not be required during the reaction. The simulation indicates that the side chains of residues involved in catalysis (i.e., His-57, Ser-195, and Asp-102) are significantly less flexible than other side chains in the protein. The 40% reduction in rms fluctuations is consistent with a comparable reduction calculated from the temperature factors obtained in the X-ray crystallographic data of gamma-chymotrypsin. The greater rigidity of active site residues seems to result from interconnected hydrogen bonding networks among the residues and between the residues and the solvent water in the active site.

Acylation↗
PubMed1993

The nature of the ion binding interactions in EF-hand peptide analogs: free energy simulation of Asp to Asn mutations.

The binding of the La3+ ion to a tridecapeptide, which is a model for the EF-hand in calcium-binding proteins, is studied in solution by free energy simulations. The calculations analyze the effect on the La3+ ion binding of the mutation of Asp to Asn for side chains that interact directly with the ion. The results are compared with the measurements of Marsden,B.J., Hodges,R.S. and Sykes,B.D. (1989) Biochemistry, 28, 8839, on the same system. They found that the Asp to Asn mutation has only a small effect on the binding; the observed differences in the free energies on changing one Asp to an Asn are between -0.3 and 1.8 kcal/mol. This result is analyzed by alchemical simulations for the tridecapeptide in the bound (loop) structure and free (extended) form. The free energy changes due to the mutation of an Asp to an Asn are large and positive for both the bound and free forms. However, since the values of the free energy changes are calculated to be similar in the two forms, the difference in the binding free energy of Asp and Asn peptides is found to be small, in agreement with experiment. By use of thermodynamic integration, the various contributions to the free energy changes are estimated. In the complexed form, the Asp to Asn mutation is favored by the reduction in the repulsive interaction with other charged residues of the peptide; it is disfavored by the reduction of the stabilization of the ion and the surrounding water has a small effect.(ABSTRACT TRUNCATED AT 250 WORDS)

Amino Acid Sequence↗
PubMed1993

Erythrocyte Li+i-Na+o countertransport in the newborn.

The present study investigates the erythrocyte Li+i-Na+o countertransport (CT) in the newborn. Erythrocytes of term infants exhibit no or very low CT activity at birth (0.025 +/- 0.011 mmol Li/l RBC.h). CT rates increase gradually and reach adult values within the 4th-7th day of life (0.225 +/- 0.018 mmol Li/l RBC.h). CT activity in healthy preterm infants at birth is not significantly different from that in term infants (0.014 +/- 0.007 mmol Li/l RBC.h). Postnatal maturation of CT activity in preterm infants is similar to that found in term infants and is independent of gestational age. However, CT rates at birth of premature infants subjected to prenatal stress are significantly higher (0.082 +/- 0.017 mmol Li/l RBC.h). It is concluded that erythrocyte CT activity in the newborn reaches adult values within the first week of life, apparently reflecting the postnatal maturation of the CT system in the cell membrane. We suggest that prenatal stress induces the synthesis of the CT protein in utero.

Antiporters↗
PubMed1992

Chromophore-protein interactions and the function of the photosynthetic reaction center: a molecular dynamics study.

The coupling between electron transfer and protein structure and dynamics in the photosynthetic reaction center of Rhodopseudomonas viridis is investigated. For this purpose molecular dynamics simulations of the essential portions (a segment of 5797 atoms) of this protein complex have been carried out. Electron transfer in the primary event is modeled by altering the charge distributions of the chromophores according to quantum chemical calculations. The simulations show (i) that fluctuations of the protein matrix, which are coupled electrostatically to electron transfer, play an important role in controlling the electron transfer rates and (ii) that the protein matrix stabilizes the separated electron pair state through rapid (200 fs) and temperature-independent dielectric relaxation. The photosynthetic reaction center resembles a polar liquid in that the internal motions of the whole protein complex, rather than only those of specific side groups, contribute to i and ii. The solvent reorganization energy is about 4.5 kcal/mol. The simulations indicate that rather small structural rearrangements and changes in motional amplitudes accompany the primary electron transfer.

Bacterial Chromatophores↗
PubMed1992

Beta-sheet coil transitions in a simple polypeptide model.

A simplified model of a polypeptide chain is used to study the dynamics of the beta-sheet-coil transition. Each amino acid residue is treated as a single quasiparticle in an effective potential that approximates the potential of mean force in solution. The model is used to study the equilibrium and dynamic aspects of the sheet-coil transition. Systems studied include ones with both strands free to move (two-strand sheet), and ones with either strand fixed in position (multistrand sheet). The equilibrium properties examined include sheet-coil equilibrium constants and their dependence on chain position. Dynamic properties are investigated by a stochastic simulation of the Brownian motion of the chain in its solvent surroundings. Time histories of the dihedral angles and residue-residue cross-strand distances are used to study the behavior of the sheet structure. Auto- and cross-correlation functions are calculated from the time histories with relaxation times of tens to hundreds of picoseconds. Sheet-coil rate constants of tens of ns-1 were found for the fixed strand cases.

Models, Molecular↗
PubMed1992

Congenital nasal stenosis in newborn infants.

Five newborn infants with evidence of nasal obstruction were shown to have congenital nasal stenosis. Conservative treatment, including temporary nasopharyngeal intubation in three, eventually resulted in symptomatic relief, so that surgery could be avoided.

Constriction, Pathologic↗