1
|
Pirnia A, Maqdisi R, Mittal S, Sener M, Singharoy A. Perspective on Integrative Simulations of Bioenergetic Domains. J Phys Chem B 2024; 128:3302-3319. [PMID: 38562105 DOI: 10.1021/acs.jpcb.3c07335] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.
Collapse
Affiliation(s)
- Adam Pirnia
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
| | - Ranel Maqdisi
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
| | - Sumit Mittal
- VIT Bhopal University, Sehore 466114, Madhya Pradesh, India
| | - Melih Sener
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
- Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Abhishek Singharoy
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1004, United States
| |
Collapse
|
2
|
Chio CC, Tse YLS. Reparameterization of Polarizable Force Fields for Studying Ion Transfer across Liquid-Liquid Interfaces. J Phys Chem B 2024; 128:1987-1999. [PMID: 38356148 DOI: 10.1021/acs.jpcb.3c07055] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2024]
Abstract
We have developed a general scheme for refining classical polarizable molecular dynamics (MD) force fields that can accurately describe the molecular interactions in systems with liquid-liquid interfaces. While ab initio MD (AIMD) simulations can naturally describe molecular interactions, they are often so computationally expensive that simulating large system sizes and/or long time scales is usually infeasible. To resolve this, we parameterized efficient and accurate classical polarizable force fields that use AIMD reference data by minimizing both the relative entropy and the root mean squared deviation in atomic forces. We utilized our new multiscale models to study chloride ion transfer across the water-dichloromethane (DCM) interface with and without the tetraethylammonium cation as the phase-transfer catalyst. Our calculated free-energy barrier for the water-DCM interface is consistent with the other reported simulation results. We further analyzed the ion-transfer process by studying the hydration shell structures around the chloride ion and the ion-pair formation to better understand the mechanism. We observed that electronic polarizability is an important factor for the studied phase-transfer catalyst to lower the free-energy barrier of the ion transfer. Using the water-benzene interface system as an additional example, we show that our parameterization scheme provides a general route for modeling different liquid-liquid interface systems even when the experimental data or force field parameters are not readily available.
Collapse
Affiliation(s)
- Chung Chi Chio
- Department of Chemistry, The Chinese University of Hong Kong, Sha Tin, New Territories, Hong Kong, China
| | - Ying-Lung Steve Tse
- Department of Chemistry, The Chinese University of Hong Kong, Sha Tin, New Territories, Hong Kong, China
| |
Collapse
|
3
|
Almeida ER, Goliatt PVZC, Dos Santos HF, Picaud F. Modeling the Cellular Uptake of Functionalized Carbon Nanohorns Loaded with Cisplatin through a Breast Cancer Cell Membrane. Mol Pharm 2024; 21:38-52. [PMID: 37646561 DOI: 10.1021/acs.molpharmaceut.3c00379] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/01/2023]
Abstract
The cisplatin encapsulation into carbon nanohorns (CNH) is a promising nanoformulation to circumvent the drug dissipation and to specifically accumulate it in tumor sites. Herein, biased molecular dynamics simulations were used to analyze the transmembrane transport of the CNH loaded with cisplatin through a breast cancer cell membrane prototype. The simulations revealed a four-stage mechanism: approach, insertion, permeation, and internalization. Despite the lowest structural disturbance of the membrane provided by the nanocarrier, the average free energy barrier for the translocation was 55.2 kcal mol-1, suggesting that the passive process is kinetically unfavorable. In contrast, the free energy profiles revealed potential wells of -6.8 kcal mol-1 along the insertion stage in the polar heads region of the membrane, which might enhance the retention of the drug in tumor sites; therefore, the most likely cisplatin delivery mechanism should involve the adsorption and retention of CNH on the surface of cancer cells, allowing the loaded cisplatin be slowly released and passively transported through the cell membrane.
Collapse
Affiliation(s)
- Eduardo R Almeida
- Núcleo de Estudos em Química Computacional (NEQC), Departamento de Química, Instituto de Ciências Exatas, Universidade Federal de Juiz de Fora (UFJF), Campus Universitário, Martelos, Juiz de Fora, Minas Gerais 36036-330, Brazil
- Laboratoire de Nanomédecine, Imagerie et Thérapeutiques, EA 4662, Centre Hospitalier Universitaire de Besançon, Université de Franche-Comté, 16 route de Gray, 25030 Besançon, Cedex, France
| | - Priscila V Z Capriles Goliatt
- Grupo de Modelagem Computacional Aplicada (GMCA), Programa de Pós-Graduação em Modelagem Computacional (PGMC), Departamento de Ciência da Computação, Instituto de Ciências Exatas, Universidade Federal de Juiz de Fora (UFJF), Campus Universitário, Martelos, Juiz de Fora, Minas Gerais 36036-330, Brazil
| | - Hélio F Dos Santos
- Núcleo de Estudos em Química Computacional (NEQC), Departamento de Química, Instituto de Ciências Exatas, Universidade Federal de Juiz de Fora (UFJF), Campus Universitário, Martelos, Juiz de Fora, Minas Gerais 36036-330, Brazil
| | - Fabien Picaud
- Laboratoire de Nanomédecine, Imagerie et Thérapeutiques, EA 4662, Centre Hospitalier Universitaire de Besançon, Université de Franche-Comté, 16 route de Gray, 25030 Besançon, Cedex, France
| |
Collapse
|
4
|
Almeida ER, Goliatt PVZC, Dos Santos HF, Picaud F. Translocation Processes of Pt(II)-Based Drugs through Human Breast Cancer Cell Membrane: In Silico Experiments. J Chem Inf Model 2023; 63:6141-6155. [PMID: 37751589 DOI: 10.1021/acs.jcim.3c00943] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/28/2023]
Abstract
Breast cancer is one of the most frequent modalities of cancer worldwide, with notable mortality. The medication based on platinum drugs (cisplatin (cddp), carboplatin (cpx), and oxaliplatin (oxa)) is a conventional chemotherapy despite severe side effects and the development of drug resistance. In order to provide a deeper molecular description of the influx and efflux processes of platinum drugs through breast cancer tissues, this study focuses on molecular dynamics (MD) simulations of the passive translocation process through a realistic plasma membrane prototype of human breast cancer cell (c_memb). The results showed that the permeation events were mainly mediated by neutral lipids (DOPC, DOPE, and cholesterol), producing a low and temporary membrane deformation. The drug insertion in the region of polar heads was the most favorable stage of the translocation mechanism, especially for cddp and oxa with potential wells of -8.6 and -9.8 kcal mol-1, respectively. However, the potentials of mean force (PMF) revealed unfavorable kinetics for the permeation of these drugs through lipid tails, with energy barriers of 28.3 (cddp), 32.2 (cpx), and 30.4 kcal mol-1 (oxa). The low permeability coefficients (P) of cpx and oxa, which were 3 and 1 orders of magnitude inferior than for cddp, resulted from the high energy barriers for their traslocation processes through the membrane. The obtained results provide a more accurate picture of the permeation of Pt(II)-based drugs through breast cancer cells, which may be relevant for the design and evaluation of new platinum complexes.
Collapse
Affiliation(s)
- Eduardo R Almeida
- Núcleo de Estudos em Química Computacional (NEQC), Departamento de Química, ICE, Universidade Federal de Juiz de Fora (UFJF), Campus Universitário, Martelos, Juiz de Fora, MG 36036-330, Brazil
- Laboratoire de Nanomédecine, Imagerie et Thérapeutiques, EA 4662, Université de Franche-Comté, Centre Hospitalier Universitaire de Besançon, 16 route de Gray, 25030 Besançon cedex, France
| | - Priscila V Z Capriles Goliatt
- Programa de Pós-graduação em Modelagem Computacional (PGMC), Departamento de Ciência da Computação, Instituto de Ciências Exatas, Universidade Federal de Juiz de Fora (UFJF), Campus Universitário, Martelos, Juiz de Fora, MG 36036-330, Brazil
| | - Hélio F Dos Santos
- Núcleo de Estudos em Química Computacional (NEQC), Departamento de Química, ICE, Universidade Federal de Juiz de Fora (UFJF), Campus Universitário, Martelos, Juiz de Fora, MG 36036-330, Brazil
| | - Fabien Picaud
- Laboratoire de Nanomédecine, Imagerie et Thérapeutiques, EA 4662, Université de Franche-Comté, Centre Hospitalier Universitaire de Besançon, 16 route de Gray, 25030 Besançon cedex, France
| |
Collapse
|
5
|
Iida S, Tomoshi K. Free energy and kinetic rate calculation via non-equilibrium molecular simulation: application to biomolecules. Biophys Rev 2022; 14:1303-1314. [PMID: 36659997 PMCID: PMC9842846 DOI: 10.1007/s12551-022-01036-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Accepted: 11/26/2022] [Indexed: 12/30/2022] Open
Abstract
Non-equilibrium molecular dynamics (NEMD) simulation has been recognized as a powerful tool for examining biomolecules and provides fruitful insights into not only non-equilibrium but also equilibrium processes. We review recent advances in NEMD simulation and relevant, fundamental results of non-equilibrium statistical mechanics. We first introduce Crooks fluctuation theorem and Jarzynski equality that relate free energy difference to work done on a physical system during a non-equilibrium process. The theorems are beneficial for the analysis of NEMD trajectories. We then describe rate theory, a framework to calculate molecular kinetics from a non-equilibrium process; this theoretical framework enables us to calculate a reaction time-mean-first passage time-from NEMD trajectories. We, in turn, present recent NEMD techniques that apply an external force to a system to enhance molecular dissociation and introduce their application to biomolecules. Lastly, we show the current status of an appropriate selection of reaction coordinates for NEMD simulation.
Collapse
Affiliation(s)
- Shinji Iida
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-Ku, Tokyo, 135-0064 Japan
| | - Kameda Tomoshi
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST), 2-4-7 Aomi, Koto-Ku, Tokyo, 135-0064 Japan
| |
Collapse
|
6
|
Chauhan M, Bhardwaj VK, Kumar A, Kumar V, Kumar P, Enayathullah MG, Thomas J, George J, Kumar BK, Purohit R, Kumar A, Kumar S. Theaflavin 3-gallate inhibits the main protease (M pro) of SARS-CoV-2 and reduces its count in vitro. Sci Rep 2022; 12:13146. [PMID: 35908093 PMCID: PMC9338964 DOI: 10.1038/s41598-022-17558-5] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2022] [Accepted: 07/27/2022] [Indexed: 12/26/2022] Open
Abstract
The main protease (Mpro) of SARS-CoV-2 has been recognized as an attractive drug target because of its central role in viral replication. Our previous preliminary molecular docking studies showed that theaflavin 3-gallate (a natural bioactive molecule derived from theaflavin and found in high abundance in black tea) exhibited better docking scores than repurposed drugs (Atazanavir, Darunavir, Lopinavir). In this study, conventional and steered MD-simulations analyses revealed stronger interactions of theaflavin 3-gallate with the active site residues of Mpro than theaflavin and a standard molecule GC373 (a known inhibitor of Mpro and novel broad-spectrum anti-viral agent). Theaflavin 3-gallate inhibited Mpro protein of SARS-CoV-2 with an IC50 value of 18.48 ± 1.29 μM. Treatment of SARS-CoV-2 (Indian/a3i clade/2020 isolate) with 200 μM of theaflavin 3-gallate in vitro using Vero cells and quantifying viral transcripts demonstrated reduction of viral count by 75% (viral particles reduced from Log106.7 to Log106.1). Overall, our findings suggest that theaflavin 3-gallate effectively targets the Mpro thus limiting the replication of the SARS-CoV-2 virus in vitro.
Collapse
Affiliation(s)
- Mahima Chauhan
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India.,Academy of Scientific and Innovative Research, Ghaziabad, Uttar Pradesh, 201002, India
| | - Vijay Kumar Bhardwaj
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India.,Academy of Scientific and Innovative Research, Ghaziabad, Uttar Pradesh, 201002, India.,Structural Bioinformatics Lab, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India
| | - Asheesh Kumar
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India.,Academy of Scientific and Innovative Research, Ghaziabad, Uttar Pradesh, 201002, India
| | - Vinod Kumar
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India
| | - Pawan Kumar
- Academy of Scientific and Innovative Research, Ghaziabad, Uttar Pradesh, 201002, India.,Chemical Technology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India
| | - M Ghalib Enayathullah
- CSIR-Center for Cellular and Molecular Biology, Annexe-II, Medical Biotechnology Complex, Uppal Road, Hyderabad, Telangana, 500007, India
| | - Jessie Thomas
- CSIR-Center for Cellular and Molecular Biology, Annexe-II, Medical Biotechnology Complex, Uppal Road, Hyderabad, Telangana, 500007, India
| | - Joel George
- CSIR-Center for Cellular and Molecular Biology, Annexe-II, Medical Biotechnology Complex, Uppal Road, Hyderabad, Telangana, 500007, India
| | - Bokara Kiran Kumar
- CSIR-Center for Cellular and Molecular Biology, Annexe-II, Medical Biotechnology Complex, Uppal Road, Hyderabad, Telangana, 500007, India.
| | - Rituraj Purohit
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India. .,Academy of Scientific and Innovative Research, Ghaziabad, Uttar Pradesh, 201002, India. .,Structural Bioinformatics Lab, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India.
| | - Arun Kumar
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India. .,Academy of Scientific and Innovative Research, Ghaziabad, Uttar Pradesh, 201002, India.
| | - Sanjay Kumar
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India
| |
Collapse
|
7
|
Kashyap P, Bhardwaj VK, Chauhan M, Chauhan V, Kumar A, Purohit R, Kumar A, Kumar S. A ricin-based peptide BRIP from Hordeum vulgare inhibits M pro of SARS-CoV-2. Sci Rep 2022; 12:12802. [PMID: 35896605 PMCID: PMC9326418 DOI: 10.1038/s41598-022-15977-y] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2022] [Accepted: 07/01/2022] [Indexed: 12/13/2022] Open
Abstract
COVID-19 pandemic caused by SARS-CoV-2 led to the research aiming to find the inhibitors of this virus. Towards this world problem, an attempt was made to identify SARS-CoV-2 main protease (Mpro) inhibitory peptides from ricin domains. The ricin-based peptide from barley (BRIP) was able to inhibit Mpro in vitro with an IC50 of 0.52 nM. Its low and no cytotoxicity upto 50 µM suggested its therapeutic potential against SARS-CoV-2. The most favorable binding site on Mpro was identified by molecular docking and steered molecular dynamics (MD) simulations. The Mpro-BRIP interactions were further investigated by evaluating the trajectories for microsecond timescale MD simulations. The structural parameters of Mpro-BRIP complex were stable, and the presence of oppositely charged surfaces on the binding interface of BRIP and Mpro complex further contributed to the overall stability of the protein-peptide complex. Among the components of thermodynamic binding free energy, Van der Waals and electrostatic contributions were most favorable for complex formation. Our findings provide novel insight into the area of inhibitor development against COVID-19.
Collapse
Affiliation(s)
- Prakriti Kashyap
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India
| | - Vijay Kumar Bhardwaj
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India.,Academy of Scientific and Innovative Research, Ghaziabad, 201002, Uttar Pradesh, India.,Structural Bioinformatics Lab, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India
| | - Mahima Chauhan
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India.,Academy of Scientific and Innovative Research, Ghaziabad, 201002, Uttar Pradesh, India
| | - Varun Chauhan
- Covid-19 Testing Facility, Dietetics & Nutrition Technology Division, Council of Scientific and Industrial Research-Institute of Himalayan Bioresource Technology (CSIR-IHBT), Palampur, H.P, India, 176061
| | - Asheesh Kumar
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India.,Academy of Scientific and Innovative Research, Ghaziabad, 201002, Uttar Pradesh, India
| | - Rituraj Purohit
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India. .,Academy of Scientific and Innovative Research, Ghaziabad, 201002, Uttar Pradesh, India. .,Structural Bioinformatics Lab, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India.
| | - Arun Kumar
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India. .,Academy of Scientific and Innovative Research, Ghaziabad, 201002, Uttar Pradesh, India.
| | - Sanjay Kumar
- Biotechnology Division, CSIR-Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh, 176061, India.
| |
Collapse
|
8
|
Almeida ER, Capriles PVSZ, Dos Santos HF. Unveiling the Releasing Processes of Pt(II)-Based Anticancer Drugs from Oxidized Carbon Nanohorn: An In Silico Study. J Phys Chem B 2022; 126:4246-4260. [PMID: 35670834 DOI: 10.1021/acs.jpcb.2c02555] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
About half of all cancer chemotherapies currently applied involve medication with the three worldwide approved Pt(II)-based drugs, cisplatin (cddp), carboplatin (cpx), and oxaliplatin (oxa), due to their notable antitumor activity for several cancers. However, this wide application is accompanied by severe side effects, such as nephrotoxicity, myelosuppression, and neurotoxicity, as a result of their low bioavailability and selectivity for cancer cells. To mitigate these drawbacks, the use of chemically functionalized carbon nanohorns (CNH) as nanocarriers represents a potential formulation since CNH has been noted for their biodegradability, biocompatibility, low toxicity, and cavities dimensionally compatible with small drugs. This work reports energetic and dynamic analyses of complexes formed by oxidized CNH (CNHox) and the cddp, cpx, and oxa drugs. Using unbiased molecular dynamics (MD) simulations, we show that the encapsulated formulations (cddp@CNHox, cpx@CNHox, and oxa@CNHox) were more stable by ∼11.0 kcal mol-1 than the adsorbed ones (cddp > CNHox, cpx > CNHox, and oxa > CNHox). This high stability, mainly governed by van der Waals interactions, was responsible for the drug confinement during the entire simulation time (200 ns). The biased MD simulations of the inclusion complexes confirmed the nonspontaneity of the drug release since the potentials of mean force (PMF) indicated the endergonic character of this process. Additionally, the releasing energy profiles pointed out that the free energy barrier (ΔΔG≠) for the escape from CNHox cavity follows the order oxa > cpx ∼ cddp, with the value for the oxa complex (21-26 kcal mol-1) found to be about 36 and 30% larger than those for cpx and cddp, respectively. While the approximate residence time (tres) of the oxa drug inside the CNHox cavity was 5.45 × 108 s, the same measure for the cddp and cpx drugs was 5.3 × 105 and 1.60 × 103 s. Simulations also revealed that the escape of oxa with the oxalate group facing the nanowindow was the most unfavorable process, giving tres = 1.09 × 109 s. Besides reinforcing and extending the nanovectorization of cddp, cpx, and oxa in CNHox for cancer chemotherapies, all features considered may provide interpretations for experimental data and encourage new investigations aiming to propose less aggressive treatments for oncological diseases.
Collapse
Affiliation(s)
- Eduardo R Almeida
- Núcleo de Estudos em Química Computacional (NEQC), Departamento de Química, Instituto de Ciências Exatas, Universidade Federal de Juiz de Fora (UFJF), Campus Universitário, Martelos, Juiz de Fora, MG 36036-330, Brazil
| | - Priscila V S Z Capriles
- Programa de Pós-graduação em Modelagem Computacional (PGMC), Departamento de Ciência da Computação, Instituto de Ciências Exatas, Universidade Federal de Juiz de Fora (UFJF), Campus Universitário, Martelos, Juiz de Fora, MG 36036-330, Brazil
| | - Hélio F Dos Santos
- Núcleo de Estudos em Química Computacional (NEQC), Departamento de Química, Instituto de Ciências Exatas, Universidade Federal de Juiz de Fora (UFJF), Campus Universitário, Martelos, Juiz de Fora, MG 36036-330, Brazil
| |
Collapse
|
9
|
Prusty S, Sarkar R, Chakraborty A, Roy S. Correlation in Domain Fluctuations Navigates Target Search of a Viral Peptide along RNA. J Phys Chem B 2021; 125:12678-12689. [PMID: 34756044 DOI: 10.1021/acs.jpcb.1c07699] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Biological macromolecules often exhibit correlations in fluctuations involving distinct domains. This study decodes their functional implications in RNA-protein recognition and target-specific binding. The target search of a peptide along RNA in a viral TAR-Tat complex is closely monitored using atomistic simulations, steered molecular dynamics simulations, free energy calculations, and a machine-learning-based clustering technique. An anticorrelated domain fluctuation is identified between the tetraloop and the bulge region in the apo form of TAR RNA that sets a hierarchy in the domain-specific fluctuations at each binding event and that directs the succeeding binding footsteps. Thus, at each binding footstep, the dynamic partner selects an RNA location for binding where it senses a higher fluctuation, which is conventionally reduced upon binding. This event stimulates an alternate domain fluctuation, which then dictates sequential binding footstep/s and thus the search progresses. Our cross-correlation maps show that the fluctuations relay from one domain to another specific domain until the anticorrelation between those interdomain fluctuations sustains. Artificial attenuation of that hierarchical domain fluctuation inhibits specific RNA binding. The binding is completed with the arrival of a few long-lived water molecules that mediate slightly distant RNA-protein sites and finally stabilize the overall complex. The study underscores the functional importance of naturally designed fluctuating RNA motifs (bulge, tetraloop) and their interplay in dictating the directionality of the search in a highly dynamic environment.
Collapse
Affiliation(s)
- Sangram Prusty
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Campus Road, Mohanpur, West Bengal 741246, India
| | - Raju Sarkar
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Campus Road, Mohanpur, West Bengal 741246, India
| | - Amrita Chakraborty
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Campus Road, Mohanpur, West Bengal 741246, India
| | - Susmita Roy
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Campus Road, Mohanpur, West Bengal 741246, India
| |
Collapse
|
10
|
Kumar Bhardwaj V, Purohit R. Taming the ringmaster of the genome (PCNA): Phytomolecules for anticancer therapy against a potential non-oncogenic target. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116437] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
11
|
Almeida ER, Dos Santos HF, Capriles PVSZ. Carbon nanohorns as nanocontainers for cisplatin: insight into their interaction with the plasma membranes of normal and breast cancer cells. Phys Chem Chem Phys 2021; 23:16376-16389. [PMID: 34318822 DOI: 10.1039/d1cp02015c] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Cisplatin (cddp)-based chemotherapy is one of the most effective therapeutic alternatives for breast cancer treatment, the most common form of cancer, despite the severe side effects related to the high toxicity and low selectivity of cddp. To circumvent these drawbacks, the encapsulation of cddp into oxidized carbon nanohorns (CNHoxs) has been shown as a promising formulation with biocompatibility and low toxicity. However, there is still a lack of studies regarding the behavior of this cddp@CNHox nanovector on the cell membranes. This study presents an in silico description of the interactions between cddp@CNHox and membrane models of cancer (C_memb) and normal (N_memb) cells referring to a typical human breast. The results revealed the interaction mechanism of the inclusion complex 3cddp@CNHox (three cddp molecules are included in the CNHox cavity) with these biomembranes, which is a multistep process including approach, landing, insertion, and penetration. The 3cddp@CNHox stability was monitored over time, and demonstrated the trapping of cddp molecules inside the CNHox cavity over all simulations. The van der Waals contribution played a primary role (∼74%) for the complex stability. Moreover, the binding free energy calculations indicated that the interaction of the 3cddp@CNHox complex with the C_memb model was slightly more favorable, on average, than with the N_memb model. Analysis of the hydrogen bonds (HBs) formed over simulations of 800 ns explains the selectivity for the C_memb model, since the total number of HBs established between the inclusion complex and the C_memb model was about three times greater than that with the N_memb model. By reinforcing the potentiality of oxidized CNHox as a nanovector of cddp, the results presented in this study may assist and drive new experimental studies with this nanomaterial, focusing on the development of less aggressive formulations for breast cancer treatment.
Collapse
Affiliation(s)
- Eduardo R Almeida
- Núcleo de Estudos em Química Computacional (NEQC), Departamento de Química, Instituto de Ciências Exatas, Universidade Federal de Juiz de Fora (UFJF), Campus Universitário, Martelos, Juiz de Fora, MG 36036-330, Brazil.
| | | | | |
Collapse
|
12
|
CityApps: A bioinformatics tool for predicting the key residues of enzymes weakly interacting with monovalent metal ions. Process Biochem 2021. [DOI: 10.1016/j.procbio.2021.03.005] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
13
|
Abstract
The properties of natural lipid bilayers are vital to the regulation of many membrane proteins. Scaffolded nanodiscs provide an in vitro lipid bilayer platform to host membrane proteins in an environment that approximates native lipid bilayers. However, the properties of scaffold-enclosed bilayers may depart significantly from those of bulk cellular membranes. Therefore, to improve the usefulness of nanodiscs it is essential to understand the properties of lipids restricted by scaffolds. We used computational molecular dynamics and modeling approaches to understand the effects of nanodisc size, scaffold type (DNA or protein), and hydrophobic modification of DNA scaffolds on bilayer stability and degree to which the properties of enclosed bilayers approximate bulk bilayers. With respect to achieving bulk bilayer behavior, we found that charge neutralization of DNA scaffolds was more important than the total hydrophobic content of their modifications: bilayer properties were better for scaffolds having a large number of short alkyl chains than those having fewer long alkyl chains. Further, complete charge neutralization of DNA scaffolds enabled better lipid binding, and more stable bilayers, as shown by steered molecular dynamics simulations that measured the force required to dislodge scaffolds from lipid bilayer patches. Considered together, our simulations provide a guide to the design of DNA-scaffolded nanodiscs suitable for studying membrane proteins.
Collapse
Affiliation(s)
- Vishal Maingi
- Department of Bioengineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Paul W K Rothemund
- Departments of Bioengineering, Computing + Mathematical Sciences, and Computation & Neural Systems, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
14
|
Bennett WFD, He S, Bilodeau CL, Jones D, Sun D, Kim H, Allen JE, Lightstone FC, Ingólfsson HI. Predicting Small Molecule Transfer Free Energies by Combining Molecular Dynamics Simulations and Deep Learning. J Chem Inf Model 2020; 60:5375-5381. [DOI: 10.1021/acs.jcim.0c00318] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- W. F. Drew Bennett
- Biochemical and Biophysical Systems Group, Biosciences and Biotechnology Division, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California, United States
| | - Stewart He
- Global Security Computing Applications, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California, United States
| | - Camille L. Bilodeau
- Biochemical and Biophysical Systems Group, Biosciences and Biotechnology Division, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California, United States
| | - Derek Jones
- Global Security Computing Applications, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California, United States
| | - Delin Sun
- Biochemical and Biophysical Systems Group, Biosciences and Biotechnology Division, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California, United States
| | - Hyojin Kim
- Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California, United States
| | - Jonathan E. Allen
- Global Security Computing Applications, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California, United States
| | - Felice C. Lightstone
- Biochemical and Biophysical Systems Group, Biosciences and Biotechnology Division, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California, United States
| | - Helgi I. Ingólfsson
- Biochemical and Biophysical Systems Group, Biosciences and Biotechnology Division, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California, United States
| |
Collapse
|
15
|
|