Implementation of Time-Averaged Restraints with UNRES Coarse-Grained Model of Polypeptide Chains.

Time-averaged restraints from nuclear magnetic resonance (NMR) measurements have been implemented in the UNRES coarse-grained model of polypeptide chains in order to develop a tool for data-assisted modeling of the conformational ensembles of multistate proteins, intrinsically disordered proteins (IDPs) and proteins with intrinsically disordered regions (IDRs), many of which are essential in cell biology. A numerically stable variant of molecular dynamics with time-averaged restraints has been introduced, in which the total energy is conserved in sections of a trajectory in microcanonical runs, the bath temperature is maintained in canonical runs, and the time-average-restraint-force components are scaled up with the length of the memory window so that the restraints affect the simulated structures. The new approach restores the conformational ensembles used to generate ensemble-averaged distances, as demonstrated with synthetic restraints.

Secondary Structure in Free and Assisted Modeling of Proteins with the Coarse-Grained UNRES Force Field.

Secondary structure is a solid scaffold on which the three-dimensional structure of a protein is built. Therefore, care must be taken to reproduce the secondary structure as accurately as possible in the simulations of protein systems. In this chapter, we summarize the physics-based energy terms that govern secondary-structure formation, the auxiliary restraints on secondary structure derived from bioinformatics and from the experimental data, and the role of those in the modeling of protein structures, dynamics, and thermodynamics with the physics-based coarse-grained UNRES force field.

Mean-Field Coupling Between Local Interactions in Proteins in Relation to Chirality, Secondary, and Supersecondary Structure Formation, and Allostery.

Coarse graining is usually considered as a tool to extend the time and size scale of simulations. However, leaving out the atomistic details to keep their fingerprints in a coarse-grained model also enables us to understand better structure formation and dynamics. In this chapter, by using our scale-consistent theory of coarse graining, we demonstrate that the coarse-grained terms corresponding to the coupling between local conformational states of amino-acid residues explain secondary-structure propagation along polypeptide backbone to stabilize -helices and -strands in proteins and direct the loops preceding and following such segments of protein structure.

Structural dynamics of the TPR domain of the peroxisomal cargo receptor Pex5 in Trypanosoma.

Peroxisomal protein import has been identified as a valid target in trypanosomiases, an important health threat in Central and South America. The importomer is built of multiple peroxins (Pex) and structural characterization of these proteins facilitates rational inhibitor development. We report crystal structures of the Trypanosoma brucei and T.

Assessment of Four Theoretical Approaches to Predict Protein Flexibility in the Crystal Phase and Solution.

In this paper, we evaluated the ability of four coarse-grained methods to predict protein flexible regions with potential biological importance, UNRES-flex, UNRES-DSSP-flex (based on the united residue model of polypeptide chains without and with secondary structure restraints, respectively), CABS-flex (based on the C-α, C-β, and side chain model), and nonlinear rigid block normal mode analysis (NOLB) with a set of 100 protein structures determined by NMR spectroscopy or X-ray crystallography, with all secondary structure types. End regions with high fluctuations were excluded from analysis. The Pearson and Spearman correlation coefficients were used to quantify the conformity between the calculated and experimental fluctuation profiles, the latter determined from NMR ensembles and X-ray -factors, respectively.

Assessment of Two Restraint Potentials for Coarse-Grained Chemical-Cross-Link-Assisted Modeling of Protein Structures.

The influence of distance restraints from chemical cross-link mass spectroscopy (XL-MS) on the quality of protein structures modeled with the coarse-grained UNRES force field was assessed by using a protocol based on multiplexed replica exchange molecular dynamics, in which both simulated and experimental cross-link restraints were employed, for 23 small proteins. Six cross-links with upper distance boundaries from 4 Å to 12 Å (azido benzoic acid succinimide (ABAS), triazidotriazine (TATA), succinimidyldiazirine (SDA), disuccinimidyl adipate (DSA), disuccinimidyl glutarate (DSG), and disuccinimidyl suberate (BS)) and two types of restraining potentials ((i) simple flat-bottom Lorentz-like potentials dependent on side chain distance (all cross-links) and (ii) distance- and orientation-dependent potentials determined based on molecular dynamics simulations of model systems (DSA, DSG, BS, and SDA)) were considered. The Lorentz-like potentials with properly set parameters were found to produce a greater number of higher-quality models compared to unrestrained simulations than the MD-based potentials, because the latter can force too long distances between side chains.

