Self-Guided Molecular Simulation to Enhance Concerted Motion.

Self-guided (SG) molecular simulation methods, namely self-guided molecular dynamics (SGMD) and self-guided Langevin dynamics (SGLD), enhance conformational search by promoting low-frequency motion. A simple local time averaging scheme is used to extract low-frequency properties with little overhead in computing costs. For molecular processes to form ordered structures like ligand binding and protein folding, it is believed that concerted motions play crucial roles.

AmberTools.

AmberTools is a free and open-source collection of programs used to set up, run, and analyze molecular simulations. The newer features contained within AmberTools23 are briefly described in this Application note.

Determination of van der Waals Parameters Using a Double Exponential Potential for Nonbonded Divalent Metal Cations in TIP3P Solvent.

A double exponential (DE) functional form for Lennard-Jones (LJ) interactions, proposed in our previous study, has many advantages over LJ potentials including a natural softcore characteristic for the convenience of the pathway-based free-energy calculations, fast convergence, and flexibility in use. In this work, we put the first step on the application of the DE functional form by identifying a DE potential, coined DE-TIP3P, for molecular simulations using the TIP3P water model. The developed DE-TIP3 potential was better than LJ potential in reproducing the experimental water properties.

Reformulation of the self-guided molecular simulation method.

Self-guided molecular/Langevin dynamics (SGMD/SGLD) simulation methods were developed to enhance conformational sampling through promoting low frequency motion of molecular systems and have been successfully applied in many simulation studies. Quantitative understanding of conformational distribution in SGLD has been achieved by separating microscopic properties according to frequency. However, a missing link between the guiding factors and conformational distributions makes it highly empirical and system dependent when choosing the values of the guiding parameters.

Membrane cholesterol dependence of cannabinoid modulation of glycine receptor.

Cannabinoids exert therapeutic effects on several diseases such as chronic pain and startle disease by targeting glycine receptors (GlyRs). Our previous studies have shown that cannabinoids target a serine residue at position 296 in the third transmembrane helix of the α1/α3 GlyR. This site is located on the outside of the ion channel protein at the lipid interface where the cholesterol concentrates.

Protonation state of the selectivity filter of bacterial voltage-gated sodium channels is modulated by ions.

The selectivity filter (SF) of bacterial voltage-gated sodium channels consists of four glutamate residues arranged in a C symmetry. The protonation state population of this tetrad is unclear. To address this question, we simulate the pore domain of bacterial voltage-gated sodium channel of Magnetococcus sp.

A double exponential potential for van der Waals interaction.

Van der Waals (vdw) interaction is an important force between atoms and molecules. Many potential functions have been proposed to model vdw interaction such as the Lennard-Jones (L-J) potential. To overcome certain drawbacks of existing function forms, this work proposes a double exponential (DE) potential that contains a repulsive exponential term and an attractive exponential term.

The homogeneity condition: A simple way to derive isotropic periodic sum potentials for efficient calculation of long-range interactions in molecular simulation.

Isotropic periodic sum (IPS) is a method to calculate long-range interactions based on the homogeneity of simulation systems. By using the isotropic periodic images of a local region to represent remote structures, long-range interactions become a function of the local conformation. This function is called the IPS potential, which folds long-ranged interactions into a short-ranged potential and can be calculated as efficiently as a cutoff method.

Comparison of the umbrella sampling and the double decoupling method in binding free energy predictions for SAMPL6 octa-acid host-guest challenges.

We calculate the absolute binding free energies of tetra-methylated octa-acids host-guest systems as a part of the SAMPL6 blind challenge (receipt ID vq30p). We employed two different free energy simulation methods, i.e.

Hydronium Ions Accompanying Buried Acidic Residues Lead to High Apparent Dielectric Constants in the Interior of Proteins.

Internal ionizable groups are known to play important roles in protein functions. A mystery that has attracted decades of extensive experimental and theoretical studies is the apparent dielectric constants experienced by buried ionizable groups, which are much higher than values expected for protein interiors. Many interpretations have been proposed, such as water penetration, conformational relaxation, local unfolding, protein intrinsic backbone fluctuations, etc.

Atomic Resolution Cryo-EM Structure of β-Galactosidase.

The advent of direct electron detectors has enabled the routine use of single-particle cryo-electron microscopy (EM) approaches to determine structures of a variety of protein complexes at near-atomic resolution. Here, we report the development of methods to account for local variations in defocus and beam-induced drift, and the implementation of a data-driven dose compensation scheme that significantly improves the extraction of high-resolution information recorded during exposure of the specimen to the electron beam. These advances enable determination of a cryo-EM density map for β-galactosidase bound to the inhibitor phenylethyl β-D-thiogalactopyranoside where the ordered regions are resolved at a level of detail seen in X-ray maps at ∼ 1.

Origin of pK Shifts of Internal Lysine Residues in SNase Studied Via Equal-Molar VMMS Simulations in Explicit Water.

Protein internal ionizable groups can exhibit large shifts in pK values. Although the environment and interaction changes have been extensively studied both experimentally and computationally, direct calculation of pK values of these internal ionizable groups in explicit water is challenging due to energy barriers in solvent interaction and in conformational transition. The virtual mixture of multiple states (VMMS) method is a new approach designed to study chemical state equilibrium.

Isotropic periodic sum for multipole interactions and a vector relation for calculation of the Cartesian multipole tensor.

