Modeling and Simulation of the NMDA Receptor at Coarse-Grained and Atomistic Levels.

N-Methyl-D-aspartate (NMDA) receptors are glutamate-gated excitatory channels that play essential roles in brain functions. While high-resolution structures were solved for an allosterically inhibited form of functional NMDA receptor, other key functional states (particularly the active open-channel state) have not yet been resolved at atomic resolutions. To decrypt the molecular mechanism of the NMDA receptor activation, structural modeling and simulation are instrumental in providing detailed information about the dynamics and energetics of the receptor in various functional states. In this chapter, we describe coarse-grained modeling of the NMDA receptor using an elastic network model and related modeling/analysis tools (e.g., normal mode analysis, flexibility and hotspot analysis, cryo-EM flexible fitting, and transition pathway modeling) based on available structures. Additionally, we show how to build an atomistic model of the active-state receptor with targeted molecular dynamics (MD) simulation and explore its energetics and dynamics with conventional MD simulation. Taken together, these modeling and simulation can offer rich structural and dynamic information which will guide experimental studies of the activation of this key receptor.

Interplay of structural preorganization and conformational sampling in UDP-glucuronic acid 4-epimerase catalysis.

Understanding enzyme catalysis as connected to protein motions is a major challenge. Here, based on temperature kinetic studies combined with isotope effect measurements, we obtain energetic description of C-H activation in NAD-dependent UDP-glucuronic acid C4 epimerase. Approach from the ensemble-averaged ground state (GS) to the transition state-like reactive conformation (TSRC) involves, alongside uptake of heat ( Δ H ‡ = 54 kJ mol-1), significant loss in entropy ( - T Δ S ‡ = 20 kJ mol-1; 298 K) and negative activation heat capacity ( Δ C p ‡ = -0.64 kJ mol-1 K-1). Thermodynamic changes suggest the requirement for restricting configurational freedom at the GS to populate the TSRC. Enzyme variants affecting the electrostatic GS preorganization reveal active-site interactions important for precise TSRC sampling and H-transfer. Collectively, our study captures thermodynamic effects associated with TSRC sampling and establishes rigid positioning for C-H activation in an enzyme active site that requires conformational flexibility in fulfillment of its natural epimerase function.

Structural basis and synergism of ATP and Na+ activation in bacterial K+ uptake system KtrAB.

The K+ uptake system KtrAB is essential for bacterial survival in low K+ environments. The activity of KtrAB is regulated by nucleotides and Na+. Previous studies proposed a putative gating mechanism of KtrB regulated by KtrA upon binding to ATP or ADP. However, how Na+ activates KtrAB and the Na+ binding site remain unknown. Here we present the cryo-EM structures of ATP- and ADP-bound KtrAB from Bacillus subtilis (BsKtrAB) both solved at 2.8 Å. A cryo-EM density at the intra-dimer interface of ATP-KtrA was identified as Na+, as supported by X-ray crystallography and ICP-MS. Thermostability assays and functional studies demonstrated that Na+ binding stabilizes the ATP-bound BsKtrAB complex and enhances its K+ flux activity. Comparing ATP- and ADP-BsKtrAB structures suggests that BsKtrB Arg417 and Phe91 serve as a channel gate. The synergism of ATP and Na+ in activating BsKtrAB is likely applicable to Na+-activated K+ channels in central nervous system.

Identification of domains in Plasmodium falciparum proteins of unknown function using DALI search on AlphaFold predictions.

Plasmodium falciparum, the causative agent of malaria, poses a significant global health challenge, yet much of its biology remains elusive. A third of the genes in the P. falciparum genome lack annotations regarding their function, impeding our understanding of the parasite's biology. In this study, we employ structure predictions and the DALI search algorithm to analyse proteins encoded by uncharacterized genes in the reference strain 3D7 of P. falciparum. By comparing AlphaFold predictions to experimentally determined protein structures in the Protein Data Bank, we found similarities to known domains in 353 proteins of unknown function, shedding light on their potential functions. The lowest-scoring 5% of similarities were additionally validated using the size-independent TM-align algorithm, confirming the detected similarities in 88% of the cases. Notably, in over 70 P. falciparum proteins the presence of domains resembling heptatricopeptide repeats, which are typically involvement in RNA binding and processing, was detected. This suggests this family, which is important in transcription in mitochondria and apicoplasts, is much larger in Plasmodium parasites than previously thought. The results of this domain search provide a resource to the malaria research community that is expected to inform and enable experimental studies.

