Structure Of Complex Research Articles

Structure prediction of protein-ligand complexes from sequence information with Umol

Protein-ligand docking is an established tool in drug discovery and development to narrow down potential therapeutics for experimental testing. However, a high-quality protein structure is required and often the protein is treated as fully or partially rigid. Here we develop an AI system that can predict the fully flexible all-atom structure of protein-ligand complexes directly from sequence information. We find that classical docking methods are still superior, but depend upon having crystal structures of the target protein. In addition to predicting flexible all-atom structures, predicted confidence metrics (plDDT) can be used to select accurate predictions as well as to distinguish between strong and weak binders. The advances presented here suggest that the goal of AI-based drug discovery is one step closer, but there is still a way to go to grasp the complexity of protein-ligand interactions fully. Umol is available at: https://github.com/patrickbryant1/Umol.

MoLPC2: Improved prediction of large protein complex structures and stoichiometry using monte carlo tree search and AlphaFold2.

Today, the prediction of structures of large protein complexes solely from their sequence information requires prior knowledge of the stoichiometry of the complex. To address this challenge, we have enhanced the Monte Carlo Tree Search algorithms in MoLPC to enable the assembly of protein complexes while simultaneously predicting their stoichiometry. In MoLPC2, we have improved the predictions by allowing sampling alternative AlphaFold predictions. Using MoLPC2, we accurately predicted the structures of 50 out of 175 non-redundant protein complexes (TM-score > = 0.8) without knowing the stoichiometry. MoLPC2 provides new opportunities for predicting protein complex structures without stoichiometry information. MoLPC2 is freely available at https://github.com/hychim/molpc2. A notebook is also available from the repository for easy use. Supplementary data are available at Bioinformatics online.

Structure and mechanism of the human CTDNEP1–NEP1R1 membrane protein phosphatase complex necessary to maintain ER membrane morphology

C-terminal Domain Nuclear Envelope Phosphatase 1 (CTDNEP1) is a noncanonical protein serine/threonine phosphatase that has a conserved role in regulating ER membrane biogenesis. Inactivating mutations in CTDNEP1 correlate with the development of medulloblastoma, an aggressive childhood cancer. The transmembrane protein Nuclear Envelope Phosphatase 1 Regulatory Subunit 1 (NEP1R1) binds CTDNEP1, but the molecular details by which NEP1R1 regulates CTDNEP1 function are unclear. Here, we find that knockdown of NEP1R1 generates identical phenotypes to reported loss of CTDNEP1 in mammalian cells, establishing CTDNEP1-NEP1R1 as an evolutionarily conserved membrane protein phosphatase complex that restricts ER expansion. Mechanistically, NEP1R1 acts as an activating regulatory subunit that directly binds and increases the phosphatase activity of CTDNEP1. By defining a minimal NEP1R1 domain sufficient to activate CTDNEP1, we determine high-resolution crystal structures of the CTDNEP1-NEP1R1 complex bound to a peptide sequence acting as a pseudosubstrate. Structurally, NEP1R1 engages CTDNEP1 at a site distant from the active site to stabilize and allosterically activate CTDNEP1. Substrate recognition is facilitated by a conserved Arg residue in CTDNEP1 that binds and orients the substrate peptide in the active site. Together, this reveals mechanisms for how NEP1R1 regulates CTDNEP1 and explains how cancer-associated mutations inactivate CTDNEP1.

Isotope-Edited Variable Temperature Infrared Spectroscopy for Measuring Transition Temperatures of Single A-T Watson-Crick Base Pairs in DNA Duplexes.

