Unveiling structural and energetic characterization of the emissive RNA alphabet anchored in the methylthieno[3,4-d]pyrimidine heterocycle core.

This study presents a comprehensive theoretical exploration of the fluorescent non-natural emissive nucleobases- mthA, mthG, mthC, and mthU derived from the methylthieno[3,4-d]pyrimidine heterocycle. Our calculations, aligning with experimental findings, reveal that these non-natural bases exert minimal influence on the geometry of classical Watson-Crick base pairs within an RNA duplex, maintaining H-bonding akin to natural bases. In terms of energy, the impact of the modified bases, but for mthG, is also found to be little significant. We delved into an in-depth analysis of the photophysical properties of these non-natural bases. This investigation unveiled a correlation between their absorption/emission peaks and the substantial impact of the modification on the energy levels of the highest unoccupied molecular orbitals (HOMO) and the lowest unoccupied molecular orbital (LUMO). Notably, this alteration in energy levels resulted in a significant reduction of the HOMO-LUMO gap, from approximately 5.4-5.5 eV in the natural bases, to roughly 3.9-4.7 eV in the modified bases. This shift led to a consequential change in absorption and emission spectra towards longer wavelengths, elucidating their bathochromic shift.

Occurrence and stability of anion-π interactions between phosphate and nucleobases in functional RNA molecules.

We present a systematic structural and energetic characterization of phosphate(OP)-nucleobase anion…π stacking interactions in RNAs. We observed OP-nucleobase stacking contacts in a variety of structural motifs other than regular helices and spanning broadly diverse sequence distances. Apart from the stacking between a phosphate and a guanine or a uracil two-residue upstream in specific U-turns, such interactions in RNA have been scarcely characterized to date. Our QM calculations showed an energy minimum at a distance between the OP atom and the nucleobase plane centroid slightly below 3 Å for all the nucleobases. By sliding the OP atom over the nucleobase plane we localized the optimal mutual positioning of the stacked moieties, corresponding to an energy minimum below -6 kcal•mol-1, for all the nucleobases, consistently with the projections of the OP atoms over the different π-rings we observed in experimental occurrences. We also found that the strength of the interaction clearly correlates with its electrostatic component, pointing to it as the most relevant contribution. Finally, as OP-uracil and OP-guanine interactions represent together 86% of the instances we detected, we also proved their stability under dynamic conditions in model systems simulated by state-of-the art DFT-MD calculations.

Discriminating physiological from non-physiological interfaces in structures of protein complexes: A community-wide study.

Reliably scoring and ranking candidate models of protein complexes and assigning their oligomeric state from the structure of the crystal lattice represent outstanding challenges. A community-wide effort was launched to tackle these challenges. The latest resources on protein complexes and interfaces were exploited to derive a benchmark dataset consisting of 1677 homodimer protein crystal structures, including a balanced mix of physiological and non-physiological complexes. The non-physiological complexes in the benchmark were selected to bury a similar or larger interface area than their physiological counterparts, making it more difficult for scoring functions to differentiate between them. Next, 252 functions for scoring protein-protein interfaces previously developed by 13 groups were collected and evaluated for their ability to discriminate between physiological and non-physiological complexes. A simple consensus score generated using the best performing score of each of the 13 groups, and a cross-validated Random Forest (RF) classifier were created. Both approaches showed excellent performance, with an area under the Receiver Operating Characteristic (ROC) curve of 0.93 and 0.94, respectively, outperforming individual scores developed by different groups. Additionally, AlphaFold2 engines recalled the physiological dimers with significantly higher accuracy than the non-physiological set, lending support to the reliability of our benchmark dataset annotations. Optimizing the combined power of interface scoring functions and evaluating it on challenging benchmark datasets appears to be a promising strategy.

Toward α-1,3/4 fucosyltransferases targeted drug discovery: In silico uncovering of promising natural inhibitors of fucosyltransferase 6.

