scholarly journals Effects of the Temperature and the pH on the Main Protease of SARS-CoV-2: A Molecular Dynamics Simulation Study

2021 ◽  
Vol 12 (6) ◽  
pp. 7239-7248

The novel coronavirus, recognized as COVID-19, is the cause of an infection outbreak in December 2019. The effect of temperature and pH changes on the main protease of SARS-CoV-2 were investigated using all-atom molecular dynamics simulation. The obtained results from the root mean square deviation (RMSD) and root mean square fluctuations (RMSF) analyses showed that at a constant temperature of 25℃ and pH=5, the conformational change of the main protease is more significant than that of pH=6 and 7. Also, by increasing temperature from 25℃ to 55℃ at constant pH=7, a remarkable change in protein structure was observed. The radial probability of water molecules around the main protease was decreased by increasing temperature and decreasing pH. The weakening of the binding energy between the main protease and water molecules due to the increasing temperature and decreasing pH has reduced the number of hydrogen bonds between the main protease and water molecules. Finding conditions that alter the conformation of the main protease could be fundamental because this change could affect the virus’s functionality and its ability to impose illness.

Molecules ◽  
2021 ◽  
Vol 26 (6) ◽  
pp. 1767
Author(s):  
Mohamed E. Abouelela ◽  
Hamdy K. Assaf ◽  
Reda A. Abdelhamid ◽  
Ehab S. Elkhyat ◽  
Ahmed M. Sayed ◽  
...  

Severe acute respiratory syndrome coronavirus (SARS-CoV-2) disease is a global rapidly spreading virus showing very high rates of complications and mortality. Till now, there is no effective specific treatment for the disease. Aloe is a rich source of isolated phytoconstituents that have an enormous range of biological activities. Since there are no available experimental techniques to examine these compounds for antiviral activity against SARS-CoV-2, we employed an in silico approach involving molecular docking, dynamics simulation, and binding free energy calculation using SARS-CoV-2 essential proteins as main protease and spike protein to identify lead compounds from Aloe that may help in novel drug discovery. Results retrieved from docking and molecular dynamics simulation suggested a number of promising inhibitors from Aloe. Root mean square deviation (RMSD) and root mean square fluctuation (RMSF) calculations indicated that compounds 132, 134, and 159 were the best scoring compounds against main protease, while compounds 115, 120, and 131 were the best scoring ones against spike glycoprotein. Compounds 120 and 131 were able to achieve significant stability and binding free energies during molecular dynamics simulation. In addition, the highest scoring compounds were investigated for their pharmacokinetic properties and drug-likeness. The Aloe compounds are promising active phytoconstituents for drug development for SARS-CoV-2.


2013 ◽  
Vol 2013 ◽  
pp. 1-12 ◽  
Author(s):  
Zhiguo Chen ◽  
Yi Fu ◽  
Wenbo Xu ◽  
Ming Li

Bacillus amyloliquefaciensribonuclease Barnase (RNase Ba) is a 12 kD (kilodalton) small extracellular ribonuclease. It has broad application prospects in agriculture, clinical medicine, pharmaceutical, and so forth. In this work, the thermal stability of Barnase has been studied using molecular dynamics simulation at different temperatures. The present study focuses on the contribution of noncovalent intramolecular interaction to protein stability and how they affect the thermal stability of the enzyme. Profiles of root mean square deviation and root mean square fluctuation identify thermostable and thermosensitive regions of Barnase. Analyses of trajectories in terms of secondary structure content, intramolecular hydrogen bonds and salt bridge interactions indicate distinct differences in different temperature simulations. In the simulations, Four three-member salt bridge networks (Asp8-Arg110-Asp12, Arg83-Asp75-Arg87, Lys66-Asp93-Arg69, and Asp54-Lys27-Glu73) have been identified as critical salt bridges for thermostability which are maintained stably at higher temperature enhancing stability of three hydrophobic cores. The study may help enlighten our knowledge of protein structural properties, noncovalent interactions which can stabilize secondary peptide structures or promote folding, and also help understand their actions better. Such an understanding is required for designing efficient enzymes with characteristics for particular applications at desired working temperatures.


2016 ◽  
Author(s):  
Yuan-Ping Pang

ABSTRACTPredicting crystallographic B-factors of a protein from a conventional molecular dynamics simulation is challenging in part because the B-factors calculated through sampling the atomic positional fluctuations in a picosecond molecular dynamics simulation are unreliable and the sampling of a longer simulation yields overly large root mean square deviations between calculated and experimental B-factors. This article reports improved B-factor prediction achieved by sampling the atomic positional fluctuations in multiple picosecond molecular dynamics simulations that use uniformly increased atomic masses by 100-fold to increase time resolution. Using the third immunoglobulin-binding domain of protein G, bovine pancreatic trypsin inhibitor, ubiquitin, and lysozyme as model systems, the B-factor root mean square deviations (mean ± standard error) of these proteins were 3.1 ± 0.2–9 ± 1 Å2for Cα and 7.3 ± 0.9–9.6 ± 0.2 Å2for Cγ, when the sampling was done, for each of these proteins, over 20 distinct, independent, and 50-picosecond high-mass molecular dynamics simulations using AMBER forcefield FF12MC or FF14SB. These results suggest that sampling the atomic positional fluctuations in multiple picosecond high-mass molecular dynamics simulations may be conducive toa prioriprediction of crystallographic B-factors of a folded globular protein.