Structural basis for the intracellular regulation of ferritin degradation.

The interaction between nuclear receptor coactivator 4 (NCOA4) and the iron storage protein ferritin is a crucial component of cellular iron homeostasis. The binding of NCOA4 to the FTH1 subunits of ferritin initiates ferritinophagy-a ferritin-specific autophagic pathway leading to the release of the iron stored inside ferritin. The dysregulation of NCOA4 is associated with several diseases, including neurodegenerative disorders and cancer, highlighting the NCOA4-ferritin interface as a prime target for drug development. Here, we present the cryo-EM structure of the NCOA4-FTH1 interface, resolving 16 amino acids of NCOA4 that are crucial for the interaction. The characterization of mutants, designed to modulate the NCOA4-FTH1 interaction, is used to validate the significance of the different features of the binding site. Our results explain the role of the large solvent-exposed hydrophobic patch found on the surface of FTH1 and pave the way for the rational development of ferritinophagy modulators.

Structural basis for EROS binding to human phagocyte NADPH oxidase NOX2.

Essential for reactive oxygen species (EROS) protein is a recently identified molecular chaperone of NOX2 (gp91phox), the catalytic subunit of phagocyte NADPH oxidase. Deficiency in EROS is a recently identified cause for chronic granulomatous disease, a genetic disorder with recurrent bacterial and fungal infections. Here, we report a cryo-EM structure of the EROS-NOX2-p22phox heterotrimeric complex at an overall resolution of 3.56Å. EROS and p22phox are situated on the opposite sides of NOX2, and there is no direct contact between them. EROS associates with NOX2 through two antiparallel transmembrane (TM) α-helices and multiple ß-strands that form hydrogen bonds with the cytoplasmic domain of NOX2. EROS binding induces a 79° upward bend of TM2 and a 48° backward rotation of the lower part of TM6 in NOX2, resulting in an increase in the distance between the two hemes and a shift of the binding site for flavin adenine dinucleotide (FAD). These conformational changes are expected to compromise superoxide production by NOX2, suggesting that the EROS-bound NOX2 is in a protected state against activation. Phorbol myristate acetate, an activator of NOX2 in vitro, is able to induce dissociation of NOX2 from EROS with concurrent increase in FAD binding and superoxide production in a transfected COS-7 model. In differentiated neutrophil-like HL-60, the majority of NOX2 on the cell surface is dissociated with EROS. Further studies are required to delineate how EROS dissociates from NOX2 during its transport to cell surface, which may be a potential mechanism for regulation of NOX2 activation.

Principles of peptide selection by the transporter associated with antigen processing.

Our ability to fight pathogens relies on major histocompatibility complex class I (MHC-I) molecules presenting diverse antigens on the surface of diseased cells. The transporter associated with antigen processing (TAP) transports nearly the entire repertoire of antigenic peptides into the endoplasmic reticulum for MHC-I loading. How TAP transports peptides specific for MHC-I is unclear. In this study, we used cryo-EM to determine a series of structures of human TAP, both in the absence and presence of peptides with various sequences and lengths. The structures revealed that peptides of eight or nine residues in length bind in a similarly extended conformation, despite having little sequence overlap. We also identified two peptide-anchoring pockets on either side of the transmembrane cavity, each engaging one end of a peptide with primarily main chain atoms. Occupation of both pockets results in a global conformational change in TAP, bringing the two halves of the transporter closer together to prime it for isomerization and ATP hydrolysis. Shorter peptides are able to bind to each pocket separately but are not long enough to bridge the cavity to bind to both simultaneously. Mutations that disrupt hydrogen bonds with the N and C termini of peptides almost abolish MHC-I surface expression. Our findings reveal that TAP functions as a molecular caliper that selects peptides according to length rather than sequence, providing antigen diversity for MHC-I presentation.

Natural history of eukaryotic DNA viruses with double jelly-roll major capsid proteins.