Sialyl Lewis X (sLex ) antigen is a fucosylated cell-surface glycan that is normally involved in cell-cell interactions. The enhanced expression of sLex on cell surface glycans, which is attributed to the upregulation of fucosyltransferase 6 (FUT6), has been implicated in facilitating metastasis in human colorectal, lung, prostate, and oral cancers. The role that the upregulated FUT6 plays in the progression of tumor to malignancy, with reduced survival rates, makes it a potential target for anticancer drugs. Unfortunately, the lack of experimental structures for FUT6 has hampered the design and development of its inhibitors. In this study, we used in silico techniques to identify potential FUT6 inhibitors. We first modeled the three-dimensional structure of human FUT6 using AlphaFold. Then, we screened the natural compound libraries from the COCONUT database to sort out potential natural products (NPs) with best affinity toward the FUT6 model. As a result of these simulations, we identified three NPs for which we predicted binding affinities and interaction patterns quite similar to those we calculated for two experimentally tested FUT6 inhibitors, that is, fucose mimetic-1 and a GDP-triazole derived compound. We also performed molecular dynamics (MD) simulations for the FUT6 complexes with identified NPs, to investigate their stability. Analysis of the MD simulations showed that the identified NPs establish stable contacts with FUT6 under dynamics conditions. On these grounds, the three screened compounds appear as promising natural alternatives to experimentally tested FUT6 synthetic inhibitors, with expected comparable binding affinity. This envisages good prospects for future experimental validation toward FUT6 inhibition.

Occurrence and stability of lone pair-π and OH-π interactions between water and nucleobases in functional RNAs.

We identified over 1000 instances of water-nucleobase stacking contacts in a variety of RNA molecules from a non-redundant set of crystal structures with resolution ≤3.0 Å. Such contacts may be of either the lone pair-π (lp-π) or the OH-π type, in nature. The distribution of the distances of the water oxygen from the nucleobase plane peaks at 3.5 Å for A, G and C, and approximately at 3.1-3.2 Å for U. Quantum mechanics (QM) calculations confirm, as expected, that the optimal energy is reached at a shorter distance for the lp-π interaction as compared to the OH-π one (3.0 versus 3.5 Å). The preference of each nucleobase for either type of interaction closely correlates with its electrostatic potential map. Furthermore, QM calculations show that for all the nucleobases a favorable interaction, of either the lp-π or the OH-π type, can be established at virtually any position of the water molecule above the nucleobase skeleton, which is consistent with the uniform projection of the OW atoms over the nucleobases ring we observed in the experimental occurrences. Finally, molecular dynamics simulations of a model system for the characterization of water-nucleobase stacking contacts confirm the stability of these interactions also under dynamic conditions.

D936Y and Other Mutations in the Fusion Core of the SARS-CoV-2 Spike Protein Heptad Repeat 1: Frequency, Geographical Distribution, and Structural Effect.

The crown of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is constituted by its spike (S) glycoprotein. S protein mediates the SARS-CoV-2 entry into the host cells. The "fusion core" of the heptad repeat 1 (HR1) on S plays a crucial role in the virus infectivity, as it is part of a key membrane fusion architecture. While SARS-CoV-2 was becoming a global threat, scientists have been accumulating data on the virus at an impressive pace, both in terms of genomic sequences and of three-dimensional structures. On 15 February 2021, from the SARS-CoV-2 genomic sequences in the GISAID resource, we collected 415,673 complete S protein sequences and identified all the mutations occurring in the HR1 fusion core. This is a 21-residue segment, which, in the post-fusion conformation of the protein, gives many strong interactions with the heptad repeat 2, bringing viral and cellular membranes in proximity for fusion. We investigated the frequency and structural effect of novel mutations accumulated over time in such a crucial region for the virus infectivity. Three mutations were quite frequent, occurring in over 0.1% of the total sequences. These were S929T, D936Y, and S949F, all in the N-terminal half of the HR1 fusion core segment and particularly spread in Europe and USA. The most frequent of them, D936Y, was present in 17% of sequences from Finland and 12% of sequences from Sweden. In the post-fusion conformation of the unmutated S protein, D936 is involved in an inter-monomer salt bridge with R1185. We investigated the effect of the D936Y mutation on the pre-fusion and post-fusion state of the protein by using molecular dynamics, showing how it especially affects the latter one.

The CASP13-CAPRI targets as case studies to illustrate a novel scoring pipeline integrating CONSRANK with clustering and interface analyses.

