scholarly journals NAPS update: network analysis of molecular dynamics data and protein–nucleic acid complexes

2019 ◽  
Vol 47 (W1) ◽  
pp. W462-W470 ◽  
Author(s):  
Broto Chakrabarty ◽  
Varun Naganathan ◽  
Kanak Garg ◽  
Yash Agarwal ◽  
Nita Parekh

Abstract Network theory is now a method of choice to gain insights in understanding protein structure, folding and function. In combination with molecular dynamics (MD) simulations, it is an invaluable tool with widespread applications such as analyzing subtle conformational changes and flexibility regions in proteins, dynamic correlation analysis across distant regions for allosteric communications, in drug design to reveal alternative binding pockets for drugs, etc. Updated version of NAPS now facilitates network analysis of the complete repertoire of these biomolecules, i.e., proteins, protein–protein/nucleic acid complexes, MD trajectories, and RNA. Various options provided for analysis of MD trajectories include individual network construction and analysis of intermediate time-steps, comparative analysis of these networks, construction and analysis of average network of the ensemble of trajectories and dynamic cross-correlations. For protein–nucleic acid complexes, networks of the whole complex as well as that of the interface can be constructed and analyzed. For analysis of proteins, protein–protein complexes and MD trajectories, network construction based on inter-residue interaction energies with realistic edge-weights obtained from standard force fields is provided to capture the atomistic details. Updated version of NAPS also provides improved visualization features, interactive plots and bulk execution. URL: http://bioinf.iiit.ac.in/NAPS/

Crystals ◽  
2019 ◽  
Vol 9 (10) ◽  
pp. 532
Author(s):  
Jonathan H. Sheehan ◽  
Jarrod A. Smith ◽  
Pradeep S. Pallan ◽  
Terry P. Lybrand ◽  
Martin Egli

The (4′→6′)-linked DNA homolog 2′,3′-dideoxy-β-D-glucopyranosyl nucleic acid (dideoxy-glucose nucleic acid or homo-DNA) exhibits stable self-pairing of the Watson–Crick and reverse-Hoogsteen types, but does not cross-pair with DNA. Molecular modeling and NMR solution studies of homo-DNA duplexes pointed to a conformation that was nearly devoid of a twist and a stacking distance in excess of 4.5 Å. By contrast, the crystal structure of the homo-DNA octamer dd(CGAATTCG) revealed a right-handed duplex with average values for helical twist and rise of ca. 15° and 3.8 Å, respectively. Other key features of the structure were strongly inclined base-pair and backbone axes in the duplex with concomitant base-pair slide and cross-strand stacking, and the formation of a dimer across a crystallographic dyad with inter-duplex base swapping. To investigate the conformational flexibility of the homo-DNA duplex and a potential influence of lattice interactions on its geometry, we used molecular dynamics (MD) simulations of the crystallographically observed dimer of duplexes and an isolated duplex in the solution state. The dimer of duplexes showed limited conformational flexibility, and key parameters such as helical rise, twist, and base-pair slide exhibited only minor fluctuations. The single duplex was clearly more flexible by comparison and underwent partial unwinding, albeit without significant lengthening. Thus, base stacking was preserved in the isolated duplex and two adenosines extruded from the stack in the dimer of duplexes were reinserted into the duplex and pair with Ts in a Hoogsteen mode. Our results confirmed that efficient stacking in homo-DNA seen in the crystal structure of a dimer of duplexes was maintained in the separate duplex. Therefore, lattice interactions did not account for the different geometries of the homo-DNA duplex in the crystal and earlier models that resembled inclined ladders with large base-pair separations that precluded efficient stacking.


Molecules ◽  
2020 ◽  
Vol 25 (4) ◽  
pp. 904 ◽  
Author(s):  
Arnold Amusengeri ◽  
Rolland Bantar Tata ◽  
Özlem Tastan Bishop

