PubMed HealthSearch

Biomedical subjects

R M Brunne

Publications and source records attributed to R M Brunne.

11 recordsLinked to original sources

Probing protein structure by solvent perturbation of NMR spectra: the surface accessibility of bovine pancreatic trypsin inhibitor.

In the absence of specific interactions, the relative attenuation of protein NMR signals due to added stable free radicals such as TEMPOL should reflect the solvent accessibility of the molecular surface. The quantitative correlation between observed attenuation and surface accessibility was investigated with a model system, i.e., the small protein bovine pancreatic trypsin inhibitor. A detailed discussion is presented on the reliability and limits of the approach, and guidelines are provided for data acquisition, treatment, and interpretation. The NMR-derived accessibilities are compared with those obtained from x-ray diffraction and molecular dynamics data. Although the time-averaged accessibilities from molecular dynamics are ideally suited to fit the NMR data, better agreement was observed between the paramagnetic attenuations of the fingerprint cross-peaks of homonuclear proton spectra and the total NH and H alpha accessibilities calculated from x-ray coordinates, than from time-averaged molecular dynamics simulations. In addition, the solvent perturbation response appears to be a promising approach for detecting the thermal conformational evolution of secondary structure elements in proteins.

Amino Acid Sequence

A pseudo-particle approach for studying protein-ligand models truncated to their active sites.

A molecular dynamics method has been developed to describe the structural and dynamic properties of protein-ligand complexes that are truncated to their active sites. The active site is comprised of the ligand and discontinuous, positionally unrestrained peptide chains. This truncated active-site complex is surrounded by big unspecific pseudo-particles representing the complete protein and the solvent. Thus, knowledge of the folding of the outer parts of the protein is not required, and the method can be applied to protein models, derived from homology modeling. The method has been tested using ligand complexes of adenylate kinase, retinol binding protein, HIV-1 protease, and human leucocyte antigen. Comparisons with their crystal structures and with results from time-demanding simulations of the whole complexes in explicit water solvent show that the ligand binding properties are conserved. Most of the hydrogen bonds between the ligand and the active-site residues are reproduced and, furthermore, the simulation time is reduced.

Amino Acid Sequence

Structure and internal dynamics of the bovine pancreatic trypsin inhibitor in aqueous solution from long-time molecular dynamics simulations.

Structural and dynamic properties of bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution are investigated using two molecular dynamics (MD) simulations: one of 1.4 ns length and one of 0.8 ns length in which atom-atom distance bounds derived from NMR spectroscopy are included in the potential energy function to make the trajectory satisfy these experimental data more closely. The simulated properties of BPTI are compared with crystal and solution structures of BPTI, and found to be in agreement with the available experimental data. The best agreement with experiment was obtained when atom-atom distance restraints were applied in a time-averaged manner in the simulation. The polypeptide segments found to be most flexible in the MD simulations coincide closely with those showing differences between the crystal and solution structures of BPTI.

Animals

Nucleotide-binding properties of adenylate kinase from Escherichia coli: a molecular dynamics study in aqueous and vacuum environments.

The complex of adenylate kinase with its transition-state inhibitor has been studied by molecular dynamics simulations in water and in vacuum environments with the GROMOS force field over a period of 300 ps. The adenylate kinase, a member of the nucleotide-binding protein family, was exemplarily chosen for the inspection of the nucleotide-binding properties in the active site. The ligand binding and the domain movements have been studied in detail over the simulation period and compared with the crystal structure. Secondary structure transitions and domain closures defined those parts of the structure which are involved in an induced-fit movement of the enzyme. The presence of more stable hydrogen bonds on the substrate side leads to the assumption that substrate binding is more specific than cosubstrate binding. Reliable results were achieved only if water was explicitly included in the stimulation.

Adenylate Kinase

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

Hydration of proteins. A comparison of experimental residence times of water molecules solvating the bovine pancreatic trypsin inhibitor with theoretical model calculations.