BACKGROUND: Properly scoring protein-protein docking models to single out the correct ones is an open challenge, also object of assessment in CAPRI (Critical Assessment of PRedicted Interactions), a community-wide blind docking experiment. We introduced in the field CONSRANK (CONSensus RANKing), the first pure consensus method. Also available as a web server, CONSRANK ranks docking models in an ensemble based on their ability to match the most frequent inter-residue contacts in it. We have been blindly testing CONSRANK in all the latest CAPRI rounds, where we showed it to perform competitively with the state-of-the-art energy and knowledge-based scoring functions. More recently, we developed Clust-CONSRANK, an algorithm introducing a contact-based clustering of the models as a preliminary step of the CONSRANK scoring process. In the latest CASP13-CAPRI joint experiment, we participated as scorers with a novel pipeline, combining both our scoring tools, CONSRANK and Clust-CONSRANK, with our interface analysis tool COCOMAPS. Selection of the 10 models for submission was guided by the strength of the emerging consensus, and their final ranking was assisted by results of the interface analysis. RESULTS: As a result of the above approach, we were by far the first scorer in the CASP13-CAPRI top-1 ranking, having high/medium quality models ranked at the top-1 position for the majority of targets (11 out of the total 19). We were also the first scorer in the top-10 ranking, on a par with another group, and the second scorer in the top-5 ranking. Further, we topped the ranking relative to the prediction of binding interfaces, among all the scorers and predictors. Using the CASP13-CAPRI targets as case studies, we illustrate here in detail the approach we adopted. CONCLUSIONS: Introducing some flexibility in the final model selection and ranking, as well as differentiating the adopted scoring approach depending on the targets were the key assets for our highly successful performance, as compared to previous CAPRI rounds. The approach we propose is entirely based on methods made available to the community and could thus be reproduced by any user.

A homozygous MRPL24 mutation causes a complex movement disorder and affects the mitoribosome assembly.

Mitochondrial ribosomal protein large 24 (MRPL24) is 1 of the 82 protein components of mitochondrial ribosomes, playing an essential role in the mitochondrial translation process. We report here on a baby girl with cerebellar atrophy, choreoathetosis of limbs and face, intellectual disability and a combined defect of complexes I and IV in muscle biopsy, caused by a homozygous missense mutation identified in MRPL24. The variant predicts a Leu91Pro substitution at an evolutionarily conserved site. Using human mutant cells and the zebrafish model, we demonstrated the pathological role of the identified variant. In fact, in fibroblasts we observed a significant reduction of MRPL24 protein and of mitochondrial respiratory chain complex I and IV subunits, as well a markedly reduced synthesis of the mtDNA-encoded peptides. In zebrafish we demonstrated that the orthologue gene is expressed in metabolically active tissues, and that gene knockdown induced locomotion impairment, structural defects and low ATP production. The motor phenotype was complemented by human WT but not mutant cRNA. Moreover, sucrose density gradient fractionation showed perturbed assembly of large subunit mitoribosomal proteins, suggesting that the mutation leads to a conformational change in MRPL24, which is expected to cause an aberrant interaction of the protein with other components of the 39S mitoribosomal subunit.

Mapping the minimum domain of the fibronectin binding site on transglutaminase 2 (TG2) and its importance in mediating signaling, adhesion, and migration in TG2-expressing cells.

The interaction between the enzyme transglutaminase 2 (TG2) and fibronectin (FN) is involved in the cell-matrix interactions that regulate cell signaling, adhesion, and migration and play central roles in pathologic conditions, particularly fibrosis and cancer. A precise definition of the exact interaction domains on both proteins could provide a tool to design novel molecules with potential therapeutic applications. Although specific residues involved in the interaction within TG2 have been analyzed, little is known regarding the TG2 binding site on FN. This site has been mapped to a large internal 45-kDa protein fragment coincident with the gelatin binding domain (GBD). With the goal of defining the minimal FN interacting domain for TG2, we produced several expression constructs encoding different portions or modules of the GBD and tested their binding and functional properties. The results demonstrate that the I8 module is necessary and sufficient for TG2-binding in vitro, but does not have functional effects on TG2-expressing cells. Modules I7 and I9 increase the strength of the binding and are required for cell adhesion. A 15-kDa fragment encompassing modules I7-9 behaves as the whole 45-kDa GBD and mediates signaling, adhesion, spreading, and migration of TG2+ cells. This study provides new insights into the mechanism for TG2 binding to FN.-Soluri, M. F., Boccafoschi, F., Cotella, D., Moro, L., Forestieri, G., Autiero, I., Cavallo, L., Oliva, R., Griffin, M., Wang, Z., Santoro, C., Sblattero, D. Mapping the minimum domain of the fibronectin binding site on transglutaminase 2 (TG2) and its importance in mediating signaling, adhesion, and migration in TG2-expressing cells.

A novel method for improved accuracy of transcription factor binding site prediction.

