Less Is More: Coarse-Grained Integrative Modeling of Large Biomolecular Assemblies with HADDOCK

Phosphorylation at the D56 residue of MtrA in Mycobacterium tuberculosis enhances its DNA binding affinity by modulating inter-domain interaction

The response regulator, MtrA, plays a major role in adaptation to the host environment, cell division, replication, and dormancy activation of Mycobacterium tuberculosis (Mtb). The phosphorylation of the response regulator MtrA alters the downstream activity, typically involving changes in DNA binding activity. However, there is a substantial knowledge gap in understanding the phosphorylation-mediated structural changes in MtrA. Additionally, the active conformation of the protein has yet to be determined. Therefore, in this study, we have investigated the phosphorylation-induced conformational changes of MtrA using all-atom molecular dynamics simulations under various phosphorylation conditions. The results from this study demonstrate that the phosphorylation at D56 (pD56-MtrA) increases the compactness of the MtrA protein by stabilizing the inter-domain interaction between the regulatory domain and DNA binding domain. Notably, the higher occupancy H-bond (over 95 %) between Arg200-Asn100 in case of the pD56-MtrA condition, which is otherwise absent in the non-phosphorylated (uMtrA) condition, suggests the importance of this interaction in the active conformation of the protein. The dynamic cross-correlation analysis reveals that phosphorylation (especially pD56-MtrA) reduces the anti-correlated motions and increases correlated motions between different domains. Moreover, the higher DNA binding affinity of pD56-MtrA compared to uMtrA supported by molecular docking and MD simulation followed by MMPBSA analysis suggests that pD56-MtrA is the possible active conformation of the MtrA protein. Overall, this investigation elucidates the key structural changes in MtrA under different phosphorylated conditions, which might help in designing novel therapeutics against tuberculosis.

Integration of molecular coarse-grained model into geometric representation learning framework for protein-protein complex property prediction

Structure-based machine learning algorithms have been utilized to predict the properties of protein-protein interaction (PPI) complexes, such as binding affinity, which is critical for understanding biological mechanisms and disease treatments. While most existing algorithms represent PPI complex graph structures at the atom-scale or residue-scale, these representations can be computationally expensive or may not sufficiently integrate finer chemical-plausible interaction details for improving predictions. Here, we introduce MCGLPPI, a geometric representation learning framework that combines graph neural networks (GNNs) with MARTINI molecular coarse-grained (CG) models to predict PPI overall properties accurately and efficiently. Extensive experiments on three types of downstream PPI property prediction tasks demonstrate that at the CG-scale, MCGLPPI achieves competitive performance compared with the counterparts at the atom- and residue-scale, but with only a third of computational resource consumption. Furthermore, CG-scale pre-training on protein domain-domain interaction structures enhances its predictive capabilities for PPI tasks. MCGLPPI offers an effective and efficient solution for PPI overall property predictions, serving as a promising tool for the large-scale analysis of biomolecular interactions.

The HADDOCK2.4 web server for integrative modeling of biomolecular complexes

