Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “MD simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64Linked to original sources
PubMed2005

Deamidation of model beta-turn cyclic peptides in the solid state.

To investigate the importance of secondary structure on peptide deamidation in the solid state, two cyclic beta-turn peptides and their linear analogs were used as models of Asn residues in structured and unstructured domains, and incorporated into poly(vinyl pyrrolidone) (PVP)-based lyophilized solids. The secondary structure of the model peptides was determined in solution and the solid state using a combination of nuclear magnetic resonance (NMR) spectroscopy, circular dichroism (CD), and Fourier transform infrared (FTIR) spectroscopy. The model beta-turn cyclic peptides were found to be type II beta-turns while the linear analogs were determined to be predominantly unstructured. Quantitatively, the cyclic peptides consisted of approximately 80% beta-turn while the linear analogs contained only 30%-35% beta-turn. To characterize the solid environment, T(g), and moisture content of the solid-state formulations were determined. Accelerated stability studies were conducted in the solid state at 37 degrees C using formulations lyophilized from solutions at pH 8.8 (0.1 M borate buffer). The effect of matrix mobility on solid-state deamidation was investigated by altering the moisture content through variation of relative humidity or the addition of a plasticizer. Cyclic peptides degraded 1.2-8 times slower than the linear analogs under all of the conditions studied. The observed rate constants, however, for all of the peptides decreased dramatically (four orders of magnitude) in the glassy solids. This suggests the greater importance of matrix mobility in solid-state degradation. Molecular dynamics (MD) simulations were also performed to explore the low energy, preferred state of the peptides, and determine the structure around the beta-turn.

Amides↗
PubMed2005

Configuration and E/Z interconversion mechanism of O(S)-allyl-S(O)-methyl-N-(acridin-9-yl) iminothiocarbonate.

The configuration and dynamic behavior of O-allyl-S-methyl-N-(acridin-9-yl)iminothiocarbonate (1) and its S-allyl-O-methyl regioisomer (2) were studied using quantum chemical calculations and by applying a novel graphical method to scatter maps obtained from MD simulations for evaluation of an NOE-weighted internuclear distance (r(NOE)). Energy calculations indicated that the Z configuration was predominant for each compound and, further, this was supported both by the calculated chemical shifts and the r(NOE). Both N-inversion- and rotation-type transition-state structures were also calculated for the E/Z isomerization process, the results indicating that the preferred interconversion mechanism for 1 is N-inversion, but contrastingly, interconversion via rotation is equally as probable as N-inversion for 2. This supports the notion that one or the other or both pathways can be active and each system needs to be assessed on a case-by-case basis.

Journal Article↗
PubMed1992

Quantitative analysis of cyclic beta-turn models.

The beta-turn is a frequently found structural unit in the conformation of globular proteins. Although the circular dichroism (CD) spectra of the alpha-helix and beta-pleated sheet are well defined, there remains some ambiguity concerning the pure component CD spectra of the different types of beta-turns. Recently, it has been reported (Hollósi, M., Kövér, K.E., Holly, S., Radics, L., & Fasman, G.D., 1987, Biopolymers 26, 1527-1572; Perczel, A., Hollósi, M., Foxman, B.M., & Fasman, G.D., 1991a, J. Am. Chem. Soc. 113, 9772-9784) that some pseudohexapeptides (e.g., the cyclo[(delta)Ava-Gly-Pro-Aaa-Gly] where Aaa = Ser, Ser(OtBu), or Gly) in many solvents adopt a conformational mixture of type I and the type II beta-turns, although the X-ray-determined conformation was an ideal type I beta-turn. In addition to these pseudohexapeptides, conformational analysis was also carried out on three pseudotetrapeptides and three pseudooctapeptides. The target of the conformation analysis reported herein was to determine whether the ring stress of the above beta-turn models has an influence on their conformational properties. Quantitative nuclear Overhauser effect (NOE) measurements yielded interproton distances. The conformational average distances so obtained were interpreted utilizing molecular dynamics (MD) simulations to yield the conformational percentages. These conformational ratios were correlated with the conformational weights obtained by quantitative CD analysis of the same compounds. The pure component CD curves of type I and type II beta-turns were also obtained, using a recently developed algorithm (Perczel, A., Tusnády, G., Hollósi, M., & Fasman, G.D., 1991b, Protein Eng. 4(6), 669-679). For the first time the results of a CD deconvolution, based on the CD spectra of 14 beta-turn models, were assigned by quantitative NOE results. The NOE experiments confirmed the ratios of the component curves found for the two major beta-turns by CD analysis. These results can now be used to enhance the conformational determination of globular proteins on the basis of their CD spectra.