Experimental methods to determine transition temperatures for individual base pair melting events in DNA duplexes are lacking despite intense interest in these thermodynamic parameters. Here, we determine the dimensions of the thymine (T) C2═O stretching vibration when it is within the DNA duplex via isotopic substitutions at other atomic positions in the structure. First, we determined that this stretching state was localized enough to specific atoms in the molecule to make submolecular scale measurements of local structure and stability in high molecular weight complexes. Next, we develop a new isotope-edited variable temperature infrared method to measure melting transitions at various locations in a DNA structure. As an initial test of this "sub-molecular scale thermometer", we applied our T13C2 difference infrared signal to measure location-dependent melting temperatures (TmL) in a DNA duplex via variable temperature attenuated total reflectance Fourier transform infrared (VT-ATR-FTIR) spectroscopy. We report that the TmL of a single Watson-Crick A-T base pair near the end of an A-T rich sequence (poly T) is ∼34.9 ± 0.7°C. This is slightly lower than the TmL of a single base pair near the middle position of the poly T sequence (TmL ∼35.6±0.2°C). In addition, we also report that the TmL of a single Watson-Crick A-T base pair near the end of a 50% G-C sequence (12-mer) is ∼52.5 ± 0.3°C, which is slightly lower than the global melting Tm of the 12-mer sequence (TmL ∼54.0±0.9°C). Our results provide direct physical evidence for end fraying in DNA sequences with our novel spectroscopic methods.

Structural insights into the cross-exon to cross-intron spliceosome switch.

Early spliceosome assembly can occur through an intron-defined pathway, whereby U1 and U2 small nuclear ribonucleoprotein particles (snRNPs) assemble across the intron1. Alternatively, it can occur through an exon-defined pathway2-5, whereby U2 binds the branch site located upstream of the defined exon and U1 snRNP interacts with the 5' splice site located directly downstream of it. The U4/U6.U5 tri-snRNP subsequently binds to produce a cross-intron (CI) or cross-exon (CE) pre-B complex, which is then converted to the spliceosomal B complex6,7. Exon definition promotes the splicing of upstream introns2,8,9 and plays a key part in alternative splicing regulation10-16. However, the three-dimensional structure of exon-defined spliceosomal complexes and the molecular mechanism of the conversion from a CE-organized to a CI-organized spliceosome, a pre-requisite for splicing catalysis, remain poorly understood. Here cryo-electron microscopy analyses of human CE pre-B complex and B-like complexes reveal extensive structural similarities with their CI counterparts. The results indicate that the CE and CI spliceosome assembly pathways converge already at the pre-B stage. Add-back experiments using purified CE pre-B complexes, coupled with cryo-electron microscopy, elucidate the order of the extensive remodelling events that accompany the formation of B complexes and B-like complexes. The molecular triggers and roles of B-specific proteins in these rearrangements are also identified. We show that CE pre-B complexes can productively bind in trans to a U1 snRNP-bound 5' splice site. Together, our studies provide new mechanistic insights into the CE to CI switch during spliceosome assembly and its effect on pre-mRNA splice site pairing at this stage.

(5-Fluoro-2,6-dioxo-1,2,3,6-tetrahydropyrimidin-1-ido-κN 1)(1,4,8,11-tetraazacyclotetradecane-κ4 N)zinc(II) perchlorate

In the structure of the title complex, [Zn(C4H2FN2O2)(C10H24N4)]ClO4, the zinc(II) ion forms coordination bonds with the four nitrogen atoms of cyclam (1,4,8,11-tetraazacyclotetradecane or [14]aneN4) as well as with the nitrogen atom of a deprotonated 5-fluorouracil ion (FU−). Cyclam adopts a trans-I type conformation within this structure. The coordination structure of the zinc(II) ion is a square pyramid with a distorted base plane formed by the four nitrogen atoms of the cyclam. FU− engages in intermolecular hydrogen bonding with neighboring FU− molecules and with the cyclam molecule.

DOX-DNA Interactions on the Nanoscale: In Situ Studies Using Tip-Enhanced Raman Scattering.

Chemotherapeutic anthracyclines, like doxorubicin (DOX), are drugs endowed with cytostatic activity and are widely used in antitumor therapy. Their molecular mechanism of action involves the formation of a stable anthracycline-DNA complex, which prevents cell division and results in cell death. It is known that elevated DOX concentrations induce DNA chain loops and overlaps. Here, for the first time, tip-enhanced Raman scattering was used to identify and localize intercalated DOX in isolated double-stranded calf thymus DNA, and the correlated near-field spectroscopic and morphologic experiments locate the DOX molecules in the DNA and provide further information regarding specific DOX-nucleobase interactions. Thus, the study provides a tool specifically for identifying intercalation markers and generally analyzing drug-DNA interactions. The structure of such complexes down to the molecular level provides mechanistic information about cytotoxicity and the development of potential anticancer drugs.

Experimental biological studies of novel tetradentate hydrazone Cu(II) complexes for potential applications in medicinal chemistry

