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 523 records · Page 29Linked to original sources
PubMed2006

Cofactor assisted gating mechanism in the active site of NADH oxidase from Thermus thermophilus.

NADH oxidase (NOX) from Thermus thermophilus is a member of a structurally homologous flavoprotein family of nitroreductases and flavin reductases. The importance of local conformational dynamics in the active site of NOX has been recently demonstrated. The enzyme activity was increased by 250% in the presence of 1 M urea with no apparent perturbation of the native structure of the protein. The present in silico results correlate with the in vitro data and suggest the possible explanation about the effect of urea on NOX activity at the molecular level. Both, X-ray structure and molecular dynamics (MD) simulations, show open conformation of the active site represented by approximately 0.9 nm distance between the indole ring of Trp47 and the isoalloxazine ring of FMN412. In this conformation, the substrate molecule can bind in the active site without sterical restraints. MD simulations also indicate more stable conformation of the active site called "closed" conformation. In this conformation, Trp47 and the isoalloxazine ring of FMN412 are so close to each other (approximately 0.5 nm) that the substrate molecule is unable to bind between them without perturbing this conformation. The open/close transition of the active site between Trp47 and the flavin ring is accompanied by release of the "tightly" bound water molecule from the active site--cofactor assisted gating mechanism. The presence of urea in aqueous solutions of NOX prohibits closing of the active site and even unlocks the closed active site because of the concomitant binding of a urea molecule in the active site cavity. The binding of urea in the active site is stabilized by formation of one/two persistent hydrogen bonds involving the carbonyl group of the urea molecule. Our report represents the first MD study of an enzyme from the novel flavoprotein family of nitroreductases and flavin reductases. The common occurrence of aromatic residues covering the active sites in homologous enzymes suggests the possibility of a general gating mechanism and the importance of local dynamics within this flavoprotein family.

Binding Sites↗
PubMed2006

Mixed atomistic and coarse-grained molecular dynamics: simulation of a membrane-bound ion channel.

The recently developed multiscale coarse-graining (MS-CG) method (Izvekov, S.; Voth, G. A. J. Phys. Chem. B 2005, 109, 2469; J. Chem. Phys. 2005, 123, 134105) is used to build a mixed all-atom and coarse-grained (AA-CG) model of the gramicidin A (gA) ion channel embedded in a dimyristoylphosphatidylcholine (DMPC) lipid bilayer and water environment. In this model, the gA peptide was described in full atomistic detail, while the lipid and water molecules were described using coarse-grained representations. The atom-CG and CG-CG interactions in the mixed AA-CG model were determined using the MS-CG method. Molecular dynamics (MD) simulations were performed using the resulting AA-CG model. The results from simulations of the AA-CG model compare very favorably to those from all-atom MD simulations of the entire system. Since the MS-CG method employs a general and systematic approach to obtain effective interactions from the underlying all-atom models, the present approach to rigorously develop mixed AA-CG models has the potential to be extended to many other systems.

Computer Simulation↗
PubMed2004

Intrinsic versus mutation dependent instability/flexibility: a comparative analysis of the structure and dynamics of wild-type transthyretin and its pathogenic variants.

Transthyretin (TTR) is one of the about 20 known human proteins associated with amyloidosis which is characterized by the accumulation of amyloid fibrils in tissues or extracellular matrix surrounding vital organs. Unlike Alzheimer's fibrils that comprise a fragment of a large precursor protein, TTR amyloid fibrils are composed of both full-length protein and fragments of the molecule. The native state of TTR is a homotetramer with eight beta-strands organized into a beta-sandwich in each monomer. To elucidate the structural reorganization mechanisms preceding amyloid formation, it is important to characterize the dynamic features of the wild-type native state as well as to reveal the influence of disease-associated mutations on the structure and dynamics. Molecular dynamics (MD) simulations complement X-ray crystallography and D-H exchange to capture the intrinsically unstable/flexible sites of the wild-type as well as the mutation dependent unstable sites of the pathogenic variants. Our results of MD simulations have shown that the Leu55-->Pro (L55P) mutation occurs in an intrinsically unstable site, leading to substantial local and global structural changes. This observation supports the early speculation that the C-strand-loop-D-strand rearrangement leads to the formation of amyloidogenic intermediates. In addition to the D strand, the alpha-helical region and the strands at the monomer-monomer interface are also intrinsically unstable. The central channel of L55P-TTR undergoes opening and closing fluctuations, which may provide an explanation for the fact that while the mutation is far from the channel, the mutant shows a substantial low binding affinity of thyroxine.