2019 ◽  
Vol 22 (8) ◽  
pp. 546-554
Author(s):  
Mohadeseh Bazoobandi ◽  
Mohammad R. Bozorgmehr ◽  
Ali Mahmoudi ◽  
Ali Morsali

Aim and Objective: The stability of the G-quadruplex structure can increase its activity in telomerase inhibiting cancer cells. In this study, a molecular dynamics simulation method was used to study the effect of three phenanthroline-based ligands on the structure of G-quadruplex at the temperatures of 20, 40, 60 and 80°C. Materials and Methods: RMSD values and frequency of calculated RMSD in the presence and absence of ligands show that ligands cause the relative stability of the G-quadruplex, particularly at low temperatures. The calculation of hydrogen bonds in Guanine-tetrads in three different quadruplex sheets shows that the effect of ligands on the sheets is not the same so that the bottom sheet of G-quadruplex is most affected by the ligands at high temperatures, and the Guaninetetrads in this sheet are far away. Conformation factor was calculated as a measure of ligands binding affinity for each of the G-quadruplex residues. Results: The results show that the studied ligands interact more with the G-quadruplex than loop areas, although with increasing temperature, the binding area also includes the G-quadruplex sheets. The contribution of each of the residues involved in the G-quadruplex binding area with ligands was also calculated. Conclusion: The calculations performed are consistent with the previous experimental observations that can help to understand the molecular mechanism of the interaction of phenanthroline and its derivatives with quadruplex.


Molecules ◽  
2021 ◽  
Vol 26 (21) ◽  
pp. 6538
Author(s):  
Md. Mominur Rahman ◽  
Md. Junaid ◽  
S. M. Zahid Hosen ◽  
Mohammad Mostafa ◽  
Lei Liu ◽  
...  

Inflammation plays an important role in different chronic diseases. Brominated indoles derived from the Australian marine mollusk Dicathais orbita (D. orbita) are of interest for their anti-inflammatory properties. This study evaluates the binding mechanism and potentiality of several brominated indoles (tyrindoxyl sulfate, tyrindoleninone, 6-bromoisatin, and 6,6′-dibromoindirubin) against inflammatory mediators cyclooxygenases-1/2 (COX-1/2) using molecular docking, followed by molecular dynamics simulation, along with physicochemical, drug-likeness, pharmacokinetic (pk), and toxicokinetic (tk) properties. Molecular docking identified that these indole compounds are anchored, with the main amino acid residues, positioned in the binding pocket of the COX-1/2, required for selective inhibition. Moreover, the molecular dynamics simulation based on root mean square deviation (RMSD), radius of gyration (Rg), solvent accessible surface area (SASA), and root mean square fluctuation (RMSF) analyses showed that these natural brominated molecules transit rapidly to a progressive constant configuration during binding with COX-1/2 and seem to accomplish a consistent dynamic behavior by maintaining conformational stability and compactness. The results were comparable to the Food and Drug Administration (FDA)-approved selective COX inhibitor, aspirin. Furthermore, the free energy of binding for the compounds assessed by molecular mechanics–Poisson–Boltzmann surface area (MM–PBSA) confirmed the binding capacity of indoles towards COX-1/2, with suitable binding energy values except for the polar precursor tyrindoxyl sulfate (with COX-1). The physicochemical and drug-likeness analysis showed zero violations of Lipinski’s rule, and the compounds are predicted to have excellent pharmacokinetic profiles. These indoles are projected to be non-mutagenic and free from hepatotoxicity, with no inhibition of human ether-a-go–go gene (hERG) I inhibitors, and the oral acute toxicity LD50 in rats is predicted to be similar or lower than aspirin. Overall, this work has identified a plausible mechanism for selective COX inhibition by natural marine indoles as potential therapeutic candidates for the mitigation of inflammation.


NANO ◽  
2011 ◽  
Vol 06 (01) ◽  
pp. 19-29 ◽  
Author(s):  
SERGEY SHITYAKOV ◽  
THOMAS DANDEKAR

Using classical all-atom molecular dynamics simulation, we investigated the molecular dynamics of palmitoyloleoylphosphatidylcholine and palmitoyloleoylphosphatidylethanolamine membrane bilayers enforced by a single-wall carbon nanotube. We postulated that an insertion of a single-wall carbon nanotube in the center of lipid membrane "strengthens" ambient lipids and prevents the whole system from further destabilization by high temperatures. We implemented root mean square deviation and root mean square fluctuation analyses of simulated structures from their initial states in order to emphasize the molecular dynamics behavior of these structures during 1000 ps simulation time at different temperatures. The data suggest that an intercalated carbon nanotube restrains the conformational freedom of adjacent lipids and hence has an impact on the membrane stabilization dynamics. On the other hand, different lipid membranes may have dissimilarities due to the differing abilities to create a bridge formation between the adherent lipid molecules. The results derived from this work may be of importance in developing stable nanosystems for construction of novel biomaterials and delivery of various biomolecules in the fields of biosensors, biomaterials, and biophysics.


Sign in / Sign up

Export Citation Format

Share Document