Closing the gap: an atomistic structural and functional perspective of S. mansoni universal stress G4LZI3 protein in complex with phenolic compounds

Author(s):  
Priscilla Masamba ◽  
Geraldene Munsamy ◽  
Abidemi Paul Kappo

Background: For decades, Praziquantel has been the undisputed drug of choice for all schistosome infections, but rising concerns due to the unelucidated mechanism of action of the drug and unavoidable reports of emerging drug resistant strains has necessitated the need for alternative treatment drug. Moreover, current apprehension has been reinforced by total dependence on the drug for treatment hence, the search for novel and effective anti-schistosomal drugs. Uses: This study made use of bioinformatic tools to determine the structural binding of the Universal G4LZI3 stress protein (USP) in complex with ten polyphenol compounds, thereby highlighting the effectiveness of these recently identified ‘lead’ molecules in the design of novel therapeutics targeted against schistosomiasis. Upregulation of the G4LZI3 USP throughout the schistosome multifaceted developmental cycle sparks interest in its potential role as a druggable target. The integration of in silico tools provides an atomistic perspective into the binding of potential inhibitors to target proteins. Conclusion: This study therefore, implemented the use of molecular dynamic (MD) simulations to provide functional and structural insight into key conformational changes upon binding of G4ZLI3 to these key phenolic compounds. Post-MD analyses revealed unique structural and conformational changes in the G4LZI3 protein in complex with curcumin and catechin respectively. These systems exhibited the highest binding energies, while the major interacting residues conserved in all the complexes provides a route map for structure-based drug design of novel compounds with enhanced inhibitory potency against the G4LZI3 protein. This study suggests an alternative approach for the development of anti-schistosomal drugs using natural compounds.

2021 ◽  
Vol 8 ◽  
Author(s):  
Edward King ◽  
Erick Aitchison ◽  
Han Li ◽  
Ray Luo

The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Raghavendar Reddy Sanganna Gari ◽  
Joel José Montalvo‐Acosta ◽  
George R. Heath ◽  
Yining Jiang ◽  
Xiaolong Gao ◽  
...  

AbstractConformational changes in ion channels lead to gating of an ion-conductive pore. Ion flux has been measured with high temporal resolution by single-channel electrophysiology for decades. However, correlation between functional and conformational dynamics remained difficult, lacking experimental techniques to monitor sub-millisecond conformational changes. Here, we use the outer membrane protein G (OmpG) as a model system where loop-6 opens and closes the β-barrel pore like a lid in a pH-dependent manner. Functionally, single-channel electrophysiology shows that while closed states are favored at acidic pH and open states are favored at physiological pH, both states coexist and rapidly interchange in all conditions. Using HS-AFM height spectroscopy (HS-AFM-HS), we monitor sub-millisecond loop-6 conformational dynamics, and compare them to the functional dynamics from single-channel recordings, while MD simulations provide atomistic details and energy landscapes of the pH-dependent loop-6 fluctuations. HS-AFM-HS offers new opportunities to analyze conformational dynamics at timescales of domain and loop fluctuations.


2015 ◽  
Vol 43 (5) ◽  
pp. 1023-1032 ◽  
Author(s):  
Thomas Stockner ◽  
Anna Mullen ◽  
Fraser MacMillan

ABC transporters are primary active transporters found in all kingdoms of life. Human multidrug resistance transporter ABCB1, or P-glycoprotein, has an extremely broad substrate spectrum and confers resistance against chemotherapy drug treatment in cancer cells. The bacterial ABC transporter MsbA is a lipid A flippase and a homolog to the human ABCB1 transporter, with which it partially shares its substrate spectrum. Crystal structures of MsbA and ABCB1 have been solved in multiple conformations, providing a glimpse into the possible conformational changes the transporter could be going through during the transport cycle. Crystal structures are inherently static, while a dynamic picture of the transporter in motion is needed for a complete understanding of transporter function. Molecular dynamics (MD) simulations and electron paramagnetic resonance (EPR) spectroscopy can provide structural information on ABC transporters, but the strength of these two methods lies in the potential to characterise the dynamic regime of these transporters. Information from the two methods is quite complementary. MD simulations provide an all atom dynamic picture of the time evolution of the molecular system, though with a narrow time window. EPR spectroscopy can probe structural, environmental and dynamic properties of the transporter in several time regimes, but only through the attachment sites of an exogenous spin label. In this review the synergistic effects that can be achieved by combining the two methods are highlighted, and a brief methodological background is also presented.