Identifying transcription factor (TF) binding sites (TFBSs) is important in the computational inference of gene regulation. Widely used computational methods of TFBS prediction based on position weight matrices (PWMs) usually have high false positive rates. Moreover, computational studies of transcription regulation in eukaryotes frequently require numerous PWM models of TFBSs due to a large number of TFs involved. To overcome these problems we developed DRAF, a novel method for TFBS prediction that requires only 14 prediction models for 232 human TFs, while at the same time significantly improves prediction accuracy. DRAF models use more features than PWM models, as they combine information from TFBS sequences and physicochemical properties of TF DNA-binding domains into machine learning models. Evaluation of DRAF on 98 human ChIP-seq datasets shows on average 1.54-, 1.96- and 5.19-fold reduction of false positives at the same sensitivities compared to models from HOCOMOCO, TRANSFAC and DeepBind, respectively. This observation suggests that one can efficiently replace the PWM models for TFBS prediction by a small number of DRAF models that significantly improve prediction accuracy. The DRAF method is implemented in a web tool and in a stand-alone software freely available at http://cbrc.kaust.edu.sa/DRAF.

Structural and Energetic Impact of Non-natural 7-Deaza-8-azaguanine, 7-Deaza-8-azaisoguanine, and Their 7-Substituted Derivatives on Hydrogen-Bond Pairing with Cytosine and Isocytosine.

The impact of 7-deaza-8-azaguanine (DAG) and 7-deaza-8-azaisoguanine (DAiG) modifications on the geometry and stability of the G:C Watson-Crick (cWW) base pair and the G:iC and iG:C reverse Watson-Crick (tWW) base pairs has been characterized theoretically. In addition, the effect on the same base pairs of seven C7-substituted DAG and DAiG derivatives, some of which have been previously experimentally characterized, has been investigated. Calculations indicate that all of these modifications have a negligible impact on the geometry of the above base pairs, and that modification of the heterocycle skeleton has a small impact on the base-pair interaction energies. Instead, base-pair interaction energies are dependent on the nature of the C7 substituent. For the 7-substituted DAG-C cWW systems, a linear correlation between the base-pair interaction energy and the Hammett constant of the 7-substituent is found, with higher interaction energies corresponding to more electron-withdrawing substituents. Therefore, the explored modifications are expected to be accommodated in both parallel and antiparallel nucleic acid duplexes without perturbing their geometry, while the strength of a base pair (and duplex) featuring a DAG modification can, in principle, be tuned by incorporating different substituents at the C7 position.

Targeting CD34+ cells of the inflamed synovial endothelium by guided nanoparticles for the treatment of rheumatoid arthritis.

Despite the advances in the treatment of rheumatoid arthritis (RA) achieved in the last few years, several patients are diagnosed late, do not respond to or have to stop therapy because of inefficacy and/or toxicity, leaving still a huge unmet need. Tissue-specific strategies have the potential to address some of these issues. The aim of the study is the development of a safe nanotechnology approach for tissue-specific delivery of drugs and diagnostic probes. CD34 + endothelial precursors were addressed in inflamed synovium using targeted biodegradable nanoparticles (tBNPs). These nanostructures were made of poly-lactic acid, poly-caprolactone, and PEG and then coated with a synovial homing peptide. Immunofluorescence analysis clearly demonstrated their capacity to selectively address CD34 + endothelial cells in synovial tissue obtained from human, mouse, and rat. Biodistribution studies in two different animal models of rheumatoid arthritis (antigen-induced arthritis/AIA and collagen-induced arthritis/CIA) confirmed the selective accumulation in inflamed joints but also evidenced the capacity of tBNP to detect early phases of the disease and the preferential liver elimination. The therapeutic effect of methotrexate (MTX)-loaded tBNPs were studied in comparison with conventional MTX doses. MTX-loaded tBNPs prevented and treated CIA and AIA at a lower dose and reduced administration frequency than MTX. Moreover, MTX-loaded tBNP showed a novel mechanism of action, in which the particles target and kill CD34 + endothelial progenitors, preventing neo-angiogenesis and, consequently, synovial inflammation. tBNPs represent a stable and safe platform to develop highly-sensitive imaging and therapeutic approaches in RA targeting specifically synovial neo-angiogenesis to reduce local inflammation.

Occurrence and stability of lone pair-π stacking interactions between ribose and nucleobases in functional RNAs.