Binding Sites↗
PubMed2006

Solvent quality changes the structure of G8 PAMAM dendrimer, a disagreement with some experimental interpretations.

We have performed approximately 20-40 ns of molecular dynamics (MD) simulations for the generation 8 PAMAM dendrimer in explicit water under varying pH conditions to study the structure of the dendrimer (approximately 156,738 atoms at low pH). This is the first report of such a long MD simulation of a larger generation PAMAM dendrimer including the effect of salt and counterions with explicit water molecules. We find that changing the pH from a high value (approximately 12) to a low value (approximately 3) changes the radius of gyration from Rg = 37.8 to 43.1 A (increasing by 13%). We also find significant back-folding of the primary amines and a large amount of water penetration inside the polymer. The increase in size with decrease in pH is consistent with our earlier studies on G3-G6 and agrees with the Monte Carlo theory by Welch and Muthukumar of G8 (Macromolecules, 1998, 31, 5892) and the experiments on G5 and G8 PAMAM dendrimer by Topp et al. (Macromolecules, 1999, 32, 7232). However, these results disagree dramatically with the interpretations of SANS experiments of G8 PAMAM dendrimers by Nisato et al. (Macromolecules, 2000, 33, 4172) who observe no change in the size of the dendrimer with variations of solution pH and ionic strength. We assume that the disagreement might arise from neglecting nonspherical shape, penetration of water and ions into the core, and aggregation, all of which might depend on pH.

Journal Article↗
PubMed2005

Anionic clays containing anti-inflammatory drug molecules: comparison of molecular dynamics simulation and measurements.

Three representative nonsteroidal anti-inflammatory drug molecules, Ibuprofen, Diclofenac, and Indomethacin, have been intercalated within the galleries of an anionic clay, Mg-Al layered double hydroxide (LDH). X-ray diffraction, IR and Raman vibrational spectroscopy and (13)C cross-polarization magic-angle spinning NMR have been used to characterize the confined drug molecules, while molecular dynamics (MD) simulations were used to probe the interlayer structure, arrangement, orientation, and geometry of the intercalated species. All three drug molecules are arranged as bilayers in the interlamellar space of the anionic clay. But while the structure of the intercalated Ibuprofen is identical to that of the molecule outside the layers, spectroscopy as well as MD simulation shows that there is a change in the geometry of Diclofenac and Indomethacin upon confinement within the galleries of the LDH. The change in geometry of Diclofenac and Indomethacin upon intercalation is shown to originate from the electrostatic interaction between the electronegative chlorine atoms on the drug molecule and the positively charged metal hydroxide sheets of the anionic clay. It is shown that these changes in the geometry of the intercalated drug molecules allow for the observed interlayer spacing to be realized without the bilayers having to interdigitate, which would otherwise have been necessary if the structure of the drug molecules had remained identical to that outside the layers. Comparisons of experimental measurements with simulation have provided a more detailed understanding of the geometry and organization of flexible drug molecules confined in the anionic clay.

Journal Article↗
PubMed1999

Comparative molecular dynamics of mesophilic and psychrophilic protein homologues studied by 1.2 ns simulations.

It is well established that the dynamic motion of proteins plays an important functional role, and that the adaptation of a protein molecule to its environment requires optimization of internal non-covalent interactions and protein-solvent interactions. Serine proteinases in general, and trypsin in particular has been used as a model system in exploring possible structural features for cold adaptation. In this study, a 500 p.s. and a 1200 p.s. molecular dynamics (MD) simulation at 300 K of both anionic salmon trypsin and cationic bovine trypsin are analyzed in terms of molecular flexibility, internal non-covalent interactions and protein-solvent interactions. The present MD simulations do not indicate any increased flexibility of the cold adapted enzyme on an overall basis. However, the apparent higher flexibility and deformability of the active site of anionic salmon trypsin may lower the activation energy for ligand binding and for catalysis, and might be a reason for the increased binding affinity and catalytic efficiency compared to cationic bovine trypsin.

