Multilayered Computational Framework for Designing Peptide Inhibitors of HVEM-LIGHT Interaction.

The herpesvirus entry mediator (HVEM) and its ligand LIGHT play crucial roles in immune system regulation, including T-cell proliferation, B-cell differentiation, and immunoglobulin secretion. However, excessive T-cell activation can lead to chronic inflammation and autoimmune diseases. Thus, inhibiting the HVEM-LIGHT interaction emerges as a promising therapeutic strategy for these conditions and in preventing adverse reactions in organ transplantation. This study focused on designing peptide inhibitors, targeting the HVEM-LIGHT interaction, using molecular dynamics (MD) simulations of 65 peptides derived from HVEM. These peptides varied in length and disulfide-bond configurations, crucial for their interaction with the LIGHT trimer. By simulating 31 HVEM domain variants, including the full-length protein, we assessed conformational changes upon LIGHT binding to understand the influence of HVEM segments and disulfide bonds on the binding mechanism. Employing multitrajectory microsecond-scale, all-atom MD simulations and molecular mechanics with generalized Born and surface area (MM-GBSA) binding energy estimation, we identified promising CRD2 domain variants with high LIGHT affinity. Notably, point mutations in these variants led to a peptide with a single disulfide bond (C58-C73) and a K54E substitution, exhibiting the highest binding affinity. The importance of the CRD2 domain and Cys58-Cys73 disulfide bond for interrupting HVEM-LIGHT interaction was further supported by analyzing truncated CRD2 variants, demonstrating similar binding strengths and mechanisms. Further investigations into the binding mechanism utilized steered MD simulations at various pulling speeds and umbrella sampling to estimate the energy profile of HVEM-based inhibitors with LIGHT. These comprehensive analyses revealed key interactions and different binding mechanisms, highlighting the increased binding affinity of selected peptide variants. Experimental circular dichroism techniques confirmed the structural properties of these variants. This study not only advances our understanding of the molecular basis of HVEM-LIGHT interactions but also provides a foundation for developing novel therapeutic strategies for immune-related disorders. Furthermore, it sets a gold standard for peptide inhibitor design in drug development due to its systematic approach.

Free-Docking and Template-Based Docking: Physics Versus Knowledge-Based Docking.

Docking methods can be used to predict the orientations of two or more molecules with respect of each other using a plethora of various algorithms, which can be based on the physics of interactions or can use information from databases and templates. The usability of these approaches depends on the type and size of the molecules, whose relative orientation will be estimated. The two most important limitations are (i) the computational cost of the prediction and (ii) the availability of the structural information for similar complexes. In general, if there is enough information about similar systems, knowledge-based and template-based methods can significantly reduce the computational cost while providing high accuracy of the prediction. However, if the information about the system topology and interactions between its partners is scarce, physics-based methods are more reliable or even the only choice. In this chapter, knowledge-, template-, and physics-based methods will be compared and briefly discussed providing examples of their usability with a special emphasis on physics-based protein-protein, protein-peptide, and protein-fullerene docking in the UNRES coarse-grained model.

A Practical Guide to All-Atom and Coarse-Grained Molecular Dynamics Simulations Using Amber and Gromacs: A Case Study of Disulfide-Bond Impact on the Intrinsically Disordered Amyloid Beta.

Intrinsically disordered proteins (IDPs) pose challenges to conventional experimental techniques due to their large-scale conformational fluctuations and transient structural elements. This work presents computational methods for studying IDPs at various resolutions using the Amber and Gromacs packages with both all-atom (Amber ff19SB with the OPC water model) and coarse-grained (Martini 3 and SIRAH) approaches. The effectiveness of these methodologies is demonstrated by examining the monomeric form of amyloid-ß (Aß42), an IDP, with and without disulfide bonds at different resolutions. Our results clearly show that the addition of a disulfide bond decreases the ß-content of Aß42; however, it increases the tendency of the monomeric Aß42 to form fibril-like conformations, explaining the various aggregation rates observed in experiments. Moreover, analysis of the monomeric Aß42 compactness, secondary structure content, and comparison between calculated and experimental chemical shifts demonstrates that all three methods provide a reasonable choice to study IDPs; however, coarse-grained approaches may lack some atomistic details, such as secondary structure recognition, due to the simplifications used. In general, this study not only explains the role of disulfide bonds in Aß42 but also provides a step-by-step protocol for setting up, conducting, and analyzing molecular dynamics (MD) simulations, which is adaptable for studying other biomacromolecules, including folded and disordered proteins and peptides.

