Assessing the Antimalarial Potentials of Phytochemicals: Virtual Screening, Molecular Dynamics and In-Vitro Investigations

2019 ◽  
Vol 16 (3) ◽  
pp. 291-300
Author(s):  
Saumya K. Patel ◽  
Mohd Athar ◽  
Prakash C. Jha ◽  
Vijay M. Khedkar ◽  
Yogesh Jasrai ◽  
...  

Background: Combined in-silico and in-vitro approaches were adopted to investigate the antiplasmodial activity of Catharanthus roseus and Tylophora indica plant extracts as well as their isolated components (vinblastine, vincristine and tylophorine). </P><P> Methods: We employed molecular docking to prioritize phytochemicals from a library of 26 compounds against Plasmodium falciparum multidrug-resistance protein 1 (PfMDR1). Furthermore, Molecular Dynamics (MD) simulations were performed for a duration of 10 ns to estimate the dynamical structural integrity of ligand-receptor complexes. </P><P> Results: The retrieved bioactive compounds viz. tylophorine, vinblastin and vincristine were found to exhibit significant interacting behaviour; as validated by in-vitro studies on chloroquine sensitive (3D7) as well as chloroquine resistant (RKL9) strain. Moreover, they also displayed stable trajectory (RMSD, RMSF) and molecular properties with consistent interaction profile in molecular dynamics simulations. </P><P> Conclusion: We anticipate that the retrieved phytochemicals can serve as the potential hits and presented findings would be helpful for the designing of malarial therapeutics.

2021 ◽  
Vol 22 (11) ◽  
pp. 5931
Author(s):  
Maria Arnittali ◽  
Anastassia N. Rissanou ◽  
Maria Amprazi ◽  
Michael Kokkinidis ◽  
Vagelis Harmandaris

In the current work we study, via molecular simulations and experiments, the folding and stability of proteins from the tertiary motif of 4-α-helical bundles, a recurrent motif consisting of four amphipathic α-helices packed in a parallel or antiparallel fashion. The focus is on the role of the loop region in the structure and the properties of the wild-type Rop (wtRop) and RM6 proteins, exploring the key factors which can affect them, through all-atom molecular dynamics (MD) simulations and supporting by experimental findings. A detailed investigation of structural and conformational properties of wtRop and its RM6 loopless mutation is presented, which display different physical characteristics even in their native states. Then, the thermal stability of both proteins is explored showing RM6 as more thermostable than wtRop through all studied measures. Deviations from native structures are detected mostly in tails and loop regions and most flexible residues are indicated. Decrease of hydrogen bonds with the increase of temperature is observed, as well as reduction of hydrophobic contacts in both proteins. Experimental data from circular dichroism spectroscopy (CD), are also presented, highlighting the effect of temperature on the structural integrity of wtRop and RM6. The central goal of this study is to explore on the atomic level how a protein mutation can cause major changes in its physical properties, like its structural stability.


2020 ◽  
Author(s):  
Sean A. Newmister ◽  
Kinshuk Raj Srivastava ◽  
Rosa V. Espinoza ◽  
Kersti Caddell Haatveit ◽  
Yogan Khatri ◽  
...  

Biocatalysis offers an expanding and powerful strategy to construct and diversify complex molecules by C-H bond functionalization. Due to their high selectivity, enzymes have become an essential tool for C-H bond functionalization and offer complementary reactivity to small-molecule catalysts. Hemoproteins, particularly cytochromes P450, have proven effective for selective oxidation of unactivated C-H bonds. Previously, we reported the in vitro characterization of an oxidative tailoring cascade in which TamI, a multifunctional P450 functions co-dependently with the TamL flavoprotein to catalyze regio- and stereoselective hydroxylations and epoxidation to yield tirandamycin A and tirandamycin B. TamI follows a defined order including 1) C10 hydroxylation, 2) C11/C12 epoxidation, and 3) C18 hydroxylation. Here we present a structural, biochemical, and computational investigation of TamI to understand the molecular basis of its substrate binding, diverse reactivity, and specific reaction sequence. The crystal structure of TamI in complex with tirandamycin C together with molecular dynamics simulations and targeted mutagenesis suggest that hydrophobic interactions with the polyene chain of its natural substrate are critical for molecular recognition. QM/MM calculations and molecular dynamics simulations of TamI with variant substrates provided detailed information on the molecular basis of sequential reactivity, and pattern of regio- and stereo-selectivity in catalyzing the three-step oxidative cascade.<br>