Animals↗
PubMed2005

Comparative molecular dynamics simulations of HIV-1 integrase and the T66I/M154I mutant: binding modes and drug resistance to a diketo acid inhibitor.

HIV-1 IN is an essential enzyme for viral replication and an interesting target for the design of new pharmaceuticals for use in multidrug therapy of AIDS. L-731,988 is one of the most active molecules of the class of beta-diketo acids. Individual and combined mutations of HIV-1 IN at residues T66, S153, and M154 confer important degrees of resistance to one or more inhibitors belonging to this class. In an effort to understand the molecular mechanism of the resistance of T66I/M154I IN to the inhibitor L-731,988 and its specific binding modes, we have carried out docking studies, explicit solvent MD simulations, and binding free energy calculations. The inhibitor was docked against different protein conformations chosen from prior MD trajectories, resulting in 2 major orientations within the active site. MD simulations have been carried out for the T66I/M154I DM IN, DM IN in complex with L-731,988 in 2 different orientations, and 1QS4 IN in complex with L-731,988. The results of these simulations show a similar dynamical behavior between T66I/M154I IN alone and in complex with L-731,988, while significant differences are observed in the mobility of the IN catalytic loop (residues 138-149). Water molecules bridging the inhibitor to residues from the active site have been identified, and residue Gln62 has been found to play an important role in the interactions between the inhibitor and the protein. This work provides information about the binding modes of L-731,988, as well as insight into the mechanism of inhibitor-resistance in HIV-1 integrase.

Amino Acid Substitution↗
PubMed2006

Structure of the virus capsid protein VP1 of enterovirus 71 predicted by some homology modeling and molecular docking studies.

The homology modeling technique has been used to construct the structure of enterovirus 71 (EV 71) capsid protein VP1. The protein is consisted of 297 amino acid residues and treated as the target. The amino acid sequence identity between the target protein and sequences of template proteins 1EAH, 1PIV, and 1D4M searched from NCBI protein BLAST and WorkBench protein tools were 38, 37, and 36%, respectively. Based on these template structures, the protein model was constructed by using the InsightII/Homology program. The protein model was briefly refined by energy minimization and molecular dynamics (MD) simulation steps. The protein model was validated using some web available servers such as ERRAT, PROCHECK, PROVE, and PROSA2003. However, an inconsistency between the docking scores and the measured activity was observed for a series of EV 71 VP1 inhibitors synthesized by Shia et al. (J Med Chem 2002, 45, 1644) and docked into the binding pocket of the protein model using the DOCK 4.0.2 program. The protein model with an EV 71 VP1 inhibitor docked and engulfed was then refined further by some MD simulation steps in the presence of water molecules. The docking scores obtained for these inhibitors after such a MD refinement were well correlated with the activities. The structure-activity relationships for the ligand-protein model system was also analyzed using the GRID-VOLSURF programs and the corresponding noncrossvalidated and crossvalidated (by leave-one-out) r2 and q2 were 0.99 and 0.61, respectively. The hydrophobic nature of the binding pocket of the protein model was also examined using the GRID21 program. The possibility of improving the potency of the current series of EV 71 VP1 inhibitors was discussed based on all the studies presented.

Amino Acid Sequence↗
PubMed2003

Theoretical study of a new DNA structure: the antiparallel Hoogsteen duplex.

