Classical molecular dynamics
Dotaz
Zobrazit nápovědu
Fluidity of lipid membranes is known to play an important role in the functioning of living organisms. The fluorescent probe Laurdan embedded in a lipid membrane is typically used to assess the fluidity state of lipid bilayers by utilizing the sensitivity of Laurdan emission to the properties of its lipid environment. In particular, Laurdan fluorescence is sensitive to gel vs liquid⁻crystalline phases of lipids, which is demonstrated in different emission of the dye in these two phases. Still, the exact mechanism of the environment effects on Laurdan emission is not understood. Herein, we utilize dipalmitoylphosphatidylcholine (DPPC) and dioleoylphosphatidylcholine (DOPC) lipid bilayers, which at room temperature represent gel and liquid⁻crystalline phases, respectively. We simulate absorption and emission spectra of Laurdan in both DOPC and DPPC bilayers with quantum chemical and classical molecular dynamics methods. We demonstrate that Laurdan is incorporated in heterogeneous fashion in both DOPC and DPPC bilayers, and that its fluorescence depends on the details of this embedding.
- MeSH
- 1,2-dipalmitoylfosfatidylcholin chemie MeSH
- 2-naftylamin analogy a deriváty chemie MeSH
- chemické modely * MeSH
- fluorescence MeSH
- fosfatidylcholiny chemie MeSH
- kvantová teorie MeSH
- laurany chemie MeSH
- lipidové dvojvrstvy chemie MeSH
- simulace molekulární dynamiky * MeSH
- Publikační typ
- časopisecké články MeSH
Abstract Aromatic stacking of nucleic acid bases is one of the key players in determining the structure and dynamics of nucleic acids. The arrangement of nucleic acid bases with extensive overlap of their aromatic rings gave rise to numerous often contradictory suggestions about the physical origins of stacking and the possible role of delocalized electrons in stacked aromatic π systems, leading to some confusion about the issue. The recent advance of computer hardware and software finally allowed the application of state of the art quantum-mechanical approaches with inclusion of electron correlation effects to study aromatic base stacking, now providing an ultimitate qualitative description of the phenomenon. Base stacking is determined by an interplay of the three most commonly encountered molecular interactions: dispersion attraction, electrostatic interaction, and short-range repulsion. Unusual (aromatic- stacking specific) energy contributions were in fact not evidenced and are not necessary to describe stacking. The currently used simple empirical potential form, relying on atom-centered constant point charges and Lennard-Jones van der Waals terms, is entirely able to reproduce the essential features of base stacking. Thus, we can conclude that base stacking is in principle one of the best described interactions in current molecular modeling and it allows to study base stacking in DNA using large-scale classical molecular dynamics simulations. Neglect of cooperativity of stacking appears to be the most serious approximation of the currently used force field form. This review summarizes recent developments in the field. It is written for an audience that is not necessarily expert in computational quantum chemistry and follows up on our previous contribution (Sponer et. al., J. Biomol. Struct. Dyn. 14, 117, (1997)). First, the applied methodology, its accuracy, and the physical nature of base stacking is briefly overviewed, including a comment on the accuracy of other molecular orbital methods and force fields. Then, base stacking is contrasted with hydrogen bonding, the other dominant force in nucleic acid structure. The sequence dependence and cooperativity of base stacking is commented on, and finally a brief introduction into recent progress in large-scale molecular dynamics simulations of nucleic acids is provided. Using four stranded DNA assemblies as an example, we demonstrate the efficacy of current molecular dynamics techniques that utilize refined and verified force fields in the study of stacking in nucleic acid molecules.
Small differences in physical and chemical properties of H2O and D2O, such as melting and boiling points or pKa, can be traced back to a slightly stronger hydrogen bonding in heavy versus normal water. In particular, deuteration reduces zero-point vibrational energies as a demonstration of nuclear quantum effects. In principle, computationally demanding quantum molecular dynamics is required to model such effects. However, as already demonstrated by Feynmann and Hibbs, zero-point vibrations can be effectively accounted for by modifying the interaction potential within classical dynamics. In the spirit of the Feymann-Hibbs approach, we develop here two water models for classical molecular dynamics by fitting experimental differences between H2O and D2O. We show that a three-site SPCE-based model accurately reproduces differences between properties of the two water isotopes, with a four-site TIP4P-2005/based model in addition capturing also the absolute values of key properties of heavy water. The present models are computationally simple enough to allow for extensive simulations of biomolecules in heavy water relevant, for example, for experimental techniques such as NMR or neutron scattering.
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Binding of muscarinic ligands, both antagonists and agonists, and their effects on the conformation of the M2 acetylcholine receptor were modeled in silico and compared to experimental data. After docking of antagonists to the M2 receptor in an inactive conformation (3UON, 5ZK3, 5ZKB, or 5ZKB) and agonists in an active conformation (4MQS), 100 ns of conventional molecular dynamics (MD) followed by 500 ns of accelerated MD was run. Conventional MD revealed ligand-specific interactions with the receptor. Antagonists stabilized the receptor in an inactive conformation during accelerated MD. The receptor in complex with various agonists attained different conformations specific to individual agonists. The magnitude of the TM6 movement correlated with agonist efficacy at the non-preferential Gs pathway. The shape of the intracellular opening where the receptor interacts with a G-protein was different for the classical agonist carbachol, super-agonist iperoxo, and Gi/o-biased partial agonists JR-6 and JR-7, being compatible with experimentally observed agonist bias at the G-protein level. Moreover, a wash-resistant binding of the unique agonist xanomeline associated with interactions with membrane lipids was formed during accelerated MD. Thus, accelerated MD is suitable for modeling of ligand-specific receptor binding and receptor conformations that is essential for the design of experiments aimed at identification of the secondary binding sites and understanding molecular mechanisms underlying receptor activation.
Nowadays, advanced computational chemistry methods offer various strategies for revealing prospective hit structures in drug development essentially through accurate binding free energy predictions. After the era of molecular docking and quantitative structure-activity relationships, much interest has been lately oriented to perturbed molecular dynamic approaches like replica exchange with solute tempering and free energy perturbation (REST/FEP) and the potential of the mean force with adaptive biasing and accelerated weight histograms (PMF/AWH). Both of these receptor-based techniques can exploit exascale CPU&GPU supercomputers to achieve high throughput performance. In this fundamental study, we have compared the predictive power of a panel of supercomputerized molecular modelling methods to distinguish the major binding modes and the corresponding binding free energies of a promising tacrine related potential antialzheimerics in human acetylcholinesterase. The binding free energies were estimated using flexible molecular docking, molecular mechanics/generalized Born surface area/Poisson-Boltzmann surface area (MM/GBSA/PBSA), transmutation REST/FEP with 12 x 5 ns/λ windows, annihilation FEP with 20 x 5 ns/λ steps, PMF with weight histogram analysis method (WHAM) and 40 x 5 ns samples, and PMF/AWH with 10 x 100 ns replicas. Confrontation of the classical approaches such as canonical molecular dynamics and molecular docking with alchemical calculations and steered molecular dynamics enabled us to show how large errors in ΔG predictions can be expected if these in silico methods are employed in the elucidation of a common case of enzyme inhibition.Communicated by Ramaswamy H. Sarma.
Since its early days in the 19th century, medicinal chemistry has concentrated its efforts on the treatment of diseases, using tools from areas such as chemistry, pharmacology, and molecular biology. The understanding of biological mechanisms and signaling pathways is crucial information for the development of potential agents for the treatment of diseases mainly because they are such complex processes. Given the limitations that the experimental approach presents, computational chemistry is a valuable alternative for the study of these systems and their behavior. Thus, classical molecular dynamics, based on Newton's laws, is considered a technique of great accuracy, when appropriated force fields are used, and provides satisfactory contributions to the scientific community. However, as many configurations are generated in a large MD simulation, methods such as Statistical Inefficiency and Optimal Wavelet Signal Compression Algorithm are great tools that can reduce the number of subsequent QM calculations. Accordingly, this review aims to briefly discuss the importance and relevance of medicinal chemistry allied to computational chemistry as well as to present a case study where, through a molecular dynamics simulation of AMPK protein (50 ns) and explicit solvent (TIP3P model), a minimum number of snapshots necessary to describe the oscillation profile of the protein behavior was proposed. For this purpose, the RMSD calculation, together with the sophisticated OWSCA method was used to propose the minimum number of snapshots.
- MeSH
- algoritmy MeSH
- farmaceutická chemie MeSH
- kvantová teorie MeSH
- lidé MeSH
- proteinkinasy aktivované AMP metabolismus chemie MeSH
- simulace molekulární dynamiky * MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- přehledy MeSH
This review provides a critical assessment of the advantages and limitations of modeling methods available for guanine quadruplex (G-DNA) molecules. We characterize the relations of simulations to the experimental techniques and explain the actual meaning and significance of the results. The following aspects are discussed: pair-additive approximation of the empirical force fields, sampling limitations stemming from the simulation time and accuracy of description of base stacking, H-bonding, sugar-phosphate backbone and ions by force fields. Several methodological approaches complementing the classical explicit solvent molecular dynamics simulations are commented on, including enhanced sampling methods, continuum solvent methods, free energy calculations and gas phase simulations. The successes and pitfalls of recent simulation studies of G-DNA are demonstrated on selected results, including studies of cation interactions and dynamics of G-DNA stems, studies of base substitutions (inosine, thioguanine and mixed tetrads), analysis of possible kinetic intermediates in folding pathway of a G-DNA stem and analysis of loop regions of G-DNA molecules.
CONTEXT: Chlordecone (CLD) and β-hexachlorocyclohexane (β-HCH) are chlorinated pesticides that coexist as persistent organic pollutants in the groundwater of several countries in the Caribbean, being an environmental issue. This work evaluates theoretically the competitive formation of host-guest complexes pesticides@cyclodextrines (CDs) as an alternative for water purification and selective separation of pesticides. METHODS: Quantum mechanical calculations based on density functional theory (DFT) and classical molecular dynamics (MD) simulations were used to achieve information on geometries, energies, structure, and dynamics of guest-host complexes in the gas phase, implicit solvent medium, and in aqueous solutions. RESULTS: DFT studies showed that interactions of both pesticides with CDs are mediated by steric factors and guided by maximization of the hydrophobic interactions either with the other pesticide or with the CD cavity's inner atoms. MD results corroborate the formation of stable complexes of both pesticides with the studied CDs. α-CD exhibited a preference for the smaller β-HCH molecule over the CLD that could not perturb the formed complex. CONCLUSIONS: The simulation of competitive formation with γ-CD illustrated that this molecule could accommodate both pesticides inside its cavity. These results suggest that CDs with smaller cavity sizes such as α-CD could be used for selective separation of β-HCH from CLD in water bodies, while γ-CD could be used for methods that aim to remove both pesticides at the same time.
- Publikační typ
- časopisecké články MeSH
Carbohydrate-protein complexes are often characterized by interactions via aromatic amino acid residues. Several mechanisms have been proposed to explain these stacking-like interactions between pyranose sugars and aromatic moieties. The physical basis of these interactions is being explained as either dispersion CH/π or hydrophobic. In order to elucidate the nature of these interactions, we performed a series of molecular dynamics simulation of hevein domain (HEV32) in complex with (β-D-GlcNAc)3. Selected OH/O and CH/π hydrogen bonds involved in carbohydrate recognition were artificially weakened in 100 ns molecular dynamics simulations. Separate weakening of either OH/O or CH/π hydrogen bonds was not sufficient to destabilize the complex. This indicates that other effects, not solely CH/π dispersion interactions, contribute significantly to the stability of the complex. Significant destabilization of complexes was reached only by simultaneous weakening of OH/O and CH/π hydrogen bonds. This also shows that classical hydrogen bonds and CH/π interactions are working in concert to stabilize this carbohydrate-protein test case.
- MeSH
- kationické antimikrobiální peptidy chemie metabolismus MeSH
- konformace sacharidů MeSH
- molekulární modely MeSH
- rostlinné lektiny chemie metabolismus MeSH
- sacharidové sekvence MeSH
- simulace molekulární dynamiky MeSH
- trisacharidy chemie MeSH
- vazebná místa MeSH
- vodík chemie MeSH
- vodíková vazba MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
A simple molecular modeling method for the characterization of polymeric drug carriers is presented. Six biodegradable polymers have been investigated as drug carriers using molecular simulations: l-polylactide, d-polylactide, chitosan, polyglycolic acid, polyethylene glycol and cellulose. Cyclosporine A has been chosen as a model drug substance. Classical molecular dynamics and docking calculations were employed to model and predict polymer-drug interactions. These interactions have been analyzed by non-bond interaction energy and interaction parameter calculated using Flory-Huggins theory. Flexibility of polymer chains has been characterized by the change of gyration radius along the molecular dynamics trajectory. The relationship between mixing energy, chain length and chain flexibility has been revealed for each polymer/drug system.
- MeSH
- cyklosporin chemie MeSH
- molekulární modely * MeSH
- nosiče léků chemie MeSH
- polymery chemie MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
