The Associative memory, Water mediated, Structure and Energy Model (AWSEM) is a coarse-grained protein force field. AWSEM contains physically motivated terms, such as hydrogen bonding, as well as a bioinformatically based local structure biasing term, which efficiently takes into account many-body effects that are modulated by the local sequence. When combined with appropriate local or global alignments to choose memories, AWSEM can be used to perform de novo protein structure prediction. Herein we present structure prediction results for a particular choice of local sequence alignment method based on short residue sequences called fragments. We demonstrate the model’s structure prediction capabilities for three levels of global homology between the target sequence and those proteins used for local structure biasing, all of which assume that the structure of the target sequence is not known. When there are no homologs in the database of structures used for local structure biasing, AWSEM calculations produce structural predictions that are somewhat improved compared with prior works using related approaches. The inclusion of a small number of structures from homologous sequences improves structure prediction only marginally but when the fragment search is restricted to only homologous sequences, AWSEM can perform high resolution structure prediction and can be used for kinetics and dynamics studies.

The Nuclear Pore Complex (NPC, ~50 MDa) is the sole passageway for the transport of macromolecules across the nuclear envelope. The NPC plays a key role in numerous critical cellular processes such as transcription, and many of its components are implicated in human diseases such as cancer. Previous work (ref 1, 2) defined the relative positions of its 456 constituent proteins (nucleoporin or Nups), based on spatial restraints derived from biophysical, electron microscopy, and proteomic data. Further elucidation of the evolutionary origin, transport mechanism, and assembly of the NPC will require higher resolution information. As part of an effort to improve upon the resolution and accuracy of the NPC structure, we set out to determine the atomic structures of the NPC components. Because it proved difficult to determine the atomic structures of whole Nups by X-ray crystallography alone, we are relying on multiple datasets that are combined computationally by our Integrative Modeling Platform (IMP) package (http://salilab.org/imp). In particular, we developed an integrative modeling approach that benefits from crystallographic structures of fragments of the protein or its homologs, Solution Small Angle X-ray Scattering (SAXS) profiles of the protein and its fragments (ref 3), NMR, and negative stain Electron Microscopy (EM) micrographs of the protein. Each dataset is converted into a set of spatial restraints on the protein structure, followed by finding a model that satisfies the restraints as well as possible using a Monte Carlo / molecular dynamics optimization procedure. The approach will be illustrated by its application to yeast Nup133.

We investigate protein-protein association using the associativememory, water-mediated, structure, and energy model (AWSEM), a coarse-grained protein folding model that has been optimized using energy-landscape theory. The potential was originally parameterized by enforcing a funneled nature for a database of dimeric interfaces but was later further optimized to create funneled folding landscapes for individual monomeric proteins. The ability of the model to predict interfaces was not tested previously. The present results show that simulated annealing of the model indeed is able to predict successfully the native interfaces of eight homodimers and four heterodimers, thus amounting to a flexible docking algorithm. We go on to address the relative importance of monomer geometry, flexibility, and nonnative intermonomeric contacts in the association process for the homodimers. Monomer surface geometry is found to be important in determining the binding interface, but it is insufficient. Using a uniform binding potential rather than the water-mediated potential results in sampling of misbound structures that are geometrically preferred but are nonetheless energetically disfavored by AWSEM, as well as in nature. Depending on the stability of the unbound monomers, nonnative contacts play different roles in the association process. For unstable monomers, thermodynamic states stabilized by nonnative interactions correspond to productive, on-pathway intermediates and can, therefore, catalyze binding through a fly-casting mechanism. For stable monomers, in contrast, states stabilized by nonnative interactions generally correspond to traps that impede binding.binding interface prediction | swapped contacts P rotein-protein interfaces encode information that is key to a molecular understanding of biological functions. The folding of proteins is well understood in the framework of energy landscape theory and its principle of minimal frustration. Are binding landscapes also funneled? Mechanistic consequences of funneled binding landscapes have been investigated using structure-based models (1-5). The agreement of these mechanisms with observation suggests that binding landscapes are generally funneled, explaining why topology is indeed a major factor in determining binding mechanisms (1). A statistical analysis of a large database of protein complexes revealed that for many of the complexes, the binding energy gap is indeed larger than expected knowing the variance of the binding energy (6), the hallmark feature of a funneled landscape (7). When further testing this idea, Papoian et al. discovered that for other complexes, to have a funneled landscape for binding, unanticipated water-mediated interactions were required. They developed a water-mediated potential encoding these interactions (8). This transferable potential was later optimized to create funneled folding landscapes that successfully predict the structure of monomeric proteins (9, 10). Therefore, there is considerable support for the idea that, like folding landsc...

Coarse-grained (CG) models of molecular systems, with fewer mechanical degrees of freedom than an all-atom model, are used extensively in chemical physics. It is generally accepted that a coarse-grained model that accurately describes equilibrium structural properties (as a result of having a well constructed CG potential energy function) does not necessarily exhibit appropriate dynamical behavior when simulated using conservative Hamiltonian dynamics for the CG degrees of freedom on the CG potential energy surface. Attempts to develop accurate CG dynamic models usually focus on replacing Hamiltonian motion by stochastic but Markovian dynamics on that surface, such as Langevin or Brownian dynamics. However, depending on the nature of the system and the extent of the coarse-graining, a Markovian dynamics for the CG degrees of freedom may not be appropriate. In this paper, we consider the problem of constructing dynamic CG models within the context of the Multi-Scale Coarse-graining (MS-CG) method of Voth and coworkers. We propose a method of converting a MS-CG model into a dynamic CG model by adding degrees of freedom to it in the form of a small number of fictitious particles that interact with the CG degrees of freedom in simple ways and that are subject to Langevin forces. The dynamic models are members of a class of nonlinear systems interacting with special heat baths that were studied by Zwanzig [J. Stat. Phys. 9, 215 (1973)]. The properties of the fictitious particles can be inferred from analysis of the dynamics of all-atom simulations of the system of interest. This is analogous to the fact that the MS-CG method generates the CG potential from analysis of equilibrium structures observed in all-atom simulation data. The dynamic models generate a non-Markovian dynamics for the CG degrees of freedom, but they can be easily simulated using standard molecular dynamics programs. We present tests of this method on a series of simple examples that demonstrate that the method provides realistic dynamical CG models that have non-Markovian or close to Markovian behavior that is consistent with the actual dynamical behavior of the all-atom system used to construct the CG model. Both the construction and the simulation of such a dynamic CG model have computational requirements that are similar to those of the corresponding MS-CG model and are good candidates for CG modeling of very large systems.

The increasing interest in the modeling of complex macromolecular systems in recent years has spurred the development of numerous coarse-graining (CG) techniques. However, many of the CG models are constructed assuming that all details beneath the resolution of CG degrees of freedom are fast and average out, which sets limits on the resolution of feasible coarse-grainings and on the range of applications of the CG models. Ultra-coarse-graining (UCG) makes it possible to construct models at any desired resolution while accounting for discrete conformational or chemical changes within the CG sites that can modulate the interactions between them. Here, we discuss the UCG methodology and its numerical implementation. We pay particular attention to the numerical mechanism for including state transitions between different conformations within CG sites because this has not been discussed previously. Using a simple example of 1,2-dichloroethane, we demonstrate the ability of the UCG model to reproduce the multiconfigurational behavior of this molecular liquid, even when each molecule is modeled with only one CG site. The methodology can also be applied to other molecular liquids and macromolecular systems with time scale separation between conformational transitions and other intramolecular motions and rotations.

Immune checkpoint inhibitors (ICIs) have obtained durable responses in many cancers, making it possible to foresee their potential in improving the health of cancer patients. However, immunotherapies are currently limited to a minority of patients and there is a need to develop a better understanding of the basic molecular mechanisms and functions of pivotal immune regulatory molecules. Immune checkpoint cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) and regulatory T (Treg) cells play pivotal roles in hindering the anticancer immunity. Treg cells suppress antigen-presenting cells (APCs) by depleting immune stimulating cytokines, producing immunosuppressive cytokines and constitutively expressing CTLA-4. CTLA-4 molecules bind to CD80 and CD86 with a higher affinity than CD28 and act as competitive inhibitors of CD28 in APCs. The purpose of this review is to summarize state-of-the-art understanding of the molecular mechanisms underlining CTLA-4 immune regulation and the correlation of the ICI response with CTLA-4 expression in Treg cells from preclinical and clinical studies for possibly improving CTLA-4-based immunotherapies, while highlighting the knowledge gap.