2013 ◽  
Vol 12 (08) ◽  
pp. 1341005 ◽  
Author(s):  
FÁTIMA PARDO-AVILA ◽  
LIN-TAI DA ◽  
YING WANG ◽  
XUHUI HUANG

RNA polymerase is the enzyme that synthesizes RNA during the transcription process. To understand its mechanism, structural studies have provided us pictures of the series of steps necessary to add a new nucleotide to the nascent RNA chain, the steps altogether known as the nucleotide addition cycle (NAC). However, these static snapshots do not provide dynamic information of these processes involved in NAC, such as the conformational changes of the protein and the atomistic details of the catalysis. Computational studies have made efforts to fill these knowledge gaps. In this review, we provide examples of different computational approaches that have improved our understanding of the transcription elongation process for RNA polymerase, such as normal mode analysis, molecular dynamic (MD) simulations, Markov state models (MSMs). We also point out some unsolved questions that could be addressed using computational tools in the future.


2020 ◽  
Vol 48 (9) ◽  
pp. e49-e49 ◽  
Author(s):  
Shreya Ghosh ◽  
Matthew J Lawless ◽  
Hanna J Brubaker ◽  
Kevin Singewald ◽  
Michael R Kurpiewski ◽  
...  

Abstract Electron paramagnetic resonance (EPR) has become an important tool to probe conformational changes in nucleic acids. An array of EPR labels for nucleic acids are available, but they often come at the cost of long tethers, are dependent on the presence of a particular nucleotide or can be placed only at the termini. Site directed incorporation of Cu2+-chelated to a ligand, 2,2′dipicolylamine (DPA) is potentially an attractive strategy for site-specific, nucleotide independent Cu2+-labelling in DNA. To fully understand the potential of this label, we undertook a systematic and detailed analysis of the Cu2+-DPA motif using EPR and molecular dynamics (MD) simulations. We used continuous wave EPR experiments to characterize Cu2+ binding to DPA as well as optimize Cu2+ loading conditions. We performed double electron-electron resonance (DEER) experiments at two frequencies to elucidate orientational selectivity effects. Furthermore, comparison of DEER and MD simulated distance distributions reveal a remarkable agreement in the most probable distances. The results illustrate the efficacy of the Cu2+-DPA in reporting on DNA backbone conformations for sufficiently long base pair separations. This labelling strategy can serve as an important tool for probing conformational changes in DNA upon interaction with other macromolecules.


2020 ◽  
Vol 21 (6) ◽  
pp. 2167
Author(s):  
Jingxuan Zhu ◽  
Congcong Li ◽  
Hengzheng Yang ◽  
Xiaoqing Guo ◽  
Tianci Huang ◽  
...  

Activation of the mitogen-activated protein kinase (MAPK) signaling pathway regulated by human MAP kinase 1 (MEK1) is associated with the carcinogenesis and progression of numerous cancers. In addition, two active mutations (P124S and E203K) have been reported to enhance the activity of MEK1, thereby eventually leading to the tumorigenesis of cancer. Trametinib is an MEK1 inhibitor for treating EML4-ALK-positive, EGFR-activated, and KRAS-mutant lung cancers. Therefore, in this study, molecular docking and molecular dynamic (MD) simulations were performed to explore the effects of inactive/active mutations (A52V/P124S and E203K) on the conformational changes of MEK1 and the changes in the interaction of MEK1 with trametinib. Moreover, steered molecular dynamic (SMD) simulations were further utilized to compare the dissociation processes of trametinib from the wild-type (WT) MEK1 and two active mutants (P124S and E203K). As a result, trametinib had stronger interactions with the non-active MEK1 (WT and A52V mutant) than the two active mutants (P124S and E203K). Moreover, two active mutants may make the allosteric channel of MEK1 wider and shorter than that of the non-active types (WT and A52V mutant). Hence, trametinib could dissociate from the active mutants (P124S and E203K) more easily compared with the WT MEK1. In summary, our theoretical results demonstrated that the active mutations may attenuate the inhibitory effects of MEK inhibitor (trametinib) on MEK1, which could be crucial clues for future anti-cancer treatment.


Molecules ◽  
2020 ◽  
Vol 25 (17) ◽  
pp. 3822 ◽  
Author(s):  
R.P. Vivek-Ananth ◽  
Abhijit Rana ◽  
Nithin Rajan ◽  
Himansu S. Biswal ◽  
Areejit Samal

