MD simulation Dotaz Zobrazit nápovědu
- MeSH
- financování organizované MeSH
- Publikační typ
- abstrakty MeSH
BACKGROUND: Guanine quadruplexes (GQs) play vital roles in many cellular processes and are of much interest as drug targets. In contrast to the availability of many structural studies, there is still limited knowledge on GQ folding. SCOPE OF REVIEW: We review recent molecular dynamics (MD) simulation studies of the folding of GQs, with an emphasis paid to the human telomeric DNA GQ. We explain the basic principles and limitations of all types of MD methods used to study unfolding and folding in a way accessible to non-specialists. We discuss the potential role of G-hairpin, G-triplex and alternative GQ intermediates in the folding process. We argue that, in general, folding of GQs is fundamentally different from funneled folding of small fast-folding proteins, and can be best described by a kinetic partitioning (KP) mechanism. KP is a competition between at least two (but often many) well-separated and structurally different conformational ensembles. MAJOR CONCLUSIONS: The KP mechanism is the only plausible way to explain experiments reporting long time-scales of GQ folding and the existence of long-lived sub-states. A significant part of the natural partitioning of the free energy landscape of GQs comes from the ability of the GQ-forming sequences to populate a large number of syn-anti patterns in their G-tracts. The extreme complexity of the KP of GQs typically prevents an appropriate description of the folding landscape using just a few order parameters or collective variables. GENERAL SIGNIFICANCE: We reconcile available computational and experimental studies of GQ folding and formulate basic principles characterizing GQ folding landscapes. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.
- MeSH
- denaturace nukleových kyselin MeSH
- DNA chemie MeSH
- G-kvadruplexy * MeSH
- guanin chemie MeSH
- kinetika MeSH
- lidé MeSH
- párování bází MeSH
- simulace molekulární dynamiky * MeSH
- telomery chemie MeSH
- vztahy mezi strukturou a aktivitou MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- přehledy MeSH
Knowledge of the structure and conformational flexibility of carbohydrates in an aqueous solvent is important to improving our understanding of how carbohydrates function in biological systems. In this study, we extend a variant of the Hamiltonian replica-exchange molecular dynamics (MD) simulation to improve the conformational sampling of saccharides in an explicit solvent. During the simulations, a biasing potential along the glycosidic-dihedral linkage between the saccharide monomer units in an oligomer is applied at various levels along the replica runs to enable effective transitions between various conformations. One reference replica runs under the control of the original force field. The method was tested on disaccharide structures and further validated on biologically relevant blood group B, Lewis X and Lewis A trisaccharides. The biasing potential-based replica-exchange molecular dynamics (BP-REMD) method provided a significantly improved sampling of relevant conformational states compared with standard continuous MD simulations, with modest computational costs. Thus, the proposed BP-REMD approach adds a new dimension to existing carbohydrate conformational sampling approaches by enhancing conformational sampling in the presence of solvent molecules explicitly at relatively low computational cost.
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
- chemie farmaceutická 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
RNA recognition motif (RRM) proteins represent an abundant class of proteins playing key roles in RNA biology. We present a joint atomistic molecular dynamics (MD) and experimental study of two RRM-containing proteins bound with their single-stranded target RNAs, namely the Fox-1 and SRSF1 complexes. The simulations are used in conjunction with NMR spectroscopy to interpret and expand the available structural data. We accumulate more than 50 μs of simulations and show that the MD method is robust enough to reliably describe the structural dynamics of the RRM-RNA complexes. The simulations predict unanticipated specific participation of Arg142 at the protein-RNA interface of the SRFS1 complex, which is subsequently confirmed by NMR and ITC measurements. Several segments of the protein-RNA interface may involve competition between dynamical local substates rather than firmly formed interactions, which is indirectly consistent with the primary NMR data. We demonstrate that the simulations can be used to interpret the NMR atomistic models and can provide qualified predictions. Finally, we propose a protocol for 'MD-adapted structure ensemble' as a way to integrate the simulation predictions and expand upon the deposited NMR structures. Unbiased μs-scale atomistic MD could become a technique routinely complementing the NMR measurements of protein-RNA complexes.
- MeSH
- konformace proteinů MeSH
- lidé MeSH
- magnetická rezonanční spektroskopie MeSH
- molekulární modely MeSH
- motiv rozpoznávající RNA genetika MeSH
- multiproteinové komplexy chemie genetika MeSH
- RNA chemie genetika MeSH
- sekvence aminokyselin genetika MeSH
- serin-arginin sestřihové faktory chemie genetika MeSH
- sestřihové faktory chemie genetika MeSH
- simulace molekulární dynamiky MeSH
- vazebná místa MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
This work extends the multi-scale computational scheme for the quantum mechanics (QM) calculations of Nuclear Magnetic Resonance (NMR) chemical shifts (CSs) in proteins that lack a well-defined 3D structure. The scheme couples the sampling of an intrinsically disordered protein (IDP) by classical molecular dynamics (MD) with protein fragmentation using the adjustable density matrix assembler (ADMA) and density functional theory (DFT) calculations. In contrast to our early investigation on IDPs (Pavlíková Přecechtělová et al., J. Chem. Theory Comput., 2019, 15, 5642-5658) and the state-of-the art NMR calculations for structured proteins, a partial re-optimization was implemented on the raw MD geometries in vibrational normal mode coordinates to enhance the accuracy of the MD/ADMA/DFT computational scheme. In addition, machine-learning based cluster analysis was performed on the scheme to explore its potential in producing protein structure ensembles (CLUSTER ensembles) that yield accurate CSs at a reduced computational cost. The performance of the cluster-based calculations is validated against results obtained with conventional structural ensembles consisting of MD snapshots extracted from the MD trajectory at regular time intervals (REGULAR ensembles). CS calculations performed with the refined MD/ADMA/DFT framework employed the 6-311++G(d,p) basis set that outperformed IGLO-III calculations with the same density functional approximation (B3LYP) and both explicit and implicit solvation. The partial geometry optimization did not universally improve the agreement of computed CSs with the experiment but substantially decreased errors associated with the ensemble averaging. A CLUSTER ensemble with 50 structures yielded ensemble averages close to those obtained with a REGULAR ensemble consisting of 500 MD frames. The cluster based calculations thus required only a fraction of the computational time.
This review focuses on the influence of oxidized phosphatidylcholines (oxPCs) on the biophysical properties of model membranes and is limited to fluorescence, EPR, and MD studies. OxPCs are divided into two classes: A) hydroxy- or hydroperoxy-dieonyl phospatidylcholines, B) phospatidylcholines with oxidized and truncated chains with either aldehyde or carboxylic group. It was shown that the presence of the investigated oxPCs in phospholipid model membranes may have the following consequences: 1) decrease of the lipid order, 2) lowering of phase transition temperatures, 3) lateral expansion and thinning of the bilayer, 4) alterations of bilayer hydration profiles, 5) increased lipid mobility, 6) augmented flip-flop, 7) influence on the lateral phase organisation, and 8) promotion of water defects and, under extreme conditions (i.e. high concentrations of class B oxPCs), disintegration of the bilayer. The effects of class A oxPCs appear to be more moderate than those observed or predicted for class B. Many of the abovementioned findings are related to the ability of the oxidized chains of certain oxPCs to reorient toward the water phase. Some of the effects appear to be moderated by the presence of cholesterol. Although those biophysical alternations are found at oxPC concentrations higher than the total oxPC concentrations found under physiological conditions, certain organelles may reach such elevated oxPC concentrations locally. It is a challenge for the future to correlate the biophysics of oxidized phospholipids to metabolic studies in order to define the significance of the findings presented herein for pathophysiology. This article is part of a Special Issue entitled: Oxidized phospholipids-their properties and interactions with proteins.
- MeSH
- biofyzikální jevy * MeSH
- elektronová paramagnetická rezonance MeSH
- fluorescence MeSH
- fosfolipidy chemie metabolismus MeSH
- lidé MeSH
- lipidové dvojvrstvy chemie metabolismus MeSH
- oxidace-redukce MeSH
- simulace molekulární dynamiky * MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- přehledy MeSH
Molecular motions of free and pheromone-bound mouse major urinary protein I, previously investigated by NMR relaxation, were simulated in 30 ns molecular dynamics (MD) runs. The backbone flexibility was described in terms of order parameters and correlation times, commonly used in the NMR relaxation analysis. Special attention was paid to the effect of conformational changes on the nanosecond time scale. Time-dependent order parameters were determined in order to separate motions occurring on different time scales. As an alternative approach, slow conformational changes were identified from the backbone torsion angle variances, and "conformationally filtered" order parameters were calculated for well-defined conformation states. A comparison of the data obtained for the free and pheromone-bound protein showed that some residues are more rigid in the bound form, but a larger portion of the protein becomes more flexible upon the pheromone binding. This finding is in general agreement with the NMR results. The higher flexibility observed on the fast (fs-ps) time scale was typically observed for the residues exhibiting higher conformational freedom on the ns time scale. An inspection of the hydrogen bond network provided a structural explanation for the flexibility differences between the free and pheromone-bound proteins in the simulations.
Hybrid QM/MM methods combine the rigor of quantum mechanical (QM) calculations with the low computational cost of empirical molecular mechanical (MM) treatment allowing to capture dynamic properties to probe critical atomistic details of enzyme reactions. Catalysis by RNA enzymes (ribozymes) has only recently begun to be addressed with QM/MM approaches and is thus still a field under development. This review surveys methodology as well as recent advances in QM/MM applications to RNA mechanisms, including those of the HDV, hairpin, and hammerhead ribozymes, as well as the ribosome. We compare and correlate QM/MM results with those from QM and/or molecular dynamics (MD) simulations, and discuss scope and limitations with a critical eye on current shortcomings in available methodologies and computer resources. We thus hope to foster mutual appreciation and facilitate collaboration between experimentalists and theorists to jointly advance our understanding of RNA catalysis at an atomistic level.
- MeSH
- biofyzika metody MeSH
- fosfáty chemie MeSH
- fosforylace MeSH
- hořčík chemie MeSH
- katalýza MeSH
- konformace nukleové kyseliny MeSH
- kvantová teorie MeSH
- lidé MeSH
- molekulární modely MeSH
- počítačová simulace MeSH
- ribozomy chemie MeSH
- RNA katalytická chemie MeSH
- RNA virová chemie MeSH
- RNA chemie MeSH
- software MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- práce podpořená grantem MeSH
- přehledy MeSH
- Research Support, N.I.H., Extramural MeSH
BACKGROUND: In the last couple of years, viral infections have been leading the globe, considered one of the most widespread and extremely damaging health problems and one of the leading causes of mortality in the modern period. Although several viral infections are discovered, such as SARS CoV-2, Langya Henipavirus, there have only been a limited number of discoveries of possible antiviral drug, and vaccine that have even received authorization for the protection of human health. Recently, another virial infection is infecting worldwide (Monkeypox, and Smallpox), which concerns pharmacists, biochemists, doctors, and healthcare providers about another epidemic. Also, currently no specific treatment is available against Monkeypox. This research gap encouraged us to develop a new molecule to fight against monkeypox and smallpox disease. So, firstly, fifty different curcumin derivatives were collected from natural sources, which are available in the PubChem database, to determine antiviral capabilities against Monkeypox and Smallpox. MATERIAL AND METHOD: Preliminarily, the molecular docking experiment of fifty different curcumin derivatives were conducted, and the majority of the substances produced the expected binding affinities. Then, twelve curcumin derivatives were picked up for further analysis based on the maximum docking score. After that, the density functional theory (DFT) was used to determine chemical characterizations such as the highest occupied molecular orbital (HOMO), lowest unoccupied molecular orbital (LUMO), softness, and hardness, etc. RESULTS: The mentioned derivatives demonstrated docking scores greater than 6.80 kcal/mol, and the most significant binding affinity was at -8.90 kcal/mol, even though 12 molecules had higher binding scores (-8.00 kcal/mol to -8.9 kcal/mol), and better than the standard medications. The molecular dynamic simulation is described by root mean square deviation (RMSD) and root-mean-square fluctuation (RMSF), demonstrating that all the compounds might be stable in the physiological system. CONCLUSION: In conclusion, each derivative of curcumin has outstanding absorption, distribution, metabolism, excretion, and toxicity (ADMET) characteristics. Hence, we recommended the aforementioned curcumin derivatives as potential antiviral agents for the treatment of Monkeypox and Smallpox virus, and more in vivo investigations are warranted to substantiate our findings.
- MeSH
- antivirové látky farmakologie MeSH
- COVID-19 * MeSH
- kurkumin * farmakologie MeSH
- lidé MeSH
- objevování léků MeSH
- opičí neštovice * MeSH
- pravé neštovice * farmakoterapie MeSH
- racionální návrh léčiv MeSH
- simulace molekulární dynamiky MeSH
- simulace molekulového dockingu MeSH
- virus varioly * MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
