Coarse-Grained Simulations of Membrane Insertion and Folding of Small Helical Proteins Using the CABS Model

Integrative modeling of diverse protein-peptide systems using CABS-dock

The CABS model can be applied to a wide range of protein-protein and protein-peptide molecular modeling tasks, such as simulating folding pathways, predicting structures, docking, and analyzing the structural dynamics of molecular complexes. In this work, we use the CABS-dock tool in two diverse modeling tasks: 1) predicting the structures of amyloid protofilaments and 2) identifying cleavage sites in the peptide substrates of proteolytic enzymes. In the first case, simulations of the simultaneous docking of amyloidogenic peptides indicated that the CABS model can accurately predict the structures of amyloid protofilaments which have an in-register parallel architecture. Scoring based on a combination of symmetry criteria and estimated interaction energy values for bound monomers enables the identification of protofilament models that closely match their experimental structures for 5 out of 6 analyzed systems. For the second task, it has been shown that CABS-dock coarse-grained docking simulations can be used to identify the positions of cleavage sites in the peptide substrates of proteolytic enzymes. The cleavage site position was correctly identified for 12 out of 15 analyzed peptides. When combined with sequence-based methods, these docking simulations may lead to an efficient way of predicting cleavage sites in degraded proteins. The method also provides the atomic structures of enzyme-substrate complexes, which can give insights into enzyme-substrate interactions that are crucial for the design of new potent inhibitors.

Hierarchical Coarse-Grained Strategy for Macromolecular Self-Assembly: Application to Hepatitis B Virus-Like Particles

Macromolecular self-assembly is at the basis of many phenomena in material and life sciences that find diverse applications in technology. One example is the formation of virus-like particles (VLPs) that act as stable empty capsids used for drug delivery or vaccine fabrication. Similarly to the capsid of a virus, VLPs are protein assemblies, but their structural formation, stability, and properties are not fully understood, especially as a function of the protein modifications. In this work, we present a data-driven modeling approach for capturing macromolecular self-assembly on scales beyond traditional molecular dynamics (MD), while preserving the chemical specificity. Each macromolecule is abstracted as an anisotropic object and high-dimensional models are formulated to describe interactions between molecules and with the solvent. For this, data-driven protein-protein interaction potentials are derived using a Kriging-based strategy, built on high-throughput MD simulations. Semi-automatic supervised learning is employed in a high performance computing environment and the resulting specialized force-fields enable a significant speed-up to the micrometer and millisecond scale, while maintaining high intermolecular detail. The reported generic framework is applied for the first time to capture the formation of hepatitis B VLPs from the smallest building unit, i.e., the dimer of the core protein HBcAg. Assembly pathways and kinetics are analyzed and compared to the available experimental observations. We demonstrate that VLP self-assembly phenomena and dependencies are now possible to be simulated. The method developed can be used for the parameterization of other macromolecules, enabling a molecular understanding of processes impossible to be attained with other theoretical models.

SPICA Force Field for Proteins and Peptides

A coarse-grained (CG) model for peptides and proteins was developed as an extension of the Surface Property fItting Coarse grAined (SPICA) force field (FF). The model was designed to examine membrane proteins that are fully compatible with the lipid membranes of the SPICA FF. A preliminary version of this protein model was created using thermodynamic properties, including the surface tension and density in the SPICA (formerly called SDK) FF. In this study, we improved the CG protein model to facilitate molecular dynamics (MD) simulations with a reproduction of multiple properties from both experiments and all-atom (AA) simulations. An elastic network model was adopted to maintain the secondary structure within a single chain. The side-chain analogues reproduced the transfer free energy profiles across the lipid membrane and demonstrated reasonable association free energy (potential of mean force) in water compared to those from AA MD. A series of peptides/proteins adsorbed onto or penetrated into the membrane simulated by the CG MD correctly predicted the penetration depths and tilt angles of peripheral and transmembrane peptides/proteins as comparable to those in the orientations of proteins in membranes (OPM) database. In addition, the dimerization free energies of several transmembrane helices within a lipid bilayer were comparable to those from experimental estimation. Application studies on a series of membrane protein assemblies, scramblases, and poliovirus capsids demonstrated the good performance of the SPICA FF.

Computational models for studying physical instabilities in high concentration biotherapeutic formulations

