PubMed Health⌕ Search

Biomedical subjects

A E Mark

Publications and source records attributed to A E Mark.

At least 19 recordsLinked to original sources
PubMed2001

Entropy calculations on a reversibly folding peptide: changes in solute free energy cannot explain folding behavior.

The configurational entropy of a beta-heptapeptide in solution at four different temperatures is calculated. The contributions of the backbone and of the side-chain atoms to the total peptide entropy are analyzed separately and the effective contribution to the entropy arising from correlations between these terms determined. The correlation between the backbone and side-chain atoms amounts to about 17% and is rather insensitive to the temperature. The correlation of motion within the backbone and within side-chains is much larger and decreases with temperature. As the peptide reversibly folds at higher temperatures, its change in entropy and enthalpy upon folding is analyzed. The change in entropy and enthalpy upon folding of the peptide alone cannot account for the observed change in free energy on folding of the peptide in solution. Enthalpic and entropic contributions of the solvent thus also play a key role. Proteins 2001;43:45-56.

Entropy↗
PubMed2000

On the temperature and pressure dependence of a range of properties of a type of water model commonly used in high-temperature protein unfolding simulations.

Molecular dynamics simulations of protein folding and unfolding are often carried out at temperatures (400-600 K) that are much higher than physiological or room temperature to speed up the (un)folding process. Use of such high temperatures changes both the protein and solvent properties considerably, compared to physiological or room temperature. Water models designed for use in conjunction with biomolecules, such as the simple point charge (SPC) model, have generally been calibrated at room temperature and pressure. To determine the distortive effect of high simulation temperatures on the behavior of such "room temperature" water models, the structural, dynamic, and thermodynamic properties of the much-used SPC water model are investigated in the temperature range from 300 to 500 K. Both constant pressure and constant volume conditions, as used in protein simulations, were analyzed. We found that all properties analyzed change markedly with increasing temperature, but no phase transition in this temperature range was observed.

Calibration↗
PubMed1999

The effect of motional averaging on the calculation of NMR-derived structural properties.

The effect of motional averaging when relating structural properties inferred from nuclear magnetic resonance (NMR) experiments to molecular dynamics simulations of peptides is considered. In particular, the effect of changing populations of conformations, the extent of sampling, and the sampling frequency on the estimation of nuclear Overhauser effect (NOE) inter-proton distances, vicinal (3)J-coupling constants, and chemical shifts are investigated. The analysis is based on 50-ns simulations of a beta-heptapeptide in methanol at 298 K, 340 K, 350 K, and 360 K. This peptide undergoes reversible folding and samples a significant proportion of the available conformational space during the simulations, with at 298 K being predominantly folded and at 360 K being predominantly unfolded. The work highlights the fact that when motional averaging is included, NMR data has only limited capacity to distinguish between a single fully folded peptide conformation and various mixtures of folded and unfolded conformations. Proteins 1999;36:542-555.

Computer Simulation↗
PubMed1999

Folding-unfolding thermodynamics of a beta-heptapeptide from equilibrium simulations.

The thermodynamics of folding and unfolding of a beta-heptapeptide in methanol solution has been studied at four different temperatures, 298 K, 340 K, 350 K, and 360 K, by molecular dynamics simulation. At each of these temperatures, the 50-ns simulations were sufficient to generate an equilibrium distribution between a relatively small number of conformations (approximately 10(2)), showing that, even above the melting temperature (approximately 340 K), the peptide does not randomly sample conformational space. The free energy of folding and the free energy difference between pairs of conformations have been calculated from their relative populations. The experimentally determined folded conformation at 298 K, a left-handed 3(1)-helix, is at each of the four temperatures the predominant conformation, with its probability and average lifetime decreasing with increasing temperature. The most common intermediates of folding and unfolding are also the same at the four temperatures. Paths and rates of interconversion between different conformations have been determined. It has been found that folding can occur through multiple pathways, not necessarily downhill in free energy, although the final step involves a reduced number of intermediates.

Models, Molecular↗
PubMed1998

Reversible peptide folding in solution by molecular dynamics simulation.

Long-standing questions on how peptides fold are addressed by the simulation at different temperatures of the reversible folding of a peptide in solution in atomic detail. Molecular dynamics simulations correctly predict the structure that is thermodynamically stable at 298 K, irrespective of the initial peptide conformation. The rate of folding and the free energy of folding at different temperatures are estimated. Although the conformational space potentially accessible to the peptide is extremely large, very few conformers (10(1) to 10(2)) are significantly populated at 20 K above the melting temperature. This implies that the search problem in peptide (or even protein) folding is surmountable using dynamics simulations.

Computer Simulation↗
PubMed1998

Molecular dynamics simulations of peptide fragments from hen lysozyme: insight into non-native protein conformations.

Molecular dynamics simulations of four peptides taken from the hen lysozyme sequence have been used to generate models for non-native protein conformations. Comparisons between the different peptides and with experimental data for denatured lysozyme and peptide fragments provides insight into the characteristics of the conformational ensembles populated in these non-native states and the dependence of their structural features on the amino acid sequence. For the denatured conformers populated local contacts dominate in determining the properties observed in the trajectories, all four peptides showing similar characteristics. These include a significant increase in the number of main-chain O(i)-NH(i+2) hydrogen bonds and hydrogen bonds involving side-chain groups, this increase compensating to a large extent for the loss of hydrogen bonds involved in helical or beta-sheet secondary structure in the native fold, and the generation of a population of collapsed states with local clusterings of hydrophobic groups. The hydrophobic clusters enable at least partial burial of many side-chains exposed by the loss of tertiary contacts on denaturation and provide models that may explain the experimentally observed protection of amides from hydrogen exchange and the existence of residual secondary structure in non-native species of lysozyme. The results suggest that this approach has an important role to play in aiding the interpretation of experimental data for conformationally disordered non-native states of proteins.