A series of mononuclear Cu(II) complexes were synthesized using a potential tetradentate hydrazone ligand obtained from the reaction between phenylhydrazine groups and 2,5-thiophenedicarboxaldehyde. The structures of ligands and complexes were elucidated through various spectroscopic techniques, confirming their composition. All complexes were found to adopt four-coordinated geometries, indicating the formation of stable structures. Spectroscopic analysis revealed that the hydrazone ligand coordinated with the Cu(II) metal ions as a dibasic tetradentate ligand by utilizing the phenolic oxygen and azomethine nitrogen atoms. The binding affinity of the complexes with calf thymus DNA (CT-DNA) was investigated using absorption and viscosity measurements, demonstrating their interaction through the intercalation mode. The binding studies showed that the Cu(II) complexes exhibited varying degrees of binding affinity, with Cu(L4) demonstrating the highest affinity, followed by Cu(L1), Cu(L2), and Cu(L3). Moreover, the DNA fragmentation properties of the Cu(II) complexes were evaluated, suggesting their potential utilization in pharmaceutical applications. The obtained results highlight the significance of these novel complexes in the field of medicinal chemistry.

Insights into coordination and ligand trends of lanthanide complexes from the Cambridge Structural Database

Understanding lanthanide coordination chemistry can help develop new ligands for more efficient separation of lanthanides for critical materials needs. The Cambridge Structural Database (CSD) contains tens of thousands of single crystal structures of lanthanide complexes that can serve as a training ground for both fundamental chemical insights and future machine learning and generative artificial intelligence models. This work aims to understand the currently available structures of lanthanide complexes in CSD by analyzing the coordination shell, donor types, and ligand types, from the perspective of rare-earth element (REE) separations. We obtain four sets of lanthanide complexes from CSD: Subset 1, all Ln-containing complexes (49472 structures); Subset 2, mononuclear Ln complexes (27858 structures); Subset 3, mononuclear Ln complexes without cyclopentadienyl ligands (Cp) (26156 structures); Subset 4, Ln complexes with at least one 1,10-phenanthroline (phen) or its derivative as a coordinating ligand (2226 structures). The subsequent analysis of lanthanide complexes in these subsets examines the trends in coordination numbers and first shell distances as well as identifies and characterizes the ligands and donor groups. In addition, examples of Ln-complexes with commercially available complexants and phen-based ligands are interrogated in detail. This systematic investigation lays the groundwork for future data-driven ligand designs for REE separations based on the structural insights into the lanthanide coordination chemistry.

Structures of the mumps virus polymerase complex via cryo-electron microscopy

The viral polymerase complex, comprising the large protein (L) and phosphoprotein (P), is crucial for both genome replication and transcription in non-segmented negative-strand RNA viruses (nsNSVs), while structures corresponding to these activities remain obscure. Here, we resolved two L–P complex conformations from the mumps virus (MuV), a typical member of nsNSVs, via cryogenic-electron microscopy. One conformation presents all five domains of L forming a continuous RNA tunnel to the methyltransferase domain (MTase), preferably as a transcription state. The other conformation has the appendage averaged out, which is inaccessible to MTase. In both conformations, parallel P tetramers are revealed around MuV L, which, together with structures of other nsNSVs, demonstrates the diverse origins of the L-binding X domain of P. Our study links varying structures of nsNSV polymerase complexes with genome replication and transcription and points to a sliding model for polymerase complexes to advance along the RNA templates.

Prediction of Paratope-Epitope Pairs Using Convolutional Neural Networks.

Antibodies play a central role in the adaptive immune response of vertebrates through the specific recognition of exogenous or endogenous antigens. The rational design of antibodies has a wide range of biotechnological and medical applications, such as in disease diagnosis and treatment. However, there are currently no reliable methods for predicting the antibodies that recognize a specific antigen region (or epitope) and, conversely, epitopes that recognize the binding region of a given antibody (or paratope). To fill this gap, we developed ImaPEp, a machine learning-based tool for predicting the binding probability of paratope-epitope pairs, where the epitope and paratope patches were simplified into interacting two-dimensional patches, which were colored according to the values of selected features, and pixelated. The specific recognition of an epitope image by a paratope image was achieved by using a convolutional neural network-based model, which was trained on a set of two-dimensional paratope-epitope images derived from experimental structures of antibody-antigen complexes. Our method achieves good performances in terms of cross-validation with a balanced accuracy of 0.8. Finally, we showcase examples of application of ImaPep, including extensive screening of large libraries to identify paratope candidates that bind to a selected epitope, and rescoring and refining antibody-antigen docking poses.

