Search PubMed⌕ Search

Biomedical subjects

A E Mark

Publications and source records attributed to A E Mark.

34 records · Page 2Linked to original sources
PubMed1995

Computational approaches to study protein unfolding: hen egg white lysozyme as a case study.

Four methods are compared to drive the unfolding of a protein: (1) high temperature (T-run), (2) high pressure (P-run), (3) by imposing a gradual increase in the mean radius of the protein using a penalty function added to the physical interaction function (F-run, radial force driven unfolding), and (4) by weak coupling of the difference between the temperature of the radially outward moving atoms and the radially inward moving atoms to an external temperature bath (K-run, kinetic energy driven unfolding). The characteristic features of the four unfolding pathways are analyzed in order to detect distortions due to the size or the type of the applied perturbation, as well as the features that are common to all of them. Hen egg white lysozyme is used as a test system. The simulations are analyzed and compared to experimental data like 1H-NMR amide proton exchange-folding competition, heat capacity, and compressibility measurements.

Amides↗
PubMed1995

Rapid non-empirical approaches for estimating relative binding free energies.

Rapid non-empirical methods for estimating binding free energies are reviewed. A novel approach based on the application of the free energy perturbation formula to a biased ensemble is presented. Preliminary results demonstrating the applicability of this approach in protein systems are shown and the potential of this method in structure-based drug design is discussed.

Models, Theoretical↗
PubMed1994

Decomposition of the free energy of a system in terms of specific interactions. Implications for theoretical and experimental studies.

Recently, a number of methods have been proposed that are designed to extract contributions to the change in free energy associated with a given perturbation or mutation of a protein originating from specific residue-residue or atom-atom interactions, both based on theoretical calculations and on experimental data. We caution here that detailed analysis based on these methods is unreliable. It is demonstrated, both from first principles using statistical mechanics and by way of example, that in a general case a meaningful decomposition of the free energy in terms of specific residue-residue or atom-atom interactions is not possible.

Azurin↗
PubMed1994

Investigation of shape variations in the antibody binding site by molecular dynamics computer simulation.

Molecular dynamics simulations have been used to investigate the flexibility and variations in the shape of the binding site of an antibody against human Rhinovirus serotype 2 (HRV2) and its complex with a 15 amino acid oligopeptide, the structure of which has been recently determined by X-ray crystallography. During the simulation of the unbound antibody the binding site, defined in terms of the hypervariable regions or complementarity determining regions (CDRs), shows significant fluctuations in shape. For the complex such variations in the shape of the binding site were reduced. The largest fluctuations in the unbound antibody occurred within the CDR-H3. The largest differences between the bound and unbound crystal structures are also associated with CDR-H3. The relative displacements of the loops have been analysed in terms of internal distortions, rigid body motions of the loops and changes with respect to the framework regions. The degree to which the motions of the loops are correlated and the variation in the volume of the binding pocket during the simulation have also been examined.

Amino Acid Sequence↗
PubMed1993

Solvent-dependent conformation and hydrogen-bonding capacity of cyclosporin A: evidence from partition coefficients and molecular dynamics simulations.

The partition coefficient of cyclosporin A (CsA) was measured in octanol/water and heptane/water by centrifugal partition chromatography. By comparison with results from model compounds, it was deduced that the hydrogen-bonding capacity of CsA changed dramatically from an apolar solvent (where it is internally H-bonded) to polar solvents (where it exposes its H-bonding groups to the solvent). Molecular dynamics simulations in water and CCl4 support the suggestion that CsA undergoes a solvent-dependent conformational changes and that the interconversion process is slow on the molecular dynamics time scale.

Amino Acid Isomerases↗
PubMed1993

An approximate but efficient method to calculate free energy trends by computer simulation: application to dihydrofolate reductase-inhibitor complexes.

