Optimal Dielectric Boundary for Binding Free Energy Estimates in the Implicit Solvent.

Accuracy of binding free energy calculations utilizing implicit solvent models is critically affected by parameters of the underlying dielectric boundary, specifically, the atomic and water probe radii. Here, a multidimensional optimization pipeline is used to find optimal atomic radii, specifically for binding calculations in the implicit solvent. To reduce overfitting, the optimization target includes separate, weighted contributions from both binding and hydration free energies.

Inclusion of Water Multipoles into the Implicit Solvation Framework Leads to Accuracy Gains.

The current practical "workhorses" of the atomistic implicit solvation─the Poisson-Boltzmann (PB) and generalized Born (GB) models─face fundamental accuracy limitations. Here, we propose a computationally efficient implicit solvation framework, the Implicit Water Multipole GB (IWM-GB) model, that systematically incorporates the effects of multipole moments of water molecules in the first hydration shell of a solute, beyond the dipole water polarization already present at the PB/GB level. The framework explicitly accounts for coupling between polar and nonpolar contributions to the total solvation energy, which is missing from many implicit solvation models.

The probability of chromatin to be at the nuclear lamina has no systematic effect on its transcription level in fruit flies.

Background: Multiple studies have demonstrated a negative correlation between gene expression and positioning of genes at the nuclear envelope (NE) lined by nuclear lamina, but the exact relationship remains unclear, especially in light of the highly stochastic, transient nature of the gene association with the NE.

Results: In this paper, we ask whether there is a causal, systematic, genome-wide relationship between the expression levels of the groups of genes in topologically associating domains (TADs) of Drosophila nuclei and the probabilities of TADs to be found at the NE. To investigate the nature of this possible relationship, we combine a coarse-grained dynamic model of the entire Drosophila nucleus with genome-wide gene expression data; we analyze the TAD averaged transcription levels of genes against the probabilities of individual TADs to be in contact with the NE in the control and lamins-depleted nuclei.

Improving the Accuracy of Physics-Based Hydration-Free Energy Predictions by Machine Learning the Remaining Error Relative to the Experiment.

The accuracy of computational models of water is key to atomistic simulations of biomolecules. We propose a computationally efficient way to improve the accuracy of the prediction of hydration-free energies (HFEs) of small molecules: the remaining errors of the physics-based models relative to the experiment are predicted and mitigated by machine learning (ML) as a postprocessing step. Specifically, the trained graph convolutional neural network attempts to identify the "blind spots" in the physics-based model predictions, where the complex physics of aqueous solvation is poorly accounted for, and partially corrects for them.

Strong interactions between highly dynamic lamina-associated domains and the nuclear envelope stabilize the 3D architecture of Drosophila interphase chromatin.

Background: Interactions among topologically associating domains (TADs), and between the nuclear envelope (NE) and lamina-associated domains (LADs) are expected to shape various aspects of three-dimensional (3D) chromatin structure and dynamics; however, relevant genome-wide experiments that may provide statistically significant conclusions remain difficult.

Results: We have developed a coarse-grained dynamical model of D. melanogaster nuclei at TAD resolution that explicitly accounts for four distinct epigenetic classes of TADs and LAD-NE interactions.

Explicit ions/implicit water generalized Born model for nucleic acids.

The ion atmosphere around highly charged nucleic acid molecules plays a significant role in their dynamics, structure, and interactions. Here we utilized the implicit solvent framework to develop a model for the explicit treatment of ions interacting with nucleic acid molecules. The proposed explicit ions/implicit water model is based on a significantly modified generalized Born (GB) model and utilizes a non-standard approach to define the solute/solvent dielectric boundary.

Spermine Condenses DNA, but Not RNA Duplexes.

Interactions between the polyamine spermine and nucleic acids drive important cellular processes. Spermine condenses DNA and some RNAs, such as poly(rA):poly(rU). A large fraction of the spermine present in cells is bound to RNA but apparently does not condense it.

Opposing Effects of Multivalent Ions on the Flexibility of DNA and RNA.

Increasing the concentration of counterions (salt) is known to reduce the bending persistence length of DNA. Here we use atomistic molecular dynamics simulations to predict that multivalent counterions have the opposite effect on double-stranded RNA, increasing its bending rigidity by at least 30%. This counterintuitive effect is observed for various tri- and tetravalent ions alike, and is robust to methodological details and the RNA sequence.

Multi-shell model of ion-induced nucleic acid condensation.

We present a semi-quantitative model of condensation of short nucleic acid (NA) duplexes induced by trivalent cobalt(iii) hexammine (CoHex) ions. The model is based on partitioning of bound counterion distribution around single NA duplex into "external" and "internal" ion binding shells distinguished by the proximity to duplex helical axis. In the aggregated phase the shells overlap, which leads to significantly increased attraction of CoHex ions in these overlaps with the neighboring duplexes.

Understanding nucleic acid structural changes by comparing wide-angle x-ray scattering (WAXS) experiments to molecular dynamics simulations.

Wide-angle x-ray scattering (WAXS) is emerging as a powerful tool for increasing the resolution of solution structure measurements of biomolecules. Compared to its better known complement, small angle x-ray scattering (SAXS), WAXS targets higher scattering angles and can enhance structural studies of molecules by accessing finer details of solution structures. Although the extension from SAXS to WAXS is easy to implement experimentally, the computational tools required to fully harness the power of WAXS are still under development.