Amino Acid Sequence↗
PubMed1996

Dynamics of ribonuclease A and ribonuclease S: computational and experimental studies.

RNase S is a complex consisting of two proteolytic fragments of RNase A: the S peptide (residues 1-20) and S protein (residues 21-124). RNase S and RNase A have very similar X-ray structures and enzymatic activities. Previous experiments have shown increased rates of hydrogen exchange and greater sensitivity to tryptic cleavage for RNase S relative to RNase A. It has therefore been asserted that the RNase S complex is considerably more dynamically flexible than RNase A. In the present study we examine the differences in the dynamics of RNase S and RNase A computationally, by MD simulations, and experimentally, using trypsin cleavage as a probe of dynamics. The fluctuations around the average solution structure during the simulation were analyzed by measuring the RMS deviation in coordinates. No significant differences between RNase S and RNase A dynamics were observed in the simulations. We were able to account for the apparent discrepancy between simulation and experiment by a simple model. According to this model, the experimentally observed differences in dynamics can be quantitatively explained by the small amounts of free S peptide and S protein that are present in equilibrium with the RNase S complex. Thus, folded RNase A and the RNase S complex have identical dynamic behavior, despite the presence of a break in polypeptide chain between residues 20 and 21 in the latter molecule. This is in contrast to what has been widely believed for over 30 years about this important fragment complementation system.

Computer Simulation↗
PubMed2003

Functional role for Tyr 31 in the catalytic cycle of chicken dihydrofolate reductase.

Despite much work, many key aspects of the mechanism of the dihydrofolate reductase (DHFR) catalyzed reduction of dihydrofolate remain unresolved. In bacterial forms of DHFR both substrate and water access to the active site are controlled by the conformation of the mobile M20 loop. In vertebrate DHFRs only one conformation of the residues corresponding to the M20 loop has been observed. Access to the active site was proposed to be controlled by residue 31. MD simulations of chicken DHFR complexed with substrates and cofactor revealed a closing of the side chain of Tyr 31 over the active site on binding of dihydrofolate. This conformational change was dependent on the presence of glutamate on the para-aminobenzoylamide moiety of dihydrofolate. In its absence, the conformation remained open. Although water could enter the active site and hydrogen bond to N5 of dihydrofolate, indicating the feasibility of water as the proton donor, this was not controlled by the conformation of Tyr 31. The water accessibility of the active site was low for both conformations of Tyr 31. However, when hydride was transferred from NADPH to C6 of dihydrofolate before protonation, the average time during which water was found in hydrogen bonding distance to N5 of dihydrofolate in the active site increased almost fivefold. These results indicated that water can serve as the Broensted acid for the protonation of N5 of dihydrofolate during the DHFR catalyzed reduction.

Animals↗
PubMed2003

Effective energy function for proteins in lipid membranes.

A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.

Amino Acids↗
PubMed2003

Characterization of the active site of DNA polymerase beta by molecular dynamics and quantum chemical calculation.

