Exploring the Properties of Curved Lipid Membranes: Comparative Analysis of Atomistic and Coarse-Grained Force Fields.

Curvature emerges as a fundamental membrane characteristic crucial for diverse biological processes, including vesicle formation, cell signaling, and membrane trafficking. Increasingly valuable insights into atomistic details governing curvature-dependent membrane properties are provided by computer simulations. Nevertheless, the underlying force field models are conventionally calibrated and tested in relation to experimentally derived parameters of planar bilayers, thereby leaving uncertainties concerning their consistency in reproducing curved lipid systems.

Dehydration of Lipid Membranes Drives Redistribution of Cholesterol Between Lateral Domains.

Cholesterol-rich lipid rafts are found to facilitate membrane fusion, central to processes like viral entry, fertilization, and neurotransmitter release. While the fusion process involves local, transient membrane dehydration, the impact of reduced hydration on cholesterol's structural organization in biological membranes remains unclear. Here, we employ confocal fluorescence microscopy and atomistic molecular dynamics simulations to investigate cholesterol behavior in phase-separated lipid bilayers under controlled hydration.

Molecular Mechanisms behind Conformational Transitions of the Influenza Virus Hemagglutinin Membrane Anchor.

Membrane fusion is a fundamental process that is exploited by enveloped viruses to enter host cells. In the case of the influenza virus, fusion is facilitated by the trimeric viral hemagglutinin protein (HA). So far, major focus has been put on its N-terminal fusion peptides, which are directly responsible for fusion initiation.

Two modes of fusogenic action for influenza virus fusion peptide.

The entry of influenza virus into the host cell requires fusion of its lipid envelope with the host membrane. It is catalysed by viral hemagglutinin protein, whose fragments called fusion peptides become inserted into the target bilayer and initiate its merging with the viral membrane. Isolated fusion peptides are already capable of inducing lipid mixing between liposomes.

Discovering highly potent antimicrobial peptides with deep generative model HydrAMP.

Antimicrobial peptides emerge as compounds that can alleviate the global health hazard of antimicrobial resistance, prompting a need for novel computational approaches to peptide generation. Here, we propose HydrAMP, a conditional variational autoencoder that learns lower-dimensional, continuous representation of peptides and captures their antimicrobial properties. The model disentangles the learnt representation of a peptide from its antimicrobial conditions and leverages parameter-controlled creativity.

Identification of CB1 Ligands among Drugs, Phytochemicals and Natural-Like Compounds: Virtual Screening and In Vitro Verification.