Derivatives of free energy differences have been calculated by molecular dynamics techniques. The systems under study were ternary complexes of Trimethoprim (TMP) with dihydrofolate reductases of E. coli and chicken liver, containing the cofactor NADPH. Derivatives are taken with respect to modification of TMP, with emphasis on altering the 3-, 4- and 5-substituents of the phenyl ring. A linear approximation allows the encompassing of a whole set of modifications in a single simulation, as opposed to a full perturbation calculation, which requires a separate simulation for each modification. In the case considered here, the proposed technique requires a factor of 1000 less computing effort than a full free energy perturbation calculation. For the linear approximation to yield a significant result, one has to find ways of choosing the perturbation evolution, such that the initial trend mirrors the full calculation. The generation of new atoms requires a careful treatment of the singular terms in the non-bonded interaction. The result can be represented by maps of the changed molecule, which indicate whether complex formation is favoured under movement of partial charges and change in atom polarizabilities. Comparison with experimental measurements of inhibition constants reveals fair agreement in the range of values covered. However, detailed comparison fails to show a significant correlation. Possible reasons for the most pronounced deviations are given.

Animals↗
PubMed1993

Can the stability of protein mutants be predicted by free energy calculations?

The use of free energy simulation techniques in the study of protein stability is critically evaluated. Results from two simulations of the thermostability mutation Asn218 to Ser218 in Subtilisin are presented. It is shown that components of the free energy change can be highly sensitive to the computational details of the simulation leading to the conclusion that free energy calculations cannot currently be used to reliably predict protein stability. The different factors that undermine the reliability are discussed.

Asparagine↗
PubMed1992

Prediction of the activity and stability effects of site-directed mutagenesis on a protein core.

When theoretical methods are used to predict the properties of a given system, such as the effects of the substitution of a specific amino acid on the activity or stability of a protein as a whole, the accuracy of the prediction is directly dependent on the validity of the underlying model. A common error, however, is to attempt to improve a basically crude model by performing one aspect of the calculation in a rigorous manner. The accuracy of the model as a whole will remain limited by the crudest approximation or weakest assumption. To demonstrate the principle that nothing can be gained by performing extensive calculations using a basically crude underlying model we compare the predictive power of three models in relation to activity and stability data for 78 triple-site sequence variants of the lambda-repressor protein. This system has recently been analysed in terms of a conceptionally simple, but computationally elaborate model for the prediction of the energy of a protein in which amino acid residues in the core of the protein have been mutated. We show that comparable, if not better agreement with the experimental data can be reached using either of two much simpler models, based on straightforward structural considerations, which do not require elaborate calculations on a computer.

Computer Simulation↗
PubMed1992

Simulation of the thermal denaturation of hen egg white lysozyme: trapping the molten globule state.

In the study of protein folding, much attention has focused on the characterization of folding intermediates. We report here molecular dynamics simulations in which the initial stages of the thermal denaturation of hen egg white lysozyme in aqueous solution are examined in detail. It is found that lysozyme unfolds in a two-stage process with the initial formation a quasi-stable state in which significant rearrangement of the secondary structure takes place. No evidence for distinct folding domains was found. The simulations suggest that the formation of well-defined secondary structure occurs after the initial collapse of the peptide chain and thus tend against the framework model of protein folding.

Animals↗
PubMed1992

On the interpretation of biochemical data by molecular dynamics computer simulation.

The application of computer simulation to molecular systems of biochemical interest is reviewed. It is shown that computer simulation is a tool complementary to experimental methods, which can be used to access atomic details inaccessible to experimental probes. Examples are given in which computer simulation augments the experimental information by providing an atomic picture of high resolution with respect to space, energy or time. The usefulness of a computer simulation largely depends on its quality. The most important factors that limit the accuracy of simulated results are discussed. The accuracy of different simulation studies can differ by orders of magnitude. The accuracy will depend on the type of biomolecular system and process studied. It will also depend on the choice of force field, the simulation set-up and the protocol that is used. A list of quality-determining factors is given, which may be useful when interpreting simulation studies appearing in the literature.

Biochemical Phenomena↗
PubMed1991

Conformational flexibility of aqueous monomeric and dimeric insulin: a molecular dynamics study.

