Insight into the Initial Stages of the Folding Process in Onconase Revealed by UNRES

Molecular Dynamics Reveal Key Steps in BAR-Related Membrane Remodeling

Endocytosis plays a complex role in pathogen-host interactions. It serves as a pathway for pathogens to enter the host cell and acts as a part of the immune defense mechanism. Endocytosis involves the formation of lipid membrane vesicles and the reshaping of the cell membrane, a task predominantly managed by proteins containing BAR (Bin1/Amphiphysin/yeast RVS167) domains. Insights into how BAR domains can remodel and reshape cell membranes provide crucial information on infections and can aid the development of treatment. Aiming at deciphering the roles of the BAR dimers in lipid membrane bending and remodeling, we conducted extensive all-atom molecular dynamics simulations and discovered that the presence of helix kinks divides the BAR monomer into two segments-the "arm segment" and the "core segment"-which exhibit distinct movement patterns. Contrary to the prior hypothesis of BAR domains working as a rigid scaffold, we found that it functions in an "Arms-Hands" mode. These findings enhance the understanding of endocytosis, potentially advancing research on pathogen-host interactions and aiding in the identification of new treatment strategies targeting BAR domains.

Self-organizing maps of unbiased ligand-target binding pathways and kinetics

The interpretation of ligand-target interactions at atomistic resolution is central to most efforts in computational drug discovery and optimization. However, the highly dynamic nature of protein targets, as well as possible induced fit effects, makes difficult to sample many interactions effectively with docking studies or even with large-scale molecular dynamics (MD) simulations. We propose a novel application of Self-Organizing Maps (SOMs) to address the sampling and dynamic mapping tasks, particularly in cases involving ligand flexibility and induced fit. The SOM approach offers a data-driven strategy to create a map of the interaction process and pathways based on unbiased MD. Furthermore, we show how the preliminary SOM mapping is complementary to kinetic analysis, with the employment of both network-based approaches and Markov state models. We demonstrate the method by comprehensively mapping a large dataset of 640 μs of unbiased trajectories sampling the recognition process between the phosphorylated YEEI peptide and its high-specificity target lck-SH2. The integration of SOM into unbiased simulation protocols significantly advances our understanding of the ligand binding mechanism. This approach serves as a potent tool for mapping intricate ligand-target interactions with unprecedented detail, thereby enhancing the characterization of kinetic properties crucial to drug design.

Deciphering the Mystery in p300 Taz2-p53 TAD2 Recognition

Intrinsically disordered proteins (IDPs) engage in various fundamental biological activities, and their behavior is of particular importance for a better understanding of the verbose but well-organized signal transduction in cells. IDPs exhibit uniquely paradoxical features with low affinity but simultaneously high specificity in recognizing their binding targets. The transcription factor p53 plays a crucial role in cancer suppression, carrying out some of its biological functions using its disordered regions, such as N-terminal transactivation domain 2 (TAD2). Exploration of the binding and unbinding processes between proteins is challenging, and the inherently disordered properties of these regions further complicate the issue. Computer simulations are a powerful tool to complement the experiments to fill gaps to explore the binding/unbinding processes between proteins. Here, we investigated the binding mechanism between p300 Taz2 and p53 TAD2 through extensive molecular dynamics (MD) simulations using the physics-based UNited RESidue (UNRES) force field with additional Go̅-like potentials. Distance restraints extracted from the NMR-resolved structures were imposed on intermolecular residue pairs to accelerate binding simulations, in which Taz2 was immobilized in a native-like conformation and disordered TAD2 was fully free. Starting from six structures with TAD2 placed at different positions around Taz2, we observed a metastable intermediate state in which the middle helical segment of TAD2 is anchored in the binding pocket, highlighting the significance of the TAD2 helix in directing protein recognition. Physics-based binding simulations show that successful binding is achieved after a series of stages, including (1) protein collisions to initiate the formation of encounter complexes, (2) partial attachment of TAD2, and finally (3) full attachment of TAD2 to the correct binding pocket of Taz2. Furthermore, machine-learning-based PathDetect-SOM was used to identify two binding pathways, the encounter complexes, and the intermediate states.

Unveiling dissociation mechanisms and binding patterns in the UHRF1-DPPA3 complex via multi-replica molecular dynamics simulations

CONTEXT