We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.

The protein p53 is a crucial tumor suppressor, often called “the guardian of the genome”; however, mutations transform p53 into a powerful cancer promoter. The oncogenic capacity of mutant p53 has been ascribed to enhanced propensity to fibrillize and recruit other cancer fighting proteins in the fibrils, yet the pathways of fibril nucleation and growth remain obscure. Here, we combine immunofluorescence three-dimensional confocal microscopy of human breast cancer cells with light scattering and transmission electron microscopy of solutions of the purified protein and molecular simulations to illuminate the mechanisms of phase transformations across multiple length scales, from cellular to molecular. We report that the p53 mutant R248Q (R, arginine; Q, glutamine) forms, both in cancer cells and in solutions, a condensate with unique properties, mesoscopic protein-rich clusters. The clusters dramatically diverge from other protein condensates. The cluster sizes are decoupled from the total cluster population volume and independent of the p53 concentration and the solution concentration at equilibrium with the clusters varies. We demonstrate that the clusters carry out a crucial biological function: they host and facilitate the nucleation of amyloid fibrils. We demonstrate that the p53 clusters are driven by structural destabilization of the core domain and not by interactions of its extensive unstructured region, in contradistinction to the dense liquids typical of disordered and partially disordered proteins. Two-step nucleation of mutant p53 amyloids suggests means to control fibrillization and the associated pathologies through modifying the cluster characteristics. Our findings exemplify interactions between distinct protein phases that activate complex physicochemical mechanisms operating in biological systems.