Entropy-based methods for formulating bottom-up ultra-coarse-grained models.

Bottom-up coarse-grained (CG) modeling is an effective means of bypassing the limited spatiotemporal scales of conventional atomistic molecular dynamics while retaining essential information from the atomistic model. A central challenge in CG modeling is the trade-off between accuracy and efficiency, as the inclusion of often pivotal many-body interaction terms in the CG force-field renders simulation markedly slower than simple pairwise models. The Ultra Coarse-Graining (UCG) method incorporates many-body terms through discrete internal state variables that modulate the CG force-field according to, e.

Data-Driven Equation-Free Dynamics Applied to Many-Protein Complexes: The Microtubule Tip Relaxation.

Microtubules (MTs) constitute the largest components of the eukaryotic cytoskeleton and play crucial roles in various cellular processes, including mitosis and intracellular transport. The property allowing MTs to cater to such diverse roles is attributed to dynamic instability, which is coupled to the hydrolysis of GTP (guanosine-5'-triphosphate) to GDP (guanosine-5'-diphosphate) within the β-tubulin monomers. Understanding the equilibrium dynamics and the structural features of both GDP- and GTP-complexed MT tips, especially at an all-atom level, remains challenging for both experimental and computational methods because of their dynamic nature and the prohibitive computational demands of simulating large, many-protein systems.

On the emergence of machine-learning methods in bottom-up coarse-graining.

Machine-learning methods have gained significant attention in the computational chemistry community as a viable approach to molecular modeling and analysis. Recent successes in utilizing neural networks to learn atomistic force-fields which 'coarse-grain' electronic structure have inspired similar applications to the thermodynamic coarse-graining of chemical and biological systems. In this review, we discuss the current viability and challenges in using machine-learning methods to represent coarse-grained force-fields, as well as the utility of machine-learning in various aspects of coarse-grained modeling.

Kinetic Implications of IP Anion Binding on the Molecular Switch of the HIV-1 Capsid Assembly.

HIV-1 capsid proteins (CA) self-assemble into a fullerene-shaped capsid, enabling cellular transport and nuclear entry of the viral genome. A structural switch comprising the Thr-Val-Gly-Gly (TVGG) motif either assumes a disordered coil or a 3 helix conformation to regulate hexamer or pentamer assembly, respectively. The cellular polyanion inositol hexakisphosphate (IP6) binds to a positively charged pore of CA capsomers rich in arginine and lysine residues mediated by electrostatic interactions.

Enhancing the Assembly Properties of Bottom-Up Coarse-Grained Phospholipids.

A plethora of key biological events occur at the cellular membrane where the large spatiotemporal scales necessitate dimensionality reduction or coarse-graining approaches over conventional all-atom molecular dynamics simulation. Constructing coarse-grained descriptions of membranes systematically from statistical mechanical principles has largely remained challenging due to the necessity of capturing amphipathic self-assembling behavior in coarse-grained models. We show that bottom-up coarse-grained lipid models can possess metastable morphological behavior and that this potential metastability has ramifications for accurate development and training.

Activation of the Influenza B M2 Proton Channel (BM2).

Influenza B viruses have cocirculated during most seasonal flu epidemics and can cause significant human morbidity and mortality due to their rapid mutation, emerging drug resistance, and severe impact on vulnerable populations. The influenza B M2 proton channel (BM2) plays an essential role in viral replication, but the mechanisms behind its symmetric proton conductance and the involvement of a second histidine (His27) cluster remain unclear. Here we performed membrane-enabled continuous constant-pH molecular dynamics simulations on wildtype BM2 and a key H27A mutant channel to explore its pH-dependent conformational switch.

Megamolecule Self-Assembly Networks: A Combined Computational and Experimental Design Strategy.

This work describes the use of computational strategies to design megamolecule building blocks for the self-assembly of lattice networks. The megamolecules are prepared by attaching four Cutinase-SnapTag fusion proteins (CS fusions) to a four-armed linker, followed by functionalizing each fusion with a terpyridine linker. This functionality is designed to participate in a metal-mediated self-assembly process to give networks.