Exploring Structural Insights of Aß42 and α-Synuclein Monomers and Heterodimer: A Comparative Study Using Implicit and Explicit Solvent Simulations.

Protein misfolding, aggregation, and fibril formation play a central role in the development of severe neurological disorders, including Alzheimer's and Parkinson's diseases. The structural stability of mature fibrils in these diseases is of great importance, as organisms struggle to effectively eliminate amyloid plaques. To address this issue, it is crucial to investigate the early stages of fibril formation when monomers aggregate into small, toxic, and soluble oligomers. However, these structures are inherently disordered, making them challenging to study through experimental approaches. Recently, it has been shown experimentally that amyloid-ß 42 (Aß42) and α-synuclein (α-Syn) can coassemble. This has motivated us to investigate the interaction between their monomers as a first step toward exploring the possibility of forming heterodimeric complexes. In particular, our study involves the utilization of various Amber and CHARMM force-fields, employing both implicit and explicit solvent models in replica exchange and conventional simulation modes. This comprehensive approach allowed us to assess the strengths and weaknesses of these solvent models and force fields in comparison to experimental and theoretical findings, ensuring the highest level of robustness. Our investigations revealed that Aß42 and α-Syn monomers can indeed form stable heterodimers, and the resulting heterodimeric model exhibits stronger interactions compared to the Aß42 dimer. The binding of α-Syn to Aß42 reduces the propensity of Aß42 to adopt fibril-prone conformations and induces significant changes in its conformational properties. Notably, in AMBER-FB15 and CHARMM36m force fields with the use of explicit solvent, the presence of Aß42 significantly increases the ß-content of α-Syn, consistent with the experiments showing that Aß42 triggers α-Syn aggregation. Our analysis clearly shows that although the use of implicit solvent resulted in too large compactness of monomeric α-Syn, structural properties of monomeric Aß42 and the heterodimer were preserved in explicit-solvent simulations. We anticipate that our study sheds light on the interaction between α-Syn and Aß42 proteins, thus providing the atom-level model required to assess the initial stage of aggregation mechanisms related to Alzheimer's and Parkinson's diseases.

Integrating Explicit and Implicit Fullerene Models into UNRES Force Field for Protein Interaction Studies.

Fullerenes, particularly C60, exhibit unique properties that make them promising candidates for various applications, including drug delivery and nanomedicine. However, their interactions with biomolecules, especially proteins, remain not fully understood. This study implements both explicit and implicit C60 models into the UNRES coarse-grained force field, enabling the investigation of fullerene-protein interactions without the need for restraints to stabilize protein structures. The UNRES force field offers computational efficiency, allowing for longer timescale simulations while maintaining accuracy. Five model proteins were studied: FK506 binding protein, HIV-1 protease, intestinal fatty acid binding protein, PCB-binding protein, and hen egg-white lysozyme. Molecular dynamics simulations were performed with and without C60 to assess protein stability and investigate the impact of fullerene interactions. Analysis of contact probabilities reveals distinct interaction patterns for each protein. FK506 binding protein (1FKF) shows specific binding sites, while intestinal fatty acid binding protein (1ICN) and uteroglobin (1UTR) exhibit more generalized interactions. The explicit C60 model shows good agreement with all-atom simulations in predicting protein flexibility, the position of C60 in the binding pocket, and the estimation of effective binding energies. The integration of explicit and implicit C60 models into the UNRES force field, coupled with recent advances in coarse-grained modeling and multiscale approaches, provides a powerful framework for investigating protein-nanoparticle interactions at biologically relevant scales without the need to use restraints stabilizing the protein, thus allowing for large conformational changes to occur. These computational tools, in synergy with experimental techniques, can aid in understanding the mechanisms and consequences of nanoparticle-biomolecule interactions, guiding the design of nanomaterials for biomedical applications.