Presently, there are no approved drugs or vaccines to treat COVID-19, which has spread to over 200 countries and at the time of writing was responsible for over 650,000 deaths worldwide. Recent studies have shown that two human proteases, TMPRSS2 and cathepsin L, play a key role in host cell entry of SARS-CoV-2. Importantly, inhibitors of these proteases were shown to block SARS-CoV-2 infection. Here, we perform virtual screening of 14,011 phytochemicals produced by Indian medicinal plants to identify natural product inhibitors of TMPRSS2 and cathepsin L. AutoDock Vina was used to perform molecular docking of phytochemicals against TMPRSS2 and cathepsin L. Potential phytochemical inhibitors were filtered by comparing their docked binding energies with those of known inhibitors of TMPRSS2 and cathepsin L. Further, the ligand binding site residues and non-covalent interactions between protein and ligand were used as an additional filter to identify phytochemical inhibitors that either bind to or form interactions with residues important for the specificity of the target proteases. This led to the identification of 96 inhibitors of TMPRSS2 and 9 inhibitors of cathepsin L among phytochemicals of Indian medicinal plants. Further, we have performed molecular dynamics (MD) simulations to analyze the stability of the protein-ligand complexes for the three top inhibitors of TMPRSS2 namely, qingdainone, edgeworoside C and adlumidine, and of cathepsin L namely, ararobinol, (+)-oxoturkiyenine and 3α,17α-cinchophylline. Interestingly, several herbal sources of identified phytochemical inhibitors have antiviral or anti-inflammatory use in traditional medicine. Further in vitro and in vivo testing is needed before clinical trials of the promising phytochemical inhibitors identified here.


2019 ◽  
Vol 20 (19) ◽  
pp. 4721 ◽  
Author(s):  
Tripathi ◽  
Vetrivel ◽  
Téletchéa ◽  
Jean ◽  
Legembre ◽  
...  

The enzyme phospholipase C gamma 1 (PLCγ1) has been identified as a potential drug target of interest for various pathological conditions such as immune disorders, systemic lupus erythematosus, and cancers. Targeting its SH3 domain has been recognized as an efficient pharmacological approach for drug discovery against PLCγ1. Therefore, for the first time, a combination of various biophysical methods has been employed to shed light on the atomistic interactions between PLCγ1 and its known binding partners. Indeed, molecular modeling of PLCγ1 with SLP76 peptide and with previously reported inhibitors (ritonavir, anethole, daunorubicin, diflunisal, and rosiglitazone) facilitated the identification of the common critical residues (Gln805, Arg806, Asp808, Glu809, Asp825, Gly827, and Trp828) as well as the quantification of their interaction through binding energies calculations. These features are in agreement with previous experimental data. Such an in depth biophysical analysis of each complex provides an opportunity to identify new inhibitors through pharmacophore mapping, molecular docking and MD simulations. From such a systematic procedure, a total of seven compounds emerged as promising inhibitors, all characterized by a strong binding with PLCγ1 and a comparable or higher binding affinity to ritonavir (∆Gbind < −25 kcal/mol), one of the most potent inhibitor reported till now.


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):  
Yangang Pan ◽  
Luda S. Shlyakhtenko ◽  
Yuri L. Lyubchenko

AbstractViral infectivity factor (Vif) is a protein that is essential for the replication of the HIV-1 virus. The key function of Vif is to disrupt the antiviral activity of APOBEC3 proteins, which mutate viral nucleic acids. Inside the cell, Vif binds to the host cell proteins Elongin-C, Elongin-B, and CBF-β, forming a four-protein complex called VCBC. The structure of VCBC in complex with the Cullin5 (Cul5) protein has been solved by X-ray crystallography, and recently, using molecular dynamic (MD) simulations, the dynamics of VCBC and VCBC-Cul5 complexes were characterized. Here, we applied time-lapse high-speed atomic force microscopy (HS-AFM) to visualize the conformational changes of the VCBC complex. We determined the three most favorable conformations of the VCBC complex, which we identified as triangle, dumbbell, and globular structures. In addition, we characterized the dynamics of each of these structures. While our data show a very dynamic behavior for all these structures, we found the triangle and dumbbell structures to be the most dynamic. These findings provide insight into the structure and dynamics of the VCBC complex and support further research into the improvement of HIV treatment, as Vif is essential for virus survival in the cell.


Sign in / Sign up

Export Citation Format

Share Document