Ubiquitin-like with PHD and RING finger domain containing protein 1 (UHRF1) is responsible for preserving the stability of genomic methylation through the recruitment of DNA methyltransferase 1 (DNMT1). However, the interaction between Developmental pluripotency associated 3 (DPPA3) and the pre-PHD-PHD (PPHD) domain of UHRF1 hinders the nuclear localization of UHRF1. This disruption has implications for potential cancer treatment strategies. Drugs that mimic the binding pattern between DPPA3 and PPHD could offer a promising approach to cancer treatment. Our study reveals that DPPA3 undergoes dissociation from the C-terminal through three different modes of helix unfolding. Furthermore, we have identified key residue pairs involved in this dissociation process and potential drug-targeting residues. These findings offer valuable insights into the dissociation mechanism of DPPA3 from PPHD and have the potential to inform the design of novel drugs targeting UHRF1 for cancer therapy.

METHODS

To comprehend the dissociation process and binding patterns of PPHD-DPPA3, we employed enhanced sampling techniques, including steered molecular dynamics (SMD) and conventional molecular dynamics (cMD). Additionally, we utilized self-organizing maps (SOM) and time-resolved force distribution analysis (TRFDA) methodologies. The Gromacs software was used for performing molecular dynamics simulations, and the AMBER FF14SB force field was applied to the protein.

Mechanism of RGD-conjugated nanodevice binding to its target protein integrin αVβ3 by atomistic molecular dynamics and machine learning

Active targeting strategies have been proposed to enhance the selective uptake of nanoparticles (NPs) by diseased cells, and recent experimental findings have proven the effectiveness of this approach. However, no mechanistic studies have yet revealed the atomistic details of the interactions between ligand-activated NPs and integrins. As a case study, here we investigate, by means of advanced molecular dynamics simulations (MD) and machine learning methods (namely equilibrium MD, binding free energy calculations and training of self-organized maps), the interaction of a cyclic-RGD-conjugated PEGylated TiO2 NP (the nanodevice) with the extracellular segment of integrin αVβ3 (the target), the latter experimentally well-known to be over-expressed in several solid tumors. Firstly, we proved that the cyclic-RGD ligand binding to the integrin pocket is established and kept stable even in the presence of the cumbersome realistic model of the nanodevice. In this respect, the unsupervised machine learning analysis allowed a detailed comparison of the ligand/integrin binding in the presence and in the absence of the nanodevice, which unveiled differences in the chemical features. Then, we discovered that unbound cyclic RGDs conjugated to the NP largely contribute to the interactions between the nanodevice and the integrin. Finally, by increasing the density of cyclic RGDs on the PEGylated TiO2 NP, we observed a proportional enhancement of the nanodevice/target binding. All these findings can be exploited to achieve an improved targeting selectivity and cellular uptake, and thus a more successful clinical outcome.

Molecular Dynamics for the Optimal Design of Functionalized Nanodevices to Target Folate Receptors on Tumor Cells

Atomistic details on the mechanism of targeting activity by biomedical nanodevices of specific receptors are still scarce in the literature, where mostly ligand/receptor pairs are modeled. Here, we use atomistic molecular dynamics (MD) simulations, free energy calculations, and machine learning approaches on the case study of spherical TiO2 nanoparticles (NPs) functionalized with folic acid (FA) as the targeting ligand of the folate receptor (FR). We consider different FA densities on the surface and different anchoring approaches, i.e., direct covalent bonding of FA γ-carboxylate or through polyethylene glycol spacers. By molecular docking, we first identify the lowest energy conformation of one FA inside the FR binding pocket from the X-ray crystal structure, which becomes the starting point of classical MD simulations in a realistic physiological environment. We estimate the binding free energy to be compared with the existing experimental data. Then, we increase complexity and go from the isolated FA to a nanosystem decorated with several FAs. Within the simulation time framework, we confirm the stability of the ligand-receptor interaction, even in the presence of the NP (with or without a spacer), and no significant modification of the protein secondary structure is observed. Our study highlights the crucial role played by the spacer, FA protonation state, and density, which are parameters that can be controlled during the nanodevice preparation step.

Metadynamics simulations for the investigation of drug loading on functionalized inorganic nanoparticles

Inorganic nanoparticles show promising properties that allow them to be efficiently used as drug carriers. The main limitation in this type of application is currently the drug loading capacity, which can be overcome with a proper functionalization of the nanoparticle surface. In this study, we present, for the first time, a computational approach based on metadynamics to estimate the binding free energy of the doxorubicin drug (DOX) to a functionalized TiO₂ nanoparticle under different pH conditions. On a thermodynamic basis, we demonstrate the robustness of our approach to capture the overall mechanism behind the pH-triggered release of DOX due to environmental pH changes. Notably, binding free energy estimations align well with what is expected for a pH-sensitive drug delivery system. Based on our results, we envision the use of metadynamics as a promising computational tool for the rational design and in silico optimization of organic ligands with improved drug carrier properties.