Serotonin Promotes Vesicular Association and Fusion by Modifying Lipid Bilayers.

The primary event in chemical neurotransmission involves the fusion of a membrane-limited vesicle at the plasma membrane and the subsequent release of its chemical neurotransmitter cargo. The cargo itself is not known to have any effect on the fusion event. However, amphiphilic monoamine neurotransmitters (e.g., serotonin and dopamine) are known to strongly interact with lipid bilayers and to affect their mechanical properties, which can in principle impact membrane-mediated processes. Here, we probe whether serotonin can enhance the association and fusion of artificial lipid vesicles in vitro. We employ fluorescence correlation spectroscopy and total internal reflection fluorescence microscopy to measure the attachment and fusion of vesicles whose lipid compositions mimic the major lipid components of synaptic vesicles. We find that the association between vesicles and supported lipid bilayers is strongly enhanced in a serotonin dose-dependent manner, and this drives an increase in the rate of spontaneous fusion. Molecular dynamics simulations and fluorescence spectroscopy data show that serotonin insertion increases the water content of the hydrophobic part of the bilayer. This suggests that the enhanced membrane association is likely driven by an energetically favorable drying transition. Other monoamines, such as dopamine and norepinephrine, but not other related species, such as tryptophan, show similar effects on membrane association. Our results reveal a lipid bilayer-mediated mechanism by which monoamines can themselves modulate vesicle fusion, potentially adding to the control toolbox for the tightly regulated process of neurotransmission in vivo.

Amyloid-ß Tetramers and Divalent Cations at the Membrane/Water Interface: Simple Models Support a Functional Role.

Charge polarization at the membrane interface is a fundamental process in biology. Despite the lower concentration compared to the abundant monovalent ions, the relative abundance of divalent cations (Ca2+, Mg2+, Zn2+, Fe2+, Cu2+) in particular spaces, such as the neuron synapse, raised many questions on the possible effects of free multivalent ions and of the required protection of membranes by the eventual defects caused by the free forms of the cations. In this work, we first applied a recent realistic model of divalent cations to a well-investigated model of a polar lipid bilayer, di-myristoyl phosphatidyl choline (DMPC). The full atomistic model allows a fairly good description of changes in the hydration of charged and polar groups upon the association of cations to lipid atoms. The lipid-bound configurations were analyzed in detail. In parallel, amyloid-ß 1-42 (Aß42) peptides assembled into tetramers were modeled at the surface of the same bilayer. Two of the protein tetramers' models were loaded with four Cu2+ ions, the latter bound as in DMPC-free Aß42 oligomers. The two Cu-bound models differ in the binding topology: one with each Cu ion binding each of the monomers in the tetramer; one with pairs of Cu ions linking two monomers into dimers, forming tetramers as dimers of dimers. The models here described provide hints on the possible role of Cu ions in synaptic plasticity and of Aß42 oligomers in storing the same ions away from lipids. The release of structurally disordered peptides in the synapse can be a mechanism to recover ion homeostasis and lipid membranes from changes in the divalent cation concentration.

Why does the herpes simplex 1 virus-encoded UL49.5 protein fail to inhibit the TAP-dependent antigen presentation?