Interactions between macromolecules, such as proteins and nucleic acids, are essential for cellular functions. Experimental methods can fail to provide all the information required to fully model biomolecular complexes at atomic resolution, particularly for large and heterogeneous assemblies. Integrative computational approaches have, therefore, gained popularity, complementing traditional experimental methods in structural biology. Here, we introduce HADDOCK2.4, an integrative modeling platform, and its updated web interface ( https://wenmr.science.uu.nl/haddock2.4 ). The platform seamlessly integrates diverse experimental and theoretical data to generate high-quality models of macromolecular complexes. The user-friendly web server offers automated parameter settings, access to distributed computing resources, and pre- and post-processing steps that enhance the user experience. To present the web server's various interfaces and features, we demonstrate two different applications: (i) we predict the structure of an antibody-antigen complex by using NMR data for the antigen and knowledge of the hypervariable loops for the antibody, and (ii) we perform coarse-grained modeling of PRC1 with a nucleosome particle guided by mutagenesis and functional data. The described protocols require some basic familiarity with molecular modeling and the Linux command shell. This new version of our widely used HADDOCK web server allows structural biologists and non-experts to explore intricate macromolecular assemblies encompassing various molecule types.

Breaking the 'don't eat me' signal: in silico design of CD47-directed peptides for cancer immunotherapy

The leading cause of death worldwide is cancer. Although there are various therapies available to treat cancer, finding a successful one can be like searching for a needle in a haystack. Immunotherapy appears to be one of those needles in the haystack of cancer treatment. Immunotherapeutic agents enhance the immune response of the patient's body to tumor cells. One of the immunotherapeutic targets, Cluster of Differentiation 47 (CD47), releases the "don't eat me" signal when it binds to its receptor, Signal Regulatory Protein (SIRPα). Tumor cells use this signal to circumvent the immune system, rendering it ineffective. To stop tumor cells from releasing the "don't eat me" signal, the CD47-SIRPα interaction is specifically targeted in this study. To do so, in silico peptides were designed based on the structural analysis of the interaction between two proteins using point mutations on the interacting residues with the other amino acids. The peptide library was designed and docked on SIRPα using computational tools. Later on, after analyzing the docked complex, the best of them was selected for MD simulation studies of 100 ns. Further analysis after MD studies was carried out to determine the possible potential anti-SIRPα peptides.

Molecular simulation reveals that pathogenic mutations in BTB/ANK domains of Arabidopsis thaliana NPR1 circumscribe the EDS1-mediated immune regulation

The NPR1 (nonexpressor of pathogenesis-related genes 1) is a key regulator of the salicylic-acid-mediated immune response caused by pathogens in Arabidopsis thaliana. Mutations C150Y and H334Y in the BTB/ANK domains of NPR1 inhibit the defense response, and transcriptional co-activity with enhanced disease susceptibility 1 (EDS1) has been revealed experimentally. This study examined the conformational changes and reduced NPR1-EDS1 interaction upon mutation using a molecular dynamics simulation. Initially, BTBC150YNPR1 and ANKH334YNPR1 were categorized as pathological mutations rather than others based on sequence conservation. A distant ortholog was used to map the common residues shared among the wild-type because the mutations were highly conserved. Overall, 179 of 373 residues were determining the secondary structures and fold versatility of conformations. In addition, the mutational hotspots Cys150, Asp152, Glu153, Cys155, His157, Cys160, His334, Arg339 and Lys370 were crucial for oligomer-to-monomer exchange. Subsequently, the atomistic simulations with free energy (MM/PB(GB)SA) calculations predicted structural displacements engaging in the N-termini α5133-178α7 linker connecting the central ANK regions (α13260-290α14 and α18320-390α22), where prominent long helices (α516) and short helices (α310) replaced with β-turns and loops disrupting hydrogen bonds and salt bridges in both mutants implicating functional regulation and activation. Furthermore, the mutation repositions the intact stability of multiple regions (L13C149-N356α20BTB/ANK-α17W301-E357α21N-ter/coiled-coil) compromising a dynamic interaction of NPR1-EDS1. By unveiling the transitions between the distinct functions of mutational perception, this study paves the way for future investigation to orchestrate additive host-adapted transcriptional reprogramming that controls defense-related regulatory mechanisms of NPR1s in plants.

Riding the wave of innovation: immunoinformatics in fish disease control

The spread of infectious illnesses has been a significant factor restricting aquaculture production. To maximise aquatic animal health, vaccination tactics are very successful and cost-efficient for protecting fish and aquaculture animals against many disease pathogens. However, due to the increasing number of immunological cases and their complexity, it is impossible to manage, analyse, visualise, and interpret such data without the assistance of advanced computational techniques. Hence, the use of immunoinformatics tools is crucial, as they not only facilitate the management of massive amounts of data but also greatly contribute to the creation of fresh hypotheses regarding immune responses. In recent years, advances in biotechnology and immunoinformatics have opened up new research avenues for generating novel vaccines and enhancing existing vaccinations against outbreaks of infectious illnesses, thereby reducing aquaculture losses. This review focuses on understanding in silico epitope-based vaccine design, the creation of multi-epitope vaccines, the molecular interaction of immunogenic vaccines, and the application of immunoinformatics in fish disease based on the frequency of their application and reliable results. It is believed that it can bridge the gap between experimental and computational approaches and reduce the need for experimental research, so that only wet laboratory testing integrated with in silico techniques may yield highly promising results and be useful for the development of vaccines for fish.

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.

Exploring the prognostic significance of PKCε variants in cervical cancer

BACKGROUND

Protein Kinase C-epsilon (PKCε) is a member of the novel subfamily of PKCs (nPKCs) that plays a role in cancer development. Studies have revealed that its elevated expression levels are associated with cervical cancer. Previously, we identified pathogenic variations in its different domains through various bioinformatics tools and molecular dynamic simulation. In the present study, the aim was to find the association of its variants rs1553369874 and rs1345511001 with cervical cancer and to determine the influence of these variants on the protein-protein interactions of PKCε, which can lead towards cancer development and poor survival rates.

METHODS

The association of the variants with cervical cancer and its clinicopathological features was determined through genotyping analysis. Odds ratio and relative risk along with Fisher exact test were calculated to evaluate variants significance and disease risk. Protein-protein docking was performed and docked complexes were subjected to molecular dynamics simulation to gauge the variants impact on PKCε's molecular interactions.

RESULTS

This study revealed that genetic variants rs1553369874 and rs1345511001 were associated with cervical cancer. Smad3 interacts with PKCε and this interaction promotes cervical cancer angiogenesis; therefore, Smad3 was selected for protein-protein docking. The analysis revealed PKCε variants promoted aberrant interactions with Smad3 that might lead to the activation of oncogenic pathways. The data obtained from this study suggested the prognostic significance of PRKCE gene variants rs1553369874 and rs1345511001.

CONCLUSION

Through further in vitro and in vivo validation, these variants can be used at the clinical level as novel prognostic markers and therapeutic targets against cervical cancer.

Hybrid computational methods combining experimental information with molecular dynamics

A goal of structural biology is to understand how macromolecules carry out their biological roles by identifying their metastable states, mechanisms of action, pathways leading to conformational changes, and the thermodynamic and kinetic relationships between those states. Integrative modeling brings structural insights into systems where traditional structure determination approaches cannot help. We focus on the synergies and challenges of integrative modeling combining experimental data with molecular dynamics simulations.

Host Cell Membrane Capture by the SARS-CoV-2 Spike Protein Fusion Intermediate

Cell entry by SARS-CoV-2 is accomplished by the S2 subunit of the spike S protein on the virion surface by capture of the host cell membrane and fusion with the viral envelope. Capture and fusion require the prefusion S2 to transit to its potent fusogenic form, the fusion intermediate (FI). However, the FI structure is unknown, detailed computational models of the FI are unavailable, and the mechanisms and timing of membrane capture and fusion are not established. Here, we constructed a full-length model of the SARS-CoV-2 FI by extrapolating from known SARS-CoV-2 pre- and postfusion structures. In atomistic and coarse-grained molecular dynamics simulations the FI was remarkably flexible and executed giant bending and extensional fluctuations due to three hinges in the C-terminal base. The simulated configurations and their giant fluctuations are quantitatively consistent with SARS-CoV-2 FI configurations measured recently using cryo-electron tomography. Simulations suggested a host cell membrane capture time of ∼2 ms. Isolated fusion peptide simulations identified an N-terminal helix that directed and maintained binding to the membrane but grossly underestimated the binding time, showing that the fusion peptide environment is radically altered when attached to its host fusion protein. The large configurational fluctuations of the FI generated a substantial exploration volume that aided capture of the target membrane, and may set the waiting time for fluctuation-triggered refolding of the FI that draws the viral envelope and host cell membrane together for fusion. These results describe the FI as machinery that uses massive configurational fluctuations for efficient membrane capture and suggest novel potential drug targets.

Machine learned coarse-grained protein force-fields: Are we there yet?

The successful recent application of machine learning methods to scientific problems includes the learning of flexible and accurate atomic-level force-fields for materials and biomolecules from quantum chemical data. In parallel, the machine learning of force-fields at coarser resolutions is rapidly gaining relevance as an efficient way to represent the higher-body interactions needed in coarse-grained force-fields to compensate for the omitted degrees of freedom. Coarse-grained models are important for the study of systems at time and length scales exceeding those of atomistic simulations. However, the development of transferable coarse-grained models via machine learning still presents significant challenges. Here, we discuss recent developments in this field and current efforts to address the remaining challenges.

Molecular determinants associated with temporal succession of SARS-CoV-2 variants in Uttar Pradesh, India

The emergence and rapid evolution of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) caused a global crisis that required a detailed characterization of the dynamics of mutational pattern of the viral genome for comprehending its epidemiology, pathogenesis and containment. We investigated the molecular evolution of the SASR-CoV-2 genome during the first, second and third waves of COVID-19 in Uttar Pradesh, India. Nanopore sequencing of the SARS-CoV-2 genome was undertaken in 544 confirmed cases of COVID-19, which included vaccinated and unvaccinated individuals. In the first wave (unvaccinated population), the 20A clade (56.32%) was superior that was replaced by 21A Delta in the second wave, which was more often seen in vaccinated individuals in comparison to unvaccinated (75.84% versus 16.17%, respectively). Subsequently, 21A delta got outcompeted by Omicron (71.8%), especially the 21L variant, in the third wave. We noticed that Q677H appeared in 20A Alpha and stayed up to Delta, D614G appeared in 20A Alpha and stayed in Delta and Omicron variants (got fixed), and several other mutations appeared in Delta and stayed in Omicron. A cross-sectional analysis of the vaccinated and unvaccinated individuals during the second wave revealed signature combinations of E156G, F157Del, L452R, T478K, D614G mutations in the Spike protein that might have facilitated vaccination breach in India. Interestingly, some of these mutation combinations were carried forward from Delta to Omicron. In silico protein docking showed that Omicron had a higher binding affinity with the host ACE2 receptor, resulting in enhanced infectivity of Omicron over the Delta variant. This work has identified the combinations of key mutations causing vaccination breach in India and provided insights into the change of [virus's] binding affinity with evolution, resulting in more virulence in Delta and more infectivity in Omicron variants of SARS-CoV-2. Our findings will help in understanding the COVID-19 disease biology and guide further surveillance of the SARS-CoV-2 genome to facilitate the development of vaccines with better efficacies.

New opportunities in integrative structural modeling

Integrative structural modeling enables structure determination of macromolecules and their complexes by integrating data from multiple sources. It has been successfully used to characterize macromolecular structures when a single structural biology technique was insufficient. Recent developments in cellular structural biology, including in-cell cryo-electron tomography and artificial intelligence-based structure prediction, have created new opportunities for integrative structural modeling. Here, we will review these opportunities along with the latest developments in integrative modeling methods and their applications. We also highlight open challenges and directions for further development.

Identification of Potential ACE2-Derived Peptide Mimetics in SARS-CoV-2 Omicron Variant Therapeutics using Computational Approaches

The COVID-19 pandemic has become a global health challenge because of the emergence of distinct variants. Omicron, a new variant, is recognized as a variant of concern (VOC) by the World Health Organization (WHO) because of its higher mutations and accelerated human infection. The infection rate is strongly dependent on the binding rate of the receptor binding domain (RBD) against human angiotensin converting enzyme-2 (ACE2human) receptor. Inhibition of protein-protein (RBDs(SARS-CoV-2/omicron)-ACE2human) interaction has been already proven to inhibit viral infection. We have systematically designed ACE2human-derived peptides and peptide mimetics that have high binding affinity toward RBDomicron. Our peptide mutational analysis indicated the influence of canonical amino acids on the peptide binding process. Herein, efforts have been made to explore the atomistic details and events of RBDs(SARS-CoV-2/omicron)-ACE2human interactions by using molecular dynamics simulation. Our studies pave a path for developing therapeutic peptidomimetics against omicron.

An in-frame deletion mutation in the degron tail of auxin coreceptor IAA2 confers resistance to the herbicide 2,4-D in Sisymbrium orientale

Synthetic auxin herbicides intersect basic plant developmental biology and applied weed management. We investigated resistance to 2,4-D in the Australian weed Sisymbrium orientale (Indian hedge mustard). We identified a mechanism involving an in-frame 27-bp deletion in the degron tail of auxin coreceptor IAA2, one member of the gene family of Aux/IAA auxin co-receptors. We show that this deletion in IAA2 is a gain-of-function mutation that confers synthetic auxin resistance. This field-evolved mechanism of resistance to synthetic auxin herbicides confirms previous biochemical studies showing the role of the Aux/IAA degron tail in regulating Aux/IAA protein degradation upon auxin perception. The deletion mutation could be generated in crops using gene-editing approaches for cross-resistance to multiple synthetic auxin herbicides.

The natural auxin indole-3-acetic acid (IAA) is a key regulator of many aspects of plant growth and development. Synthetic auxin herbicides such as 2,4-D mimic the effects of IAA by inducing strong auxinic-signaling responses in plants. To determine the mechanism of 2,4-D resistance in a Sisymbrium orientale (Indian hedge mustard) weed population, we performed a transcriptome analysis of 2,4-D-resistant (R) and -susceptible (S) genotypes that revealed an in-frame 27-nucleotide deletion removing nine amino acids in the degron tail (DT) of the auxin coreceptor Aux/IAA2 (SoIAA2). The deletion allele cosegregated with 2,4-D resistance in recombinant inbred lines. Further, this deletion was also detected in several 2,4-D-resistant field populations of this species. Arabidopsis transgenic lines expressing the SoIAA2 mutant allele were resistant to 2,4-D and dicamba. The IAA2-DT deletion reduced binding to TIR1 in vitro with both natural and synthetic auxins, causing reduced association and increased dissociation rates. This mechanism of synthetic auxin herbicide resistance assigns an in planta function to the DT region of this Aux/IAA coreceptor for its role in synthetic auxin binding kinetics and reveals a potential biotechnological approach to produce synthetic auxin-resistant crops using gene-editing.

Interface refinement of low- to medium-resolution Cryo-EM complexes using HADDOCK2.4

A wide range of cellular processes requires the formation of multimeric protein complexes. The rise of cryo-electron microscopy (cryo-EM) has enabled the structural characterization of these protein assemblies. The density maps produced can, however, still suffer from limited resolution, impeding the process of resolving structures at atomic resolution. In order to solve this issue, monomers can be fitted into low- to medium-resolution maps. Unfortunately, the models produced frequently contain atomic clashes at the protein-protein interfaces (PPIs), as intermolecular interactions are typically not considered during monomer fitting. Here, we present a refinement approach based on HADDOCK2.4 to remove intermolecular clashes and optimize PPIs. A dataset of 14 cryo-EM complexes was used to test eight protocols. The best-performing protocol, consisting of a semi-flexible simulated annealing refinement with centroid restraints on the monomers, was able to decrease intermolecular atomic clashes by 98% without significantly deteriorating the quality of the cryo-EM density fit.

Unified Nanotechnology Format: One Way to Store Them All

The domains of DNA and RNA nanotechnology are steadily gaining in popularity while proving their value with various successful results, including biosensing robots and drug delivery cages. Nowadays, the nanotechnology design pipeline usually relies on computer-based design (CAD) approaches to design and simulate the desired structure before the wet lab assembly. To aid with these tasks, various software tools exist and are often used in conjunction. However, their interoperability is hindered by a lack of a common file format that is fully descriptive of the many design paradigms. Therefore, in this paper, we propose a Unified Nanotechnology Format (UNF) designed specifically for the biomimetic nanotechnology field. UNF allows storage of both design and simulation data in a single file, including free-form and lattice-based DNA structures. By defining a logical and versatile format, we hope it will become a widely accepted and used file format for the nucleic acid nanotechnology community, facilitating the future work of researchers and software developers. Together with the format description and publicly available documentation, we provide a set of converters from existing file formats to simplify the transition. Finally, we present several use cases visualizing example structures stored in UNF, showcasing the various types of data UNF can handle.

Leveraging crosslinking mass spectrometry in structural and cell biology

Crosslinking mass spectrometry (crosslinking-MS) is a versatile tool providing structural insights into protein conformation and protein-protein interactions. Its medium-resolution residue-residue distance restraints have been used to validate protein structures proposed by other methods and have helped derive models of protein complexes by integrative structural biology approaches. The use of crosslinking-MS in integrative approaches is underpinned by progress in estimating error rates in crosslinking-MS data and in combining these data with other information. The flexible and high-throughput nature of crosslinking-MS has allowed it to complement the ongoing resolution revolution in electron microscopy by providing system-wide residue-residue distance restraints, especially for flexible regions or systems. Here, we review how crosslinking-MS information has been leveraged in structural model validation and integrative modeling. Crosslinking-MS has also been a key technology for cell biology studies and structural systems biology where, in conjunction with cryoelectron tomography, it can provide structural and mechanistic insights directly in situ.

Protein-Protein Docking: Past, Present, and Future

The biological significance of proteins attracted the scientific community in exploring their characteristics. The studies shed light on the interaction patterns and functions of proteins in a living body. Due to their practical difficulties, reliable experimental techniques pave the way for introducing computational methods in the interaction prediction. Automated methods reduced the difficulties but could not yet replace experimental studies as the field is still evolving. Interaction prediction problem being critical needs highly accurate results, but none of the existing methods could offer reliable performance that can parallel with experimental results yet. This article aims to assess the existing computational docking algorithms, their challenges, and future scope. Blind docking techniques are quite helpful when no information other than the individual structures are available. As more and more complex structures are being added to different databases, information-driven approaches can be a good alternative. Artificial intelligence, ruling over the major fields, is expected to take over this domain very shortly.

Accurate Estimation of Solvent Accessible Surface Area for Coarse-Grained Biomolecular Structures with Deep Learning

Coarse-grained (CG) models of biomolecules have been widely used in protein/ribonucleic acid (RNA) three-dimensional structure prediction, docking, drug design, and molecular simulations due to their superiority in computational efficiency. Most of these applications strongly depend on the reasonable estimation of solvation free energy, which requires the accurate calculation of solvent accessible surface area (SASA). Although algorithms for SASA calculations with all-atom protein and RNA structures have been well-established, accurately estimating the SASA based on CG structures is extremely challenging. In this work, we developed a deep learning-based SASA estimator (DeepCGSA), which can provide almost perfect SASA estimation based on CG structures of protein and RNA molecules. Extensive testing analysis showed that for three types of widely used CG protein models, including the Cα-based, Cα-Cβ, and Martini models, the correlation coefficients between the predicted values and the reference values can be as high as 0.95-0.99, which perform dramatically better than available methods. In addition, the new method can be used for CG RNA structures and unfolded protein structures with much improved accuracy. We anticipate that DeepCGSA will be highly useful in the protein/RNA structure prediction, drug design, and other applications, in which accurate estimations of SASA for CG biomolecular structures are critically important.

Protein-Protein Docking with Large-Scale Backbone Flexibility Using Coarse-Grained Monte-Carlo Simulations

Most of the protein-protein docking methods treat proteins as almost rigid objects. Only the side-chains flexibility is usually taken into account. The few approaches enabling docking with a flexible backbone typically work in two steps, in which the search for protein-protein orientations and structure flexibility are simulated separately. In this work, we propose a new straightforward approach for docking sampling. It consists of a single simulation step during which a protein undergoes large-scale backbone rearrangements, rotations, and translations. Simultaneously, the other protein exhibits small backbone fluctuations. Such extensive sampling was possible using the CABS coarse-grained protein model and Replica Exchange Monte Carlo dynamics at a reasonable computational cost. In our proof-of-concept simulations of 62 protein-protein complexes, we obtained acceptable quality models for a significant number of cases.

The milk-derived lactoferrin inhibits V-ATPase activity by targeting its V1 domain

Lactoferrin (Lf), a bioactive milk protein, exhibits strong anticancer and antifungal activities. The search for Lf targets and mechanisms of action is of utmost importance to enhance its effective applications. A common feature among Lf-treated cancer and fungal cells is the inhibition of a proton pump called V-ATPase. Lf-driven V-ATPase inhibition leads to cytosolic acidification, ultimately causing cell death of cancer and fungal cells. Given that a detailed elucidation of how Lf and V-ATPase interact is still missing, herein we aimed to fill this gap by employing a five-stage computational approach. Molecular dynamics simulations of both proteins were performed to obtain a robust sampling of their conformational landscape, followed by clustering, which allowed retrieving representative structures, to then perform protein-protein docking. Subsequently, molecular dynamics simulations of the docked complexes and free binding energy calculations were carried out to evaluate the dynamic binding process and build a final ranking based on the binding affinities. Detailed atomist analysis of the top ranked complexes clearly indicates that Lf binds to the V₁ cytosolic domain of V-ATPase. Particularly, our data suggest that Lf binds to the interfaces between A/B subunits, where the ATP hydrolysis occurs, thus inhibiting this process. The free energy decomposition analysis further identified key binding residues that will certainly aid in the rational design of follow-up experimental studies, hence bridging computational and experimental biochemistry.

Perspectives on High-Throughput Ligand/Protein Docking With Martini MD Simulations

Molecular docking is central to rational drug design. Current docking techniques suffer, however, from limitations in protein flexibility and solvation models and by the use of simplified scoring functions. All-atom molecular dynamics simulations, on the other hand, feature a realistic representation of protein flexibility and solvent, but require knowledge of the binding site. Recently we showed that coarse-grained molecular dynamics simulations, based on the most recent version of the Martini force field, can be used to predict protein/ligand binding sites and pathways, without requiring any a priori information, and offer a level of accuracy approaching all-atom simulations. Given the excellent computational efficiency of Martini, this opens the way to high-throughput drug screening based on dynamic docking pipelines. In this opinion article, we sketch the roadmap to achieve this goal.

Integrative modeling of membrane-associated protein assemblies

Membrane proteins are among the most challenging systems to study with experimental structural biology techniques. The increased number of deposited structures of membrane proteins has opened the route to modeling their complexes by methods such as docking. Here, we present an integrative computational protocol for the modeling of membrane-associated protein assemblies. The information encoded by the membrane is represented by artificial beads, which allow targeting of the docking toward the binding-competent regions. It combines efficient, artificial intelligence-based rigid-body docking by LightDock with a flexible final refinement with HADDOCK to remove potential clashes at the interface. We demonstrate the performance of this protocol on eighteen membrane-associated complexes, whose interface lies between the membrane and either the cytosolic or periplasmic regions. In addition, we provide a comparison to another state-of-the-art docking software, ZDOCK. This protocol should shed light on the still dark fraction of the interactome consisting of membrane proteins.

Computational structure modeling for diverse categories of macromolecular interactions

Computational protein-protein docking is one of the most intensively studied topics in structural bioinformatics. The field has made substantial progress through over three decades of development. The development began with methods for rigid-body docking of two proteins, which have now been extended in different directions to cover the various macromolecular interactions observed in a cell. Here, we overview the recent developments of the variations of docking methods, including multiple protein docking, peptide-protein docking, and disordered protein docking methods.

Protein-ligand binding with the coarse-grained Martini model

The detailed understanding of the binding of small molecules to proteins is the key for the development of novel drugs or to increase the acceptance of substrates by enzymes. Nowadays, computer-aided design of protein–ligand binding is an important tool to accomplish this task. Current approaches typically rely on high-throughput docking essays or computationally expensive atomistic molecular dynamics simulations. Here, we present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein–ligand interactions of small drug-like molecules. Remarkably, we achieve high accuracy without the need of any a priori knowledge of binding pockets or pathways. Our approach is applied to a range of systems from the well-characterized T4 lysozyme over members of the GPCR family and nuclear receptors to a variety of enzymes. The presented results open the way to high-throughput screening of ligand libraries or protein mutations using the coarse-grained Martini model.

Computer-aided design of protein-ligand binding is important for the development of novel drugs. Here authors present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein-ligand binding interactions of small drug-like molecules.

Towards gaining sight of multiscale events: utilizing network models and normal modes in hybrid methods

With the explosion of normal mode analyses (NMAs) based on elastic network models (ENMs) in the last decade, and the proven precision of MD simulations for visualizing interactions at atomic scale, many hybrid methods have been proposed in recent years. These aim at exploiting the best of both worlds: the atomic precision of MD that often fall short of exploring time and length scales of biological interest, and the capability of ENM-NMA to predict the cooperative and often functional rearrangements of large structures and assemblies, albeit at low resolution. We present an overview of recent progress in the field with examples of successful applications highlighting the utility of such hybrid methods and pointing to emerging future directions guided by advances in experimental characterization of biomolecular systems structure and dynamics.

Coarse-grained (hybrid) integrative modeling of biomolecular interactions

The computational modeling field has vastly evolved over the past decades. The early developments of simplified protein systems represented a stepping stone towards establishing more efficient approaches to sample intricated conformational landscapes. Downscaling the level of resolution of biomolecules to coarser representations allows for studying protein structure, dynamics and interactions that are not accessible by classical atomistic approaches. The combination of different resolutions, namely hybrid modeling, has also been proved as an alternative when mixed levels of details are required. In this review, we provide an overview of coarse-grained/hybrid models focusing on their applicability in the modeling of biomolecular interactions. We give a detailed list of ready-to-use modeling software for studying biomolecular interactions allowing various levels of coarse-graining and provide examples of complexes determined by integrative coarse-grained/hybrid approaches in combination with experimental information.

An overview of data-driven HADDOCK strategies in CAPRI rounds 38-45

Our information-driven docking approach HADDOCK has demonstrated a sustained performance since the start of its participation to CAPRI. This is due, in part, to its ability to integrate data into the modeling process, and to the robustness of its scoring function. We participated in CAPRI both as server and manual predictors. In CAPRI rounds 38-45, we have used various strategies depending on the available information. These ranged from imposing restraints to a few residues identified from literature as being important for the interaction, to binding pockets identified from homologous complexes or template-based refinement/CA-CA restraint-guided docking from identified templates. When relevant, symmetry restraints were used to limit the conformational sampling. We also tested for a large decamer target a new implementation of the MARTINI coarse-grained force field in HADDOCK. Overall, we obtained acceptable or better predictions for 13 and 11 server and manual submissions, respectively, out of the 22 interfaces. Our server performance (acceptable or higher-quality models when considering the top 10) was better (59%) than the manual (50%) one, in which we typically experiment with various combinations of protocols and data sources. Again, our simple scoring function based on a linear combination of intermolecular van der Waals and electrostatic energies and an empirical desolvation term demonstrated a good performance in the scoring experiment with a 63% success rate across all 22 interfaces. An analysis of model quality indicates that, while we are consistently performing well in generating acceptable models, there is room for improvement for generating/identifying higher quality models.

MARTINI-Based Protein-DNA Coarse-Grained HADDOCKing

Modeling biomolecular assemblies is an important field in computational structural biology. The inherent complexity of their energy landscape and the computational cost associated with modeling large and complex assemblies are major drawbacks for integrative modeling approaches. The so-called coarse-graining approaches, which reduce the degrees of freedom of the system by grouping several atoms into larger “pseudo-atoms,” have been shown to alleviate some of those limitations, facilitating the identification of the global energy minima assumed to correspond to the native state of the complex, while making the calculations more efficient. Here, we describe and assess the implementation of the MARTINI force field for DNA into HADDOCK, our integrative modeling platform. We combine it with our previous implementation for protein-protein coarse-grained docking, enabling coarse-grained modeling of protein-nucleic acid complexes. The system is modeled using MARTINI topologies and interaction parameters during the rigid body docking and semi-flexible refinement stages of HADDOCK, and the resulting models are then converted back to atomistic resolution by an atom-to-bead distance restraints-guided protocol. We first demonstrate the performance of this protocol using 44 complexes from the protein-DNA docking benchmark, which shows an overall ~6-fold speed increase and maintains similar accuracy as compared to standard atomistic calculations. As a proof of concept, we then model the interaction between the PRC1 and the nucleosome (a former CAPRI target in round 31), using the same information available at the time the target was offered, and compare all-atom and coarse-grained models.