scholarly journals Thermodynamics and free energy landscape of BAR-domain dimerization from molecular simulations

2020 ◽  
Author(s):  
Adip Jhaveri ◽  
Dhruw Maisuria ◽  
Matthew Varga ◽  
Dariush Mohammadyani ◽  
Margaret E Johnson
Keyword(s):  
Free Energy ◽  
Energy Surface ◽  
Md Simulations ◽  
Coarse Grained ◽  
Enhanced Sampling ◽  
Free Energies ◽  
Free Energy Surface ◽  
Bar Domain ◽  
Bar Domains ◽  
Protein Models

AbstractNearly all proteins interact specifically with other proteins, often forming reversible bound structures whose stability is critical to function. Proteins with BAR domains function to bind to, bend, and remodel biological membranes, where the dimerization of BAR domains is a key step in this function. Here we characterize the binding thermodynamics of homodimerization between the LSP1 BAR domain proteins in solution, using Molecular Dynamics (MD) simulations. By combining the MARTINI coarse-grained protein models with enhanced sampling through metadynamics, we construct a two-dimensional free energy surface quantifying the bound versus unbound ensembles as a function of two distance variables. Our simulations portray a heterogeneous and extraordinarily stable bound ensemble for these modeled LSP1 proteins. The proper crystal structure dimer has a large hydrophobic interface that is part of a stable minima on the free energy surface, which is enthalpically the minima of all bound structures. However, we also find several other stable nonspecific dimers with comparable free energies to the specific dimer. Through structure-based clustering of these bound structures, we find that some of these ‘nonspecific’ contacts involve extended tail regions that help stabilize the higher-order oligomers formed by BAR-domains, contacts that are separated from the homodimer interface. We find that the known membrane-binding residues of the LSP1 proteins rarely participate in any of the bound interfaces, but that both patches of residues are aligned to interact with the membrane in the specific dimer. Hence, we would expect a strong selection of the specific dimer in binding to the membrane. The effect of a 100mM NaCl buffer reduces the relative stability of nonspecific dimers compared to the specific dimer, indicating that it would help prevent aggregation of the proteins. With these results, we provide the first free energy characterization of interaction pathways in this important class of membrane sculpting domains, revealing a variety of interfacial contacts outside of the specific dimer that may help stabilize its oligomeric assemblies.

scholarly journals THE STRUCTURE OF THE FREE ENERGY SURFACE OF COARSE-GRAINED OFF-LATTICE PROTEIN MODELS

2007 ◽  
Vol 18 (01) ◽  
pp. 99-106 ◽  
Author(s):  
ETHEM AKTÜRK ◽  
HANDAN ARKIN ◽  
TARIK ÇELİK
Keyword(s):  
Energy Surface ◽  
Energy Landscape ◽  
Low Energy ◽  
Coarse Grained ◽  
Free Energy Surface ◽  
Protein Model ◽  
Energy Part ◽  
Folding Pathways ◽  
Protein Models ◽  
Lattice Protein Models

We have performed multicanonical simulations of hydrophobic-hydrophilic heteropolymers with a simple effective, coarse-grained off-lattice model to study the structure and the topology of the energy surface. The multicanonical method samples the whole rugged energy landscape, in particular the low-energy part, and enables one to better understand the critical behaviors and visualize the folding pathways of the considered protein model.

scholarly journals Locating landmarks on high-dimensional free energy surfaces

2015 ◽  
Vol 112 (11) ◽  
pp. 3235-3240 ◽  
Author(s):  
Ming Chen ◽  
Tang-Qing Yu ◽  
Mark E. Tuckerman
Keyword(s):  
Free Energy ◽  
Degrees Of Freedom ◽  
Energy Surface ◽  
Saddle Points ◽  
Coarse Graining ◽  
Graph Representation ◽  
Coarse Grained ◽  
High Dimensional ◽  
Free Energy Surface ◽  
Energy Surfaces

Coarse graining of complex systems possessing many degrees of freedom can often be a useful approach for analyzing and understanding key features of these systems in terms of just a few variables. The relevant energy landscape in a coarse-grained description is the free energy surface as a function of the coarse-grained variables, which, despite the dimensional reduction, can still be an object of high dimension. Consequently, navigating and exploring this high-dimensional free energy surface is a nontrivial task. In this paper, we use techniques from multiscale modeling, stochastic optimization, and machine learning to devise a strategy for locating minima and saddle points (termed “landmarks”) on a high-dimensional free energy surface “on the fly” and without requiring prior knowledge of or an explicit form for the surface. In addition, we propose a compact graph representation of the landmarks and connections between them, and we show that the graph nodes can be subsequently analyzed and clustered based on key attributes that elucidate important properties of the system. Finally, we show that knowledge of landmark locations allows for the efficient determination of their relative free energies via enhanced sampling techniques.