Herpes simplex virus 1 (HSV-1) is a well-studied herpesvirus that causes various human diseases. Like other herpesviruses, HSV-1 produces the transmembrane glycoprotein N (gN/UL49.5 protein), which has been extensively studied, but its function in HSV-1 remains largely unknown. The amino-acid sequences and lengths of UL49.5 proteins differ between herpesvirus species. It is, therefore, crucial to determine whether and to what extent the spatial structure of UL49.5 orthologs that are transporter associated with antigen processing (TAP) inhibitors (i.e., of bovine herpesvirus 1; BoHV-1) differ from that of non-TAP inhibitors (i.e., of HSV-1). Our study aimed to examine the 3D structure of the HSV-1-encoded UL49.5 protein in an advanced model of the endoplasmic reticulum (ER) membrane using circular dichroism, 2D nuclear magnetic resonance, and multiple-microsecond all-atom molecular dynamics simulations in an ER membrane mimetic environment. According to our findings, the N-terminus of the HSV-1-encoded UL49.5 adopts a highly flexible, unordered structure in the extracellular part due to the presence of a large number of proline and glycine residues. In contrast to the BoHV-1-encoded homolog, the transmembrane region of the HSV-1-encoded UL49.5 is formed by a single long transmembrane α-helix, rather than two helices oriented perpendicularly, while the cytoplasmic part of the protein (C-terminus) has a short unordered structure. Our findings provide valuable experimental structural information on the HSV-1-encoded UL49.5 protein and offer, based on the obtained structure, insight into its lack of biological activity in inhibiting the TAP-dependent antigen presentation pathway.

Unusual Robustness of Neurotransmitter Vesicle Membranes against Serotonin-Induced Perturbations.

Nature confines hundreds of millimolar of amphiphilic neurotransmitters, such as serotonin, in synaptic vesicles. This appears to be a puzzle, as the mechanical properties of lipid bilayer membranes of individual major polar lipid constituents of synaptic vesicles [phosphatidylcholine (PC), phosphatidylethanolamine (PE), and phosphatidylserine (PS)] are significantly affected by serotonin, sometimes even at few millimolar concentrations. These properties are measured by atomic force microscopy, and their results are corroborated by molecular dynamics simulations. Complementary 2H solid-state NMR measurements also show that the lipid acyl chain order parameters are strongly affected by serotonin. The resolution of the puzzle lies in the remarkably different properties displayed by the mixture of these lipids, at molar ratios mimicking those of natural vesicles (PC:PE:PS:Cholesterol = 3:5:2:5). Bilayers constituting of these lipids are minimally perturbed by serotonin, and show only a graded response at physiological concentrations (>100 mM). Significantly, the cholesterol (up to 33% molar ratio) plays only a minor role in dictating these mechanical perturbations, with PC:PE:PS:Cholesterol = 3:5:2:5 and 3:5:2:0 showing similar perturbations. We infer that nature uses an emergent mechanical property of a specific mixture of lipids, all individually vulnerable to serotonin, to appropriately respond to physiological serotonin levels.

The N-terminal Proline Hinge Motif Controls the Structure of Bovine Herpesvirus 1-encoded Inhibitor of the Transporter Associated with Antigen Processing Required for its Immunomodulatory Function.

Due to unique features, proline residues may control protein structure and function. Here, we investigated the role of 52PPQ54 residues, indicated by the recently established experimental 3D structure of bovine herpesvirus 1-encoded UL49.5 protein as forming a characteristic proline hinge motif in its N-terminal domain. UL49.5 acts as a potent inhibitor of the transporter associated with antigen processing (TAP), which alters the antiviral immune response. Mechanisms employed by UL49.5 to affect TAP remain undetermined on a molecular level. We found that mutations in the 52PPQ54 region had a vast impact on its immunomodulatory function, increasing cell surface MHC class I expression, TAP levels, and peptide transport efficiency. This inhibitory effect was specific for UL49.5 activity towards TAP but not towards the viral glycoprotein M. To get an insight into the impact of proline hinge modifications on structure and dynamics, we performed all-atom and coarse-grained molecular dynamics studies on the native protein and PPQ mutants. The results demonstrated that the proline hinge sequence with its highly rigid conformation served as an anchor into the membrane. This anchor was responsible for the structural and dynamical behavior of the whole protein, constraining the mobility of the C-terminus, increasing the mobility of the transmembrane region, and controlling the accessibility of the C-terminal residues to the cytoplasmic environment. Those features appear crucial for TAP binding and inhibition. Our findings significantly advance the structural understanding of the UL49.5 protein and its functional regions and support the importance of proline motifs for the protein structure.

Mechanical Stability of Ribonuclease A Heavily Depends on the Redox Environment.