Impact of AlphaFold on structure prediction of protein complexes: The CASP15-CAPRI experiment.

We present the results for CAPRI Round 54, the 5th joint CASP-CAPRI protein assembly prediction challenge. The Round offered 37 targets, including 14 homodimers, 3 homo-trimers, 13 heterodimers including 3 antibody-antigen complexes, and 7 large assemblies. On average ~70 CASP and CAPRI predictor groups, including more than 20 automatics servers, submitted models for each target.

Toward Consistent Physics-Based Modeling of Local Backbone Structures and Chirality Change of Proteins in Coarse-Grained Approaches.

A reliable representation of local interactions is critical for the accuracy of modeling protein structure and dynamics at both the all-atom and coarse-grained levels. The development of local (mainly torsional) potentials was focused on careful parametrization of the predetermined (usually Fourier) formulas rather than on their physics-based derivation. In this Perspective we discuss the state-of-the-art methods for modeling local interactions, including the scale-consistent theory developed in our laboratory, which implies that the coarse-grained torsional potentials inseparably depend on the virtual-bond angles adjacent to a given dihedral and that multitorsional terms should be considered.

Pragmatic Coarse-Graining of Proteins: Models and Applications.

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail.

Long-Time Dynamics of Selected Molecular-Motor Components Using a Physics-Based Coarse-Grained Approach.

Molecular motors are essential for the movement and transportation of macromolecules in living organisms. Among them, rotatory motors are particularly efficient. In this study, we investigated the long-term dynamics of the designed left-handed alpha/alpha toroid (PDB: 4YY2), the RBM2 flagellum protein ring from (PDB: 6SD5), and the V-type Na+-ATPase rotor in (PDB: 2BL2) using microcanonical and canonical molecular dynamics simulations with the coarse-grained UNRES force field, including a lipid-membrane model, on a millisecond laboratory time scale.

UNRES-GPU for physics-based coarse-grained simulations of protein systems at biological time- and size-scales.

Summary: The UNited RESisdue (UNRES) package for coarse-grained simulations, which has recently been optimized to treat large protein systems, has been implemented on Graphical Processor Units (GPUs). An over 100-time speed-up of the GPU code (run on an NVIDIA A100) with respect to the sequential code and an 8.5 speed-up with respect to the parallel Open Multi-Processing (OpenMP) code (run on 32 cores of 2 AMD EPYC 7313 Central Processor Units (CPUs)) has been achieved for large proteins (with size over 10 000 residues).

Long-time scale simulations of virus-like particles from three human-norovirus strains.

The dynamics of the virus like particles (VLPs) corresponding to the GII.4 Houston, GII.2 SMV, and GI.

Improvements and new functionalities of UNRES server for coarse-grained modeling of protein structure, dynamics, and interactions.