scholarly journals Exploration of Free Energy Surface and Thermal Effects on Relative Population and Infrared Spectrum of the Be6B11− Fluxional Cluster

Materials ◽  
2020 ◽  
Vol 14 (1) ◽  
pp. 112
Author(s):  
Carlos Emiliano Buelna-Garcia ◽  
José Luis Cabellos ◽  
Jesus Manuel Quiroz-Castillo ◽  
Gerardo Martinez-Guajardo ◽  
Cesar Castillo-Quevedo ◽  
...  
Keyword(s):  
Free Energy ◽  
Infrared Spectrum ◽  
Infrared Spectra ◽  
Thermal Effects ◽  
Global Minimum ◽  
Energy Surface ◽  
Weighted Sums ◽  
Free Energy Surface ◽  
Temperature Dependent ◽  
Relative Population

The starting point to understanding cluster properties is the putative global minimum and all the nearby local energy minima; however, locating them is computationally expensive and difficult. The relative populations and spectroscopic properties that are a function of temperature can be approximately computed by employing statistical thermodynamics. Here, we investigate entropy-driven isomers distribution on Be6B11− clusters and the effect of temperature on their infrared spectroscopy and relative populations. We identify the vibration modes possessed by the cluster that significantly contribute to the zero-point energy. A couple of steps are considered for computing the temperature-dependent relative population: First, using a genetic algorithm coupled to density functional theory, we performed an extensive and systematic exploration of the potential/free energy surface of Be6B11− clusters to locate the putative global minimum and elucidate the low-energy structures. Second, the relative populations’ temperature effects are determined by considering the thermodynamic properties and Boltzmann factors. The temperature-dependent relative populations show that the entropies and temperature are essential for determining the global minimum. We compute the temperature-dependent total infrared spectra employing the Boltzmann factor weighted sums of each isomer’s infrared spectrum and find that at finite temperature, the total infrared spectrum is composed of an admixture of infrared spectra that corresponds to the spectra of the lowest-energy structure and its isomers located at higher energies. The methodology and results describe the thermal effects in the relative population and the infrared spectra.

scholarly journals The Influence of Starch Origin on the Properties of Starch Films: Packaging Performance

Materials ◽  
2021 ◽  
Vol 14 (5) ◽  
pp. 1146
Author(s):  
Zuzanna Żołek-Tryznowska ◽  
Alicja Kałuża
Keyword(s):  
Tensile Strength ◽  
Free Energy ◽  
Surface Free Energy ◽  
Food Packaging ◽  
Energy Surface ◽  
Free Energy Surface ◽  
Casting Method ◽  
Starch Solution ◽  
Film Properties ◽  
Starch Films

Starch films can be used as materials for food packaging purposes. The goal of this study is to compare how the starch origin influence the selected starch film properties. The films were made from various starches such as that from maize, potato, oat, rice, and tapioca using 50%w of glycerine as a plasticizer. The obtained starch-based films were made using the well-known casting method from a starch solution in water. The properties of the films that were evaluated were tensile strength, water vapour transition rate, moisture content, wettability, and their surface free energy. Surface free energy (SFE) and its polar and dispersive components were calculated using the Owens-Wendt-Rabel-Kaelbe approach. The values of SFE in the range of 51.64 to 70.81 mJ∙m−2 for the oat starch-based film and the maize starch-based film. The films revealed worse mechanical properties than those of conventional plastics for packaging purposes. The results indicated that the poorest tensile strength was exhibited by the starch-based films made from oat (0.36 MPa) and tapioca (0.78 MPa) and the greatest tensile strength (1.49 MPa) from potato.

Efficient approach to include molecular polarizations using charge and atom dipole response kernels to calculate free energy gradients in the QM/MM scheme

2015 ◽  
Vol 17 (40) ◽  
pp. 26955-26968 ◽  
Author(s):  
Toshio Asada ◽  
Kanta Ando ◽  
Koji Sakurai ◽  
Shiro Koseki ◽  
Masataka Nagaoka
Keyword(s):  
Free Energy ◽  
Energy Surface ◽  
Free Energy Surface ◽  
Efficient Approach ◽  
Reaction Processes ◽  
Dipole Response