The phylum Preplasmiviricota (kingdom Bamfordvirae, realm Varidnaviria) is a broad assemblage of diverse viruses with comparatively short double-stranded DNA genomes (<50 kbp) that produce icosahedral capsids built from double jelly-roll major capsid proteins. Preplasmiviricots infect hosts from all cellular domains, testifying to their ancient origin, and, in particular, are associated with six of the seven supergroups of eukaryotes. Preplasmiviricots comprise four major groups of viruses, namely, polintons, polinton-like viruses (PLVs), virophages, and adenovirids. We used protein structure modeling and analysis to show that protein-primed DNA polymerases (pPolBs) of polintons, virophages, and cytoplasmic linear plasmids encompass an N-terminal domain homologous to the terminal proteins (TPs) of prokaryotic PRD1-like tectivirids and eukaryotic adenovirids that are involved in protein-primed replication initiation, followed by a viral ovarian tumor-like cysteine deubiquitinylase (vOTU) domain. The vOTU domain is likely responsible for the cleavage of the TP from the large pPolB polypeptide and is inactivated in adenovirids, in which TP is a separate protein. Many PLVs and transpovirons encode a distinct derivative of polinton-like pPolB that retains the TP, vOTU, and pPolB polymerization palm domains but lacks the exonuclease domain and instead contains a superfamily 1 helicase domain. Analysis of the presence/absence and inactivation of the vOTU domains and replacement of pPolB with other DNA polymerases in eukaryotic preplasmiviricots enabled us to outline a complete scenario for their origin and evolution.

Dynamics underlie the drug recognition mechanism by the efflux transporter EmrE.

The multidrug efflux transporter EmrE from Escherichia coli requires anionic residues in the substrate binding pocket for coupling drug transport with the proton motive force. Here, we show how protonation of a single membrane embedded glutamate residue (Glu14) within the homodimer of EmrE modulates the structure and dynamics in an allosteric manner using NMR spectroscopy. The structure of EmrE in the Glu14 protonated state displays a partially occluded conformation that is inaccessible for drug binding by the presence of aromatic residues in the binding pocket. Deprotonation of a single Glu14 residue in one monomer induces an equilibrium shift toward the open state by altering its side chain position and that of a nearby tryptophan residue. This structural change promotes an open conformation that facilitates drug binding through a conformational selection mechanism and increases the binding affinity by approximately 2000-fold. The prevalence of proton-coupled exchange in efflux systems suggests a mechanism that may be shared in other antiporters where acid/base chemistry modulates access of drugs to the substrate binding pocket.

Structure and distinct supramolecular organization of a PSII-ACPII dimer from a cryptophyte alga Chroomonas placoidea.

Cryptophyte algae are an evolutionarily distinct and ecologically important group of photosynthetic unicellular eukaryotes. Photosystem II (PSII) of cryptophyte algae associates with alloxanthin chlorophyll a/c-binding proteins (ACPs) to act as the peripheral light-harvesting system, whose supramolecular organization is unknown. Here, we purify the PSII-ACPII supercomplex from a cryptophyte alga Chroomonas placoidea (C. placoidea), and analyze its structure at a resolution of 2.47 Å using cryo-electron microscopy. This structure reveals a dimeric organization of PSII-ACPII containing two PSII core monomers flanked by six symmetrically arranged ACPII subunits. The PSII core is conserved whereas the organization of ACPII subunits exhibits a distinct pattern, different from those observed so far in PSII of other algae and higher plants. Furthermore, we find a Chl a-binding antenna subunit, CCPII-S, which mediates interaction of ACPII with the PSII core. These results provide a structural basis for the assembly of antennas within the supercomplex and possible excitation energy transfer pathways in cryptophyte algal PSII, shedding light on the diversity of supramolecular organization of photosynthetic machinery.

Regiodivergent biosynthesis of bridged bicyclononanes.

Medicinal compounds from plants include bicyclo[3.3.1]nonane derivatives, the majority of which are polycyclic polyprenylated acylphloroglucinols (PPAPs). Prototype molecules are hyperforin, the antidepressant constituent of St. John's wort, and garcinol, a potential anticancer compound. Their complex structures have inspired innovative chemical syntheses, however, their biosynthesis in plants is still enigmatic. PPAPs are divided into two subclasses, named type A and B. Here we identify both types in Hypericum sampsonii plants and isolate two enzymes that regiodivergently convert a common precursor to pivotal type A and B products. Molecular modelling and substrate docking studies reveal inverted substrate binding modes in the two active site cavities. We identify amino acids that stabilize these alternative binding scenarios and use reciprocal mutagenesis to interconvert the enzymatic activities. Our studies elucidate the unique biochemistry that yields type A and B bicyclo[3.3.1]nonane cores in plants, thereby providing key building blocks for biotechnological efforts to sustainably produce these complex compounds for preclinical development.

Binding energies and hydrogen bonds effects on DNA-cisplatin interactions: a DFT-xTB study.

