GAP positions catalytic H-Ras residue Q61 for GTP hydrolysis in molecular dynamics simulations, complicating chemical rescue of Ras deactivation.

Functional interaction of Ras signaling proteins with upstream, negative regulatory GTPase activating proteins (GAPs) represents a crucial step in cellular decision making related to growth and survival. Key components of the catalytic transition state for Ras deactivation by GAP-accelerated hydrolysis of Ras-bound guanosine triphosphate (GTP) are thought to include an arginine residue from the GAP (the arginine finger), a glutamine residue from Ras (Q61), and a water molecule that is likely coordinated by Q61 to engage in nucleophilic attack on GTP. Here, we use in-vitro fluorescence experiments to show that 0.1-100 mM concentrations of free arginine, imidazole, and other small nitrogenous molecule fail to accelerate GTP hydrolysis, even in the presence of the catalytic domain of a mutant GAP lacking its arginine finger (R1276A NF1). This result is surprising given that imidazole can chemically rescue enzyme activity in arginine-to-alanine mutant protein tyrosine kinases (PTKs) that share many active site components with Ras/GAP complexes. Complementary all-atom molecular dynamics (MD) simulations reveal that an arginine finger GAP mutant still functions to enhance Ras Q61-GTP interaction, though less extensively than wild-type GAP. This increased Q61-GTP proximity may promote more frequent fluctuations into configurations that enable GTP hydrolysis as a component of the mechanism by which GAPs accelerate Ras deactivation in the face of arginine finger mutations. The failure of small molecule analogs of arginine to chemically rescue catalytic deactivation of Ras is consistent with the idea that the influence of the GAP goes beyond the simple provision of its arginine finger. However, the failure of chemical rescue in the presence of R1276A NF1 suggests that the GAPs arginine finger is either unsusceptible to rescue due to exquisite positioning or that it is involved in complex multivalent interactions. Therefore, in the context of oncogenic Ras proteins with mutations at codons 12 or 13 that inhibit arginine finger penetration toward GTP, drug-based chemical rescue of GTP hydrolysis may have bifunctional chemical/geometric requirements that are more difficult to satisfy than those that result from arginine-to-alanine mutations in other enzymes for which chemical rescue has been demonstrated.

Inferring Pathways of Oxidative Folding from Prefolding Free Energy Landscapes of Disulfide-Rich Toxins.