Changing Your Martini Can Still Give You a Hangover.

The Martini 3.0 coarse-grained force field, which was parametrized to better capture transferability in top-down coarse-grained models, is analyzed to assess its accuracy in representing thermodynamic and structural properties with respect to the underlying atomistic representation of the system. These results are compared to those obtained following the principles of statistical mechanics that start from the same underlying atomistic system.

Lipid organization by the Caveolin-1 complex.

Caveolins are lipid-binding proteins that can organize membrane remodeling and oligomerize into the 8S complex. The CAV1-8S complex comprises a disk-like structure, about 15 nm in diameter, with a central beta barrel. Further oligomerization of 8S complexes remodels the membrane into caveolae vessels, with a dependence on cholesterol concentration.

Quantitative insights into the mechanism of proton conduction and selectivity for the human voltage-gated proton channel Hv1.

Human voltage-gated proton (hHv1) channels are crucial for regulating essential biological processes such as immune cell respiratory burst, sperm capacitation, and cancer cell migration. Despite the significant concentration difference between protons and other ions in physiological conditions, hHv1 demonstrates remarkable proton selectivity. Our calculations of single-proton, cation, and anion permeation free energy profiles quantitatively demonstrate that the proton selectivity of the wild-type channel originates from its strong proton affinity via the titration of the key residues D112 and D174, although the channel imposes similar kinetic blocking effects for protons compared to other ions.

Activation of the influenza B M2 proton channel (BM2).

Influenza B viruses have co-circulated during most seasonal flu epidemics and can cause significant human morbidity and mortality due to their rapid mutation, emerging drug resistance, and severe impact on vulnerable populations. The influenza B M2 proton channel (BM2) plays an essential role in viral replication, but the mechanisms behind its symmetric proton conductance and the involvement of a second histidine (His27) cluster remain unclear. Here we perform the membrane-enabled continuous constant-pH molecular dynamics simulations on wildtype BM2 and a key H27A mutant to explore its pH-dependent conformational switch.

Lipid Organization by the Caveolin-1 Complex.

Caveolins are lipid-binding proteins that can organize membrane remodeling and oligomerize into the 8S-complex. The CAV1 8S-complex comprises a disk-like structure, about 15nm in diameter, with a central beta barrel. Further oligomerization of 8S-complexes remodels the membrane into caveolae vessels, with a dependence on cholesterol concentration.

QM/CG-MM: Systematic Embedding of Quantum Mechanical Systems in a Coarse-Grained Environment with Accurate Electrostatics.

Quantum Mechanics/Molecular Mechanics (QM/MM) can describe chemical reactions in molecular dynamics (MD) simulations at a much lower cost than MD. Still, it is prohibitively expensive for many systems of interest because such systems usually require long simulations for sufficient statistical sampling. Additional MM degrees of freedom are often slow and numerous but secondary in interest.

Understanding dynamics in coarse-grained models. IV. Connection of fine-grained and coarse-grained dynamics with the Stokes-Einstein and Stokes-Einstein-Debye relations.

Applying an excess entropy scaling formalism to the coarse-grained (CG) dynamics of liquids, we discovered that missing rotational motions during the CG process are responsible for artificially accelerated CG dynamics. In the context of the dynamic representability between the fine-grained (FG) and CG dynamics, this work introduces the well-known Stokes-Einstein and Stokes-Einstein-Debye relations to unravel the rotational dynamics underlying FG trajectories, thereby allowing for an indirect evaluation of the effective rotations based only on the translational information at the reduced CG resolution. Since the representability issue in CG modeling limits a direct evaluation of the shear stress appearing in the Stokes-Einstein and Stokes-Einstein-Debye relations, we introduce a translational relaxation time as a proxy to employ these relations, and we demonstrate that these relations hold for the ambient conditions studied in our series of work.

Mechanism of phosphate release from actin filaments.