Predicting Emission of Heteroleptic Iridium Complexes using Artificial Chemical Intelligence.

We report a deep learning-based approach to accurately predict the emission spectraof phosphorescent heteroleptic [Ir(C^N)2(NN)]+complexes, enabling the rapid discovery of novel Ir(III) chromophores for diverse applications including organic light-emitting diodes and solar fuel cells. The deep learning models utilize graph neural networks and other chemical features in architectures that reflect the inherent structure of the heteroleptic complexes, composed of C^Nand N^Nligands, and are thus geared towards efficient training over the dataset. By leveraging experimental emission data, our models reliably predict the full emission spectra of these complexes across various emission profiles, surpassing the accuracy of conventional DFT and correlated wavefunction methods, while simultaneously achieving robustness to the presence of imperfect (noisy, low-quality) training spectra. We showcase the potential applications for these and related models for \insilico\ prediction of complexes with tailored emission properties, as well as in "design of experiment'' contexts to reduce the synthetic burden of high-throughput screening. In the latter case, we demonstrate that the models allow to exploit a limited amount of experimental data to explore a wide range of chemical space, thus leveraging a modest synthetic effort.

SACRED ARCHITECTURE OF JAPAN DURING THE SPREAD OF BUDDHISM: LOCAL TRADITIONS AND CHINESE INFLUENCES

The article examines the sacred complexes of Buddhism and Shinto from the time of the beginning of the spread of Buddhism in Japan in the Nara period (645-783) in order to identify the nature of the mutual influence of two architectural traditions – the continental one, represented by the Buddhist architecture of China, and the local Japanese one, represented by Shinto architecture. As architectural examples, two authentic objects of both religious movements built in the same time period in the area of Nara were chosen, namely: the Toshodai-ji monastery complex, founded in 759 by the Buddhist monk Jianzhen, and the Kasuga shrine complex, built in the early eighth century by representatives of the powerful Fujiwara clan, the protector of imperial power. The analysis of the selected objects is carried out in four aspects: the nature of the organization of space, the nature of the interaction of structural elements, the nature of their accessibility and semantic content. The architecture of the Toshodai-ji monastery complex is dominated by continental influences, which are reflected in the organization of the complex, the structure and interior of its main pavilion kondo, and its semantics. The architecture of the complex shows the beginning of the process of transformation of local architectural traditions into Buddhist architecture, which is mainly reflected in the the structure and organization of the interior space of the kondo. The Kasuga shrine complex, on the contrary, although influenced by continental techniques, remains largely within the framework of local influences, which prevail mainly at the level of organization of the complex and at the level of its semantics. Continental influences are manifested in attempts to organize the structure of the shrine complex following the examples of Buddhist architectural complexes, as well as in the architectural form of the hondo pavilions. The identification of the characteristic features of the two architectural traditions was carried out in parallel for each of the selected objects, and based on the results obtained, the nature of the mutual influence of the two architectural concepts was compared in quantitative terms.

Structural Diversity and Rare η1 Cu-C Interactions in CuI Complexes of 1,2,3-Triazole-Functionalized Bisphosphines.