The specific folding pattern and function of RNA molecules lies in various weak interactions, in addition to the strong base-base pairing and stacking. One of these relatively weak interactions, characterized by the stacking of the O4' atom of a ribose on top of the heterocycle ring of a nucleobase, has been known to occur but has largely been ignored in the description of RNA structures. We identified 2015 ribose-base stacking interactions in a high-resolution set of non-redundant RNA crystal structures. They are widespread in structured RNA molecules and are located in structural motifs other than regular stems. Over 50% of them involve an adenine, as we found ribose-adenine contacts to be recurring elements in A-minor motifs. Fewer than 50% of the interactions involve a ribose and a base of neighboring residues, while approximately 30% of them involve a ribose and a nucleobase at least four residues apart. Some of them establish inter-domain or inter-molecular contacts and often implicate functionally relevant nucleotides. In vacuo ribose-nucleobase stacking interaction energies were calculated by quantum mechanics methods. Finally, we found that lone pair-π stacking interactions also occur between ribose and aromatic amino acids in RNA-protein complexes.

Energetics and dynamics of the non-natural fluorescent 4AP:DAP base pair.

The fluorescent non-natural 4-aminophthalimide (4AP) base, when paired to the complementary 2,4-diaminopyrimidine (DAP) nucleobase, is accommodated in a B-DNA duplex being efficiently recognized and incorporated by DNA polymerases. To complement the experimental studies and rationalize the impact of the above non-natural bases on the structure, stability and dynamics of nucleic acid structures, we performed quantum mechanics (QM) calculations along with classical molecular dynamics (MD) simulations. QM calculations were initially focused on the geometry and energetics of the 4AP:DAP non-natural pair and of H-bonded base pairs between 4AP and all the natural bases in their classical Watson-Crick geometries. The QM calculations indicate that the 4AP:DAP pair, despite the fact that it can form 3 H-bonds in a classic Watson-Crick geometry, has a stability comparable to the A:T pair. Then, we extended the study to reverse Watson-Crick geometries, characteristic of parallel strands. MD simulations were carried out on two 13-mer DNA duplexes, featuring a central 4AP:DAP or A:T pair, respectively. No major structural deformation of the duplex was observed during the MD simulation. Snapshots from the MD simulations were subjected to QM calculations to investigate the 4AP:DAP interaction energy when embedded into a duplex structure, and to investigate the impact of the two non-natural bases on the stacking interactions with adjacent bases in the DNA duplex. We found a slight increase in stacking interactions involving the 4AP:DAP pair, counterbalanced by a moderate decrease in H-bonding interactions of the 4AP:DAP and of the adjacent base pairs in the duplex. The results of our study are in agreement with experimental data and complement them by providing an insight into which factors contribute positively and which factors contribute negatively to the structural compatibility of the fluorescent 4AP:DAP pair with a B-DNA structure.

Theoretical characterization of sulfur-to-selenium substitution in an emissive RNA alphabet: impact on H-bonding potential and photophysical properties.

We employ density functional theory (DFT) and time-dependent DFT (TDDFT) calculations to investigate the structural, energetic and optical properties of a new computationally designed RNA alphabet, where the nucleobases, tsA, tsG, tsC, and tsU (ts-bases), have been derived by replacing sulfur with selenium in the previously reported tz-bases, based on the isothiazolo[4,3-d]pyrimidine heterocycle core. We find out that the modeled non-natural bases have minimal impact on the geometry and energetics of the classical Watson-Crick base pairs, thus potentially mimicking the natural bases in a RNA duplex in terms of H-bonding. In contrast, our calculations indicate that H-bonded base pairs involving the Hoogsteen edge of purines are destabilized as compared to their natural counterparts. We also focus on the photophysical properties of the non-natural bases and correlate their absorption/emission peaks to the strong impact of the modification on the energy of the lowest unoccupied molecular orbital. It is indeed stabilized by roughly 1.1-1.6 eV as compared to the natural analogues, resulting in a reduction of the gap between the highest occupied and the lowest unoccupied molecular orbital from 5.3-5.5 eV in the natural bases to 3.9-4.2 eV in the modified ones, with a consequent bathochromic shift in the absorption and emission spectra. Overall, our analysis clearly indicates that the newly modelled ts-bases are expected to exhibit better fluorescent properties as compared to the previously reported tz-bases, while retaining similar H-bonding properties. In addition, we show that a new RNA alphabet based on size-extended benzo-homologated ts-bases can also form stable Watson-Crick base pairs with the natural complementary nucleobases.

Structural Basis for Mutations of Human Aquaporins Associated to Genetic Diseases.