Animals↗
PubMed1997

Solvent structure at a hydrophobic protein surface.

The impact of an extensive, uniform and hydrophobic protein surface on the behavior of the surrounding solvent is investigated. In particular, focus is placed on the possible enhancement of the structure of water at the interface, one model for the hydrophobic effect. Solvent residence times and radial distribution functions are analyzed around three types of atomic sites (methyl, polar, and positively charged sites) in 1 ns molecular dynamics simulations of the alpha-helical polypeptide SP-C in water, in methanol and in chloroform. For comparison, water residence times at positively and negatively charged sites are obtained from a simulation of a highly charged alpha-helical polypeptide from the protein titin in water. In the simulations the structure of water is not enhanced at the hydrophobic protein surface, but instead is disrupted and devoid of positional correlation beyond the first solvation sphere. Comparing solvents of different polarity, no clear trend toward the most polar solvent being more ordered is found. In addition, comparison of the water residence times at nonpolar, polar, positively charged, or negatively charged sites on the surface of SP-C or titin does not reveal pronounced or definite differences. It is shown, however, that the local environment may considerably affect solvent residence times. The implications of this work for the interpretation of the hydrophobic effect are discussed.

Amino Acid Sequence↗
PubMed1995

Fluctuation and cross-correlation analysis of protein motions observed in nanosecond molecular dynamics simulations.

Nanosecond molecular dynamics simulations of bovine pancreatic trypsin inhibitor and lysozyme in water are analyzed in terms of backbone atomic positional fluctuations and dynamical cross-correlations. It is found that although the molecular systems are stable, B-factors calculated over a time period as long as 500 ps are not representative for the motions within the proteins. This is especially true for the most mobile residues. On a nanosecond time-scale, the B-factors calculated from the simulations of the proteins in solution are considerably larger than those obtained by structure refinement of the proteins in crystals, based on X-ray data. The time evolution of the atomic fluctuations shows that for large portions of the proteins under study, atomic positional fluctuations are not yet converged after a nanosecond. Cross-correlations do not converge faster than the fluctuations themselves. Most display very erratic behavior if the sampling covers less than about 200 ps. It is also shown that inclusion of mobile atoms into the procedure used to remove rigid-body motion from the simulation can lead to spurious correlations between the motions of the atoms at the surface of the protein.

Animals↗
PubMed1995

Comparison of MD simulations and NMR experiments for hen lysozyme. Analysis of local fluctuations, cooperative motions, and global changes.

Three 1000 ps molecular dynamics simulations of hen lysozyme have been compared with a range of experimental NMR parameters in order to gain insight into the dynamical properties of the protein and to assess the significance of the motional events observed in the simulations. The simulations, one in vacuum and two in water, were used to estimate interproton distances (for comparison with NOE data), 3JHN alpha and 3J alpha beta coupling constants and 1H-15N order parameters. Comparison of these values with experimental data, particularly NOEs, enabled force field-induced changes to the structure during the simulations to be recognized. It has been shown, however, that these changes can be largely eliminated by slight modifications to the force field. Using a simulation performed in water with this modified force field, it has been found that 1H-15N order parameters calculated for side chain groups in particular correlate well with experimental values and reflect the substantial dependence of these motional properties on the environment, particularly surface exposure, in which the side chain is found. In this case, the simulation then provides models for the motional processes giving rise to the observed experimental data. The results indicate that the order parameter values reflect primarily the number of torsion angles about which rotameric interchange occurs. In addition to local motions, the two different domains of lysozyme have been found to behave differently in the simulations. Possible implications of these differences for the interpretation of unfolding simulations and experimental observations of folding intermediates for lysozyme are discussed.

Animals↗
PubMed1995

Investigation of protein unfolding and stability by computer simulation.

Structural, dynamic and energetic properties of proteins in solution can be studied in atomic detail by molecular dynamics computer simulation. Protein unfolding can be caused by a variety of driving forces induced in different ways: increased temperature or pressure, change of solvent composition, or protein amino acid mutation. The stability and unfolding of four different proteins (bovine pancreatic trypsin inhibitor, hen egg white lysozyme, the surfactant protein C and the DNA-binding domain of the 434 repressor) have been studied by applying the afore-mentioned driving forces and also to some artificial forces. The results give a picture of protein (in)stability and possible unfolding pathways, and are compared to experimental data where possible.

Amino Acid Sequence↗
PubMed1995

The effect of environment on the stability of an integral membrane helix: molecular dynamics simulations of surfactant protein C in chloroform, methanol and water.

A series of three molecular dynamics simulations at 300 K in explicit solvent environments of chloroform, methanol and water has been performed on the pulmonary surfactant lipoprotein, SP-C, comprising several consecutive valine residues in order to investigate the stability of the alpha-helical conformation. Two additional simulations were performed on truncated SP-C with a five-residue N-terminal deletion at 300 K and 500 K in water, the high temperature run in order to increase the rate of peptide denaturation. Indications of destabilization appear in chloroform during 1 ns while the SP-C alpha-helix is remarkably stable during 1 ns in methanol and water. In particular the polyvalyl part comprising residues Val15 to Val21 remains intact even at elevated temperature, and the valines do not disrupt the alpha-helical conformation. The valyl-rotamer sampling is partly restricted. Unfolding takes place successively along the primary sequence starting from the C-terminal end. Factors affecting polypeptide stability in molecular dynamics simulations are addressed. The intrinsic helix-forming tendency of valine residues and its dependence on the sequence context, and the role of the solvent environment in stabilizing or destabilizing an alpha-helical fold, are discussed.

Amino Acid Sequence↗
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↗