2020 ◽  
Vol 16 (6) ◽  
pp. 784-795
Author(s):  
Krisnna M.A. Alves ◽  
Fábio José Bonfim Cardoso ◽  
Kathia M. Honorio ◽  
Fábio A. de Molfetta

Background:: Leishmaniosis is a neglected tropical disease and glyceraldehyde 3- phosphate dehydrogenase (GAPDH) is a key enzyme in the design of new drugs to fight this disease. Objective:: The present study aimed to evaluate potential inhibitors of GAPDH enzyme found in Leishmania mexicana (L. mexicana). Methods: A search for novel antileishmanial molecules was carried out based on similarities from the pharmacophoric point of view related to the binding site of the crystallographic enzyme using the ZINCPharmer server. The molecules selected in this screening were subjected to molecular docking and molecular dynamics simulations. Results:: Consensual analysis of the docking energy values was performed, resulting in the selection of ten compounds. These ligand-receptor complexes were visually inspected in order to analyze the main interactions and subjected to toxicophoric evaluation, culminating in the selection of three compounds, which were subsequently submitted to molecular dynamics simulations. The docking results showed that the selected compounds interacted with GAPDH from L. mexicana, especially by hydrogen bonds with Cys166, Arg249, His194, Thr167, and Thr226. From the results obtained from molecular dynamics, it was observed that one of the loop regions, corresponding to the residues 195-222, can be related to the fitting of the substrate at the binding site, assisting in the positioning and the molecular recognition via residues responsible for the catalytic activity. Conclusion:: he use of molecular modeling techniques enabled the identification of promising compounds as inhibitors of the GAPDH enzyme from L. mexicana, and the results obtained here can serve as a starting point to design new and more effective compounds than those currently available.


Molecules ◽  
2021 ◽  
Vol 26 (6) ◽  
pp. 1711
Author(s):  
Mohamed Ahmed Khaireh ◽  
Marie Angot ◽  
Clara Cilindre ◽  
Gérard Liger-Belair ◽  
David A. Bonhommeau

The diffusion of carbon dioxide (CO2) and ethanol (EtOH) is a fundamental transport process behind the formation and growth of CO2 bubbles in sparkling beverages and the release of organoleptic compounds at the liquid free surface. In the present study, CO2 and EtOH diffusion coefficients are computed from molecular dynamics (MD) simulations and compared with experimental values derived from the Stokes-Einstein (SE) relation on the basis of viscometry experiments and hydrodynamic radii deduced from former nuclear magnetic resonance (NMR) measurements. These diffusion coefficients steadily increase with temperature and decrease as the concentration of ethanol rises. The agreement between theory and experiment is suitable for CO2. Theoretical EtOH diffusion coefficients tend to overestimate slightly experimental values, although the agreement can be improved by changing the hydrodynamic radius used to evaluate experimental diffusion coefficients. This apparent disagreement should not rely on limitations of the MD simulations nor on the approximations made to evaluate theoretical diffusion coefficients. Improvement of the molecular models, as well as additional NMR measurements on sparkling beverages at several temperatures and ethanol concentrations, would help solve this issue.


CrystEngComm ◽  
2021 ◽  
Author(s):  
Andrey Sarikov ◽  
Anna Marzegalli ◽  
Luca Barbisan ◽  
Massimo Zimbone ◽  
Corrado Bongiorno ◽  
...  

In this work, annihilation mechanism of stacking faults (SFs) in epitaxial 3C-SiC layers grown on Si(001) substrates is studied by molecular dynamics (MD) simulations. The evolution of SFs located in...


Polymers ◽  
2021 ◽  
Vol 13 (3) ◽  
pp. 347
Author(s):  
Wenlin Zhang ◽  
Lingyi Zou

We apply molecular dynamics (MD) simulations to investigate crystal nucleation in incompatible polymer blends under deep supercooling conditions. Simulations of isothermal nucleation are performed for phase-separated blends with different degrees of incompatibility. In weakly segregated blends, slow and incompatible chains in crystallizable polymer domains can significantly hinder the crystal nucleation and growth. When a crystallizable polymer is blended with a more mobile species in interfacial regions, enhanced molecular mobility leads to the fast growth of crystalline order. However, the incubation time remains the same as that in pure samples. By inducing anisotropic alignment near the interfaces of strongly segregated blends, phase separation also promotes crystalline order to grow near interfaces between different polymer domains.