CONTEXT: A systematic study of hydrogen bonds in base pairs and the interaction of cisplatin with DNA fragments was carried out. Structure, binding energies, and electron density were analyzed. xTB has proven to be an accurate method for obtaining structures and binding energies in DNA structures. Our xTB values for DNA base binding energy were in the same order and in some cases better than CAM-B3LYP values compared to experimental values. Double-stranded DNA-cisplatin structures have been calculated and the hydrogen bonds of water molecules are a decisive factor contributing to the preference for the cisplatin-Guanine interaction. Higher values of the water hydrogen bonding energies were obtained in cisplatin-Guanine structures. Furthermore, the electrostatic potential was used to investigate and improve the analysis of DNA-cisplatin structures. METHODS: We applied the xTB method and the CAM-B3LYP functional combined with def2-SVP basis set to perform and analyze of the bonding energies of the cisplatin interaction and the effects of the hydrogen bonds. Results were calculated employing the xTB and the ORCA software.

Quantitative entropy-enthalpy compensation in intraprotein interactions from model compound data.

Many small globular proteins exist in only two states-the physiologically relevant folded state and an inactive unfolded state. The active state is stabilized by numerous weak attractive contacts, including hydrogen bonds, other polar interactions, and the hydrophobic effect. Knowledge of these interactions is key to understanding the fundamental equilibrium thermodynamics of protein folding and stability. We focus on one such interaction, that between amide and aromatic groups. We provide a statistically convincing case for quantitative, linear entropy-enthalpy compensation in forming aromatic-amide interactions using published model compound transfer-free energy data.

TRAF1 from a Structural Perspective.

Tumor necrosis factor receptor-associated factor (TRAF) proteins play pivotal roles in a multitude of cellular signaling pathways, encompassing immune response, cell fate determination, development, and thrombosis. Their involvement in these processes hinges largely on their ability to interact directly with diverse receptors via the TRAF domain. Given the limited binding interface, understanding how specific TRAF domains engage with various receptors and how structurally similar binding interfaces of TRAF family members adapt their distinct binding partners has been the subject of extensive structural investigations over several decades. This review presents an in-depth exploration of the current insights into the structural and molecular diversity exhibited by the TRAF domain and TRAF-binding motifs across a range of receptors, with a specific focus on TRAF1.

A Survey of Deep Learning Methods for Estimating the Accuracy of Protein Quaternary Structure Models.

The quality prediction of quaternary structure models of a protein complex, in the absence of its true structure, is known as the Estimation of Model Accuracy (EMA). EMA is useful for ranking predicted protein complex structures and using them appropriately in biomedical research, such as protein-protein interaction studies, protein design, and drug discovery. With the advent of more accurate protein complex (multimer) prediction tools, such as AlphaFold2-Multimer and ESMFold, the estimation of the accuracy of protein complex structures has attracted increasing attention. Many deep learning methods have been developed to tackle this problem; however, there is a noticeable absence of a comprehensive overview of these methods to facilitate future development. Addressing this gap, we present a review of deep learning EMA methods for protein complex structures developed in the past several years, analyzing their methodologies, data and feature construction. We also provide a prospective summary of some potential new developments for further improving the accuracy of the EMA methods.

Installation of LYRM proteins in early eukaryotes to regulate the metabolic capacity of the emerging mitochondrion.

Core mitochondrial processes such as the electron transport chain, protein translation and the formation of Fe-S clusters (ISC) are of prokaryotic origin and were present in the bacterial ancestor of mitochondria. In animal and fungal models, a family of small Leu-Tyr-Arg motif-containing proteins (LYRMs) uniformly regulates the function of mitochondrial complexes involved in these processes. The action of LYRMs is contingent upon their binding to the acylated form of acyl carrier protein (ACP). This study demonstrates that LYRMs are structurally and evolutionarily related proteins characterized by a core triplet of α-helices. Their widespread distribution across eukaryotes suggests that 12 specialized LYRMs were likely present in the last eukaryotic common ancestor to regulate the assembly and folding of the subunits that are conserved in bacteria but that lack LYRM homologues. The secondary reduction of mitochondria to anoxic environments has rendered the function of LYRMs and their interaction with acylated ACP dispensable. Consequently, these findings strongly suggest that early eukaryotes installed LYRMs in aerobic mitochondria as orchestrated switches, essential for regulating core metabolism and ATP production.

Evaluation of DFT methods for predicting geometries and NMR spectra of Bi(III) dithiocarbamate complexes with antitumor properties.