This manuscript describes the synthesis of copper complexes of 1,2,3-triazolyl-phosphines: o-Ph2P(C6H4){1,2,3-N3CC6H5)C(PPh2)} (L1), (C6H5){1,2,3-N3C(C6H4(o-PPh2))-C(PPh2)} (L2), 3-Ph2P(C5H3N){1,2,3-N3C(C6H5)C(PPh2)} (L3), o-Ph2P(C6H4){1,2,3-N3C(C5H5N)C(PPh2)} (L4), and {(3,5-Ph2PC6H4-o)21,2,3-N3CCH} (L5). The reactions of L1-L3 with CuI salts afforded dimeric complexes having the general formula [Cu2(μ -X)2L2] (L = L1, X = Cl, Br and I: 1 - 3; L= L2, X = Cl, Br and I: 4- 6; L = L3; X = Cl, Br, and I: 7-9). The reaction of L4 with CuI in a 1:2 molar ratio afforded 1D-coordination polymer [{(CuI)2{o-Ph2P(C6H4){1,2,3-N3C(C5H5N)C(PPh2)}-μ-((k1-N)(k2-P,P))}}]n (10). The reaction of L5 with cuprous halides (CuX) (X = Br, I) yielded mononuclear complexes [CuX{(3,5-Ph2PC6H4-o)21,2,3-N3CCH}-κ2P,P] (X = Br, 12; I, 13). Crystal structures of complexes 12 and 13 showed close interactions between CuI and the triazole C7 carbon. These relatively short Cu···C7 separations may be due to the η1-C interaction (dπ-pπ bond) between the triazolic carbon atom (via pz orbital) and CuI or three-centered two-electron interaction between CuI and the triazolic C-H bond. The existence of the Cu···C interaction was further evinced by the QTAIM analysis in compounds 12 and 13.

Unconventional PDZ Recognition Revealed in α7 nAChR-PICK1 Complexes.

PDZ domains are modular domains that conventionally bind to C terminal or internal motifs of target proteins to control cellular functions through the regulation of protein complex assemblies. Almost all reported structures of PDZ-target protein complexes rely on fragments or peptides as target proteins. No intact target protein complexed with PDZ was structurally characterized. In this study, we used NMR spectroscopy and other biochemistry and biophysics tools to uncover insights into structural coupling between the PDZ domain of protein interacting with C-kinase 1 (PICK1) and α7 nicotinic acetylcholine receptors (α7 nAChR). Notably, the intracellular domains of both α7 nAChR and PICK1 PDZ exhibit a high degree of plasticity in their coupling. Specifically, the MA helix of α7 nAChR interacts with residues lining the canonical binding site of the PICK1 PDZ, while flexible loops also engage in protein-protein interactions. Both hydrophobic and electrostatic interactions mediate the coupling. Overall, the resulting structure of the α7 nAChR-PICK1 complex reveals an unconventional PDZ binding mode, significantly expanding the repertoire of functionally important PDZ interactions.

Molecular Mechanistic Insights into Dipeptidyl Peptidase-IV Inhibitory Peptides to Decipher the Structural Basis of Activity.

Dipeptidyl peptidase-IV (DPP-IV) inhibiting peptides have attracted increased attention because of their possible beneficial effects on glycemic homeostasis. However, the structural basis underpinning their activities has not been well understood. This study combined computational and in vitro investigations to explore the structural basis of DPP-IV inhibitory peptides. We first superimposed the Xaa-Pro-type peptide-like structures from several crystal structures of DPP-IV ligand-protein complexes to analyze the recognition interactions of DPP-IV to peptides. Thereafter, a small set of Xaa-Pro-type peptides was designed to explore the effect of key interactions on inhibitory activity. The intramolecular interaction of Xaa-Pro-type peptides at the first and third positions from the N-terminus was pivotal to their inhibitory activities. Residue interactions between DPP-IV and residues of the peptides at the fourth and fifth positions of the N-terminus contributed significantly to the inhibitory effect of Xaa-Pro-type tetrapeptides and pentapeptides. Based on the interaction descriptors, quantitative structure-activity relationship (QSAR) studies with the DPP-IV inhibitory peptides resulted in valid models with high R2 values (0.90 for tripeptides; 0.91 for tetrapeptides and pentapeptides) and Q2 values (0.33 for tripeptides; 0.68 for tetrapeptides and pentapeptides). Taken together, the structural information on DPP-IV and peptides in this study facilitated the development of novel DPP-IV inhibitory peptides.

Differential Studies on the Structure of Lignin-Carbohydrate Complexes (LCC) in Alkali-Extracted Plant Hemicelluloses.