2016 ◽  
Vol 18 (37) ◽  
pp. 25806-25816 ◽  
Author(s):  
Carlos Navarro-Retamal ◽  
Anne Bremer ◽  
Jans Alzate-Morales ◽  
Julio Caballero ◽  
Dirk K. Hincha ◽  
...  

Unfolding of intrinsically unstructured full-length LEA proteins in a differentially crowded environment can be modeled by 30 ns MD simulations in accordance with experimental data.


Symmetry ◽  
2021 ◽  
Vol 13 (7) ◽  
pp. 1259
Author(s):  
Maksymilian Dziura ◽  
Basel Mansour ◽  
Mitchell DiPasquale ◽  
P. Charukeshi Chandrasekera ◽  
James W. Gauld ◽  
...  

In this review, we delve into the topic of the pulmonary surfactant (PS) system, which is present in the respiratory system. The total composition of the PS has been presented and explored, from the types of cells involved in its synthesis and secretion, down to the specific building blocks used, such as the various lipid and protein components. The lipid and protein composition varies across species and between individuals, but ultimately produces a PS monolayer with the same role. As such, the composition has been investigated for the ways in which it imposes function and confers peculiar biophysical characteristics to the system as a whole. Moreover, a couple of theories/models that are associated with the functions of PS have been addressed. Finally, molecular dynamic (MD) simulations of pulmonary surfactant have been emphasized to not only showcase various group’s findings, but also to demonstrate the validity and importance that MD simulations can have in future research exploring the PS monolayer system.


2005 ◽  
Vol 1 (4) ◽  
pp. 204-209
Author(s):  
O.B. Malcıoğlu ◽  
Ş. Erkoç

The minimum energy structures of CmTin microclusters and nanoparticles have been investigated theoretically by performing molecular–dynamics (MD) simulations. Selected crystalline and completely random initial geometries are considered. The potential energy function (PEF) used in the calculations includes two– and three–body atomic interactions for C-Ti binary systems. Molecular–dynamics simulations have been performed at 1 K and 300 K. It has been found that initial geometry has a very strong influence on relaxed geometry


2013 ◽  
Vol 9 ◽  
pp. 118-134 ◽  
Author(s):  
Jutta Erika Helga Köhler ◽  
Nicole Grczelschak-Mick

Four highly ordered hydrogen-bonded models of β-cyclodextrin (β-CD) and its inclusion complex with benzene were investigated by three different theoretical methods: classical quantum mechanics (QM) on AM1 and on the BP/TZVP-DISP3 level of approximation, and thirdly by classical molecular dynamics simulations (MD) at different temperatures (120 K and 273 to 300 K). The hydrogen bonds at the larger O2/O3 rim of empty β-CDs prefer the right-hand orientation, e.g., O3-H…O2-H in the same glucose unit and bifurcated towards …O4 and O3 of the next glucose unit on the right side. On AM1 level the complex energy was −2.75 kcal mol−1 when the benzene molecule was located parallel inside the β-CD cavity and −2.46 kcal mol−1 when it was positioned vertically. The AM1 HOMO/LUMO gap of the empty β-CD with about 12 eV is lowered to about 10 eV in the complex, in agreement with data from the literature. AM1 IR spectra displayed a splitting of the O–H frequencies of cyclodextrin upon complex formation. At the BP/TZVP-DISP3 level the parallel and vertical positions from the starting structures converged to a structure where benzene assumes a more oblique position (−20.16 kcal mol−1 and −20.22 kcal mol−1, resp.) as was reported in the literature. The character of the COSMO-RS σ-surface of β-CD was much more hydrophobic on its O6 rim than on its O2/O3 side when all hydrogen bonds were arranged in a concerted mode. This static QM picture of the β-CD/benzene complex at 0 K was extended by MD simulations. At 120 K benzene was mobile but always stayed inside the cavity of β-CD. The trajectories at 273, 280, 290 and 300 K certainly no longer displayed the highly ordered hydrogen bonds of β-CD and benzene occupied many different positions inside the cavity, before it left the β-CD finally at its O2/O3 side.


Sign in / Sign up

Export Citation Format

Share Document