Predicting the rate of nonfacilitated permeation of solutes across lipid bilayers is important to drug design, toxicology, and signaling. These rates can be estimated using molecular dynamics simulations combined with the inhomogeneous solubility-diffusion model, which requires calculation of the potential of mean force and position-dependent diffusivity of the solute along the transmembrane axis. In this paper, we assess the efficiency and accuracy of several methods for the calculation of the permeability of a model DMPC bilayer to urea, benzoic acid, and codeine. We compare umbrella sampling, replica exchange umbrella sampling, adaptive biasing forces, and multiple-walker adaptive biasing forces for the calculation of the transmembrane PMF. No definitive advantage for any of these methods in their ability to predict the permeability coefficient Pm was found, provided that a sufficiently long equilibration is performed. For diffusivities, a Bayesian inference method was compared to a Generalized Langevin method, both being sensitive to chosen parameters and the slow relaxation of membrane defects. Agreement within 1.5 log units of the computed Pm with experiment is found for all permeants and methods. Remaining discrepancies can likely be attributed to limitations of the force field as well as slowly relaxing collective movements within the lipid environment. Numerical calculations based on model profiles show that Pm can be reliably estimated from only a few data points, leading to recommendations for calculating Pm from simulations.

Liquid state Overhauser Effect Dynamic Nuclear Polarization (ODNP) has experienced a recent resurgence of interest. The ODNP technique described here relies on the double resonance of electron spin resonance (ESR) at the most common, i.e. X-band (~ 10 GHz), frequency and 1H nuclear magnetic resonance (NMR) at ~ 15 MHz. It requires only a standard continuous wave (cw) ESR spectrometer with an NMR probe inserted or built into an X-band cavity. Our focus lies on reviewing a new and powerful manifestation of ODNP as a high frequency NMR relaxometry tool that probes dipolar cross relaxation between the electron spins and the 1H nuclear spins at X-band frequencies. This technique selectively measures the translational mobility of water within a volume extending 0.5–1.5 nm outward from a nitroxide radical spin probe that is attached to a targeted site of a macromolecule. This method has been applied to study the dynamics of water that hydrates or permeates the surface or interior of proteins, polymers, and lipid membrane vesicles. We begin by reviewing the recent advances that have helped develop ODNP into a tool for mapping the dynamic landscape of hydration water with sub-nanometer locality. In order to bind this work coherently together, and to place it in the context of the extensive body of research in the field of NMR relaxometry, we then rephrase the analytical model and extend the description of the ODNP-derived NMR signal enhancements. This extended model highlights several aspects of ODNP data analysis, including the importance of considering all possible effects of microwave sample heating, the need to consider the error associated with various relaxation rates, and the unique ability of ODNP to probe the electron–1H cross-relaxation process, which is uniquely sensitive to fast (tens of ps) dynamical processes. By implementing the relevant corrections in a stepwise fashion, this paper draws a consensus result from previous ODNP procedures, and then shows how such data can be further corrected to yield clear and reproducible saturation of the NMR hyperpolarization process. Finally, drawing on these results, we broadly survey the previous ODNP dynamics literature. We find that the vast number of published, empirical hydration dynamics data can be reproducibly classified into regimes of surface, interfacial, vs buried water dynamics.

We present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.

aSilica oligomerization is the key reaction in zeolite synthesis. NaOH is a common additive in the zeolite synthesis that decreases the reaction rate of smaller silica oligomers and also affects the final structure of the zeolite. Here we report a study of the role of sodium in the initial stages of silica oligomerization. We performed ab initio molecular dynamics simulations using ecplicit aqueous solution in order to obtain the free energy profile and study the behavior of sodium during the reaction. Our study confirms that sodium decreases the reaction rates of oligomerization for smaller silica chains.Analysis of the molecular dynamics trajectories shows that sodium does not increase the reaction barriers by direct coordination to the silica. However, sodium is often present in the second solvation shell of the reactive atoms. Correlation between sodium presence in the first or the second shell of the reactive oxygen and a decrease in hydrogen bonding for that oxygen was found for the first reaction step. Therefore, the presence of sodium could contribute to an increase in reaction barriers for silica oligomerization by some rearrangement of the hydrogen bond network of water solution around the reactants.

Protein aggregation plays a critical role in the pathogenesis of neurodegenerative diseases, and the mechanism of its progression is poorly understood. Here, we examine the structural and dynamic characteristics of transiently evolving protein aggregates under ambient conditions by directly probing protein surface water diffusivity, local protein segment dynamics, and interprotein packing as a function of aggregation time, along the third repeat domain and C terminus of Δtau187 spanning residues 255-441 of the longest isoform of human tau. These measurements were achieved with a set of highly sensitive magnetic resonance tools that rely on sitespecific electron spin labeling of Δtau187. Within minutes of initiated aggregation, the majority of Δtau187 that is initially homogeneously hydrated undergoes structural transformations to form partially structured aggregation intermediates. This is reflected in the dispersion of surface water dynamics that is distinct around the third repeat domain, found to be embedded in an intertau interface, from that of the solvent-exposed C terminus. Over the course of hours and in a rate-limiting process, a majority of these aggregation intermediates proceed to convert into stable β-sheet structured species and maintain their stacking order without exchanging their subunits. The population of β-sheet structured species is >5% within 5 min of aggregation and gradually grows to 50-70% within the early stages of fibril formation, while they mostly anneal block-wisely to form elongated fibrils. Our findings suggest that the formation of dynamic aggregation intermediates constitutes a major event occurring in the earliest stages of tau aggregation that precedes, and likely facilitates, fibril formation and growth.biological water | soluble oligomers | amyloid formation | site-directed spin labeling | Alzheimer's disease

The time step of atomistic molecular dynamics (MD) simulations is determined by the fastest motions in the system and is typically limited to 2 fs. An increasingly popular approach is to increase the mass of the hydrogen atoms to ∼3 amu and decrease the mass of the parent atom by an equivalent amount. This approach, known as hydrogen-mass repartitioning (HMR), permits time steps up to 4 fs with reasonable simulation stability. While HMR has been applied in many published studies to date, it has not been extensively tested for membrane-containing systems. Here, we compare the results of simulations of a variety of membranes and membrane–protein systems run using a 2 fs time step and a 4 fs time step with HMR. For pure membrane systems, we find almost no difference in structural properties, such as area-per-lipid, electron density profiles, and order parameters, although there are differences in kinetic properties such as the diffusion constant. Conductance through a porin in an applied field, partitioning of a small peptide, hydrogen-bond dynamics, and membrane mixing show very little dependence on HMR and the time step. We also tested a 9 Å cutoff as compared to the standard CHARMM cutoff of 12 Å, finding significant deviations in many properties tested. We conclude that HMR is a valid approach for membrane systems, but a 9 Å cutoff is not.

No abstract

No abstract