Aquaporins (AQPs) are among the best structural-characterized membrane proteins, fulfilling the role of allowing water flux across cellular membranes. Thus far, 34 single amino acid polymorphisms have been reported in HUMSAVAR for human aquaporins as disease-related. They affect AQP2, AQP5 and AQP8, where they are associated with nephrogenic diabetes insipidus, keratoderma and colorectal cancer, respectively. For half of these mutations, although they are mostly experimentally characterized in their dysfunctional phenotypes, a structural characterization at a molecular level is still missing. In this work, we focus on such mutations and discuss what the structural defects are that they appear to cause. To achieve this aim, we built a 3D molecular model for each mutant and explored the effect of the mutation on all of their structural features. Based on these analyses, we could collect the structural defects of all the pathogenic mutations (here or previously analysed) under few main categories, that we found to nicely correlate with the experimental phenotypes reported for several of the analysed mutants. Some of the structural analyses we present here provide a rationale for previously experimentally observed phenotypes. Furthermore, our comprehensive overview can be used as a reference frame for the interpretation, on a structural basis, of defective phenotypes of other aquaporin pathogenic mutants.

Novel mutation in mitochondrial Elongation Factor EF-Tu associated to dysplastic leukoencephalopathy and defective mitochondrial DNA translation.

The mitochondrial Elongation Factor Tu (EF-Tu), encoded by the TUFM gene, is a highly conserved GTPase, which is part of the mitochondrial protein translation machinery. In its activated form it delivers the aminoacyl-tRNAs to the A site of the mitochondrial ribosome. We report here on a baby girl with severe infantile macrocystic leukodystrophy with micropolygyria and a combined defect of complexes I and IV in muscle biopsy, caused by a novel mutation identified in TUFM. Using human mutant cells and the yeast model, we demonstrate the pathological role of the novel variant. Moreover, results of a molecular modeling study suggest that the mutant is inactive in mitochondrial polypeptide chain elongation, probably as a consequence of its reduced ability to bind mitochondrial aa-tRNAs. Four patients have so far been described with mutations in TUFM, and, following the first description of the disease in a single patient, we describe similar clinical and neuroradiological features in an additional patient.

An atlas of RNA base pairs involving modified nucleobases with optimal geometries and accurate energies.

Posttranscriptional modifications greatly enhance the chemical information of RNA molecules, contributing to explain the diversity of their structures and functions. A significant fraction of RNA experimental structures available to date present modified nucleobases, with half of them being involved in H-bonding interactions with other bases, i.e. 'modified base pairs'. Herein we present a systematic investigation of modified base pairs, in the context of experimental RNA structures. To this end, we first compiled an atlas of experimentally observed modified base pairs, for which we recorded occurrences and structural context. Then, for each base pair, we selected a representative for subsequent quantum mechanics calculations, to find out its optimal geometry and interaction energy. Our structural analyses show that most of the modified base pairs are non Watson-Crick like and are involved in RNA tertiary structure motifs. In addition, quantum mechanics calculations quantify and provide a rationale for the impact of the different modifications on the geometry and stability of the base pairs they participate in.

Analysis of the interface variability in NMR structure ensembles of protein-protein complexes.

NMR structures consist in ensembles of conformers, all satisfying the experimental restraints, which exhibit a certain degree of structural variability. We analyzed here the interface in NMR ensembles of protein-protein heterodimeric complexes and found it to span a wide range of different conservations. The different exhibited conservations do not simply correlate with the size of the systems/interfaces, and are most probably the result of an interplay between different factors, including the quality of experimental data and the intrinsic complex flexibility. In any case, this information is not to be missed when NMR structures of protein-protein complexes are analyzed; especially considering that, as we also show here, the first NMR conformer is usually not the one which best reflects the overall interface. To quantify the interface conservation and to analyze it, we used an approach originally conceived for the analysis and ranking of ensembles of docking models, which has now been extended to directly deal with NMR ensembles. We propose this approach, based on the conservation of the inter-residue contacts at the interface, both for the analysis of the interface in whole ensembles of NMR complexes and for the possible selection of a single conformer as the best representative of the overall interface. In order to make the analyses automatic and fast, we made the protocol available as a web tool at: https://www.molnac.unisa.it/BioTools/consrank/consrank-nmr.html.

CONSRANK: a server for the analysis, comparison and ranking of docking models based on inter-residue contacts.

SUMMARY: Herein, we present CONSRANK, a web tool for analyzing, comparing and ranking protein-protein and protein-nucleic acid docking models, based on the conservation of inter-residue contacts and its visualization in 2D and 3D interactive contact maps. AVAILABILITY AND IMPLEMENTATION: CONSRANK is accessible as a public web tool at https://www.molnac.unisa.it/BioTools/consrank/. CONTACT: romina.oliva@uniparthenope.it.