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,189 records · Page 66Linked to original sources

Three dimensional model of severe acute respiratory syndrome coronavirus helicase ATPase catalytic domain and molecular design of severe acute respiratory syndrome coronavirus helicase inhibitors.

The modeling of the severe acute respiratory syndrome coronavirus helicase ATPase catalytic domain was performed using the protein structure prediction Meta Server and the 3D Jury method for model selection, which resulted in the identification of 1JPR, 1UAA and 1W36 PDB structures as suitable templates for creating a full atom 3D model. This model was further utilized to design small molecules that are expected to block an ATPase catalytic pocket thus inhibit the enzymatic activity. Binding sites for various functional groups were identified in a series of molecular dynamics calculation. Their positions in the catalytic pocket were used as constraints in the Cambridge structural database search for molecules having the pharmacophores that interacted most strongly with the enzyme in a desired position. The subsequent MD simulations followed by calculations of binding energies of the designed molecules were compared to ATP identifying the most successful candidates, for likely inhibitors - molecules possessing two phosphonic acid moieties at distal ends of the molecule.

Amino Acid Sequence↗

Impact of static and dynamic A-form heterogeneity on the determination of RNA global structural dynamics using NMR residual dipolar couplings.

We examined how static and dynamic deviations from the idealized A-form helix propagate into errors in the principal order tensor parameters determined using residual dipolar couplings (rdcs). A 20-ns molecular dynamics (MD) simulation of the HIV-1 transactivation response element (TAR) RNA together with a survey of spin relaxation studies of RNA dynamics reveals that pico-to-nanosecond local motions in non-terminal Watson-Crick base-pairs will uniformly attenuate base and sugar one bond rdcs by approximately 7%. Gaussian distributions were generated for base and sugar torsion angles through statistical comparison of 40 RNA X-ray structures solved to <3.0 A resolution. For a typical number (>or=11) of one bond C-H base and sugar rdcs, these structural deviations together with rdc uncertainty (1.5 Hz) lead to average errors in the magnitude and orientation of the principal axis of order that are <9% and <4 degrees, respectively. The errors decrease to <5% and <4 degrees for >or=17 rdcs. A protocol that allows for estimation of error in A-form order tensors due to both angular deviations and rdc uncertainty (Aform-RDC) is validated using theoretical simulations and used to analyze rdcs measured previously in TAR in the free state and bound to four distinct ligands. Results confirm earlier findings that the two TAR helices undergo large changes in both their mean relative orientation and dynamics upon binding to different targets.

Base Pairing↗

Effect of constraints by threonine on proline containing alpha-helix--a molecular dynamics approach.

Proline plays an important role in the secondary structure of proteins. In the pursuit of understanding its structural role, Proline containing helices with constraints have been studied by employing molecular dynamics (MD) technique. In the present study, the constraint introduced is a threonine residue, whose sidechain has intramolecular hydrogen bond interaction with the backbone oxygen atom. The three systems that have been chosen for characterization are: (1) Ace-(Ala)12-Thr-Pro-(Ala)10-NHMe, (2) Ace-(Ala)13-Pro-Ala-Thr-(Ala)8-NHMe and (3) Ace-(Ala)13-Pro-(Ala)3-Thr-(Ala)6-NHMe. The equilibrium structures and structural transitions have been identified by monitoring the backbone dihedral angles, bend related parameters and the hydrogen bond interactions. The MD averages and root mean square (r.m.s.) fluctuations are compared and discussed. Energy minimization has been carried out on selected MD simulated points in order to analyze the characteristics of different conformations.

Amino Acid Sequence↗

Adjusting force distributions in functional site of scorpion toxin BMK M1 by cooperative effect of disulfide bonds.

We decided to investigate the influence of the presence of disulfide bonds on the pattern of force distribution in the functional site of scorpion toxin BMK M1 in its functional state. Therefore, a series of short time molecular dynamics (MD) simulations were performed on this toxin in the native state and disulfide bond broken states. The comparison of disulfide bond broken states with the native state showed that the electrostatic potential energy of important functional residues in the reverse turn and C-terminal regions were modulated by the cooperative effect of all disulfide bonds in the molecule. Furthermore, our results revealed that disulfide bonds also play a cooperative role in modulating (1) the amplitude of the fluctuations of the functional segments and (2) the correlation of motions between important functional residue pairs in this toxin. Therefore, we can conclude that the disulfide bonds have cooperation to adjust the pattern of force distribution in the functional site of this toxin in its functional state.

Animals↗

Intrinsic activity and comparative molecular dynamics of buspirone analogues at the 5-HT(1A) receptors.