After ATP-actin monomers assemble filaments, the ATP's [Formula: see text]-phosphate is hydrolyzedwithin seconds and dissociates over minutes. We used all-atom molecular dynamics simulations to sample the release of phosphate from filaments and study residues that gate release. Dissociation of phosphate from Mg is rate limiting and associated with an energy barrier of 20 kcal/mol, consistent with experimental rates of phosphate release.

Cracked actin filaments as mechanosensitive receptors.

Actin filament networks are exposed to mechanical stimuli, but the effect of strain on actin filament structure has not been well established in molecular detail. This is a critical gap in understanding because the activity of a variety of actin-binding proteins has recently been determined to be altered by actin filament strain. We therefore used all-atom molecular dynamics simulations to apply tensile strains to actin filaments and find that changes in actin subunit organization are minimal in mechanically strained, but intact, actin filaments.

The SARS-CoV-2 nucleoprotein associates with anionic lipid membranes.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a lipid-enveloped virus that acquires its lipid bilayer from the host cell it infects. SARS-CoV-2 can spread from cell to cell or from patient to patient by undergoing assembly and budding to form new virions. The assembly and budding of SARS-CoV-2 is mediated by several structural proteins known as envelope (E), membrane (M), nucleoprotein (N), and spike (S), which can form virus-like particles (VLPs) when co-expressed in mammalian cells.

Interaction of MRI Contrast Agent [Gd(DOTA)] with Lipid Membranes: A Molecular Dynamics Study.

Contrast agents are important imaging probes in clinical MRI, allowing the identification of anatomic changes that otherwise would not be possible. Intensive research on the development of new contrast agents is being made to image specific pathological markers or sense local biochemical changes. The most widely used MRI contrast agents are based on gadolinium(III) complexes.

Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package.

Simulating chemically reactive phenomena such as proton transport on nanosecond to microsecond and beyond time scales is a challenging task. methods are unable to currently access these time scales routinely, and traditional molecular dynamics methods feature fixed bonding arrangements that cannot account for changes in the system's bonding topology. The Multiscale Reactive Molecular Dynamics (MS-RMD) method, as implemented in the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) software package for the LAMMPS molecular dynamics code, offers a method to routinely sample longer time scale reactive simulation data with statistical precision.

Making the cut: Multiscale simulation of membrane remodeling.

Biological membranes are dynamic heterogeneous materials, and their shape and organization are tightly coupled to the properties of the proteins in and around them. However, the length scales of lipid and protein dynamics are far below the size of membrane-bound organelles, much less an entire cell. Therefore, multiscale modeling approaches are often necessary to build a comprehensive picture of the interplay of these factors, and have provided critical insights into our understanding of membrane dynamics.

Prediction of the essential intermolecular contacts for side-binding of VASP on F-actin.

Vasodilator-stimulated phosphoprotein (VASP) family proteins play a crucial role in mediating the actin network architecture in the cytoskeleton. The Ena/VASP homology 2 (EVH2) domain in each of the four identical arms of the tetrameric VASP consists of a loading poly-Pro region, a G-actin-binding domain (GAB), and an F-actin-binding domain (FAB). Together, the poly-Pro, GAB, and FAB domains allow VASP to bind to sides of actin filaments in a bundle, and recruit profilin-G-actin to processively elongate the filaments.

Molecular Tuning of Reactivity of Zeolite Protons in HZSM-5.

In acidic HZSM-5 zeolite, the reactivity of a methanol molecule interacting with the zeolite proton is amenable to modification via coadsorbing a stochiometric amount of an electron density donor to form the [()(CHOH)(HZ)] complex. The rate of the methanol in this complex undergoing dehydration to dimethyl ether was determined for a series of with proton affinity (PA) ranging from 659 kJ mol for CF to 825 kJ mol for CHO and was found to follow the expression: Ln(Rate) - Ln(Rate) = β(PA - PA), where = N is the reference and β and γ are constants. This trend is probably due to the increased stability of the solvated proton in the [()(CHOH)(HZ)] complex with increasing PA.