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. cruzi tetratricopeptide repeat domain (TPR) of the cytoplasmic peroxisomal targeting signal 1 (PTS1) receptor Pex5. The structure of the TPR domain of TbPex5 represents an apo-form of the receptor which, together with the previously determined structure of the complex of TbPex5 TPR and PTS1 demonstrate significant receptor dynamics associated with signal peptide recognition. The structure of the complex of TPR domain of TcPex5 with PTS1 provided in this study details the molecular interactions that guide signal peptide recognition at the atomic level in the pathogenic species currently perceived as the most relevant among Trypanosoma. Small - angle X - ray scattering (SAXS) data obtained in solution supports the crystallographic findings on the compaction of the TPR domains of TbPex5 and TcPex5 upon interaction with the cargo.

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 (BS3)) 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, BS3, 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. Therefore, the flat-bottom Lorentz-like potentials are recommended to represent cross-link restraints. It was also found that significant improvement of model quality upon the introduction of cross-link restraints is obtained when the sum of differences of indices of cross-linked residues exceeds 150.

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. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

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. We extend the treatment to split the residue-based torsional potentials into the site-based regular and improper torsional potentials. These considerations are illustrated with the revised torsional potentials and improper-torsional potentials involving the l-alanine residue and the improper-torsional potential corresponding to serine-residue enantiomerization. Applications of the new approach in coarse-grained modeling and revising all-atom force fields are discussed.

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 Salmonella (PDB: 6SD5), and the V-type Na+-ATPase rotor in Enterococcus hirae (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. Our results demonstrate that rotational motion can occur with zero total angular momentum in the microcanonical regime and that thermal motions can be converted into net rotation in the canonical regime, as previously observed in simulations of smaller cyclic molecules. For 6SD5 and 2BL2, net rotation (with a ratcheting pattern) occurring only about the pivot of the respective system was observed in canonical simulations. The extent and direction of the rotation depended on the initial conditions. This result suggests that rotatory molecular motors can convert thermal oscillations into net rotational motion. The energy from ATP hydrolysis is required probably to set the direction and extent of rotation. Our findings highlight the importance of molecular-motor structures in facilitating movement and transportation within living organisms.

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). Due to the averaging over the fine-grain degrees of freedom, 1 time unit of UNRES simulations is equivalent to about 1000 time units of laboratory time; therefore, millisecond time scale of large protein systems can be reached with the UNRES-GPU code. AVAILABILITY AND IMPLEMENTATION: The source code of UNRES-GPU along with the benchmarks used for tests is available at https://projects.task.gda.pl/eurohpcpl-public/unres.

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.1 Norwalk strains of human noroviruses (HuNoV) that cause gastroenteritis was investigated by means of long-time (about 30 µs in the laboratory timescale) molecular dynamics simulations with the coarse-grained UNRES force field. The main motion of VLP units turned out to be the bending at the junction between the P1 subdomain (that sits in the VLP shell) and the P2 subdomain (that protrudes outside) of the major VP1 protein, this resulting in a correlated wagging motion of the P2 subdomains with respect to the VLP surface. The fluctuations of the P2 subdomain were found to be more pronounced and the P2 domain made a greater angle with the normal to the VLP surface for the GII.2 strain, which could explain the inability of this strain to bind the histo-blood group antigens (HBGAs).

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. Parallel performance of the code and comparison of its performance with that of AMBER, GROMACS and MARTINI 3 is presented.

Exploring the inhibitory potential of in silico-designed small peptides on Helicobacter pylori Hp0231 (DsbK), a periplasmic oxidoreductase involved in disulfide bond formation.