Computational prediction of the behavior of concentrated protein solutions is particularly advantageous in early development stages of biotherapeutics when material availability is limited and a large set of formulation conditions needs to be explored. This review provides an overview of the different computational paradigms that have been successfully used in modeling undesirable physical behaviors of protein solutions with a particular emphasis on high-concentration drug formulations. This includes models ranging from all-atom simulations, coarse-grained representations to macro-scale mathematical descriptions used to study physical instability phenomena of protein solutions such as aggregation, elevated viscosity, and phase separation. These models are compared and summarized in the context of the physical processes and their underlying assumptions and limitations. A detailed analysis is also given for identifying protein interaction processes that are explicitly or implicitly considered in the different modeling approaches and particularly their relations to various formulation parameters. Lastly, many of the shortcomings of existing computational models are discussed, providing perspectives and possible directions toward an efficient computational framework for designing effective protein formulations.

Multiscale Modeling of Amyloid Fibrils Formed by Aggregating Peptides Derived from the Amyloidogenic Fragment of the A-Chain of Insulin

Computational prediction of molecular structures of amyloid fibrils remains an exceedingly challenging task. In this work, we propose a multi-scale modeling procedure for the structure prediction of amyloid fibrils formed by the association of ACC_1-13 aggregation-prone peptides derived from the N-terminal region of insulin’s A-chain. First, a large number of protofilament models composed of five copies of interacting ACC_1-13 peptides were predicted by application of CABS-dock coarse-grained (CG) docking simulations. Next, the models were reconstructed to all-atom (AA) representations and refined during molecular dynamics (MD) simulations in explicit solvent. The top-scored protofilament models, selected using symmetry criteria, were used for the assembly of long fibril structures. Finally, the amyloid fibril models resulting from the AA MD simulations were compared with atomic force microscopy (AFM) imaging experimental data. The obtained results indicate that the proposed multi-scale modeling procedure is capable of predicting protofilaments with high accuracy and may be applied for structure prediction and analysis of other amyloid fibrils.

Sampling of the conformational landscape of small proteins with Monte Carlo methods

Computer simulation provides an increasingly realistic picture of large-scale conformational change of proteins, but investigations remain fundamentally constrained by the femtosecond timestep of molecular dynamics simulations. For this reason, many biologically interesting questions cannot be addressed using accessible state-of-the-art computational resources. Here, we report the development of an all-atom Monte Carlo approach that permits the modelling of the large-scale conformational change of proteins using standard off-the-shelf computational hardware and standard all-atom force fields. We demonstrate extensive thermodynamic characterization of the folding process of the α-helical Trp-cage, the Villin headpiece and the β-sheet WW-domain. We fully characterize the free energy landscape, transition states, energy barriers between different states, and the per-residue stability of individual amino acids over a wide temperature range. We demonstrate that a state-of-the-art intramolecular force field can be combined with an implicit solvent model to obtain a high quality of the folded structures and also discuss limitations that still remain.

Flexible docking of peptides to proteins using CABS-dock

Molecular docking of peptides to proteins can be a useful tool in the exploration of the possible peptide binding sites and poses. CABS-dock is a method for protein-peptide docking that features significant conformational flexibility of both the peptide and the protein molecules during the peptide search for a binding site. The CABS-dock has been made available as a web server and a standalone package. The web server is an easy to use tool with a simple web interface. The standalone package is a command-line program dedicated to professional users. It offers a number of advanced features, analysis tools and support for large-sized systems. In this article, we outline the current status of the CABS-dock method, its recent developments, applications, and challenges ahead.

Extension of the UNRES Coarse-Grained Force Field to Membrane Proteins in the Lipid Bilayer

The physics-based UNRES coarse-grained force field for the simulations of protein structure and dynamics has been extended to treat membrane proteins. The lipid bilayer has been modeled by introducing a continuous nonpolar phase with the water-interface region of appropriate thickness. The potentials for average electrostatic and correlation interactions of the peptide groups have been rescaled to account for the reduction of the dielectric permittivity compared to the water phase and new potentials for protein side-chain-side-chain interactions inside and across the lipid phase have been introduced. The model was implemented in the UNRES package for coarse-grained simulations of proteins, and the package with the new functionality was tested for total energy conservation and thermostat behavior in microcanonical and canonical molecular dynamics simulations runs, respectively. The method was validated by running unrestricted ab initio blind-prediction tests of 10 short α-helical membrane proteins, all runs started from the extended structures. The modified UNRES force field was able to predict correctly the overall folds of the membrane proteins studied.

Modeling of Disordered Protein Structures Using Monte Carlo Simulations and Knowledge-Based Statistical Force Fields