Short, cysteine-rich peptides can exist in stable or metastable structural ensembles due to the number of possible patterns of formation of their disulfide bonds. One interesting subset of this peptide group is the conotoxins, which are produced by aquatic snails in the family Conidae. The µ conotoxins, which are antagonists and blockers of the voltage-gated sodium channel, exist in a folding spectrum: on one end of the spectrum are more hirudin-like folders, which form disulfide bonds and then reshuffle them, leading to an ensemble of kinetically trapped isomers, and on the other end are more BPTI-like folders, which form the native disulfide bonds one by one in a particular order, leading to a preponderance of conformations existing in a single stable state. In this Article, we employ the composite diffusion map approach to study the unified free energy surface of prefolding µ-conotoxin equilibrium. We identify the two most important nonlinear collective modes of the unified folding landscape and demonstrate that in the absence of their disulfides, the conotoxins can be thought of as largely disordered polymers. A small increase in the number of hydrophobic residues in the protein shifts the free energy landscape toward hydrophobically collapsed coil conformations responsible for cysteine proximity in hirudin-like folders, compared to semiextended coil conformations with more distal cysteines in BPTI-like folders. Overall, this work sheds important light on the folding processes and free energy landscapes of cysteine-rich peptides and demonstrates the extent to which sequence and length contribute to these landscapes.

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework.

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. Here, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (∼98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ∼1.5 million atoms is discussed in the context of MuMMI.

Drug Discovery by Automated Adaptation of Chemical Structure and Identity.

Computer-aided drug design offers the potential to dramatically reduce the cost and effort required for drug discovery. While screening-based methods are valuable in the early stages of hit identification, they are frequently succeeded by iterative, hypothesis-driven computations that require recurrent investment of human time and intuition. To increase automation, we introduce a computational method for lead refinement that combines concerted dynamics of the ligand/protein complex via molecular dynamics simulations with integrated Monte Carlo-based changes in the chemical formula of the ligand. This approach, which we refer to as ligand-exchange Monte Carlo molecular dynamics, accounts for solvent- and entropy-based contributions to competitive binding free energies by coupling the energetics of bound and unbound states during the ligand-exchange attempt. Quantitative comparison of relative binding free energies to reference values from free energy perturbation, conducted in vacuum, indicates that ligand-exchange Monte Carlo molecular dynamics simulations sample relevant conformational ensembles and are capable of identifying strongly binding compounds. Additional simulations demonstrate the use of an implicit solvent model. We speculate that the use of chemical graphs in which exchanges are only permitted between ligands with sufficient similarity may enable an automated search to capture some of the benefits provided by human intuition during hypothesis-guided lead refinement.

NaCl aggregation in water at elevated temperatures and pressures: Comparison of classical force fields.

The properties of water vary dramatically with temperature and density. This can be exploited to control its effectiveness as a solvent. Thus, supercritical water is of keen interest as solvent in many extraction processes. The low solubility of salts in lower density supercritical water has even been suggested as a means of desalination. The high temperatures and pressures required to reach supercritical conditions can present experimental challenges during collection of required physical property and phase equilibria data, especially in salt-containing systems. Molecular simulations have the potential to be a valuable tool for examining the behavior of solvated ions at these high temperatures and pressures. However, the accuracy of classical force fields under these conditions is unclear. We have, therefore, undertaken a parametric study of NaCl in water, comparing several salt and water models at 200 bar-600 bar and 450 K-750 K for a range of salt concentrations. We report a comparison of structural properties including ion aggregation, hydrogen bonding, density, and static dielectric constants. All of the force fields qualitatively reproduce the trends in the liquid phase density. An increase in ion aggregation with decreasing density holds true for all of the force fields. The propensity to aggregate is primarily determined by the salt force field rather than the water force field. This coincides with a decrease in the water static dielectric constant and reduced charge screening. While a decrease in the static dielectric constant with increasing NaCl concentration is consistent across all model combinations, the salt force fields that exhibit more ionic aggregation yield a slightly smaller dielectric decrement.

The RIT1 C-terminus associates with lipid bilayers via charge complementarity.

RIT1 is a member of the Ras superfamily of small GTPases involved in regulation of cellular signaling. Mutations to RIT1 are involved in cancer and developmental disorders. Like many Ras subfamily members, RIT1 is localized to the plasma membrane. However, RIT1 lacks the C-terminal prenylation that helps many other subfamily members adhere to cellular membranes. We used molecular dynamics simulations to examine the mechanisms by which the C-terminal peptide (CTP) of RIT1 associates with lipid bilayers. We show that the CTP is unstructured and that its membrane interactions depend on lipid composition. While a 12-residue region of the CTP binds strongly to anionic bilayers containing phosphatidylserine lipids, the CTP termini fray from the membrane allowing for accommodation of the RIT1 globular domain at the membrane-water interface.

Thermodynamics, dynamics, and structure of supercritical water at extreme conditions.

Molecular dynamics (MD) simulations to understand the thermodynamic, dynamic, and structural changes in supercritical water across the Frenkel line and the melting line have been performed. The two-phase thermodynamic model [J. Phys. Chem. B, 2010, 114(24), 8191-8198] and the velocity autocorrelation functions are used to locate the Frenkel line and to calculate the thermodynamic and dynamic properties. The Frenkel lines obtained from the two-phase thermodynamic model and the velocity autocorrelation criterion do not agree with each other. Structural characteristics and the translational diffusion dynamics of water suggest that this inconsistency could arise from the two oscillatory modes in water, which are associated with the bending of hydrogen bonds and intermolecular collisions inside the first coordination shell. The overall results lead us to conclude that the universality of the Frenkel line as a dynamic crossover line from rigid to nonrigid fluids is preserved in water.

Electrical conductivity, ion pairing, and ion self-diffusion in aqueous NaCl solutions at elevated temperatures and pressures.

We have performed classical molecular dynamics (MD) simulations of aqueous sodium chloride (NaCl) solutions from 298 to 674 K at 200 bars to understand the influence of ion pairing and ion self-diffusion on electrical conductivity in high-temperature/high-pressure salt solutions. Conductivity data obtained from the MD simulation highlight an apparent anomaly, namely, a conductivity maximum as temperature increases along an isobar, which has been also observed in experimental studies. By examining both velocity autocorrelation and cross-correlation terms of the Green-Kubo integral, we quantitatively demonstrate that the conductivity anomaly arises mainly from a competition between the single-ion self-diffusion and the contact ion pair formation. The velocity autocorrelation function in conjunction with structural analysis suggests that diffusive motion of ions is suppressed at high temperatures due to the persistence of an inner hydration shell. The contribution of velocity cross-correlation functions between oppositely charged ions becomes significant at the onset of the conductivity decrease. Structural analysis based on Voronoi tessellation and pair correlation functions indicates that the fraction of contact ion pairs increases as temperature increases. Spatial decomposition of the electrical conductivity also indicates that the formation of contact ion pairs significantly decreases the electrical conductivity compared to Nernst-Einstein conductivity, but the contribution of distant opposite charges cannot be ignored except at the highest temperature due to unscreened long-range interactions.

Simulations of NaCl Aggregation from Solution: Solvent Determines Topography of Free Energy Landscape.

The partition-enabled analysis of cluster histograms (PEACH) method is used to calculate the free energy surface of NaCl aggregation using cluster statistics from MD simulations of small systems (40-90 ions plus solvent) in four solvents. In all cases (pure methanol, pure water, and two methanol/water mixtures) NaCl clusters show a transition from amorphous to rocksalt structure with increasing cluster size. The crossover sizes, and the apparent kinetic barrier to ordering, increase with increasing water content. Implications for the proposed two-step mechanism of NaCl crystal nucleation (in which the ordered structure emerges from a large disordered cluster), and how this mechanism might depend on solvent and on degree of supersaturation, are discussed. In pure water, nonideal crowding effects that promote clustering are identified from systematic concentration-dependent deviations between simulation results and the PEACH model fit. In contrast, the ability of PEACH to fit aggregation statistics in mixed solvents is consistent with negligible interactions between ions in different clusters. © 2018 Wiley Periodicals, Inc.

Cluster Free Energies from Simple Simulations of Small Numbers of Aggregants: Nucleation of Liquid MTBE from Vapor and Aqueous Phases.

We introduce a global fitting analysis method to obtain free energies of association of noncovalent molecular clusters using equilibrated cluster size distributions from unbiased constant-temperature molecular dynamics (MD) simulations. Because the systems simulated are small enough that the law of mass action does not describe the aggregation statistics, the method relies on iteratively determining a set of cluster free energies that, using appropriately weighted sums over all possible partitions of N monomers into clusters, produces the best-fit size distribution. The quality of these fits can be used as an objective measure of self-consistency to optimize the cutoff distance that determines how clusters are defined. To showcase the method, we have simulated a united-atom model of methyl tert-butyl ether (MTBE) in the vapor phase and in explicit water solution over a range of system sizes (up to 95 MTBE in the vapor phase and 60 MTBE in the aqueous phase) and concentrations at 273 K. The resulting size-dependent cluster free energy functions follow a form derived from classical nucleation theory (CNT) quite well over the full range of cluster sizes, although deviations are more pronounced for small cluster sizes. The CNT fit to cluster free energies yielded surface tensions that were in both cases lower than those for the simulated planar interfaces. We use a simple model to derive a condition for minimizing non-ideal effects on cluster size distributions and show that the cutoff distance that yields the best global fit is consistent with this condition.

Coarse-grained molecular simulations of the melting kinetics of small unilamellar vesicles.

Simulations of small unilamellar lipid bilayer vesicles have been performed to model their response to an instantaneous rise in temperature, starting from an initial low-temperature structure, to temperatures near or above the main chain transition temperature. The MARTINI coarse-grained force-field was used to construct slabs of gel-phase DPPC bilayers, which were assembled into truncated icosahedral structures containing 13,165 or 31,021 lipids. Equilibration at 280 K produced structures with several (5-8) domains, characterized by facets of lipids packed in the gel phase connected by disordered ridges. Instantaneous heating to final temperatures ranging from 290 K to 310 K led to partial or total melting over 500 ns trajectories, accompanied by changes in vesicle shape and the sizes and arrangements of remaining gel-phase domains. At temperatures that produced partial melting, the gel-phase lipid content of the vesicles followed an exponential decay, similar in form and timescale to the sub-microsecond phase of melting kinetics observed in recent ultrafast IR temperature-jump experiments. The changing rate of melting appears to be the outcome of a number of competing contributions, but changes in curvature stress arising from the expansion of the bilayer area upon melting are a major factor. The simulations give a more detailed picture of the changes that occur in frozen vesicles following a temperature jump, which will be of use for the interpretation of temperature-jump experiments on vesicles.