The structure of a new form of duplex DNA, the antiparallel Hoogsteen duplex, is studied in polyd(AT) sequences by means of state-of-the-art molecular dynamics simulations in aqueous solution. The structure, which was found to be stable in all of the simulations, has many similarities with the standard Watson-Crick duplex in terms of general structure, flexibility, and molecular recognition patterns. Accurate MM-PB/SA (and MM-GB/SA) analysis shows that the new structure has an effective energy similar to that of the B-type duplex, while it is slightly disfavored by intramolecular entropic considerations. Overall, MD simulations strongly suggest that the antiparallel Hoogsteen duplex is an accessible structure for a polyd(AT) sequence, which might compete under proper experimental conditions with normal B-DNA. MD simulations also suggest that chimeras containing Watson-Crick duplex and Hoogsteen antiparallel helices might coexist in a common structure, but with the differential characteristics of both type of structures preserved.

DNA↗
PubMed1995

Molecular dynamics simulation of anesthetic-phospholipid bilayer interactions.

To probe the hypothesis of a lipid-mediated mechanism of general anesthetic action on a molecular level, and to help elucidate the nature of the interactions of bioactive compounds with membranes, the effects of trichloroethylene (TCE), an inhalational general anesthetic, on a dioleoylphosphatidylcholine (DOPC) lipid bilayer have been investigated by molecular dynamics (MD) simulations at 37 degrees C and 1 atm and the results compared with 31P and 2H NMR experimental studies (Ref 1). The model used included a single TCE molecule embedded in a lipid bilayer consisting of 24 DOPC molecules and an 8 A layer of explicit water of solvation in each polar head group region of the bilayer, together with constant-pressure periodic boundary conditions in three dimensions. A comparison of the bilayer properties calculated in the presence and absence of the anesthetic led to the detection of three major perturbations of the bilayer caused by the anesthetic at 1 atm: i) an increase in the ratio of the effective areas of hydrocarbon tails and the head group per lipid, predicting the tendency of lipids near the anesthetic site of action to form a hexagonal phase (HII); ii) a slight increase in the frequency of chain dihedral angles found in the gauche conformation; and iii) a significant increase in the lateral mean-square displacement of lipid molecules, an indication of increased lipid lateral diffusion and membrane fluidity. The pressure antagonism of these effects was also studied by MD simulations at pressures of 200 and 400 atm. The study of the pressure reversibility of these effects at 200 and 400 atm indicated that they were partially prevented at 200 atm and essentially blocked at 400 atm, suggesting their probable relevance to the pressure reversal effect seen with general anesthesia. These results may thus provide insights into the interaction between general anesthetics and similar small organic molecules with membranes.

Computer Graphics↗
PubMed2006

Including quantum effects in the dynamics of complex (i.e., large) molecular systems.

The development in the 1950s and 1960s of crossed molecular beam methods for studying chemical reactions at the single-collision molecular level stimulated the need and desire for theoretical methods to describe these and other dynamical processes in molecular systems. Chemical dynamics theory has made great strides in the ensuing decades, so that methods are now available for treating the quantum dynamics of small molecular systems essentially completely. For the large molecular systems that are of so much interest nowadays (e.g., chemical reactions in solution, in clusters, in nanostructures, in biological systems, etc.), however, the only generally available theoretical approach is classical molecular dynamics (MD) simulations. Much effort is currently being devoted to the development of approaches for describing the quantum dynamics of these complex systems. This paper reviews some of these approaches, especially the use of semiclassical approximations for adding quantum effects to classical MD simulations, also showing some new versions that should make these semiclassical approaches even more practical and accurate.

Journal Article↗
PubMed2003

Sequence specific DNA binding of Ets-1 transcription factor: molecular dynamics study on the Ets domain--DNA complexes.