Disulfide bonds are covalent bonds that connect nonlocal fragments of proteins, and they are unique post-translational modifications of proteins. They require the oxidizing environment to be stable, which occurs for example during oxidative stress; however, in a cell the reductive environment is maintained, lowering their stability. Despite many years of research on disulfide bonds, their role in the protein life cycle is not fully understood and seems to strictly depend on a system or process in which they are involved. In this article, coarse-grained UNited RESidue (UNRES), and all-atom Assisted Model Building with Energy Refinement (AMBER) force fields were applied to run a series of steered molecular dynamics (SMD) simulations of one of the most studied, but still not fully understood, proteins─ribonuclease A (RNase A). SMD simulations were performed to study the mechanical stability of RNase A in different oxidative-reductive environments. As disulfide bonds (and any other covalent bonds) cannot break/form in any classical all-atom force field, we applied additional restraints between sulfur atoms of reduced cysteines which were able to mimic the breaking of the disulfide bonds. On the other hand, the coarse-grained UNRES force field enables us to study the breaking/formation of the disulfide bonds and control the reducing/oxidizing environment owing to the presence of the designed distance/orientation-dependent potential. This study reveals that disulfide bonds have a strong influence on the mechanical stability of RNase A only in a highly oxidative environment. However, the local stability of the secondary structure seems to play a major factor in the overall stability of the protein. Both our thermal unfolding and mechanical stretching studies show that the most stable disulfide bond is Cys65-Cys72. The breaking of disulfide bonds Cys26-Cys84 and Cys58-Cys110 is associated with large force peaks. They are structural bridges, which are mostly responsible for stabilizing the RNase A conformation, while the presence of the remaining two bonds (Cys65-Cys72 and Cys40-Cys95) is most likely connected with the enzymatic activity rather than the structural stability of RNase A in the cytoplasm. Our results prove that disulfide bonds are indeed stabilizing fragments of the proteins, but their role is strongly redox environment-dependent.

Amyloid ß Dodecamer Disrupts the Neuronal Membrane More Strongly than the Mature Fibril: Understanding the Role of Oligomers in Neurotoxicity.

The amyloid cascade hypothesis states that senile plaques, composed of amyloid ß (Aß) fibrils, play a key role in Alzheimer's disease (AD). However, recent experiments have shown that Aß oligomers are more toxic to neurons than highly ordered fibrils. The molecular mechanism underlying this observation remains largely unknown. One of the possible scenarios for neurotoxicity is that Aß peptides create pores in the lipid membrane that allow Ca2+ ions to enter cells, resulting in a signal of cell apoptosis. Hence, one might think that oligomers are more toxic due to their higher ability to create ion channels than fibrils. In this work, we study the effect of Aß42 dodecamer and fibrils on a neuronal membrane, which is similar to that observed in AD patients, using all-atom molecular dynamics simulations. Due to short simulation times, we cannot observe the formation of pores, but useful insight on the early events of this process has been obtained. Namely, we showed that dodecamer distorts the lipid membrane to a greater extent than fibrils, which may indicate that ion channels can be more easily formed in the presence of oligomers. Based on this result, we anticipate that oligomers are more toxic than mature fibrils, as observed experimentally. Moreover, the Aß-membrane interaction was found to be governed by the repulsive electrostatic interaction between Aß and the ganglioside GM1 lipid. We calculated the bending and compressibility modulus of the membrane in the absence of Aß and obtained good agreement with the experiment. We predict that the dodecamer will increase the compressibility modulus but has little effect on the bending modulus. Due to the weak interaction with the membrane, fibrils insignificantly change the membrane elastic properties.

Computational Models for the Study of Protein Aggregation.

Protein aggregation has been studied by many groups around the world for many years because it can be the cause of a number of neurodegenerative diseases that have no effective treatment. Obtaining the structure of related fibrils and toxic oligomers, as well as describing the pathways and main factors that govern the self-organization process, is of paramount importance, but it is also very difficult. To solve this problem, experimental and computational methods are often combined to get the most out of each method. The effectiveness of the computational approach largely depends on the construction of a reasonable molecular model. Here we discussed different versions of the four most popular all-atom force fields AMBER, CHARMM, GROMOS, and OPLS, which have been developed for folded and intrinsically disordered proteins, or both. Continuous and discrete coarse-grained models, which were mainly used to study the kinetics of aggregation, are also summarized.

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.

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 CD28-CD86 protein complex structure using different level of resolution approach.