Cannabinoid receptor type 1 (CB1) is an important modulator of many key physiological functions and thus a compelling molecular target. However, safe CB1 targeting is a non-trivial task. In recent years, there has been a surge of data indicating that drugs successfully used in the clinic for years (e.

Granger Causality Analysis of Chignolin Folding.

Constantly advancing computer simulations of biomolecules provide huge amounts of data that are difficult to interpret. In particular, obtaining insights into functional aspects of macromolecular dynamics, often related to cascades of transient events, calls for methodologies that depart from the well-grounded framework of equilibrium statistical physics. One of the approaches toward the analysis of complex temporal data which has found applications in the fields of neuroscience and econometrics is Granger causality analysis.

Membrane-Bound Configuration and Lipid Perturbing Effects of Hemagglutinin Subunit 2 N-Terminus Investigated by Computer Simulations.

Hemagglutinin (HA) mediated fusion of influenza virus envelope with host lipid membrane is a critical step warrantying virus entry to the cell. Despite tremendous advances in structural biology methods, the knowledge concerning the details of HA2 subunit insertion into the target membrane and its subsequent bilayer perturbing effect is still rather limited. Herein, based on a set of molecular dynamics simulations, we investigate the structure and interaction with lipid membrane of the N-terminal HA2 region comprising a trimer of fusion peptides (HAfps) tethered by flexible linkers to a fragment of coiled-coil stem structure.

Entropy-based distance cutoff for protein internal contact networks.

Protein structure networks (PSNs) have long been used to provide a coarse yet meaningful representation of protein structure, dynamics, and internal communication pathways. An important question is what criteria should be applied to construct the network so that to include relevant interresidue contacts while avoiding unnecessary connections. To address this issue, we systematically considered varying residue distance cutoff length and the probability threshold for contact formation to construct PSNs based on atomistic molecular dynamics in order to assess the amount of mutual information within the resulting representations.

GridSolvate: A Web Server for the Prediction of Biomolecular Hydration Properties.

We present a novel web server, named gridSolvate, dedicated to the prediction of biomolecular hydration properties. Given a solute in atomic representation, such as a protein or protein-ligand complex, the server determines positions and excess chemical potential of buried and first hydration shell water molecules. Calculations are based on our semiexplicit hydration model that provides computational efficiency close to implicit solvent approaches, yet captures a number of physical effects unique to explicit solvent representation.

Conserved internal hydration motifs in protein kinases.

Crystal structures of diverse protein kinase catalytic subunits reveal a number of water molecules whose positions within the protein core are better preserved than amino acid types in many functionally important locations. It remains unknown whether they play any particular role, and whether their removal, disturbing local interaction patterns to no smaller degree than amino acid mutations, can affect kinase stability and function. In this study, we apply an array of computational approaches to characterize hydration of kinase catalytic subunits.

Restriction of S-adenosylmethionine conformational freedom by knotted protein binding sites.

S-adenosylmethionine (SAM) is one of the most important enzyme substrates. It is vital for the function of various proteins, including large group of methyltransferases (MTs). Intriguingly, some bacterial and eukaryotic MTs, while catalysing the same reaction, possess significantly different topologies, with the former being a knotted one.

Water-mediated conformational preselection mechanism in substrate binding cooperativity to protein kinase A.

Substrate binding cooperativity in protein kinase A (PKA) seems to involve allosteric coupling between the two binding sites. It received significant attention, but its molecular basis still remains not entirely clear. Based on long molecular dynamics of PKA and its complexes, we characterized an allosteric pathway that links ATP binding to the redistribution of states adopted by a protein substrate positioning segment in favor of those that warrant correct binding.

Explicit Solvent Hydration Benchmark for Proteins with Application to the PBSA Method.

Explicit and implicit solvent models have a proven record of delivering hydration free energies of small, druglike solutes in reasonable agreement with experiment. Hydration of macromolecules, such as proteins, is to a large extent uncharted territory, with few results shedding light on quantitative consistency between different solvent models, let alone their ability to reproduce real water. In this work, based on extensive explicit solvent simulations employing TIP3P and SPC/E water models we analyze hydration free energy changes between fixed conformations of 5 diverse proteins, including large multidomain structures.

Principles for Tuning Hydrophobic Ligand-Receptor Binding Kinetics.

We investigate how to tune the rate of hydrophobic ligand-receptor association due to the role of solvent in adjustable receptor pockets by explicit-water molecular dynamics (MD) simulations. Our model considers the binding of a spherical ligand (key/guest) to a concave surface recess in a nonpolar wall as receptor (lock/host). We systematically modify the receptor's physicochemical properties in terms of geometry and dispersion attraction which, in turn, alter the water occupancy and fluctuations within the pocket.

Solvent Fluctuations Induce Non-Markovian Kinetics in Hydrophobic Pocket-Ligand Binding.

We investigate the impact of water fluctuations on the key-lock association kinetics of a hydrophobic ligand (key) binding to a hydrophobic pocket (lock) by means of a minimalistic stochastic model system. The latter describes the collective hydration behavior of the pocket by bimodal fluctuations of the water-pocket interface that dynamically couples to the diffusive motion of the approaching ligand via the hydrophobic interaction. This leads to a Markovian set of overdamped stochastic equations in 2D-coordinate-space spanned by the interface position and the ligand position.

Prediction of Water Binding to Protein Hydration Sites with a Discrete, Semiexplicit Solvent Model.

Buried water molecules are ubiquitous in protein structures and are found at the interface of most protein-ligand complexes. Determining their distribution and thermodynamic effect is a challenging yet important task, of great of practical value for the modeling of biomolecular structures and their interactions. In this study, we present a novel method aimed at the prediction of buried water molecules in protein structures and estimation of their binding free energies.

Hydration in discrete water (II): from neutral to charged solutes.

In our previous work, we introduced a solvation model based on discrete solvent representation and demonstrated its ability to estimate hydration free energies for neutral solutes. Here, we present modifications extending the applicability of the model to charged solutes. They include improvements in the representation of the first hydration shell and systematic treatment of long-range interactions.

Elastic Network Models of Nucleic Acids Flexibility.

Elastic network models (ENMs) are a useful tool for describing large scale motions in protein systems. While they are well validated in the context of proteins, relatively little is known about their applicability to nucleic acids, whose different architecture does not necessarily warrant comparable performance. In this study we thoroughly evaluate and optimize the efficiency of popular ENMs for capturing RNA and DNA flexibility.