Accurate evaluation of charge asymmetry in aqueous solvation.

Charge hydration asymmetry (CHA)-a characteristic dependence of hydration free energy on the sign of the solute charge-quantifies the asymmetric response of water to electric field at microscopic level. Accurate estimates of CHA are critical for understanding hydration effects ubiquitous in chemistry and biology. However, measuring hydration energies of charged species is fraught with significant difficulties, which lead to unacceptably large (up to 300%) variation in the available estimates of the CHA effect.

Why double-stranded RNA resists condensation.

The addition of small amounts of multivalent cations to solutions containing double-stranded DNA leads to inter-DNA attraction and eventual condensation. Surprisingly, the condensation is suppressed in double-stranded RNA, which carries the same negative charge as DNA, but assumes a different double helical form. Here, we combine experiment and atomistic simulations to propose a mechanism that explains the variations in condensation of short (25 base-pairs) nucleic acid (NA) duplexes, from B-like form of homopolymeric DNA, to mixed sequence DNA, to DNA:RNA hybrid, to A-like RNA.

Introducing Charge Hydration Asymmetry into the Generalized Born Model.

The effect of charge hydration asymmetry (CHA)-non-invariance of solvation free energy upon solute charge inversion-is missing from the standard linear response continuum electrostatics. The proposed charge hydration asymmetric-generalized Born (CHA-GB) approximation introduces this effect into the popular generalized Born (GB) model. The CHA is added to the GB equation via an analytical correction that quantifies the specific propensity of CHA of a given water model; the latter is determined by the charge distribution within the water model.

Charge hydration asymmetry: the basic principle and how to use it to test and improve water models.

Charge hydration asymmetry (CHA) manifests itself in the experimentally observed strong dependence of free energy of ion hydration on the sign of the ion charge. This asymmetry is not consistently accounted for by popular models of solvation; its magnitude varies greatly between the models. While it is clear that CHA is somehow related to charge distribution within a water molecule, the exact nature of this relationship is unknown.

Binding free energy and counterion release for adsorption of the antimicrobial peptide lactoferricin B on a POPG membrane.

Molecular dynamics (MD) simulations are used to study the interaction of an anionic palmitoyl-oleoyl-phosphatidylglycerol (POPG) bilayer with the cationic antimicrobial peptide bovine lactoferricin (LFCinB) in a 100 mM NaCl solution at 310 K. The interaction of LFCinB with a POPG bilayer is employed as a model system for studying the details of membrane adsorption selectivity of cationic antimicrobial peptides. Seventy eight 4 ns MD production run trajectories of the equilibrated system, with six restrained orientations of LFCinB at 13 different separations from the POPG membrane, are generated to determine the free energy profile for the peptide as a function of the distance between LFCinB and the membrane surface.

Prediction of binding free energy for adsorption of antimicrobial peptide lactoferricin B on a POPC membrane.

Molecular dynamics (MD) simulations are used to study the interaction of a zwitterionic palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayer with the cationic antimicrobial peptide bovine lactoferricin (LFCinB) in a 100 mM NaCl solution at 310 K. The interaction of LFCinB with POPC is used as a model system for studying the details of membrane-peptide interactions, with the peptide selected because of its antimicrobial nature. Seventy-two 3 ns MD simulations, with six orientations of LFCinB at 12 different distances from a POPC membrane, are carried out to determine the potential of mean force (PMF) or free energy profile for the peptide as a function of the distance between LFCinB and the membrane surface.

Nonequilibrium molecular dynamics calculation of the conductance of the KcsA potassium ion channel.

We have developed an efficient method for the atomistic determination of the conductance of a biological ion channel model by applying an external field to the conducting ions only. The underlying theory is discussed and demonstrated in a simple test system consisting of two ions in a box of water. Finally, the theory is applied to the experimentally determined structure of the KcsA potassium channel from which a conductance in reasonable agreement with the experimental result is predicted.

Unified modeling of conductance kinetics for low- and high-conductance potassium ion channels.

A kinetics model is proposed for the description of ion conductance of low- and high-conductance potassium ion channels. The model describes ion permeation through the selectivity filter, which is assumed to be the only conductance determining part of the open channel. The filter occupancy can vary from zero to three ions, affecting the ion entry and exit rates.

Conduction through the inward rectifier potassium channel, Kir2.1, is increased by negatively charged extracellular residues.

Ion channel conductance can be influenced by electrostatic effects originating from fixed "surface" charges that are remote from the selectivity filter. To explore whether surface charges contribute to the conductance properties of Kir2.1 channels, unitary conductance was measured in cell-attached recordings of Chinese hamster ovary (CHO) cells transfected with Kir2.

Non-Michaelis-Menten kinetics model for conductance of low-conductance potassium ion channels.

A reduced kinetics model is proposed for ion permeation in low-conductance potassium ion channels with zero net electrical charge in the selectivity filter region. The selectivity filter is assumed to be the only conductance-determining part of the channel. Ion entry and exit rate constants depend on the occupancy of the filter due to ion-ion interactions.