Hemicellulose extracted by alkali treatment is of interest because of the advantages of its intact sugar structure and high degree of polymerization. However, the hemicellulose extracted by alkali treatment contained more lignin fragments and the presence of a lignin-carbohydrate complex (LCC), which affected the isolation and purification of hemicellulose and its comprehensive utilization. Therefore, the evaluation of the LCC structure of different types of lignocellulosic resources is of great significance. In this study, the LCC structures of hardwoods and Gramineae were enriched in alkaline systems. Information on the composition, structural proportions, and connection patterns of LCC samples was discussed. The similarities and differences between the LCC structures of different units of raw materials were comparatively studied. The results indicated that the monosaccharide fractions were higher in the LCC of Gramineae compared to hardwoods. The composition of the lignin fraction was dominated by G and S units. The phenyl glycosidic (PhGlc) bond is the predominant LCC linkage under alkali-stabilized conditions. In addition, Gramineae PhGlc types are more numerous compared to hardwoods. The results of the study provide insights into the differences in the chemical composition and structural features of LCC in different plants and provide important guidance for the optimization of the process of purifying hemicellulose.

Cryo-EM structures of the human Elongator complex at work

tRNA modifications affect ribosomal elongation speed and co-translational folding dynamics. The Elongator complex is responsible for introducing 5-carboxymethyl at wobble uridine bases (cm5U34) in eukaryotic tRNAs. However, the structure and function of human Elongator remain poorly understood. In this study, we present a series of cryo-EM structures of human ELP123 in complex with tRNA and cofactors at four different stages of the reaction. The structures at resolutions of up to 2.9 Å together with complementary functional analyses reveal the molecular mechanism of the modification reaction. Our results show that tRNA binding exposes a universally conserved uridine at position 33 (U33), which triggers acetyl-CoA hydrolysis. We identify a series of conserved residues that are crucial for the radical-based acetylation of U34 and profile the molecular effects of patient-derived mutations. Together, we provide the high-resolution view of human Elongator and reveal its detailed mechanism of action.

Leveraging the Fragment Molecular Orbital and MM-GBSA Methods in Virtual Screening for the Discovery of Novel Non-Covalent Inhibitors Targeting the TEAD Lipid Binding Pocket.

The Hippo pathway controls organ size and homeostasis and is linked to numerous diseases, including cancer. The transcriptional enhanced associate domain (TEAD) family of transcription factors acts as a receptor for downstream effectors, namely yes-associated protein (YAP) and transcriptional co-activator with PDZ-binding motif (TAZ), which binds to various transcription factors and is essential for stimulated gene transcription. YAP/TAZ-TEAD facilitates the upregulation of multiple genes involved in evolutionary cell proliferation and survival. TEAD1-4 overexpression has been observed in different cancers in various tissues, making TEAD an attractive target for drug development. The central drug-accessible pocket of TEAD is crucial because it undergoes a post-translational modification called auto-palmitoylation. Crystal structures of the C-terminal TEAD complex with small molecules are available in the Protein Data Bank, aiding structure-based drug design. In this study, we utilized the fragment molecular orbital (FMO) method, molecular dynamics (MD) simulations, shape-based screening, and molecular mechanics-generalized Born surface area (MM-GBSA) calculations for virtual screening, and we identified a novel non-covalent inhibitor-BC-001-with IC50 = 3.7 μM in a reporter assay. Subsequently, we optimized several analogs of BC-001 and found that the optimized compound BC-011 exhibited an IC50 of 72.43 nM. These findings can be used to design effective TEAD modulators with anticancer therapeutic implications.

Dynamics of Protein-RNA Interfaces Using All-Atom Molecular Dynamics Simulations.

Facing the current challenges posed by human health diseases requires the understanding of cell machinery at a molecular level. The interplay between proteins and RNA is key for any physiological phenomenon, as well protein-RNA interactions. To understand these interactions, many experimental techniques have been developed, spanning a very wide range of spatial and temporal resolutions. In particular, the knowledge of tridimensional structures of protein-RNA complexes provides structural, mechanical, and dynamical pieces of information essential to understand their functions. To get insights into the dynamics of protein-RNA complexes, we carried out all-atom molecular dynamics simulations in explicit solvent on nine different protein-RNA complexes with different functions and interface size by taking into account the bound and unbound forms. First, we characterized structural changes upon binding and, for the RNA part, the change in the puckering. Second, we extensively analyzed the interfaces, their dynamics and structural properties, and the structural waters involved in the binding, as well as the contacts mediated by them. Based on our analysis, the interfaces rearranged during the simulation time showing alternative and stable residue-residue contacts with respect to the experimental structure.