CONTEXT: Bismuth complexes with dithiocarbamate ligands have attracted attention because of their biological applications, such as antimicrobial, antileishmanial, and anticancer properties. These complexes have high cytotoxic activity against cancer cells, being more active than the standard drugs cisplatin, doxorubicin, and tamoxifen. In the present study, we investigated the ability of some DFT methods to reproduce the geometries and NMR spectra of the Bi(III) dithiocarbamate complexes, selected based on their proven antitumor activity. Our investigation revealed that the M06-L/def2-TZVP/ECP/CPCM method presented good accuracy in predicting geometries, while the TPSSh/def2-SVP/ECP/CPCM method proved effective in analyzing the 13C NMR spectra of these molecules. In general, all examined methods exhibited comparable performance in predicting 1H NMR signals. METHODS: Calculations were performed with the Gaussian 09 program using the def2-SVP and def2-TZVP basis sets, employing relativistic effective core potential (ECP) for Bi and using the CPCM solvent model. The exchange-correlation functionals BP86, PBE, OLYP, M06-L, B3LYP, B3LYP-D3, M06-2X, TPSSh, CAM-B3LYP, and ωB97XD were used in the study. Geometry optimizations were started from crystallographic structures available at the Cambridge Structural Database. The theoretical results were compared with experimental data using the mean root-mean-square deviation (RMSD), mean absolute deviations (MAD), and linear correlation coefficient (R2).

New tetrahydroisoquinoline-4-carbonitrile derivatives as potent agents against cyclin-dependent kinases, crystal structures, and computational studies.

The synthesis of two new hexahydroisoquinoline-4-carbonitrile derivatives (3a and 3b) is reported along with spectroscopic data and their crystal structures. In compound 3a, the intramolecular O-H···O hydrogen bond constraints the acetyl and hydroxyl groups to be syn. In the crystal, inversion dimers are generated by C-H···O hydrogen bonds and are connected into layers parallel to (10-1) by additional C-H···O hydrogen bonds. The layers are stacked with Cl···S contacts 0.17 Å less than the sum of the respective van der Waals radii. The conformation of the compound 3b is partially determined by the intramolecular O-H···O hydrogen bond. A puckering analysis of the tetrahydroisoquinoline unit was performed. In the crystal, O-H···O and C-H···O hydrogen bonds together with C-H···π(ring) interactions form layers parallel to (01-1) which pack with normal van der Waals interactions. To understand the binding efficiency and stability of the title molecules, molecular docking, and 100 ns dynamic simulation analyses were performed with CDK5A1. To rationalize their structure-activity relationship(s), a DFT study at the B3LYP/6-311++G** theoretical level was also done. The 3D Hirshfled surfaces were also taken to investigate the crystal packings of both compounds. In addition, their ADMET properties were explored.Communicated by Ramaswamy H. Sarma.

On the role of native contact cooperativity in protein folding.

The consistency of energy landscape theory predictions with available experimental data, as well as direct evidence from molecular simulations, have shown that protein folding mechanisms are largely determined by the contacts present in the native structure. As expected, native contacts are generally energetically favorable. However, there are usually at least as many energetically favorable nonnative pairs owing to the greater number of possible nonnative interactions. This apparent frustration must therefore be reduced by the greater cooperativity of native interactions. In this work, we analyze the statistics of contacts in the unbiased all-atom folding trajectories obtained by Shaw and coworkers, focusing on the unfolded state. By computing mutual cooperativities between contacts formed in the unfolded state, we show that native contacts form the most cooperative pairs, while cooperativities among nonnative or between native and nonnative contacts are typically much less favorable or even anticooperative. Furthermore, we show that the largest network of cooperative interactions observed in the unfolded state consists mainly of native contacts, suggesting that this set of mutually reinforcing interactions has evolved to stabilize the native state.

[Advances in the application of AlphaFold2: a protein structure prediction model].

Protein structure prediction is an important research field in life sciences and medicine, and it is also a key application scenario of artificial intelligence in scientific research. AlphaFold2 is a protein structure prediction system developed by DeepMind based on deep learning, capable of efficiently generating the atomic-scale spatial structure of a protein from the amino acid sequence. It has demonstrated superior performance in the prediction of protein structures since its inception, thus attracting much attention and research. This paper introduces the model architecture, highlights, limitations, and application progress of AlphaFold2. Furthermore, it briefs the capabilities, highlights, and limitations of several other types of protein structure prediction models and prospects the future development direction in this field.