A 1 ns trajectory from a molecular dynamics study of 1.4 ns total length was used for a detailed analysis of the residence times of water molecules located near 227 selected bovine pancreatic trypsin inhibitor (BPTI) atoms. The simulation was performed using the GROMOS force field, with apolar hydrogen atoms treated as united atoms, and the SPC/E water model. The system consisted of 568 BPTI atoms and 2371 water molecules. The theoretical results are in good agreement with experimental data available from nuclear magnetic resonance spectroscopy. The residence times of individual water molecules coming near a given BPTI atom, as obtained from the simulation, vary greatly and range between 10 and 500 ps. The effective residence time, calculated using a correlation function technique from the presence of all individual water molecules visiting the hydration shell of a given BPTI atom, never exceeds 200 ps. The average residence time near backbone and side-chain atoms is approximately 39 ps and 24 ps, respectively. The shortest residence times, on average, are found near charged atoms (19 ps), whereas near non-polar and polar side-chain atoms the residence times are 25 ps and 36 ps, respectively. There is no apparent correlation between the residence times of the hydration water molecules of solvent-accessible residues and their location in different regular or non-regular secondary structures.

Amino Acid Sequence

Dynamical properties of bovine pancreatic trypsin inhibitor from a molecular dynamics simulation at 5000 atm.

Molecular dynamics simulations of bovine pancreatic trypsin inhibitor in water have been performed with coupling to pressure baths at 1 atm and at 5000 atm. The positional fluctuations of atoms in the alpha- and gamma-positions are slightly decreased at 5000 atm. The mobility of the backbone phi- and psi-angles is not affected with respect to the root mean square fluctuations and the rate of torsional angle transitions. The amplitude of libration of sidechain chi-angles remains nearly the same for both pressures, but the rate of torsional angle transitions decreases on average by 30% when increasing the hydrostatic pressure to 5000 atm.

Amino Acid Sequence

Structure refinement using time-averaged J-coupling constant restraints.

We describe a new penalty function for use in restrained molecular dynamics simulations which allows experimental J-coupling information to be enforced as a time-averaged, rather than instantaneous, quantity. The pseudo-energy term has been formulated in terms of a calculated J value (a measured quantity) rather than the relevant dihedral angle (a derived quantity). This accounts for the distinct non-linearity of the coupling constant with respect to either Cartesian coordinates or dihedral angles. Example simulations of the cyclic decapeptide antamanide show the procedure's ability to enforce experimental restraints while exploring a large region of conformational space and producing a relatively small disturbance of the physical force field.

Amino Acid Sequence

Comparative conformational analysis and in vitro pharmacological evaluation of three cyclic hexapeptide NK-2 antagonists.

The conformational analysis of three cyclic hexapeptides is presented. Cyclo-(-Gln6-Trp7-Phe8-Gly9-Leu10-D-Met11-) (1) and cyclo-(-Gln6-Trp7-Phe8-Gly9-Leu10-Met11-) (2) are NK-2 antagonists in the hamster trachea assay, whereas cyclo-(-Gln6-Trp7-Phe8-(R)-Gly9-[ANC-2]Leu10-Met11+ ++-) (3), where Gly9[ANC-2]Leu10 represents (2S)-2-((3R)-3-amino-2-oxo-1-pyrrolidinyl)-4-methylpentanoyl, is inactive as agonist and antagonist in this assay. In DMSO, the NMR results cannot be interpreted as being consistent with a single conformation. However, the combined interpretation of results from NMR spectroscopy, restrained molecular dynamics simulations with application of proton-proton distance information from ROESY spectra, and pharmacological results leads to a reduced number of conformational domains for each peptide, which can be compared with each other and may be classified as responsible for their biological activity. Trying to match the conformational domains approximately with regular beta- and gamma-turns, we find a gamma n-turn at the position of the methionine occurring in all peptides. For the active peptides 1 and 2 we arrive at an inverse gamma i-turn at Phe8, and beta I'- or beta II-turns with Gly9 and Leu10 at the corner positions, these beta-turns having a similar topology with respect to the linking peptide unit. Other conformational domains common to only 1 and 2 support their classification as responsible for the biological activity.

Amino Acid Sequence

Conformational regions of Boc-Ala-Aib-Ala-OMe. Sampling with molecular dynamics simulations using time averaging of distance restraints.

The method of time averaging of distance restraints in molecular dynamics simulations is applied to Boc-Ala-Aib-Ala-OMe in order to demonstrate the improved sampling properties of this method compared to conventional distance restraining. Two conformational regions, beta-turn type II and gamma-turn, are seen during MD runs at a simulation temperature of 500 K, while in simulations with conventional distance restraining, no conformational transitions could be observed for temperatures up to 1000 K.

Amino Acid Sequence