A series of molecular dynamics simulations have been used to investigate the nature of monomeric and dimeric insulin in aqueous solution. It is shown that in the absence of crystal contacts both monomeric and dimeric insulin have a high degree of intrinsic flexibility. Neither of the two monomer conformations of 2Zn crystalline insulin appears to be favored in solution nor is the asymmetry of the crystal dimer reduced in the absence of crystal contacts. A shift is observed in the relative positions of molecules 1 and 2 in the dimer compared with that found in the crystal, which may have consequences for the prediction of the effects of mutants in the monomer-monomer interface designed to alter the self-association properties of insulin.

Amino Acid Sequence↗
PubMed1990

The self-association of zinc-free bovine insulin. Four model patterns and their significance.

Data fitting procedures both with and without correction for non-ideality are applied to experimentally measured concentration distributions of zinc-free insulin obtained over a wide range of experimental conditions (pH 2, 7, and 10, ionic strengths 0.05 and 0.1, t = 25 degrees C and 37 degrees C) using four different model self-association patterns based on known physico-chemical properties of insulin in solution and patterns already in the literature. It is shown that three of these must be considered satisfactory descriptions of the insulin system of equilibria in aqueous solution as judged by critical curve-fitting criteria. The significant differences, as well as common features of these are assessed by comparing the distribution of monomeric and polymeric forms at three insulin concentrations of practical utility, "serum", "pharmacological", and "physicochemical". The value of analysing protein self-association with the aid of explicit equations formulated for a specific model, even when it cannot be demonstrated to be unique, is discussed with particular reference to osmotic pressure measurements made by others on insulins modified by recombinant DNA techniques aimed at making them essentially monomeric.

Animals↗
PubMed1990

The reversible cross-linking of receptors by ligands: theory for the prediction of binding responses.

The general concept of receptor aggregation through the action of a cross-linking ligand is considered, three models being examined in detail. In two, ligand self-interaction leads to the formation of receptor cross-links whereas, in the third, receptor cross-links are formed through a single ligand bridge. Binding equations in closed form are formulated for each case and are shown to predict both concave and convex Scatchard plots in binding studies conducted with such systems. The significant point which emerges is that while each of the systems gives rise to very different types of binding responses, the form of the binding response is always strongly dependent on the total concentration of receptor in the system.

Algorithms↗
PubMed1988

The binding of an indefinitely associating ligand to acceptor: consideration of monovalent ligand species binding to a multivalent acceptor.

Currently available binding theory is extended to incorporate the concept of indefinite self-association of the ligand. Binding equations are formulated in closed form for the case of the binding to a multivalent acceptor of a ligand capable of isodesmically indefinitely self-associating in a "head-to-tail" mode such that each ligand state bears one site capable of interacting with the acceptor. It is shown both mathematically and by way of numerical example that this system will give rise exclusively to binding curves convex to the r-axis in Scatchard format. Thus, the system provides another example of a binding mechanism capable of generating an apparent negatively co-operative binding response.

Binding Sites↗
PubMed1987

The self-association of zinc-free bovine insulin. A single model based on interactions in the crystal that describes the association pattern in solution at pH 2, 7 and 10.

Sedimentation equilibrium studies are used to establish that a new pattern for the self-association of zinc-free insulin in solution is applicable over a wide range of conditions of pH, ionic strength and temperature. In this pattern, which is based on information from the existing literature on the X-ray crystal structure of insulin, the insulin monomer is viewed as having two distinct faces both capable of self-interaction. Sedimentation equilibrium experiments were analysed using expressions formulated for this association pattern that describe the dependence of weight average molecular weight and monomer concentration on total protein concentration. It has thereby been possible to obtain values for the two association constants which govern the system for each set of conditions studied, due allowance having been made for composition dependent non-ideality effects. Furthermore, by relating the pH, temperature and ionic strength dependence of the association constants with properties of various amino acid residues on the surface of the insulin monomer, it has also been possible to assign tentatively each constant to a particular reaction domain.

Animals↗