Добавил:
Sekretar
kiopkiopkiop18@yandex.ru
t.me/Prokururor I Вовсе не секретарь, но почту проверяю
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Ординатура / Хирургия / Библиотека им академика М.И. Перельмана / Книга_5435_Библиотеки_им_академика_М_И_Перельмана
.pdf
8.1.8.4 Nucleic acids
MD simulations have been instrumental in understanding the dynamics and struc-
tural properties of nucleic acids such as DNA and RNA. Simulations have elucidated
DNA replication, transcription, and repair mechanisms. They have provided insights
into RNA folding, ribosome function, and RNA–protein interactions. MD simulations
have also contributed to studying nucleic acid-based drug delivery systems and gene
editing techniques.
8.1.8.5 Molecular recognition and binding
MD simulations have shed light on the mechanisms of molecular recognition and
binding events in biological systems. They have elucidated the principles of ligand
binding to r eceptors, protein–ligand recognition, and the dynamics of binding pock-
ets. Simulations have provided insights into conformational changes upon binding
and the role of water molecules in binding events.
8.1.8.6 Drug resistance and mechanisms of action
MD simulations have helped in understanding drug resistance mechanisms and the
structural basis of drug action. Simulations have revealed the structural changes in
drug targets that confer resistance. They have also elucidated the effects of mutations
on drug binding and the mechanisms by which drugs inhibit target proteins.
8.1.8.7 The role of molecular dynamics simulations in antibody designing
Antibodies have emerged as a highly desirable class of biotherapeutics with less toxic-
ity and their innate ability to selectively recognize and bind to a diverse range of tar-
gets [22]. In the antibody designing process, researchers can analyze the molecular
interactions between the antibody and antigen and identify the regions that contrib-
ute most significantly to the binding affinity by employing MD simulations and other
computational methods [23].
In a recent study, Bekker et al. [24] utilized MD simulations to forecast the tem-
perature stability of single-domain antibodies and compared their predictions to ex-
perimental results. As a result, they were able to identify the crucial residues causing
the instability and effectively create a mutant antibody with increased heat stability
[24, 25].
164 Disha Tewari et al.
https://t.me/med1917
8.1.8.8 To evaluate the mobility or flexibility of biomolecules
X-ray crystallography and cryo-EM are experimental methods often used to determine
the average structure of a biomolecule. However, it is possible to quantitatively ana-
lyze the degree of movement and structural fluctuations in various sites of a molecule
under equilibrium conditions by simulating this structure. Additionally, simulations
can shed light on the dynamic behavior of salt ions and water molecules, which are
necessary for the functioning of protein and ligand binding [26–28].
8.1.8.9 Simulation in drug discovery process
In the drug development process, target validation is followed by lead compound iden-
tification and optimization of the identified compounds as a potential drug candidate.
Common methods for estimating the best ligand binding sites and binding energies are
docking and scoring. However, the scoring features of existing docking systems have
numerous flaws that lead to incorrect estimates of binding affinities [29].
With MD simulati ons, we can overcome these limitations and provides a more
accurate assessment of compound binding affinity as the simulations take into ac-
count the impact of water and the dynamics of the binding partners [30].
Overall, MD simulations have had a transformative impact on our understanding
of biological systems. They have provided atomistic-level details, dynamic insights,
and predictive capabilities that complement experimental techniques. MD simulations
have accelerated drug discovery efforts, facilitated the design of novel therapeutics,
and advanced our knowledge of biological processes, contributing to numerous scien-
tific breakthroughs in the field of biology and medicine.
8.1.9 Exploring molecular dynamics simulation with GROMACS
There are multiple software available for running MD simulations. One of the most
well-known programs, open-source GROMACS, is the main focus of this chapter.
For classical MD simulations, there is a software package called Gromacs. It is one
of the most widely used and the fastest MD simulation programs in the world, along
with NAMD, Amber, and LAMMPS. The term GROMACS stands for Groningen Machine
for Chemical Simulations [11]. It is used to carry out MD simulations of both nonbio-
logical and biological systems including polymers, proteins, nucleic acids, lipids and
other macromolecules. GROMACS calculates the nonbonding interactions very quickly
and effectively. Some of the key features of Gromacs are as follows:
8 Design of target hit molecules using molecular dynamic simulations 165
https://t.me/med1917
– It is useful to research how systems (such as protein–ligand complexes) behave
under different pressures, temperatures, and pH levels.
– GROMACS utilizes various types of force fields that can be used for energy mini-
mization of different systems.
– It facilitates the analysis of the dynamic behavior of atoms and molecules by en-
abling the monitoring of their physical movements at regular time intervals.
– Using GROMACS, it is simple to study the simulation of bonding interactions be-
tween two protein structures.
Currently, there are seven tutorials available for this program:
1. Lysozyme in water: This tutorial aims to provide new users with a fundamental
introduction to the tools needed to set up, run, and evaluate basic analysis on a
“typical” system using GROMACS.
2. KALP15 in DPPC: This specific tutorial is intended for experienced users who have
used MD simulations in the past and are interested in investigating membrane pro-
teins and understanding more about force field construction and modification.
3. Umbrella sampling: It helps to estimate the potential of mean force along a sin-
gle, linear degree of freedom. This technique is used to calculate free energy pro-
files along a specific reaction coordinate. It involves restraining the system at
different positions along the reaction coordinate and running multiple MD simu-
lations under these restraints. The resulting data is then used to reconstruct the
free energy landscape using methods such as weighted histogram analysis.
4. Biphasic systems: It helps in the formation of a water-cyclohexane biphasic
system.
5. Protein–ligand complex: In this program, users are guided through the process
of working with a protein–ligand system, emphasizing the importance of accurate
ligand parameterization and topology management.
6. Free energy of solvation: The tutorial explains the steps involved in performing
a basic free energy calculation, specifically targeting the removal of van der
Waals interactions between a straightforward molecule (methane) and water.
7. Virtual sites: With the help of this tutorial, users can manually create virtual
sites for the relatively simple CO
2
molecule, which is a linear, triatomic molecule.
GROMACS is command-line-driven software, and its usage involves writing input files,
executing simulations using terminal commands, and analyzing results. The GRO-
MACS website (www.gromacs.org) offers comprehensive documentation, tutorials,
and user forums to assist with learning and troubleshooting. Additionally, there are
graphical user interfaces (GUIs) ava ilable, such as GROMACS-GUI and VMD, which
provide a more intuitive and visual way to set up simulations and analyze results, if
you prefer a graphical workflow.
Remember, performing MD simulations requires computational resources, especially
for larger systems or longer simulation times. HPC clusters or specialized computing re-
166 Disha Tewari et al.
https://t.me/med1917
sources are often employed to carry out MD simulations efficiently. Overall, exploring
MD simulation with GROMACS can be a rewarding experience, allowing you to gain in-
sights into the dynamic behavior of molecular systems and study phenomena that are
otherwise challenging to access.
8.2 The general steps involved in designing hit
molecules using molecular dynamics (MD)
simulations
1. Target and ligand selection: The initial step is to choose a suitable target protein
for the MD simulations. This may involve selecting a protein target that is known
to be involved in a particular disease or condition.
Once the target protein has been selected, the next step is to choose a ligand
that will work well in the simulation. The ligand that we choose may be a known
inhibitor or activator of the target protein or a compound with a similar chemical
structure.
2. Preparation of the system: The next step is preparation of the system for the MD
simulations. To do this, a 3D model of the target protein and ligand must be con-
structed, and the system can then be solvated in a suitable solvent (usually water).
Gromacs supplies a number of tools for system setup, including pdb2gmx for cre-
ating topologies and solvates for introducing water molecules to the system.
3. System equilibratio n: Before beginning the production of MD simulations, the
system must be stabilized by equilibrating it. This involves minimizing the energy
of the system, gradually increasing the temperature, and allowing the solvent and
ions to equilibrate around the protein–ligand complex. Gromacs provides various
equilibration protocols such as the energy minimization, NVT and NPT equilibra-
tion, and position-restrained MD simulations.
4. MD simulations production: Once the system has reached equilibrium, produc-
tion MD simulations can be run. This involves running the simulation for an ex-
tended period, typically on the order of tens to hundreds of nanoseconds. During
the simulation, the ligand may be allowed to move freely or be restrained in a
particular position.
5. Data analysis: After the production of MD simulations, the resulting trajectories
are analyzed to identify potential hit molecules. This may involve calculating dif-
ferent properties, such as binding affinity, binding free energy, and binding
modes, and comparing the results to experimental data or other computational
methods. Gromacs provides various analysis tools such as gmx energy for calcu-
lating energy terms, gmx rms for calculating root mean square deviation, and
gmx trjconv for extracting snapshots of the simulation trajectory.
8 Design of target hit molecules using molecular dynamic simulations 167
https://t.me/med1917
Thus, Gromacs provides several key aspects, such as efficient simulation performance,
support for enhanced sampling methods, and a range of analysis tools, which make it
a popular choice for researchers in the field of drug discovery and design (Figure 8.1).
8.3 Gromacs protocol applied for molecular dynamic
simulation targeting as an example of COVID-19
MD simulations can be applied to various fields. In this chapter, we target COVID-19
as an example to explain the results of GROMACS. Here is an example protocol for
using GROMACS in a MD simulation targeting COVID-19.
The COVID-19 pandemic caused by SARS-CoV-2 has resulted in a global health crisis,
highlighting the urgent need for effective antiviral treatments. In this study, we investi-
gate the potential of phytochemical compounds derived from Kigelia Africana, a medici-
nal plant indigenous to Africa, as candidates for targeting COVID-19. We employ MD
simulations using the GROMACS software package to explore the interactions between
the identified phytochemicals and the key viral proteins involved in the infection cycle
of SARS-CoV-2. First, we collect and analyze the phytochemical composition of Kigelia
Africana extracts through experimental techniques and literature surveys. We identify
several bioactive compounds present in the plant, such as flavonoids, saponins, and ter-
penoids, which have demonstrated antiviral properties in previous studies. Next, using
Figure 8.1: Systematic representation of GROMACS used to perform molecular dynamics simulation.
168 Disha Tewari et al.
https://t.me/med1917
the three-dimensional structures of the relevant viral proteins, including the main pro-
tease (Mpro) obtained from the Protein Data Bank, we prepare the protein systems for
MD simulations. We parameterize the Kigelia Africana phytochemicals using appropri-
ate force fields and perform system optimization and equilibration steps. Subsequently,
we conduct extensive MD simulations to examine the binding modes, stability, and dy-
namics of the phytochemicals when interacting with the target viral proteins. We ana-
lyze the trajectories to determine the key molecular interactions, including hydrogen
bonding, hydrophobic interactions, and electrostatic interactions, which contribute to
the stability of the protein–ligand complexes.
Table 8.2 provides an overview of some commonly used GROMACS commands and
their respective functionalities. These commands cover tasks such as system setup, simu-
lation execution, trajectory analysis, energy calculation, and structural analysis. Remem-
bertorefertotheGROMACSdocumentationorspecificcommandhelpfordetailedusage
instructions and additional options available for each command.
Figure 8.2: Snap shot of molecular dynamic simulation command by using GROMACS.
8 Design of target hit molecules using molecular dynamic simulations 169
https://t.me/med1917
8.3.1 Results
8.3.1.1 Root mean square deviation (RMSD)
The structural stability of the Mpro-ligand complex was studied using the RMSD (root
mean square deviation) during the 100 ns simulation (Table 8.3 and Fig 8.3A). The av-
erage RMSD for protein was 0.18 ± 0.03 nm, whereas Mpro-X 77 and Mpro-Luteolin7-O-
glucoside protein–ligand complex system had RMSD values of 0.17 ± 0.02 nm and 0.18
± 0.01 nm, respectively.
8.3.1.2 Root mean square fluctuation (RMSF)
Root mean square fluctuations (RMSF), a crucial structural measurement, is used in
order to measure the flexibility and rigidity of the protein–drug complexes (Table 8.3
and Fig 8.3B). A lower RMSF number indicates the system’s higher stability, while a
higher RMSF value signifies greater flexibility during the MD simulation. For native pro-
tein, reference molecule and Mpro complexes, the average RMSF values were 0.09 ±
0.02 nm, 0.13 ± 0.04 nm, and 0.10 ± 0.03 nm, respectively.
Table 8.2: Commonly used GROMACS commands and their functionalities.
Command Description
gmx pdbgmx Converts a PDB file to a GROMACS topology file
gmx editconf Modifies the box size and shape of the system
gmx solvate Adds water molecules to the system
gmx grompp Prepares the input files for MD simulation
gmx mdrun Performs the MD simulation
gmx trjconv Converts, manipulates, and analyzes trajectories
gmx energy Calculates various energy components of the system
gmx rms Calculates root-mean-square deviation (RMSD) of atomic positions
gmx rmsf Calculates root-mean-square fluctuation (RMSF) of atomic positions
gmx hbond Analyzes hydrogen bonds in the system
gmx gyrate Calculates the radius of gyration of a group of atoms
gmx mindist Determines the minimum distance between two groups of atoms
gmx density Computes the density profile of the system
gmx potential Calculates the potential energy distribution
gmx trjorder Analyzes molecular order parameters in the system
gmx cluster Performs clustering analysis of structures
gmx free energy Calculates free energy differences using various methods
170 Disha Tewari et al.
https://t.me/med1917
Figure 8.3: GROMACS analysis results: (A) RMSD, (B) RMSF, (C) RG, and (D) HBOND.
8 Design of target hit molecules using molecular dynamic simulations 171
https://t.me/med1917
8.3.1.3 Radius of gyration (R
g
)
The degree of compactness in the protein structure caused by the presence or absence
of ligands is determined by the radius of gyration (R
g
)[31].TheaverageR
g
value
(Table 8.3) for native protein Mpro was found to be around 1.88 ± 0.05. However, R
g
values of the Mpro-Luteolin7-O-glucoside complex and reference molecule are 1.88 ±
0.05 nm and 1.88 ± 0.08 nm, respectively.
8.3.1.4 Hydrogen bond analysis
Hydrogen bonds have a major role in stabilizing the binding of ligands. The higher
the number of hydrogen bonds (Table 8.3 and Fig 8.3D), the greater the binding affin-
ity toward p rotein [32]. The hydrogen bond analysis between Mpro and the ligand
was calculated, which depicts the H-bond interaction between the protein and the li-
gand during the 100 ns MD simulations.
Based on our simulations, we identify specific phytochemical compounds from Kige-
lia Africana that exhibit strong binding affini ty and favorable interactions with the
viral proteins. These compounds potentially interfere with essential viral p rocesses,
such as proteolytic cleavage and viral entry, thereby inhibiting viral r eplication and
spread. Our findi ngs highlight the potential of Kigelia Africana phytochemicals as
promising candidates for the development of novel therapeutics against COVID-19.
Further experimental studies and validation are necessary to confirm the antiviral
activity of these phytochemicals and explore their efficacy in inhibiting SARS-CoV-2
infection. The insights gained from this study contribute to the ongoing efforts in
drug discovery and provide a basis for future investigations into natural compounds
as potential treatments for COVID-19.
Table 8.3: The average values of different parameters, RMSD, RMSF, R
g
, and H-bond.
S.
No.
Phytochemicals Average Rmsd
(nm)
Average Rmsf
(nm)
Average
R
g
(nm)
H-bond
(Å)
. Native protein (Mpro) . ± . . ± . . ± . –
. Mpro-reference complex . ± . . ± . . ± .
. Mpro-Luteolin-O-
glucoside complex
. ± . . ± . . ± .
172 Disha Tewari et al.
https://t.me/med1917
References
[1] McCammon, J. A., Gelin, B. R., & Karplus, M. Dynamics of folded proteins. Nature, 1977, 267,
585–590. https://doi.org/10.1038/267585a0.
[2] Warshel, A., & Levitt, M. Theoretical studies of enzymic reactions: Dielectric, electrostatic and steric
stabilization of the carbonium ion in the reaction of lysozyme. Journal of Molecular Biology, 1976,
103, 227–249. https://doi.org/10.1016/0022-2836(76)90311-9.
[3] Roccatano, D., Barthel, A., & Zacharias, M. Structural flexibility of the nucleosome core particle at
atomic resolution studied by molecular dynamics simulation. Biopolymers, 2007, 85, 407–421.
https://doi.org/10.1002/bip.20690.
[4] Sharma, S., Ding, F., & Dokholyan, N. V. Multiscale modeling of nucleosome dynamics. Biophysical
Journal, 2007, 92, 1457–1470. https://doi.org/10.1529/biophysj.106.094805.
[5] Brandman, R., Brandman, Y., & Pande, V. S. A-site residues move independently from P-site residues
in all-atom molecular dynamics simulations of the 70S bacterial ribosome. PLoS ONE, 2012, 7,
e29377. https://doi.org/10.1371/journal.pone.0029377.
[6] Tinoco, I., & Wen, J. D. Simulation and analysis of single-ribosome translation. Physical Biology,
2009, 6, 025006. https://doi.org/10.1088/1478-3975/6/2/025006.
[7] Karplus, M., & McCammon, J. A. Molecular dynamics simulations of biomolecules. Nature Structural
Biology, 2002, 9, 646–652. https://doi.org/10.1038/nsb0902-646.
[8] Rahman, A. Correlations in the motion of atoms in liquid argon. Physical Review, 1964, 136,
A405–A411. https://doi.org/10.1103/PhysRev.136.A405.
[9] Karâa, W. B. A., & Dey, N. (eds) Biomedical Image Analysis and Mining Techniques for Improved
Health Outcomes. 2016, IGI Global. https://doi.org/10.4018/978-1-4666-8811-7.
[10] Grubmüller, H., Heller, H., Windemuth, A., & Schulten, K. Generalized verlet algorithm for efficient
molecular dynamics simulations with long-range interactions. Molecular Simulation, 1991, 6,
121–142. https://doi.org/10.1080/08927029108022142.
[11] Abraham, M. J., Murtola, T., Schulz, R., Páll, S., Smith, J. C., Hess, B., & Lindahl, E. GROMACS: High
performance molecular simulations through multi-level parallelism from laptops to
supercomputers. SoftwareX, 2015, 1–2, 19–25. https://doi.org/10.1016/j.softx.2015.06.001.
[12] Harder, E., Damm, W., Maple, J., Wu, C., Reboul, M., Xiang, J. Y., Wang, L., Lupyan, D., Dahlgren, M. K.,
Knight, J. L., Kaus, J. W., Cerutti, D. S., Krilov, G., Jorgensen, W. L., Abel, R., & Friesner, R. A. OPLS3: A
force field providing broad coverage of drug-like small molecules and proteins. Journal of Chemical
Theory and Computation, 2016, 12, 281–296. https://doi.org/10.1021/acs.jctc.5b00864.
[13] Huang, J., Rauscher, S., Nawrocki, G., Ran, T., Feig, M., de Groot, B. L., Grubmüller, H., & MacKerell,
A. D. CHARMM36m: An improved force field for folded and intrinsically disordered proteins. Nature
Methods, 2017, 14, 71–73. https://doi.org/10.1038/nmeth.4067.
[14] Robustelli, P., Piana, S., & Shaw, D. E. Developing a molecular dynamics force field for both folded
and disordered protein states. Proceedings of the National Academy of Sciences, 2018, 115.
https://doi.org/10.1073/pnas.1800690115.
[15] Klauda, J. B., Venable, R. M., Freites, J. A., O’Connor, J. W., Tobias, D. J., Mondragon-Ramirez, C.,
Vorobyov, I., MacKerell, A. D., & Pastor, R. W. Update of the CHARMM all-atom additive force field
for lipids: Validation on six lipid types. The Journal of Physical Chemistry B, 2010, 114, 7830–7843.
https://doi.org/10.1021/jp101759q.
[16] Vanommeslaeghe, K., & MacKerell, A. D. Automation of the CHARMM general force field (CGenFF)
I: Bond perception and atom typing. Journal of Chemical Information and Modeling, 2012, 52,
3144–3154. https://doi.org/10.1021/ci300363c.
[17] Cornell, W. D., Cieplak, P., Bayly, C. I., Gould, I. R., Merz, K. M., Ferguson, D. M., Spellmeyer, D. C.,
Fox, T., Caldwell, J. W., & Kollman, P. A. A second generation force field for the simulation of
8 Design of target hit molecules using molecular dynamic simulations 173
https://t.me/med1917
Соседние файлы в папке Библиотека им академика М.И. Перельмана