Molecular dynamics (MD) simulations for Ets-1 ETS domain-DNA complexes were performed to investigate the mechanism of sequence-specific recognition of the GGAA DNA core by the ETS domain. Employing the crystal structure of the Ets-1 ETS domain-DNA complex as a starting structure we carried out MD simulations of: (i). the complex between Ets-1 ETS domain and a 14 base-pair DNA containing GGAA core sequence (ETS-GGAA); (ii). the complex between the ETS domain and a DNA having single base-pair mutation, GGAG sequence (ETS-GGAG); and (iii). the 14 base-pair DNA alone (GGAA). Comparative analyses of the MD structures of ETS-GGAA and ETS-GGAG reveal that the DNA bending angles and the ETS domain-DNA phosphate interactions are similar in these complexes. These results support that the GGAA core sequence is distinguished from the mutated GGAG sequence by a direct readout mechanism in the Ets-1 ETS domain-DNA complex. Further analyses of the direct contacts in the interface between the helix-3 region of Ets-1 and the major groove of the core DNA sequence clearly show that the highly conserved arginine residues, Arg391 and Arg394, play a critical role in binding to the GGAA core sequence. These arginine residues make bidentate contacts with the nucleobases of GG dinucleotides in GGAA core sequence. In ETS-GGAA, the hydroxyl group of Tyr395 is hydrogen bonded to N7 nitrogen of A(3) (the third adenosine in the GGAA core), while the hydroxyl group makes a contact with N4 nitrogen of C(4') (the complementary nucleotide of the fourth guanosine G(4) in the GGAG sequence) in the ETS-GGAG complex. We have found that this difference in behavior of Tyr395 results in the relatively large motion of helix-3 in the ETS-GGAG complex, causing the collapse of bidentate contacts between Arg391/Arg394 and the GG dinucleotides in the GGAG sequence.

Amino Acid Sequence↗
PubMed2003

Rigid-body dynamics in the isothermal-isobaric ensemble: a test on the accuracy and computational efficiency.

We have developed a time-reversible rigid-body (rRB) molecular dynamics algorithm in the isothermal-isobaric (NPT) ensemble. The algorithm is an extension of rigid-body dynamics [Matubayasi and Nakahara, J Chem Phys 1999, 110, 3291] to the NPT ensemble on the basis of non-Hamiltonian statistical mechanics [Martyna, G. J. et al., J Chem Phys 1994, 101, 4177]. A series of MD simulations of water as well as fully hydrated lipid bilayer systems have been undertaken to investigate the accuracy and efficiency of the algorithm. The rRB algorithm was shown to be superior to the state-of-the-art constraint-dynamics algorithm SHAKE/RATTLE/ROLL, with respect to computational efficiency. However, it was revealed that both algorithms produced accurate trajectories of molecules in the NPT as well as NVT ensembles, as long as a reasonably short time step was used. A couple of multiple time-step (MTS) integration schemes were also examined. The advantage of the rRB algorithm for computational efficiency increased when the MD simulation was carried out using MTS on parallel processing computer systems; total computer time for MTS-MD of a lipid bilayer using 64 processors was reduced by about 40% using rRB instead of SHAKE/RATTLE/ROLL.

Journal Article↗
PubMed1994

A molecular dynamics approach for the generation of complete protein structures from limited coordinate data.

Generation of full protein coordinates from limited information, e.g., the C alpha coordinates, is an important step in protein homology modeling and structure determination, and molecular dynamics (MD) simulations may prove to be important in this task. We describe a new method, in which the protein backbone is built quickly in a rather crude way and then refined by minimization techniques. Subsequently, the side chains are positioned using extensive MD calculations. The method is tested on two proteins, and results compared to proteins constructed using two other MD-based methods. In the first method, we supplemented an existing backbone building method with a new procedure to add side chains. The second one largely consists of available methodology. The constructed proteins are compared to the corresponding X-ray structures, which became available during this study, and they are in good agreement (backbone RMS values of 0.5-0.7 A, and all-atom RMS values of 1.5-1.9 A). This comparative study indicates that extensive MD simulations are able, to some extent, to generate details of the native protein structure, and may contribute to the development of a standardized methodology to predict reliably (parts of) protein structures when only partial coordinate data are available.

Amino Acids↗
PubMed2003

The role of the dielectric barrier in narrow biological channels: a novel composite approach to modeling single-channel currents.

A composite continuum theory for calculating ion current through a protein channel of known structure is proposed, which incorporates information about the channel dynamics. The approach is utilized to predict current through the Gramicidin A ion channel, a narrow pore in which the applicability of conventional continuum theories is questionable. The proposed approach utilizes a modified version of Poisson-Nernst-Planck (PNP) theory, termed Potential-of-Mean-Force-Poisson-Nernst-Planck theory (PMFPNP), to compute ion currents. As in standard PNP, ion permeation is modeled as a continuum drift-diffusion process in a self-consistent electrostatic potential. In PMFPNP, however, information about the dynamic relaxation of the protein and the surrounding medium is incorporated into the model of ion permeation by including the free energy of inserting a single ion into the channel, i.e., the potential of mean force along the permeation pathway. In this way the dynamic flexibility of the channel environment is approximately accounted for. The PMF profile of the ion along the Gramicidin A channel is obtained by combining an equilibrium molecular dynamics (MD) simulation that samples dynamic protein configurations when an ion resides at a particular location in the channel with a continuum electrostatics calculation of the free energy. The diffusion coefficient of a potassium ion within the channel is also calculated using the MD trajectory. Therefore, except for a reasonable choice of dielectric constants, no direct fitting parameters enter into this model. The results of our study reveal that the channel response to the permeating ion produces significant electrostatic stabilization of the ion inside the channel. The dielectric self-energy of the ion remains essentially unchanged in the course of the MD simulation, indicating that no substantial changes in the protein geometry occur as the ion passes through it. Also, the model accounts for the experimentally observed saturation of ion current with increase of the electrolyte concentration, in contrast to the predictions of standard PNP theory.

Cell Membrane Permeability↗
PubMed2001

The role of cation-pi interactions in biomolecular association. Design of peptides favoring interactions between cationic and aromatic amino acid side chains.

Cation-pi interactions between amino acid side chains are increasingly being recognized as important structural and functional features of proteins and other biomolecules. Although these interactions have been found in static protein structures, they have not yet been detected in dynamic biomolecular systems. We determined, by (1)H NMR spectroscopic titrations, the energies of cation-pi interactions of the amino acid derivative AcLysOMe (1) with AcPheOEt (2) and with AcTyrOEt (3) in aqueous and three organic solvents. The interaction energy is substantial; it ranges from -2.1 to -3.4 kcal/mol and depends only slightly on the dielectric constant of the solvent. To assess the effects of auxiliary interactions and structural preorganization on formation of cation-pi interactions, we studied these interactions in the association of pentapeptides. Upon binding of the positively-charged peptide AcLysLysLysLysLysNH(2) (5) to the negatively-charged partner AcAspAspXAspAspNH(2) (6), in which X is Leu (6a), Tyr (6b), and Phe (6c), multiple interactions occur. Association of the two pentapeptides is dynamic. Free peptides and their complex are in fast exchange on the NMR time-scale, and 2D (1)H ROESY spectra of the complex of the two pentapeptides do not show intermolecular ROESY peaks. Perturbations of the chemical shifts indicated that the aromatic groups in peptides 6b and 6c were affected by the association with 5. The association constants K(A) for 5 with 6a and with 6b are nearly equal, (4.0 +/- 0.7) x 10(3) and (5.0 +/- 1.0) x 10(3) M(-)(1), respectively, while K(A) for 5 with 6c is larger, (8.3 +/- 1.3) x 10(3) M(-)(1). Molecular-dynamics (MD) simulations of the pentapeptide pairs confirmed that their association is dynamic and showed that cation-pi contacts between the two peptides are stereochemically possible. A transient complex between 5 and 6 with a prominent cation-pi interaction, obtained from MD simulations, was used as a template to design cyclic peptides C(X) featuring persistent cation-pi interactions. The cyclic peptide C(X) had a sequence in which X is Tyr, Phe, and Leu. The first two peptides do, but the third does not, contain the aromatic residue capable of interacting with a cationic Lys residue. This covalent construct offered conformational stability over the noncovalent complexes and allowed thorough studies by 2D NMR spectroscopy. Multiple conformations of the cyclic peptides C(Tyr) and C(Phe) are in slow exchange on the NMR time-scale. In one of these conformations, cation-pi interaction between Lys3 and Tyr9/Phe9 is clearly evident. Multiple NOEs between the side chains of residues 3 and 9 are observed; chemical-shift changes are consistent with the placement of the side chain of Lys3 over the aromatic ring. In contrast, the cyclic peptide C(Leu) showed no evidence for close approach of the side chains of Lys3 and Leu9. The cation-pi interaction persists in both DMSO and aqueous solvents. When the disulfide bond in the cyclic peptide C(Phe) was removed, the cation-pi interaction in the acyclic peptide AC(Phe) remained. To test the reliability of the pK(a) criterion for the existence of cation-pi interactions, we determined residue-specific pK(a) values of all four Lys side chains in all three cyclic peptides C(X). While NOE cross-peaks and perturbations of the chemical shifts clearly show the existence of the cation-pi interaction, pK(a) values of Lys3 in C(Tyr) and in C(Phe) differ only marginally from those values of other lysines in these dynamic peptides. Our experimental results with dynamic peptide systems highlight the role of cation-pi interactions in both intermolecular recognition at the protein-protein interface and intramolecular processes such as protein folding.

Amino Acid Sequence↗
PubMed2005

Insights into the induced fit mechanism in antithrombin-heparin interaction using molecular dynamics simulations.

Heparin was isolated in the beginning of the 20th century and until today remains as one of the most important drugs able to interfere with the haemostatic process. Due to the side effects produced by heparin therapy, new promising drugs have been developed, as the synthetic pentasaccharide (synthetically derived from the sequence GlcN-GlcA-GlcN-IdoA-GlcN). The anticoagulant activity of this compound is based on potentiation of antithrombin (AT) inhibitory activity upon serine proteinases of clotting cascade, a mechanism based on the conformational modification of AT. In this context, we present here a molecular dynamics (MD) study of the interaction between the synthetic pentasaccharide and AT. The obtained data correctly predicted an induced fit mechanism in AT-pentasaccharide interaction, showing a solvent-exposed P1 residue instead of a hided conformation. Also, the specific contribution of important amino acid residues to the overall process was also characterized, both in (2)S(0) and (1)C(4) conformations of IdoA residue, suggesting that there is no conformational requirement to the interaction of this residue with AT. Altogether, the results show that MD simulations could be used to characterize and quantify the interaction of synthetic compounds with AT, predicting its specific capacity to induce conformational changes in AT structure. Thus, MD simulations of heparin (and heparin-derived)-AT interactions are proposed here as a powerful tool to assist and support drug design of new antithrombotic agents.

Allosteric Regulation↗
PubMed1995

Interaction of an amphiphilic peptide with a phospholipid bilayer surface by molecular dynamics simulation study.

Corticotropin-releasing factor (CRF) is the principal neuroregulator of adrenocorticotropic hormone (ACTH) secretion. Previous experiments have demonstrated that CRF binds avidly to the surface of single egg phosphatidylcholine vesicles and its amphiphilic secondary structure might play an important role in the function. In this study, the interaction of the residues 13-41 in human CRF with the surface of a DOPC bilayer was investigated by molecular dynamics (MD) simulation in order to understand the role of the membrane surface in the formation of the amphiphilic alpha helix as well as to determine the effects of the peptide on the lipid bilayer. The model used included 60 DOPC molecules, 1 helical peptide (CRF13-41) on the bilayer surface, and explicit waters of solvation in the lipid polar head group regions, together with constant-volume periodic boundary conditions in three dimensions. The MD simulation was carried out for 510 ps. In addition, CRF13-41, initially in a helical form, was simulated in vacuo as a control. The results indicate that while it was completely unstable in vacuo, the peptide helical form was generally maintained on the bilayer surface, but with distortions near the terminal ends. The peptide was confined to the bilayer headgroup/water region, similar to that reported from neutron diffraction measurement of tripeptides bound to the phosphatidylcholine bilayer surface (Ref 1). The amphiphilicity of the peptide matched that of the bilayer headgroup environment, with the hydrophilic side oriented toward water and the hydrophobic side making contact with the bilayer hydrocarbon core. These results support the hypothesis that the amphiphilic environment of a membrane surface is important in the induction of peptide amphiphilic alpha-helical secondary structure. Two major effects of the peptide on the lipids were found: the first CH2 segment in the lipid chains was significantly disordered and the lipid headgroup distribution was broadened towards the water region.

Amino Acid Sequence↗