In this era of precision medicine, insights into the resistance mechanism of drugs are integral for the development of potent therapeutics. Here, we sought to understand the contribution of four point mutations (N51I, C59R, S108N, and I164L) within the active site of the malaria parasite enzyme dihydrofolate reductase (DHFR) towards the resistance of the antimalarial drug pyrimethamine. Homology modeling was used to obtain full-length models of wild type (WT) and mutant DHFR. Molecular docking was employed to dock pyrimethamine onto the generated structures. Subsequent all-atom molecular dynamics (MD) simulations and binding free-energy computations highlighted that pyrimethamine’s stability and affinity inversely relates to the number of mutations within its binding site and, hence, resistance severity. Generally, mutations led to reduced binding affinity to pyrimethamine and increased conformational plasticity of DHFR. Next, dynamic residue network analysis (DRN) was applied to determine the impact of mutations and pyrimethamine binding on communication dispositions of DHFR residues. DRN revealed residues with distinctive communication profiles, distinguishing WT from drug-resistant mutants as well as pyrimethamine-bound from pyrimethamine-free models. Our results provide a new perspective on the understanding of mutation-induced drug resistance.


2014 ◽  
Vol 2014 ◽  
pp. 1-19 ◽  
Author(s):  
Kuan-Chung Chen ◽  
Calvin Yu-Chian Chen

The peroxisome proliferator-activated receptors (PPARs) related to regulation of lipid metabolism, inflammation, cell proliferation, differentiation, and glucose homeostasis by controlling the related ligand-dependent transcription of networks of genes. They are used to be served as therapeutic targets against metabolic disorder, such as obesity, dyslipidemia, and diabetes; especially, PPAR-γis the most extensively investigated isoform for the treatment of dyslipidemic type 2 diabetes. In this study, we filter compounds of traditional Chinese medicine (TCM) using bioactivities predicted by three distinct prediction models before the virtual screening. For the top candidates, the molecular dynamics (MD) simulations were also utilized to investigate the stability of interactions between ligand and PPAR-γprotein. The top two TCM candidates, 5-hydroxy-L-tryptophan and abrine, have an indole ring and carboxyl group to form the H-bonds with the key residues of PPAR-γprotein, such as residues Ser289 and Lys367. The secondary amine group of abrine also stabilized an H-bond with residue Ser289. From the figures of root mean square fluctuations (RMSFs), the key residues were stabilized in protein complexes with 5-Hydroxy-L-tryptophan and abrine as control. Hence, we propose 5-hydroxy-L-tryptophan and abrine as potential lead compounds for further study in drug development process with the PPAR-γprotein.


2020 ◽  
Author(s):  
Jordi Juárez-Jiménez ◽  
Philip Tew ◽  
Michael o'connor ◽  
Salome Llabres ◽  
Rebecca Sage ◽  
...  

<p>Molecular dynamics (MD) simulations are increasingly used to elucidate relationships between protein structure, dynamics and their biological function. Currently it is extremely challenging to perform MD simulations of large-scale structural rearrangements in proteins that occur on millisecond timescales or beyond, as this requires very significant computational resources, or the use of cumbersome ‘collective variable’ enhanced sampling protocols. Here we describe a framework that combines ensemble MD simulations and virtual-reality visualization (eMD-VR) to enable users to interactively generate realistic descriptions of large amplitude, millisecond timescale protein conformational changes in proteins. Detailed tests demonstrate that eMD-VR substantially decreases the computational cost of folding simulations of a WW domain, without the need to define collective variables <i>a priori</i>. We further show that eMD-VR generated pathways can be combined with Markov State Models to describe the thermodynamics and kinetics of large-scale loop motions in the enzyme cyclophilin A. Our results suggest eMD-VR is a powerful tool for exploring protein energy landscapes in bioengineering efforts. </p>


2020 ◽  
Author(s):  
Jordi Juárez-Jiménez ◽  
Philip Tew ◽  
Michael o'connor ◽  
Salome Llabres ◽  
Rebecca Sage ◽  
...  