Immune system plays essential role in functioning of higher organisms. Its hyperactivity can lead to autoimmune diseases or even anaphylactic shock while hypoactivity leads to proneness to infections or even cancer. T-cells play crucial role in immunity mechanisms and their activation and inhibition is strictly controlled by the regulatory proteins, such as CD28 and CTLA-4. Activity of these proteins is controlled by a pair of ligands, named CD80 and CD86, which can non-covalently bound to their receptors. While structure of human CTLA-4-CD86 complex in known, there is still no available structure for the CD28-CD86 system. To obtain the reliable structure of CD28-CD86 complex we first validated our methodology on the CTLA-4-CD86 system. Then coarse-grained UNRES-dock molecular docking simulation was performed followed by all-atom molecular dynamics simulations. As a result, we obtained a complete CD28-CD86 complex structure on atomistic level, in which interaction interface is consistent with available data. We also determined the kinetic properties for CTLA4-CD86 and CD28-CD86 complexes with use of coarse-grained model and determined the key residues for complex formation with use of Robetta, PPCheck and HawkDock servers. Our results not only verify high accuracy of the UNRES-dock method, but also provide a highly reliable model of the CD28-CD86 complex, which can be used in further studies and drug design.

Plant isoquinoline alkaloids as potential neurodrugs: A comparative study of the effects of benzo[c]phenanthridine and berberine-based compounds on ß-amyloid aggregation.

Herein we present a comparative study of the effects of isoquinoline alkaloids belonging to benzo[c]phenanthridine and berberine families on ß-amyloid aggregation. Results obtained using a Thioflavine T (ThT) fluorescence assay and circular dichroism (CD) spectroscopy suggested that the benzo[c]phenanthridine nucleus, present in both sanguinarine and chelerythrine molecules, was directly involved in an inhibitory effect of Aß1-42 aggregation. Conversely, coralyne, that contains the isomeric berberine nucleus, significantly increased propensity for Aß1-42 to aggregate. Surface Plasmon Resonance (SPR) experiments provided quantitative estimation of these interactions: coralyne bound to Aß1-42 with an affinity (KD = 11.6 µM) higher than benzo[c]phenanthridines. Molecular docking studies confirmed that all three compounds are able to recognize Aß1-42 in different aggregation forms suggesting their effective capacity to modulate the Aß1-42 self-recognition mechanism. Molecular dynamics simulations indicated that coralyne increased the ß-content of Aß1-42, in early stages of aggregation, consistent with fluorescence-based promotion of the Aß1-42 self-recognition mechanism by this alkaloid. At the same time, sanguinarine induced Aß1-42 helical conformation corroborating its ability to delay aggregation as experimentally proved in vitro. The investigated compounds were shown to interfere with aggregation of Aß1-42 demonstrating their potential as starting leads for the development of therapeutic strategies in neurodegenerative diseases.

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.

Computational Model to Unravel the Function of Amyloid-ß Peptides in Contact with a Phospholipid Membrane.

Divalent cations have a strong impact on the properties of phospholipid membranes, where amyloid-ß peptides exert effects related to possible functional or pathological roles. In this work, we use an atomistic computational model of dimyristoyl-phosphatidylcholine (DMPC) membrane bilayers. We perturb this model with a simple model of divalent cations (Mg2+) and with a single amyloid-ß (Aß) peptide of 42 residues, both with and without a single Cu2+ ion bound to the N-terminus. In agreement with the experimental results reported in the literature, the model confirms that divalent cations locally destabilize the DMPC membrane bilayer and, for the first time, that the monomeric form of Aß helps in avoiding the interactions between divalent cations and DMPC, preventing significant effects on the DMPC bilayer properties. These results are discussed in the frame of a protective role of the diluted Aß peptide floating around phospholipid membranes.

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.