Introduction: Helicobacter pylori is a bacterium that colonizes the gastric epithelium, which affects millions of people worldwide. H. pylori infection can lead to various gastrointestinal diseases, including gastric adenocarcinoma and mucosa-associated lymphoid tissue lymphoma. Conventional antibiotic therapies face challenges due to increasing antibiotic resistance and patient non-compliance, necessitating the exploration of alternative treatment approaches. In this study, we focused on Hp0231 (DsbK), an essential component of the H. pylori Dsb (disulfide bond) oxidative pathway, and investigated peptide-based inhibition as a potential therapeutic strategy. Methods: Three inhibitory peptides designed by computational modeling were evaluated for their effectiveness using a time-resolved fluorescence assay. We also examined the binding affinity between Hp0231 and the peptides using microscale thermophoresis. Results and discussion: Our findings demonstrate that in silico-designed synthetic peptides can effectively inhibit Hp0231-mediated peptide oxidation. Targeting Hp0231 oxidase activity could attenuate H. pylori virulence without compromising bacterial viability. Therefore, peptide-based inhibitors of Hp0231 could be candidates for the development of new targeted strategy, which does not influence the composition of the natural human microbiome, but deprive the bacterium of its pathogenic properties.

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. The remaining seven were synthesized, and their properties were measured in this work. Multiplexed replica-exchange simulations of eight-chain systems were conducted at 12 temperatures from 260 to 370 K at concentrations from 0.421 to 5.78 mM, corresponding to the experimental conditions. The temperature profiles of the fractions of monomers and octamers showed a clear transition corresponding to aggregate dissociation. Low simulated transition temperatures were obtained for the peptides, which did not precipitate after incubation, as well as for the H-GNNQQNY-NH2 prion-protein fragment, which forms small fibrils. A substantial amount of inter-strand ß-sheets was found in most of the systems. The results suggest that UNRES simulations can be used to assess peptide aggregation except for glutamine- and asparagine-rich peptides, for which a revision of the UNRES sidechain-sidechain interaction potentials appears necessary.

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. In this manuscript, we combine knowledge of the infection mechanism with theoretical methods to answer the question of the virus's selectivity. We proposed a two-stage infection mechanism. In the first step, the virus interacts with the ACE2 receptor, with the "proper strength". When the interaction is too strong, the virus will remain in an "improper position"; if the interaction is too weak, the virus will "run away" from the cell. We also indicated three residues (positions 30, 31, and 353) located on the ACE2 protein-binding interface, which seems to be crucial for successful infection. Our results indicate that these residues are necessary for the initiation of the infection process.

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.0 program and the GAFF force field. The potentials of mean force of hydrophobic dimer formation, as well as the excess free energy, excess energy, excess entropy, and excess heat capacity corresponding to the formation of the contact minimum, were determined and analyzed. For both systems, the depth of the contact minimum in the potential of mean force was found to increase with both temperature and ionic strength. The excess heat capacity of the association at the contact minimum and T = 298 K was found to be negative and to decrease, while the excess entropy and energy were found to be positive and to increase for both systems, the changes being more pronounced for the hexane dimer. The excess heat capacity is also greater in absolute value for the hexane dimer.

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.edu.pl) for physics-based simulations with the coarse-grained UNRES model of polypeptide chains. The improvements include the replacement of the old code with the recently optimized one and adding the recent scale-consistent variant of the UNRES force field, which performs better in the modeling of proteins with the ß and the α+ß structures. The scope of applications of the package was extended to data-assisted simulations with restraints from nuclear magnetic resonance (NMR) and chemical crosslink mass-spectroscopy (XL-MS) measurements. NMR restraints can be input in the NMR Exchange Format (NEF), which has become a standard. Ambiguous NMR restraints are handled without expert intervention owing to a specially designed penalty function. The server can be used to run smaller jobs directly or to prepare input data to run larger production jobs by using standalone installations of UNRES.

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. The force field has been derived and parameterized based on the principles of statistical-mechanics, which makes it independent of structural databases and applicable to treat nonstandard situations such as, proteins that contain D-amino-acid residues. Powered by Langevin dynamics and its replica-exchange extensions, UNRES has found a variety of applications, including ab initio and database-assisted protein-structure prediction, simulating protein-folding pathways, exploring protein free-energy landscapes, and solving biological problems. This chapter provides a summary of UNRES and a guide for potential users regarding the application of the UNRES package in a variety of research tasks.

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. Consequently, Langevin dynamics is a natural means of describing the motion of a system at the coarse-grained level and the potential of mean force is the physical basis of the coarse-grained force fields. Moreover, the choice of coarse-grained variables and the fact that coarse-grained sites often do not have spherical symmetry implies a non-diagonal inertia tensor. We describe selected coarse-grained models used in molecular dynamics simulations, including the most popular MARTINI model developed by Marrink's group and the UNICORN model of biological macromolecules developed in our laboratory. We conclude by discussing examples of the application of coarse-grained molecular dynamics to study biologically important processes.

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. Our consensus-fragment procedure was used to run template-assisted predictions. Each group also processed the Nuclear Magnetic Resonance (NMR)- and Small Angle X-Ray Scattering (SAXS)-data assisted targets. The average Global Distance Test Total Score (GDT_TS) of the 'Model 1' predictions were 29.17, 39.32, and 56.37 for the UNRES, UNRES-contact, and UNRES-template predictions, respectively, increasing by 0.53, 2.24, and 3.76, respectively, compared to CASP13. It was also found that the GDT_TS of the UNRES models obtained in ab initio mode and in the contact-assisted mode decreases with the square root of chain length, while the exponent in this relationship is 0.20 for the UNRES-template group models and 0.11 for the best performing AlphaFold2 models, which suggests that incorporation of database information, which stems from protein evolution, brings in long-range correlations, thus enabling the correction of force-field inaccuracies.

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). Eight were difficult targets for which only distantly related templates were found for the individual subunits. Twenty-five CAPRI groups including eight automatic servers submitted ~1250 models per target. Twenty groups including six servers participated in the CAPRI scoring challenge submitted ~190 models per target. The accuracy of the predicted models was evaluated using the classical CAPRI criteria. The prediction performance was measured by a weighted scoring scheme that takes into account the number of models of acceptable quality or higher submitted by each group as part of their five top-ranking models. Compared to the previous CASP-CAPRI challenge, top performing groups submitted such models for a larger fraction (70-75%) of the targets in this Round, but fewer of these models were of high accuracy. Scorer groups achieved stronger performance with more groups submitting correct models for 70-80% of the targets or achieving high accuracy predictions. Servers performed less well in general, except for the MDOCKPP and LZERD servers, who performed on par with human groups. In addition to these results, major advances in methodology are discussed, providing an informative overview of where the prediction of protein assemblies currently stands.

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. RESULTS: We present a new protein-protein and peptide-protein docking functionality of the UNRES package, which includes a variable degree of conformational flexibility. UNRES-Dock protocol was tested on a set of 55 complexes with size from 43 to 587 amino-acid residues, showing that structures of the complexes can be predicted with good quality, if the sampling of the conformational space is sufficient, especially for flexible peptide-protein systems. The developed automatized protocol has been implemented in the standalone UNRES package and in the UNRES server. AVAILABILITY AND IMPLEMENTATION: UNRES server: http://unres-server.chem.ug.edu.pl; UNRES package and data used in testing of UNRES-Dock: http://unres.pl. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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. Collectively, the conformational dynamics enables these proteins to fulfill specific functions. Thus, most experimental observables are averaged over the conformations that constitute an ensemble. In this article, we review the recent developments in the concept and methods for the determination of the dynamic structures of flexible peptides and proteins. In particular, we describe ways to extract information from nuclear magnetic resonance small-angle X-ray scattering (SAXS), and chemical cross-linking coupled with mass spectroscopy (XL-MS) measurements. All these techniques can be used to obtain ensemble-averaged restraints or to re-weight the simulated conformational ensembles.

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. Moreover, analytical formulas for the multibody (correlation) terms, which appear to be crucial for coarse-grained modeling of many of the regular structures such as, e.g., protein α-helices and ß-sheets, can be derived in a systematic way. Implementation of the developed theory to the UNIfied COarse-gRaiNed (UNICORN) model of biological macromolecules, which consists of the UNRES (for proteins), NARES-2P (for nucleic acids), and SUGRES-1P (for polysaccharides) components, and is being developed in our laboratory is described. Successful applications of UNICORN to the prediction of protein structure, simulating the folding and stability of proteins and nucleic acids, and solving biological problems are discussed.