Isotropic periodic sum (IPS) is a method to calculate long-range interactions based on the homogeneity of simulation systems. By using the isotropic periodic images of a local region to represent remote structures, long-range interactions become a function of the local conformation. This function is called the IPS potential; it folds long-ranged interactions into a short-ranged potential and can be calculated as efficiently as a cutoff method.

The ubiquitin ligase Ubr4 controls stability of podocin/MEC-2 supercomplexes.

The PHB-domain protein podocin maintains the renal filtration barrier and its mutation is an important cause of hereditary nephrotic syndrome. Podocin and its Caenorhabditis elegans orthologue MEC-2 have emerged as key components of mechanosensitive membrane protein signalling complexes. Whereas podocin resides at a specialized cell junction at the podocyte slit diaphragm, MEC-2 is found in neurons required for touch sensitivity.

A Virtual Mixture Approach to the Study of Multistate Equilibrium: Application to Constant pH Simulation in Explicit Water.

Chemical and thermodynamic equilibrium of multiple states is a fundamental phenomenon in biology systems and has been the focus of many experimental and computational studies. This work presents a simulation method to directly study the equilibrium of multiple states. This method constructs a virtual mixture of multiple states (VMMS) to sample the conformational space of all chemical states simultaneously.

Self-guided Langevin dynamics via generalized Langevin equation.

Self-guided Langevin dynamics (SGLD) is a molecular simulation method that enhances conformational search and sampling via acceleration of the low frequency motions of the system. This acceleration is produced via introduction of a guiding force which breaks down the detailed-balance property of the dynamics, implying that some reweighting is necessary to perform equilibrium sampling. Here, we eliminate the need of reweighing and show that the NVT and NPT ensembles are sampled exactly by a new version of self-guided motion involving a generalized Langevin equation (GLE) in which the random force is modified so as to restore detailed-balance.

2.2 Å resolution cryo-EM structure of β-galactosidase in complex with a cell-permeant inhibitor.

Cryo-electron microscopy (cryo-EM) is rapidly emerging as a powerful tool for protein structure determination at high resolution. Here we report the structure of a complex between Escherichia coli β-galactosidase and the cell-permeant inhibitor phenylethyl β-D-thiogalactopyranoside (PETG), determined by cryo-EM at an average resolution of ~2.2 angstroms (Å).

Maturation of the HIV-1 core by a non-diffusional phase transition.

The formation of the HIV-1 core is the final step in the viral maturation pathway, resulting in the formation of infectious virus. Most current models for HIV-1 core formation suggest that, upon proteolytic cleavage from the immature Gag, capsid (CA) dissociates into the viral interior before reforming into the core. Here we present evidence for an alternate view of core formation by taking advantage of our serendipitous observation of large membrane-enclosed structures in HIV-1 supernatants from infected cells.

Comparative phosphoproteomic analysis of mammalian glomeruli reveals conserved podocin C-terminal phosphorylation as a determinant of slit diaphragm complex architecture.

Glomerular biology is dependent on tightly controlled signal transduction networks that control phosphorylation of signaling proteins such as cytoskeletal regulators or slit diaphragm proteins of kidney podocytes. Cross-species comparison of phosphorylation events is a powerful mean to functionally prioritize and identify physiologically meaningful phosphorylation sites. Here, we present the result of phosphoproteomic analyses of cow and rat glomeruli to allow cross-species comparisons.

Phosphoproteomic analysis reveals regulatory mechanisms at the kidney filtration barrier.

Diseases of the kidney filtration barrier are a leading cause of ESRD. Most disorders affect the podocytes, polarized cells with a limited capacity for self-renewal that require tightly controlled signaling to maintain their integrity, viability, and function. Here, we provide an atlas of in vivo phosphorylated, glomerulus-expressed proteins, including podocyte-specific gene products, identified in an unbiased tandem mass spectrometry-based approach.

Targeted conformational search with map-restrained self-guided Langevin dynamics: application to flexible fitting into electron microscopic density maps.

We present a map-restrained self-guided Langevin dynamics (MapSGLD) simulation method for efficient targeted conformational search. The targeted conformational search represents simulations under restraints defined by experimental observations and/or by user specified structural requirements. Through map-restraints, this method provides an efficient way to maintain substructures and to set structure targets during conformational searching.

Accurate high-throughput structure mapping and prediction with transition metal ion FRET.

Mapping the landscape of a protein's conformational space is essential to understanding its functions and regulation. The limitations of many structural methods have made this process challenging for most proteins. Here, we report that transition metal ion FRET (tmFRET) can be used in a rapid, highly parallel screen, to determine distances from multiple locations within a protein at extremely low concentrations.

Enhanced Sampling in Free Energy Calculations: Combining SGLD with the Bennett's Acceptance Ratio and Enveloping Distribution Sampling Methods.

One of the key requirements for the accurate calculation of free energy differences is proper sampling of conformational space. Especially in biological applications, molecular dynamics simulations are often confronted with rugged energy surfaces and high energy barriers, leading to insufficient sampling and, in turn, poor convergence of the free energy results. In this work, we address this problem by employing enhanced sampling methods.

Maintain rigid structures in Verlet based cartesian molecular dynamics simulations.

An algorithm is presented to maintain rigid structures in Verlet based cartesian molecular dynamics (MD) simulations. After each unconstrained MD step, the coordinates of selected particles are corrected to maintain rigid structures through an iterative procedure of rotation matrix computation. This algorithm, named as SHAPE and implemented in CHARMM program suite, avoids the calculations of Lagrange multipliers, so that the complexity of computation does not increase with the number of particles in a rigid structure.