The description of protein disordered states is important for understanding protein folding mechanisms and their functions. In this short review, we briefly describe a simulation approach to modeling protein interactions, which involve disordered peptide partners or intrinsically disordered protein regions, and unfolded states of globular proteins. It is based on the CABS coarse-grained protein model that uses a Monte Carlo (MC) sampling scheme and a knowledge-based statistical force field. We review several case studies showing that description of protein disordered states resulting from CABS simulations is consistent with experimental data. The case studies comprise investigations of protein⁻peptide binding and protein folding processes. The CABS model has been recently made available as the simulation engine of multiscale modeling tools enabling studies of protein⁻peptide docking and protein flexibility. Those tools offer customization of the modeling process, driving the conformational search using distance restraints, reconstruction of selected models to all-atom resolution, and simulation of large protein systems in a reasonable computational time. Therefore, CABS can be combined in integrative modeling pipelines incorporating experimental data and other modeling tools of various resolution.

Modeling of Protein Structural Flexibility and Large-Scale Dynamics: Coarse-Grained Simulations and Elastic Network Models

Fluctuations of protein three-dimensional structures and large-scale conformational transitions are crucial for the biological function of proteins and their complexes. Experimental studies of such phenomena remain very challenging and therefore molecular modeling can be a good alternative or a valuable supporting tool for the investigation of large molecular systems and long-time events. In this minireview, we present two alternative approaches to the coarse-grained (CG) modeling of dynamic properties of protein systems. We discuss two CG representations of polypeptide chains used for Monte Carlo dynamics simulations of protein local dynamics and conformational transitions, and highly simplified structure-based elastic network models of protein flexibility. In contrast to classical all-atom molecular dynamics, the modeling strategies discussed here allow the quite accurate modeling of much larger systems and longer-time dynamic phenomena. We briefly describe the main features of these models and outline some of their applications, including modeling of near-native structure fluctuations, sampling of large regions of the protein conformational space, or possible support for the structure prediction of large proteins and their complexes.

Helix formation and stability in membranes

In this article we review current understanding of basic principles for the folding of membrane proteins, focusing on the more abundant alpha-helical class. Membrane proteins, vital to many biological functions and implicated in numerous diseases, fold into their active conformations in the complex environment of the cell bilayer membrane. While many membrane proteins rely on the translocon and chaperone proteins to fold correctly, others can achieve their functional form in the absence of any translation apparatus or other aides. Nevertheless, the spontaneous folding process is not well understood at the molecular level. Recent findings suggest that helix fraying and loop formation may be important for overall structure, dynamics and regulation of function. Several types of membrane helices with ionizable amino acids change their topology with pH. Additionally we note that some peptides, including many that are rich in arginine, and a particular analogue of gramicidin, are able passively to translocate across cell membranes. The findings indicate that a final protein structure in a lipid-bilayer membrane is sequence-based, with lipids contributing to stability and regulation. While much progress has been made toward understanding the folding process for alpha-helical membrane proteins, it remains a work in progress. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Refining amino acid hydrophobicity for dynamics simulation of membrane proteins

Coarse-grained (CG) models have been successful in simulating the chemical properties of lipid bilayers, but accurate treatment of membrane proteins and lipid-protein molecular interactions remains a challenge. The CgProt force field, original developed with the multiscale coarse graining method, is assessed by comparing the potentials of mean force for sidechain insertion in a DOPC bilayer to results reported for atomistic molecular dynamics simulations. Reassignment of select CG sidechain sites from the apolar to polar site type was found to improve the attractive interfacial behavior of tyrosine, phenylalanine and asparagine as well as charged lysine and arginine residues. The solvation energy at membrane depths of 0, 1.3 and 1.7 nm correlates with experimental partition coefficients in aqueous mixtures of cyclohexane, octanol and POPC, respectively, for sidechain analogs and Wimley-White peptides. These experimental values serve as important anchor points in choosing between alternate CG models based on their observed permeation profiles, particularly for Arg, Lys and Gln residues where the all-atom OPLS solvation energy does not agree well with experiment. Available partitioning data was also used to reparameterize the representation of the peptide backbone, which needed to be made less attractive for the bilayer hydrophobic core region. The newly developed force field, CgProt 2.4, correctly predicts the global energy minimum in the potentials of mean force for insertion of the uncharged membrane-associated peptides LS3 and WALP23. CgProt will find application in studies of lipid-protein interactions and the conformational properties of diverse membrane protein systems.