In CNS, the 5-hydroxytryptamine(1A) (5-HT(1A)) receptors exist in two different populations with different behavioural and physiological effects: (1) somatodendritic autoreceptors located pre-synaptically of 5-HT containing neurons and (2) receptors located post-synaptic to 5-HT containing neurons. Clinical studies have shown that 5-HT(1A) partial agonists have anxiolytic properties, while antagonists of pre-synaptical autoreceptors shorten the onset time of selective serotonin reuptake inhibitors (SSRIs). In the present study, the pre- and post-synaptic activity of structural analogues of buspirone was evaluated in animal models. A three dimensional model of the 5-HT(1A) receptor was used to study their interaction modes and helical displacements upon receptor binding. The predicted receptor-ligand interactions indicated similarities in the receptor binding modes for all buspirone analogues, and no clear relationship between receptor contact residues and activity at pre- and post-synaptic receptors. Comparative molecular dynamics (MD) simulations for 650ps indicated that pre-synaptic antagonistic behaviour is connected to large displacements of transmembrane helix (TMH) 7 upon binding, while pre-synaptic agonistic behaviour is connected to large displacements of TMH2 and small displacements of TMH7. Post-synaptic partial agonist behaviour is connected to large displacements of TMH4 and TMH5 upon binding, while post-synaptic antagonists only slightly displace these helices.

Animals↗

Discovery of NAT-6-321056 as a novel modulator of VEGFR2 signaling to suppress tumor angiogenesis.