It is well established that the fully formed polymerase active site of the DNA repair enzyme, polymerase beta (pol beta), including two bound Mg2+ cations and the nucleoside triphosphate (dNTP) substrate, exists at only one point in the catalytic cycle just prior to the chemical nucleotidyl transfer step. The structure of the active conformation has been the subject of much interest as it relates to the mechanism of the chemical step and also to the question of fidelity assurance. Although crystal structures of ternary pol beta-(primer-template) DNA-dNTP complexes have provided the main structural features of the active site, they are necessarily incomplete due to intentional alterations (e.g., removal of the 3'OH groups from primer and substrate) needed to obtain a structure from midcycle. Working from the crystal structure closest to the fully formed active site [Protein Data Bank (PDB) code: 1bpy], two molecular dynamics (MD) simulations of the solvated ternary complex were performed: one with the missing 3'OHs restored, via modeling, to the primer and substrate, and the other without restoration of the 3'OHs. The results of the simulations, together with ab initio optimizations on simplified active-site models, indicate that the missing primer 3'OH in the crystal structure is responsible for a significant perturbation in the coordination sphere of the catalytic cation and allow us to suggest several corrections and additions to the active-site structure as observed by crystallography. In addition, the calculations help to resolve questions raised regarding the protonation states of coordinating ligands.

Binding Sites↗
PubMed2004

Immunoglobulin kappa light chain and its amyloidogenic mutants: a molecular dynamics study.

AL amyloidosis and LCDD are pathological conditions caused by extracellural deposition of monoclonal Ig light chain variable domains. In the former case, deposits have a form of amyloid fibrils, in the latter, amorphous aggregates. 1REI kappa light chain variable domain and its two point mutants, R61N and D82I, were chosen for the analysis in this work. Wild 1REI does not create deposits in vitro, while R61N aggregates as amyloid fibrils and D82I creates amorphous aggregates. Both mutated residues create a conserved salt bridge; thus, substitution of any of them should decrease V(L) domain stability. For these three proteins, 5 ns MD simulations were conducted in temperatures of 300 K and 400 K, with protonated and unprotonated acidic residues, mimicking acidic and neutral experimental pH conditions (3 sets: N300, N400, and A400). The analysis of trajectories focused on characterization of changes in conformational behavior and stability of Ig kappa light chain variable domain caused by single aminoacid substitutions that were experimentally proved to enhance aggregation propensity, both in the form of amyloid and amorphous aggregates. Residue D82 turns out to be involved not only in R61-D82 but also in K45-D82 interaction, which was not observed in the X-ray structure, but frequently populated simulations of 1REI. The substitution D82I excludes both interactions, resulting in substantial destabilization (i.e., easier aggregation). Examination of behavior of edge regions of V(L) beta-sandwich reveals significant alterations in D82I mutant compared to wild 1REI, while relatively small changes occur in R61N. This suggests that mild and slow destabilization is the reason of the conversion of V(L) to partially folded amyloidogenic intermediate structure.

Amino Acids, Acidic↗
PubMed2001

Helix-coil transition of PrP106-126: molecular dynamic study.

A set of 34 molecular dynamic (MD) simulations totaling 305 ns of simulation time of the prion protein-derived peptide PrP106-126 was performed with both explicit and implicit solvent models. The objective of these simulations is to investigate the relative stability of the alpha-helical conformation of the peptide and the mechanism for conversion from the helix to a random-coil structure. At neutral pH, the wild-type peptide was found to lose its initial helical structure very fast, within a few nanoseconds (ns) from the beginning of the simulations. The helix breaks up in the middle and then unwinds to the termini. The spontaneous transition into the random coil structure is governed by the hydrophobic interaction between His(111) and Val(122). The A117V mutation, which is linked to GSS disease, was found to destabilize the helix conformation of the peptide significantly, leading to a complete loss of helicity approximately 1 ns faster than in the wild-type. Furthermore, the A117V mutant exhibits a different mechanism for helix-coil conversion, wherein the helix begins to break up at the C-terminus and then gradually to unwind towards the N-terminus. In most simulations, the mutation was found to speed up the conversion through an additional hydrophobic interaction between Met(112) and the mutated residue Val(117), an interaction that did not exist in the wild-type peptide. Finally, the beta-sheet conformation of the wild-type peptide was found to be less stable at acidic pH due to a destabilization of the His(111)-Val(122), since at acidic pH this histidine is protonated and is unlikely to participate in hydrophobic interaction.

Amino Acid Substitution↗
PubMed2004

Analysis of CDK2 active-site hydration: a method to design new inhibitors.

The interactions between the protein and the solvent were analyzed, and protein regions with a high density of water molecules, as well as structural water molecules, were determined by using molecular dynamics (MD) simulations. A number of water molecules that were in contact with the protein for the whole trajectory were determined. Their interaction energies and hydrogen bonds with protein residues were analyzed. Altogether, 39, 27, 49, and 32 water molecules bound to the protein were found for trajectories of the free CDK2, CDK2/ATP, CDK2/roscovitine, and CDK2/isopentenyladenine complexes, respectively. Positions of observed water molecules were compared with X-ray crystallography data. Special attention was paid to water molecules in the active site of the enzyme, and especially to the deep pocket, where the N9 roscovitine side-chain is buried. Exchange of active-site water molecules with bulk water through the tunnel from the pocket was observed. In the CDK2/isopentenyladenine complex simulation, two water molecules that arrange interaction between the inhibitor and the enzyme via an H-bond were observed. Two stable water molecules in the trajectory of the free CDK2 were found that occupy the same position as the nitrogens N3 and N9 of the isopentenyladenine or N1 and N6 nitrogens of the adenosine triphosphate (ATP). The positions of structural water molecules were compared with the positions of substrate polar groups and crystallographic water molecules found in the Brookhaven Protein Data Bank for various CDK2 complexes. It was concluded that tracing tightly bound water molecules may substantially help in designing new inhibitors.

Adenine↗
PubMed2004

Exploring protein native states and large-scale conformational changes with a modified generalized born model.

Implicit solvation models provide, for many applications, a reasonably accurate and computationally effective way to describe the electrostatics of aqueous solvation. Here, a popular analytical Generalized Born (GB) solvation model is modified to improve its accuracy in calculating the solvent polarization part of free energy changes in large-scale conformational transitions, such as protein folding. In contrast to an earlier GB model (implemented in the AMBER-6 program), the improved version does not overstabilize the native structures relative to the finite-difference Poisson-Boltzmann continuum treatment. In addition to improving the energy balance between folded and unfolded conformers, the algorithm (available in the AMBER-7 and NAB molecular modeling packages) is shown to perform well in more than 50 ns of native-state molecular dynamics (MD) simulations of thioredoxin, protein-A, and ubiquitin, as well as in a simulation of Barnase/Barstar complex formation. For thioredoxin, various combinations of input parameters have been explored, such as the underlying gas-phase force fields and the atomic radii. The best performance is achieved with a previously proposed modification to the torsional potential in the Amber ff99 force field, which yields stable native trajectories for all of the tested proteins, with backbone root-mean-square deviations from the native structures being approximately 1.5 A after 6 ns of simulation time. The structure of Barnase/Barstar complex is regenerated, starting from an unbound state, to within 1.9 A relative to the crystal structure of the complex.

Bacterial Proteins↗
PubMed2004

Accurate and efficient description of protein vibrational dynamics: comparing molecular dynamics and Gaussian models.

Current all-atom potential based molecular dynamics (MD) allows the identification of a protein's functional motions on a wide-range of timescales, up to few tens of nanoseconds. However, functional, large-scale motions of proteins may occur on a timescale currently not accessible by all-atom potential based MD. To avoid the massive computational effort required by this approach, several simplified schemes have been introduced. One of the most satisfactory is the Gaussian network approach based on the energy expansion in terms of the deviation of the protein backbone from its native configuration. Here, we consider an extension of this model that captures in a more realistic way the distribution of native interactions due to the introduction of effective side-chain centroids. Since their location is entirely determined by the protein backbone, the model is amenable to the same exact and computationally efficient treatment as previous simpler models. The ability of the model to describe the correlated motion of protein residues in thermodynamic equilibrium is established through a series of successful comparisons with an extensive (14 ns) MD simulation based on the AMBER potential of HIV-1 protease in complex with a peptide substrate. Thus, the model presented here emerges as a powerful tool to provide preliminary, fast yet accurate characterizations of protein near-native motion.

Computer Simulation↗
PubMed2004

Is histoaspartic protease a serine protease with a pepsin-like fold?

The primary structure of the so-called histoaspartic protease from Plasmodium falciparum has a very high percentage of identity and homology with the pepsin-like enzyme plasmepsin II. A homology modeling approach was used to calculate the three-dimensional structure of the enzyme. Molecular dynamics (MD) simulations were applied to find those structural properties of the histoaspartic protease that had a tendency to remain stable during all runs. The results have shown that hydrogen-bonded residues Ser37-His34-Asp214 are arranged without any strain, in a manner that resembles the active site of a serine protease, while Ser38 and Asn39 take up positions appropriate to formation of an oxyanion hole. Although there are several important differences between the enzyme and plasmepsin II, all of the structural features associated with a typical pepsin-like aspartic protease are present in the final model of the histoaspartic protease. A possibility that this enzyme may function as a serine protease is discussed.

Animals↗
PubMed2005

Theoretical investigations of prostatic acid phosphatase.

The phosphotyrosyl protein phosphatase activity of prostatic acid phosphatase (PAP) has been well established. It has also been suggested that PAP partly regulates the activity of growth factor receptors by dephosphorylating the autophosphorylysable tyrosines in them. We studied the binding of the peptides from epidermal growth factor receptor (EGFR) and its homolog (ErbB-2), corresponding to their autophosphorylation sites, to PAP using theoretical modeling and molecular dynamics (MD) simulation methods. Nine different peptides, each with a phosphotyrosine residue, were docked on human PAP. The binding energies of these peptide-PAP complexes were calculated theoretically and compared to experimentally obtained affinities. The peptide Ace--DNLpYYWD--NH2 from ErbB-2(1197-1203) showed the most favorable free energy of binding when estimated theoretically. The results demonstrate that the presence of another tyrosine residue proximate to C-terminal of autophosphorylysable Tyr enhances the binding affinity considerably. The presence of a bulky group instead prevents the binding, as is observed in case of peptide Ace--NLYpYWDQ--NH2 which failed to bind, both in theoretical calculations and experiments. Thus we demonstarted that PAP could potentially bind to EGFR and Erbb-2 and dephosphorylate them. Thus it could be involved in the regulation of the function of such receptors. In addition, complexes of a peptide from AngiotensinII and phosphotyrosine(pY) with human PAP were also modeled. The effects of different protonation states of the titratable active site residues on ligand (pY) binding have also been investigated. For a favorable binding His12 and Asp258 should be neutral, His257 should be positively charged and the phosphate group of the ligand should be in PO(4) (3-) state. Furthermore, the analysis of protein motion as observed during simulations suggests the loop-loop contact in the PAP dimer to be of importance in cooperativity.

Acid Phosphatase↗
PubMed2005

Analysis of correlated domain motions in IgG light chain reveals possible mechanisms of immunological signal transduction.

It was shown experimentally that binding of a micelle composed of Congo red molecules to immunological complexes leads to the enhanced stability of the latter, and simultaneously prevents binding of a complement molecule (C1q). The dye binds in a cavity created by the removal of N-terminal polypeptide chain, as observed experimentally in a model system-immunoglobulin G (IgG) light chain dimer. Molecular Dynamics (MD) simulations of three forms of IgG light chain dimer, with and without the dye, were performed to investigate the role of N-terminal fragment and self-assembled ligand in coupling between V and C domains. Root-mean-square distance (RMSD) time profiles show that removal of N-terminal fragment leads to destabilization of V domain. A micelle composed of four self-assembled dye molecules stabilizes and fixes the domain. Analysis of root-mean-square fluctuation (RMSF) values and dynamic cross-correlation matrices (DCCM) reveals that removal of N-terminal fragment results in complete decoupling between V and C domains. Binding of self-assembled Congo red molecules improves the coupling, albeit slightly. The disruption of a small beta-sheet composed of N- and C-terminal fragments of the domain (NC sheet) is the most likely reason for the decoupling. Self-assembled ligand, bound in the place originally occupied by N-terminal fragment, is not able to take over the function of the beta-sheet. Lack of correlation of motions between residues in V and C domains denotes that light chain-Congo red complexes have hampered ability to transmit conformational changes between domains. This is a likely explanation of the lack of complement binding by immunological complexes, which bind Congo red, and supports the idea that the NC sheet is the key structural fragment taking part in immunological signal transduction.

Computer Simulation↗
PubMed2005

The betaalphabetaalphabeta elementary supersecondary structure of the Rossmann fold from porcine lactate dehydrogenase exhibits characteristics of a molten globule.

Protein classifications show that the Rossmann fold, which consists of two betaalphabetaalphabeta motifs (BABAB) related by a rough twofold axis, is the most populated alphabeta fold, and that the betaalphabeta submotif (BAB) is a widespread elementary structural arrangement. Herein, we report MD simulations, circular dichroism and NMR analyses on BAB and BABAB from porcine lactate dehydrogenase to evaluate their intrinsic stability. Our results demonstrate that BAB is not stable in solution and is not a folding nucleus. We also find that BABAB, despite its appearance of a functional and structural unit, is not an independent and thermodynamically stable folding unit. Rather, we show that BABAB retains most native secondary structure but very little tertiary structure, thus displaying characteristics of a molten globule.

Animals↗
PubMed2005

Flexibility of prolyl oligopeptidase: molecular dynamics and molecular framework analysis of the potential substrate pathways.

The flexibility of prolyl oligopeptidase has been investigated using molecular dynamics (MD) and molecular framework approaches to delineate the route of the substrate to the active site. The selectivity of the enzyme is mediated by a seven-bladed beta-propeller that in the crystal structure does not indicate the possible passage for the substrate to the catalytic center. Its open topology however, could allow the blades to move apart and let the substrate into the large central cavity. Flexibility analysis of prolyl oligopeptidase structure using the FIRST (Floppy Inclusion and Rigid Substructure Topology) approach and the atomic fluctuations derived from MD simulations demonstrated the rigidity of the propeller domain, which does not permit the substrate to approach the active site through this domain. Instead, a smaller tunnel at the inter-domain region comprising the highly flexible N-terminal segment of the peptidase domain and a facing hydrophilic loop from the propeller (residues 192-205) was identified by cross-correlation analysis and essential dynamics as the only potential pathway for the substrate. The functional importance of the flexible loop has been also verified by kinetic analysis of the enzyme with a split loop. Catalytic effect of engineered disulfide bridges was rationalized by characterizing the concerted motions of the two domains.

Animals↗
PubMed2006

Peptide fragment studies on the folding elements of dihydrofolate reductase from Escherichia coli.

One of the necessary conditions for a protein to be foldable is the presence of a complete set of "folding elements" (FEs) that are short, contiguous peptide segments distributed over an amino acid sequence. The FE-assembly model of protein folding has been proposed, in which the FEs play a role in guiding structure formation through FE-FE interactions early in folding. However, two major issues remain to be clarified regarding the roles of the FEs in determining protein foldability. Are the FEs AFUs that can form nativelike structures in isolation? Is the presence of only the FEs without mutual connections a sufficient condition for a protein to be foldable? Here, we address these questions using peptide fragments corresponding to the FEs of DHFR from Escherichia coli. We show by CD measurement that the FE peptides are unfolded under the native conditions, and some of them have the propensities toward non-native helices. MD simulations also show the non-native helical propensities of the peptides, and the helix contents estimated from the simulations are well correlated with those estimated from the CD in TFE. Thus, the FEs of DHFR are not AFUs, suggesting the importance of the FEs in nonlocal interactions. We also show that equimolar mixtures of the FE peptides do not induce any structural formation. Therefore, mutual connections between the FEs, which should strengthen the nonlocal FE-FE interactions, are also one of the necessary conditions for a protein to be foldable.

Circular Dichroism↗