<p>Molecular dynamics (MD) simulations are increasingly used to elucidate relationships between protein structure, dynamics and their biological function. Currently it is extremely challenging to perform MD simulations of large-scale structural rearrangements in proteins that occur on millisecond timescales or beyond, as this requires very significant computational resources, or the use of cumbersome ‘collective variable’ enhanced sampling protocols. Here we describe a framework that combines ensemble MD simulations and virtual-reality visualization (eMD-VR) to enable users to interactively generate realistic descriptions of large amplitude, millisecond timescale protein conformational changes in proteins. Detailed tests demonstrate that eMD-VR substantially decreases the computational cost of folding simulations of a WW domain, without the need to define collective variables <i>a priori</i>. We further show that eMD-VR generated pathways can be combined with Markov State Models to describe the thermodynamics and kinetics of large-scale loop motions in the enzyme cyclophilin A. Our results suggest eMD-VR is a powerful tool for exploring protein energy landscapes in bioengineering efforts. </p>


Blood ◽  
2007 ◽  
Vol 110 (11) ◽  
pp. 2917-2917
Author(s):  
Tai-Sung Lee ◽  
Steven Potts ◽  
Hagop Kantarjian ◽  
Jorge Cortes ◽  
Francis Giles ◽  
...  

Abstract Molecular dynamics (MD) simulations on the complex of imatinib with the wild-type, T315I, and other 10 P-loop mutants of the tyrosine kinase Bcr-Abl have been performed to study the imatinib resistance mechanism at the atomic level. MD simulations show that large scale computational simulations could offer insight information that a static structure or simple homology modeling methods cannot provide for studying the Bcr-Abl imatinib resistance problem, especially in the case of conformational changes due to remote mutations. By utilizing the Molecular Mechanics/Poisson-Boltzmann surface area (MM-PBSA) techniques and analyzing the interactions between imatinib and individual residues, imatinib resistance mechanisms not previously thought have been revealed. Non-directly contacted P-loop mutations either unfavorably change the direct electrostatic interactions with imatinib, or cause the conformational changes influencing the contact energies between imatinib and other non-P-loop residues. We demonstrate that imatinib resistance of T315I mainly comes from the breakdown of the interactions between imatinib and E286 and M290, contradictory to previously suggested that the missing hydrogen bonding is the main contribution. We also demonstrate that except for the mutations of the direct contact residues, such as L248 and Y253, the unfavorable electrostatic interaction between P-loop and imatinib is the main reason for resistance for the P-loop mutations. Furthermore, in Y255H, protonation of the histidin is essential for rendering this mutation resistant to Gleevec. Our results demonstrate that MD is a powerful way to verify and predict clinical response or resistance to imatinib and other potential drugs.


2002 ◽  
Vol 82 (1) ◽  
pp. 93-98 ◽  
Author(s):  
Thomas W. Lynch ◽  
Dorina Kosztin ◽  
Mark A. McLean ◽  
Klaus Schulten ◽  
Stephen G. Sligar

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Balint Dudas ◽  
Daniel Toth ◽  
David Perahia ◽  
Arnaud B. Nicot ◽  
Erika Balog ◽  
...  

AbstractSulfotransferases (SULTs) are phase II drug-metabolizing enzymes catalyzing the sulfoconjugation from the co-factor 3′-phosphoadenosine 5′-phosphosulfate (PAPS) to a substrate. It has been previously suggested that a considerable shift of SULT structure caused by PAPS binding could control the capability of SULT to bind large substrates. We employed molecular dynamics (MD) simulations and the recently developed approach of MD with excited normal modes (MDeNM) to elucidate molecular mechanisms guiding the recognition of diverse substrates and inhibitors by SULT1A1. MDeNM allowed exploring an extended conformational space of PAPS-bound SULT1A1, which has not been achieved up to now by using classical MD. The generated ensembles combined with docking of 132 SULT1A1 ligands shed new light on substrate and inhibitor binding mechanisms. Unexpectedly, our simulations and analyses on binding of the substrates estradiol and fulvestrant demonstrated that large conformational changes of the PAPS-bound SULT1A1 could occur independently of the co-factor movements that could be sufficient to accommodate large substrates as fulvestrant. Such structural displacements detected by the MDeNM simulations in the presence of the co-factor suggest that a wider range of drugs could be recognized by PAPS-bound SULT1A1 and highlight the utility of including MDeNM in protein–ligand interactions studies where major rearrangements are expected.


Sign in / Sign up

Export Citation Format

Share Document