Vascular endothelial growth factor receptor 2 (VEGFR2) is a master regulator of angiogenesis and cancer progression. However, current VEGFR2 modulators face significant challenges, including off-target toxicity and acquired resistance, underscoring the urgent need for novel therapeutic agents with improved efficacy and safety profiles. Here, we reported that virtual screening of 39,442 natural products from the ZINC natural products-derived library, coupled with molecular docking and molecular dynamics (MD) simulations to evaluate the binding stability of candidate compounds, identified NAT-6-321056 as a highly promising modulator of VEGFR2 signaling. Biological evaluations demonstrated that NAT-6-321056 exerted potent inhibition on the growth of a broad spectrum of cancer cells, including both solid tumors and hematological malignancies. In EA.hy 926 endothelial cells and SK-N-DZ neuroblast cells, the compound significantly suppressed proliferation, migration, and invasion. Microscale thermophoresis (MST) confirmed direct binding of NAT-6-321056 to VEGFR2 with favorable affinity. Kinase profiling against a panel of 33 kinases indicated that NAT-6-321056 exhibited a multi-kinase modulation profile. Mechanistic studies revealed that NAT-6-321056 suppressed the expression of hypoxia-inducible factor 1-alpha (HIF-1&#x3b1;) and was associated with reduced VEGFR2 phosphorylation and attenuation of the downstream ERK/JNK/AKT signaling pathways. Moreover, NAT-6-321056 exhibited robust in vivo anti-angiogenic effects in both the chick chorioallantoic membrane (CAM) assay and transgenic zebrafish vascular fluorescence imaging models. Computational absorption, distribution, metabolism, excretion, and toxicity (ADMET) prediction suggested acceptable drug-like properties. Collectively, these findings demonstrated that NAT-6-321056 is a promising modulator of VEGFR2 signaling with potent anti-angiogenic activity and represents a viable candidate for cancer therapy.

Vascular Endothelial Growth Factor Receptor-2↗

The conserved residue Phe273(282) of PPARalpha(gamma), beyond the ligand-binding site, functions in binding affinity through solvation effect.

Peroxisome proliferator-activated receptors (PPARs) belong to the members of the nuclear receptor superfamily, and play important roles in lipid and glucose homeostasis. Residue Phe282 in PPARgamma (Phe273 in PPARalpha), beyond the ligand-binding site, is a conserved amino acid across several nuclear receptors and in all PPAR subfamily. In this work, we firstly investigated the influence of Phe282(273)Ala mutation on the binding affinity of PPARgamma(alpha) against a series of agonists by use of surface plasmon resonance (SPR) technique and cellular transcriptional activation analysis. Phe282(273)Ala mutation decreases the binding affinities of the ligands to the receptors in certain degrees, from several to 1000-folds. Phe282Ala mutation dramatically reduced the binding affinity of PPARgamma to GI262570, however, this mutation did not affect PPARalpha binding to this ligand, thereby suggesting that the Phe282 and Phe273 are associated with the selectivity of GI262570 binding to PPARgamma and PPARalpha. The mutation reduced the transcriptional activation activities of the receptors induced by the ligand binding, and the decrease degree is generally in agreement with the binding affinities of the ligands to the receptors. The 5 ns MD simulations for the wild-type and mutated PPARgamma showed that the mutation did not influence the flexibility of the receptor. There is no repulsion between Phe282 and the proceeding loop of AF2. However, substitution of Phe282 by alanine enlarged the entrance of the binding pocket and abolished the repulsive interaction between solvent water molecules and this hydrophobic residue, thus more water molecules can enter into the binding pocket. It needs more energy to exclude the extra water molecules for a ligand binding to the mutated receptor. In addition, the extra water molecules abolish some of H-bonds between the ligand and receptors. Therefore, solvent effect may be concluded as the major source of the decrease of binding affinity for the mutated receptors to the ligands, and thereby of the decrease of their transcriptional activation activities.

Amino Acid Sequence↗

In situ fabrication of self-transformable and hydrophilic poly(ethylene glycol) derivative-modified polysulfone membranes.

A self-transformable sulfonated poly(ethylene glycol) acrylate diblock copolymer (PEG-SO3A/OA) entrapped into polysulfone membrane was studied. The asymmetric membrane structure was prepared by a phase inversion process. The induced hydrophilicity by reorientation of diblock copolymer at the interface was evaluated by contact angle measurement, platelet adhesion test, and electron spectroscopy for chemical analysis (ESCA) depth profiling with ion sputtering. Molecular dynamic (MD) simulations as a function of copolymer density were also performed to obtain optimum interfacial structure information. The dependency of water clustering behavior as a hydrophilicity parameter was described in terms of an atom-atom radial distribution function (RDF). The results showed that the sulfonated diblock copolymer enhanced the hydrophilicity and long-term stability more than the copolymer having no hydrophobic block. Also, according to the ESCA, oxygen composition significantly began to decrease along the membrane depth, indicating the reorientation of diblock chains. The copolymer-entrapped surfaces significantly induced the degree of water clustering, and the resulting equilibrium rearrangement of interfacial structure was distinctly dependent upon the density of the copolymer.

Biocompatible Materials↗

Essential structural profile of a dual functional inhibitor against cyclooxygenase-2 (COX-2) and 5-lipoxygenase (5-LOX): molecular docking and 3D-QSAR analyses on DHDMBF analogues.

It is recently proposed that compounds with equal capabilities of inhibiting COX and 5-LOX, both are key enzymes involved in the arachidonic acid (AA) cascade, are expected to be safer non-steroidal anti-inflammatory drugs (NSAIDs). To dig out helpful information in designing dual functional inhibitors against the two enzymes, homology modeling, molecular dynamics (MD) simulations, automated docking, and 3D-QSAR analyses were performed in this study on 21 COX-2/5-LOX dual inhibitors, namely, 7-tert-butyl-2,3-dihydro-3,3-dimethylbenzofuran (DHDMBF) analogues. A 3D-model of 5-LOX was built based on the high-resolution X-ray structure of rabbit reticulocyte 15-lipoxygenase. Molecular docking was then applied to locate the binding orientations and conformations of DHDMBF analogues with COX-2 and 5-LOX, respectively, leading to highly predictive CoMFA models constructed on the basis of the binding conformations with q2 values of 0.782 and 0.634 for COX-2 and 5-LOX, respectively. In addition, CoMFA field distributions were found in good agreement with the structural characteristics of the corresponding binding sites. Both the docking simulations and QSAR analyses suggest that new potent dual inhibitors should share a structural feature with a moderately bulky group at R2 position and a rather negatively charged group around the position of the carbonyl group of DHDMBFs. Therefore, the final 3D-QSAR models and the information of the inhibitor-enzyme interaction should be useful in developing new NSAIDs as anti-inflammation drugs with favorable safety profile.

Amino Acid Sequence↗

Identifying inhibitors of the SARS coronavirus proteinase.

The Severe Acute Respiratory Syndrome (SARS) is a serious respiratory illness that has recently been reported in parts of Asia and Canada. In this study, we use molecular dynamics (MD) simulations and docking techniques to screen 29 approved and experimental drugs against the theoretical model of the SARS CoV proteinase as well as the experimental structure of the transmissible gastroenteritis virus (TGEV) proteinase. Our predictions indicate that existing HIV-1 protease inhibitors, L-700,417 for instance, have high binding affinities and may provide good starting points for designing SARS CoV proteinase inhibitors.

Antiviral Agents↗

Structures and stability of salt-bridge in aqueous solution.

Structures and stability of salt-bridges in aqueous solutions were investigated using a complex formed from the guanidinium (Gdm+) and formate (FmO-) ions as a model system. The Test-particle model (T-model) potentials to describe the interactions in the Gdm+-H2O, FmO(-)-H2O and Gdm+-FmO- complexes were constructed, tested and applied in molecular dynamics (MD) simulations of the aqueous solutions at 298 K. The three-dimensional structures and energetic of the hydrogen bond (H-bond) networks of water in the first hydration shells of the Gdm+ and FmO- ions, as well as the Gdm+-FmO- complex, were visualized and analyzed using various probability distribution (PD) maps. The structures of the average potential energy landscapes at the H-bond networks were employed to characterize the stability and dynamic behavior of water molecules in the first hydration shells of the solutes. It was observed that water molecules in the first hydration shell of the close-contact Gdm+-FmO- complex form associated H-bond networks, which introduce a net stabilization effect to the ion-pair, whereas those in the interstitial H-bond network destabilize and break the solvent-separated Gdm+-FmO- complex. The present results showed that, in order to provide complete insights into the structures and stability of ion-pairs in aqueous solutions, explicit water molecules have to be included in the model calculations.

Crystallography, X-Ray↗

Bead-like passage of chloride ions through ClC chloride channels.

The ClC chloride channels control the ionic composition of the cytoplasm and the volume of cells, and regulate electrical excitability. Recently, it has been proposed that prokaryotic ClC channels are H+-Cl- exchange transporter. Although X-ray and molecular dynamics (MD) studies of bacterial ClC channels have investigated the filter open-close and ion permeation mechanism of channels, details have remained unclear. We performed MD simulations of ClC channels involving H+, Na+, K+, or H3O+ in the intracellular region to elucidate the open-close mechanism, and to clarify the role of H+ ion an H+-Cl- exchange transporter. Our simulations revealed that H+ and Na+ caused channel opening and the passage of Cl- ions. Na+ induced a bead-like string of Cl- -Na+-Cl--Na+-Cl- ions to form and permeate through ClC channels to the intracellular side with the widening of the channel pathway.

Bacterial Proteins↗

Ion selectivity in potassium channels.

Potassium channels are tetrameric membrane-spanning proteins that provide a selective pore for the conduction of K(+) across the cell membranes. One of the main physiological functions of potassium channels is efficient and very selective transport of K(+) ions through the membrane to the cell. Classical views of ion selectivity are summarized within a historical perspective, and contrasted with the molecular dynamics (MD) simulations free energy perturbation (FEP) performed on the basis of the crystallographic structure of the KcsA phospholipid membrane. The results show that the KcsA channel does not select for K(+) ions by providing a binding site of an appropriate (fixed) cavity size. Rather, selectivity for K(+) arises directly from the intrinsic local physical properties of the ligands coordinating the cation in the binding site, and is a robust feature of a pore symmetrically lined by backbone carbonyl groups. Further analysis reveals that it is the interplay between the attractive ion-ligand (favoring smaller cation) and repulsive ligand-ligand interactions (favoring larger cations) that is the basic element governing Na(+)/K(+) selectivity in flexible protein binding sites. Because the number and the type of ligands coordinating an ion directly modulate such local interactions, this provides a potent molecular mechanism to achieve and maintain a high selectivity in protein binding sites despite a significant conformational flexibility.

Bacterial Proteins↗

Exploration of the conformational space of myosin recovery stroke via molecular dynamics.

Muscle contractions are driven by cyclic conformational changes of myosin, whose molecular mechanisms of operation are being elucidated by recent advances in crystallographic studies and single molecule experiments. To complement such structural studies and consider the energetics of the conformational changes of myosin head, umbrella sampling molecular dynamics (MD) simulations were performed with the all-atom model of the scallop myosin sub-fragment 1 (S1) with a bound ATP in solution in explicit water using the crystallographic near-rigor and transition state conformations as two references. The constraints on RMSD reaction coordinates used for the umbrella sampling were found to steer the conformational changes efficiently, and relatively close correlations have been observed between the set of characteristic structural changes including the lever arm rotation and the closing of the nucleotide binding pocket. The lever arm angle and key residue interaction distances in the nucleotide binding pocket and the relay helix show gradual changes along the recovery stroke reaction coordinate, consistent with previous crystallographic and computational minimum energy studies. Thermal fluctuations, however, appear to make the switch-2 coordination of ATP more flexible than suggested by crystal structures. The local solvation environment of the fluorescence probe, Trp 507 (scallop numbering), also appears highly mobile in the presence of thermal fluctuations.

Animals↗

Interaction energy decomposition in protein-protein association: a quantum mechanical study of barnase-barstar complex.

Protein-protein interactions are very important in the function of a cell. Computational studies of these interactions have been of interest, but often they have utilized classical modelling techniques. In recent years, quantum mechanical (QM) treatment of entire proteins has emerged as a powerful approach to study biomolecular systems. Herein, we apply a semi-empirical divide and conquer (DC) methodology coupled with a dielectric continuum model for the solvent, to explore the contribution of electrostatics, polarization and charge transfer to the interaction energy between barnase and barstar in their complex form. Molecular dynamic (MD) simulation was performed to account for the dynamic behavior of the complex. The results show that electrostatics, charge transfer and polarization favor the formation of the complex. Our study shows that electrostatics dominates the interaction between barnase and barstar ( approximately 73%), while charge transfer and polarization are approximately 21% and approximately 6%, respectively. Close inspection of the polarization and charge-transfer effects on the charge distribution of the complex reveals the existence of two, well localized, regions in barstar. The first region includes the residues between P27 and Y47 and the second region is between N65 and D83. Since no such regions could be detected in barnase clearly suggests that barstar is well optimized for efficiently binding barnase. Furthermore, using our interaction energy decomposition scheme, we were able to identify all residues that have been experimentally determined to be important for the complex formation and to suggest other residues never have been investigated. This suggests that our approach will be useful as an aid in further understanding protein-protein contacts for the ultimate goal to produce successful inhibitors for protein complexes.

Bacterial Proteins↗

The conformational free-energy map for solvated neocarrabiose.

A Ramachandran map of the conformational potential of mean force (pmf) for neocarrabiose in water was obtained using molecular dynamics (MD) simulations with umbrella sampling. The potential energy map calculated in a previous study for this molecule in vacuum exhibited a global minimum located at (phi = 81 degrees, psi = -141 degrees). However, the global minimum on the new pmf map in aqueous solution is located in an area centered around (phi = 175 degrees, psi = 180 degrees), indicating a considerable solvent shift. This new global minimum-energy solution conformation was found to correspond to the experimental value obtained from NMR-NOE measurements, and is also consistent with the experimental crystal structure for neocarrabiose and the fiber diffraction conformation for iota-carrageenan. The global minimum of the solution pmf and its local topology were found to be approximately reproduced by quick vacuum conformational energy mapping using several approximations that mimic solvation effects by de-emphasizing intramolecular hydrogen bonding.

Carrageenan↗

Molecular dynamics and atomic charge calculations in the study of heparin conformation in aqueous solution.

HF/6-31G** and molecular dynamics (MD) simulations were used to evaluate the performance of different atomic charge basis sets (i.e., Mulliken, Lowdin, and Electrostatic Potential Derived Charges--ESP) in heparin simulations. HF/3-21 G calculations were also used to study the NMR conformation of the IdoA residue. The results thus obtained indicated that ESP and Lowdin charges gave the better results in heparin simulations, followed by Mulliken charges, and that the minimum-energy conformation of IdoA can be different from that observed by NMR spectroscopy by less than 1 Angstrom. However, it was found that this small conformational modification is capable of inducing a change of almost 200 kJ/mol in the interactions of heparin with the surrounding environment, which is a meaningful amount of energy in the context of ligand-receptor interactions. This information can be potentially of great relevance in the design of heparin-derived antithrombotic compounds.

Carbohydrate Conformation↗

A distance ruler for RNA using EPR and site-directed spin labeling.

As a basic model study for measuring distances in RNA molecules using continuous wave (CW) EPR spectroscopy, site-directed spin-labeled 10-mer RNA duplexes and HIV-1 TAR RNA motifs with various interspin distances were examined. The spin labels were attached to the 2'-NH2 positions of appropriately placed uridines in the duplexes, and interspin distances were measured from both molecular dynamics simulations (MD) and Fourier deconvolution methods (FD). The 10-mer duplexes have interspin distances ranging from 10 A to 30 A based on MD; however, dipolar line broadening of the CW EPR spectrum is only observed for the RNAs for predicted interspin distances of 10-21 A and not for distances over 25 A. The conformational changes in TAR (transactivating responsive region) RNA in the presence and in the absence of different divalent metal ions were monitored by measuring distances between two nucleotides in the bulge region. The predicted interspin distances obtained from the FD method and those from MD calculations match well for both the model RNA duplexes and the structural changes predicted for TAR RNA. These results demonstrate that distance measurement using EPR spectroscopy is a potentially powerful method to help predict the structures of RNA molecules.

Base Sequence↗