In this paper we report the improvements and extensions of the UNRES server (https://unres-server.chem.ug.

Optimization of parallel implementation of UNRES package for coarse-grained simulations to treat large proteins.

We report major algorithmic improvements of the UNRES package for physics-based coarse-grained simulations of proteins. These include (i) introduction of interaction lists to optimize computations, (ii) transforming the inertia matrix to a pentadiagonal form to reduce computing and memory requirements, (iii) removing explicit angles and dihedral angles from energy expressions and recoding the most time-consuming energy/force terms to minimize the number of operations and to improve numerical stability, (iv) using OpenMP to parallelize those sections of the code for which distributed-memory parallelization involves unfavorable computing/communication time ratio, and (v) careful memory management to minimize simultaneous access of distant memory sections. The new code enables us to run molecular dynamics simulations of protein systems with size exceeding 100,000 amino-acid residues, reaching over 1 ns/day (1 μs/day in all-atom timescale) with 24 cores for proteins of this size.

Prediction of Aggregation of Biologically-Active Peptides with the UNRES Coarse-Grained Model.

The UNited RESidue (UNRES) model of polypeptide chains was applied to study the association of 20 peptides with sizes ranging from 6 to 32 amino-acid residues. Twelve of those were potentially aggregating hexa- or heptapeptides excised from larger proteins, while the remaining eight contained potentially aggregating sequences, functionalized by attaching larger ends rich in charged residues. For 13 peptides, the experimental data of aggregation were used.

Theoretical Investigation of the Coronavirus SARS-CoV-2 (COVID-19) Infection Mechanism and Selectivity.

The SARS-CoV-2 virus, commonly known as COVID-19, first occurred in December 2019 in Wuhan, Hubei Province, China. Since then, it has become a tremendous threat to human health. With a pandemic threat, it is in the significant interest of the scientific world to establish its method of infection.

Influence of Temperature and Salt Concentration on the Hydrophobic Interactions of Adamantane and Hexane.

One of the definitions of hydrophobic interactions is the aggregation of nonpolar particles in a polar solvent, such as water. While this phenomenon appears to be very simple, it is crucial for many complex processes, such as protein folding, to take place. In this work, the hydrophobic association of adamantane and hexane at various temperatures and ionic strengths was studied using molecular dynamics simulations with the AMBER 16.

Recent Developments in Data-Assisted Modeling of Flexible Proteins.

Many proteins can fold into well-defined conformations. However, intrinsically-disordered proteins (IDPs) do not possess a defined structure. Moreover, folded multi-domain proteins often digress into alternative conformations.

Modeling the Structure, Dynamics, and Transformations of Proteins with the UNRES Force Field.

The physics-based united-residue (UNRES) model of proteins ( www.unres.pl ) has been designed to carry out large-scale simulations of protein folding.

Theory and Practice of Coarse-Grained Molecular Dynamics of Biologically Important Systems.

Molecular dynamics with coarse-grained models is nowadays extensively used to simulate biomolecular systems at large time and size scales, compared to those accessible to all-atom molecular dynamics. In this review article, we describe the physical basis of coarse-grained molecular dynamics, the coarse-grained force fields, the equations of motion and the respective numerical integration algorithms, and selected practical applications of coarse-grained molecular dynamics. We demonstrate that the motion of coarse-grained sites is governed by the potential of mean force and the friction and stochastic forces, resulting from integrating out the secondary degrees of freedom.

Prediction of protein assemblies, the next frontier: The CASP14-CAPRI experiment.

We present the results for CAPRI Round 50, the fourth joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of twelve targets, including six dimers, three trimers, and three higher-order oligomers. Four of these were easy targets, for which good structural templates were available either for the full assembly, or for the main interfaces (of the higher-order oligomers).

Modeling protein structures with the coarse-grained UNRES force field in the CASP14 experiment.

The UNited RESidue (UNRES) force field was tested in the 14th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP14), in which larger oligomeric and multimeric targets were present compared to previous editions. Three prediction modes were tested (i) ab initio (the UNRES group), (ii) contact-assisted (the UNRES-contact group), and (iii) template-assisted (the UNRES-template group). For most of the targets, the contact restraints were derived from the server models top-ranked by the DeepQA method, while the DNCON2 method was used for 11 targets.

UNRES-Dock-protein-protein and peptide-protein docking by coarse-grained replica-exchange MD simulations.

Motivation: The majority of the proteins in living organisms occur as homo- or hetero-multimeric structures. Although there are many tools to predict the structures of single-chain proteins or protein complexes with small ligands, peptide-protein and protein-protein docking is more challenging. In this work, we utilized multiplexed replica-exchange molecular dynamics (MREMD) simulations with the physics-based heavily coarse-grained UNRES model, which provides more than a 1000-fold simulation speed-up compared with all-atom approaches to predict structures of protein complexes.

Scale-consistent approach to the derivation of coarse-grained force fields for simulating structure, dynamics, and thermodynamics of biopolymers.

In this chapter the scale-consistent approach to the derivation of coarse-grained force fields developed in our laboratory is presented, in which the effective energy function originates from the potential of mean force of the system under consideration and embeds atomistically detailed interactions in the resulting energy terms through use of Kubo's cluster-cumulant expansion, appropriate selection of the major degrees of freedom to be averaged out in the derivation of analytical approximations to the energy terms, and appropriate expression of the interaction energies at the all-atom level in these degrees of freedom. Our approach enables the developers to find correct functional forms of the effective coarse-grained energy terms, without having to import them from all-atom force fields or deriving them on a heuristic basis. In particular, the energy terms derived in such a way exhibit correct dependence on coarse-grained geometry, in particular on site orientation.