An efficient approach to evaluate free energy gradients within the QM/MM framework has been proposed to clarify reaction processes on the free energy surface.

scholarly journals Free energy surface of initial cap formation in carbon nanotube growth

2021 ◽  
Author(s):  
Satoru Fukuhara ◽  
Yasushi Shibuta
Keyword(s):  
Carbon Nanotubes ◽  
Free Energy ◽  
Carbon Nanotube ◽  
Energy Surface ◽  
Catalyst Surface ◽  
Important Process ◽  
Free Energy Surface ◽  
Carbon Nanotube Growth ◽  
Nanotube Growth ◽  
Carbon Network

Initial cap formation is an important process of carbon nanotubes (CNTs) growth where hexagonal carbon network is lifted off from the catalyst surface. In this study, free energy surface (FES)...

scholarly journals Potential Drug Candidates for SARS-CoV-2 Using Computational Screening and Enhanced Sampling Methods

2020 ◽  
Author(s):  
Sadanandam Namsani ◽  
Debabrata Pramanik ◽  
Mohd Aamir Khan ◽  
Sudip Roy ◽  
Jayant Singh
Keyword(s):  
Free Energy ◽  
High Performance ◽  
Search Space ◽  
Computational Method ◽  
Enhanced Sampling ◽  
Free Energy Surface ◽  
Drug Candidates ◽  
Drug Molecules ◽  
Computational Screening ◽  
Graphical Processing

<div><div><div><p>Here we report new chemical entities that are highly specific in binding towards the 3-chymotrypsin- like cysteine protease (3CLpro) protein present in the novel SARS-CoV2 virus. The viral 3CLpro</p><p>protein controls coronavirus replication. Therefore, 3CLpro is identified as a target for drug molecules. We have implemented an enhanced sampling method in combination with molecular dynamics and docking to bring down the computational screening search space to four molecules that could be synthesised and tested against COVID-19. Our computational method is much more robust than any other method available for drug screening e.g., docking, because of sampling of the free energy surface of the binding site of the protein (including the ligand) and use of explicit solvent. We have considered all possible interactions between all the atoms present in the protein, ligands, and water. Using high performance computing with graphical processing units we are able to perform large number of simulations within a month's time and converge to 4 most strongly bound ligands (by free energy and other scores) from a set of 17 ligands with lower docking scores. Based on our results and analysis, we claim with high confidence, that we have identified four potential ligands. Out of those, one particular ligand is the most promising candidate, based on free energy data, for further synthesis and testing against SARS-CoV-2 and might be effective for the cure of COVID-19.</p></div></div></div>

scholarly journals Residue-Level Contact Reveals Modular Domain Interactions of PICK1 Are Driven by Both Electrostatic and Hydrophobic Forces

2021 ◽  
Vol 7 ◽  
Author(s):  
Amy O. Stevens ◽  
Yi He
Keyword(s):  
Hydrophobic Interactions ◽  
Pdz Domain ◽  
Md Simulations ◽  
Intramolecular Interactions ◽  
Coarse Grained ◽  
Scaffolding Protein ◽  
Concave Surface ◽  
Stable Complex ◽  
Bar Domain ◽  
C Terminus

PICK1 is a multi-domain scaffolding protein that is uniquely comprised of both a PDZ domain and a BAR domain. While previous experiments have shown that the PDZ domain and the linker positively regulate the BAR domain and the C-terminus negatively regulates the BAR domain, the details of internal regulation mechanisms are unknown. Molecular dynamics (MD) simulations have been proven to be a useful tool in revealing the intramolecular interactions at atomic-level resolution. PICK1 performs its biological functions in a dimeric form which is extremely computationally demanding to simulate with an all-atom force field. Here, we use coarse-grained MD simulations to expose the key residues and driving forces in the internal regulations of PICK1. While the PDZ and BAR domains do not form a stable complex, our simulations show the PDZ domain preferentially interacting with the concave surface of the BAR domain over other BAR domain regions. Furthermore, our simulations show that the short helix in the linker region can form interactions with the PDZ domain. Our results reveal that the surface of the βB-βC loop, βC strand, and αA-βD loop of the PDZ domain can form a group of hydrophobic interactions surrounding the linker helix. These interactions are driven by hydrophobic forces. In contrast, our simulations reveal a very dynamic C-terminus that most often resides on the convex surface of the BAR domain rather than the previously suspected concave surface. These interactions are driven by a combination of electrostatic and hydrophobic interactions.

Sign in / Sign up

Export Citation Format

Share Document