Side-chain Interactions Research Articles

The Generation of ROS by Exposure to Trihalomethanes Promotes the IκBα/NF-κB/p65 Complex Dissociation in Human Lung Fibroblast.

Background: Disinfection by-products used to obtain drinking water, including halomethanes (HMs) such as CH2Cl2, CHCl3, and BrCHCl2, induce cytotoxicity and hyperproliferation in human lung fibroblasts (MRC-5). Enzymes such as superoxide dismutase (SOD), catalase (CAT), and glutathione peroxidase (GPx) modulate these damages through their biotransformation processes, potentially generating toxic metabolites. However, the role of the oxidative stress response in cellular hyperproliferation, modulated by nuclear factor-kappa B (NF-κB), remains unclear. Methods: In this study, MRC-5 cells were treated with these compounds to evaluate reactive oxygen species (ROS) production, lipid peroxidation, phospho-NF-κB/p65 (Ser536) levels, and the activities of SOD, CAT, and GPx. Additionally, the interactions between HMs and ROS with the IκBα/NF-κB/p65 complex were analyzed using molecular docking. Results: Correlation analysis among biomarkers revealed positive relationships between pro-oxidant damage and antioxidant responses, particularly in cells treated with CH2Cl2 and BrCHCl2. Conversely, negative relationships were observed between ROS levels and NF-κB/p65 levels in cells treated with CH2Cl2 and CHCl3. The estimated relative free energy of binding using thermodynamic integration with the p65 subunit of NF-κB was -3.3 kcal/mol for BrCHCl2, -3.5 kcal/mol for both CHCl3 and O2•, and -3.6 kcal/mol for H2O2. Conclusions: Chloride and bromide atoms were found in close contact with IPT domain residues, particularly in the RHD region involved in DNA binding. Ser281 is located within this domain, facilitating the phosphorylation of this protein. Similarly, both ROS interacted with the IPT domain in the RHD region, with H2O2 forming a side-chain oxygen interaction with Leu280 adjacent to the phosphorylation site of p65. However, the negative correlation between ROS and phospho-NF-κB/p65 suggests that steric hindrance by ROS on the C-terminal domain of NF-κB/p65 may play a role in the antioxidant response.

Residue-Specific Epitope Mapping of the PD-1/Nivolumab Interaction Using X-ray Footprinting Mass Spectrometry.

X-ray footprinting coupled with mass spectrometry (XFMS) presents a novel approach in structural biology, offering insights into protein conformation and dynamics in the solution state. The interaction of the cancer-immunotherapy monoclonal antibody nivolumab with its antigen target PD-1 was used to showcase the utility of XFMS against the previously published crystal structure of the complex. Changes in side-chain solvent accessibility, as determined by the oxidative footprint of free PD-1 versus PD-1 bound to nivolumab, agree with the binding interface side-chain interactions reported from the crystal structure of the complex. The N-linked glycosylation sites of PD-1 were confirmed through an LC-MS/MS-based deglycosylation analysis of asparagine deamidation. In addition, subtle changes in side-chain solvent accessibility were observed in the C'D loop region of PD-1 upon complex formation with nivolumab.

Geometric descriptors for beta turns.

Beta turns, in which the protein backbone abruptly changes direction over four amino acid residues, are the most common type of protein secondary structure after alpha helices and beta sheets and play key structural and functional roles. Previous work has produced classification systems for turn geometry at multiple levels of precision, but these operate in backbone dihedral-angle (Ramachandran) space, and the absence of a local Euclidean-space coordinate system and structural alignment for turns, or of any systematic Euclidean-space characterization of turn backbone shape, presents challenges for the visualization, comparison and analysis of the wide range of turn conformations and the design of turns and the structures that incorporate them. This work derives a turn-local coordinate system that implicitly aligns turns, together with a set of geometric descriptors that characterize the bulk BB shapes of turns and describe modes of structural variation not explicitly captured by existing systems. These modes are shown to be meaningful by the demonstration of clear relationships between descriptor values and the electrostatic energy of the beta-turn H-bond, the overrepresentations of key side-chain motifs, and the structural contexts of turns. Geometric turn descriptors complement Ramachandran-space classifications, and they can be used to select turn structures for compatibility with particular side-chain interactions or contexts. Potential applications include protein design and other tasks in which an enhanced Euclidean-space characterization of turns may improve understanding or performance. The web-based tools ExploreTurns, MapTurns, and ProfileTurn, available at www.betaturn.com, incorporate turn-local coordinates and turn descriptors and demonstrate their utility.

Molecular basis of hyper-thermostability in the thermophilic archaeal aldolase MfnB

Methanogenic archaea are chemolithotrophic prokaryotes that can reduce carbon dioxide with hydrogen gas to form methane. These microorganisms make a significant contribution to the global carbon cycle, with methanogenic archaea from anoxic environments estimated to contribute > 500 million tons of global methane annually. Archaeal methanogenesis is dependent on the methanofurans; aminomethylfuran containing coenzymes that act as the primary C1 acceptor molecule during carbon dioxide fixation. Although the biosynthetic pathway to the methanofurans has been elucidated, structural adaptations which confer thermotolerance to Mfn enzymes from extremophilic archaea are yet to be investigated. Here we focus on the methanofuran biosynthetic enzyme MfnB, which catalyses the condensation of two molecules of glyceralde-3-phosphate to form 4‑(hydroxymethyl)-2-furancarboxaldehyde-phosphate. In this study, MfnB enzymes from the hyperthermophile Methanocaldococcus jannaschii and the mesophile Methanococcus maripaludis have been recombinantly overexpressed and purified to homogeneity. Thermal unfolding studies, together with steady-state kinetic assays, demonstrate thermoadaptation in the M. jannaschii enzyme. Molecular dynamics simulations have been used to provide a structural explanation for the observed properties. These reveal a greater number of side chain interactions in the M. jannaschii enzyme, which may confer protection from heating effects by enforcing spatial residue constraints.

Assignment of a Physical Energy Scale for the Dimensionless Interaction Energies within the PRIME20 Peptide Model.

We present a calibration scheme to determine the conversion factors from a coarse-grained stochastic approximation Monte Carlo approach using the PRIME20 peptide interaction model into atomistic force-field interaction energies at full explicit aqueous solvation. The conversion from coarse-grained to atomistic structures was done according to our previously established inverse coarse-graining protocol. We provide a physical energy scale for both the backbone hydrogen bonding interactions and the sidechain interactions by correlating the dimensionless energy descriptors of the PRIME20 model with the energies averaged over molecular dynamics simulations. The conversion factor for these interactions turns out to be around 2kJ/mol for the backbone interactions, and zero for the sidechain interactions. We discuss these surprisingly small values in terms of their molecular interpretation.

The interaction of the azetidine thiazole side chain with the active site loop (ASL) 3 drives the evolution of IMP metallo-β-lactamase against tebipenem.

Imipenemase (IMP) metallo-β-lactamases (MBLs) hydrolyze almost all available β-lactams including carbapenems and are not inhibited by any commercially available β-lactamase inhibitor. Tebipenem (TP) pivoxil is the first orally available carbapenem and possesses a unique bicyclic azetidine thiazole moiety located at the R2 position. TP has potent in vitro activity against Enterobacterales producing extended-spectrum and/or AmpC β-lactamases. Thus far, the activity of TP against IMP-producing strains is understudied. To address this knowledge gap, we explored the structure activity relationships of IMP MBLs by investigating whether IMP-6, IMP-10, IMP-25, and IMP-78 [MBLs with expanded hydrolytic activity against meropenem (MEM)] would demonstrate enhanced activity against TP. Most of the Escherichia coli DH10B strains expressing IMP-1 variants displayed a ≥twofold MIC difference between TP and MEM, while those expressing VIM or NDM variants demonstrated comparable MICs. Catalytic efficiency (kcat/KM) values for the TP hydrolysis by IMP-1, IMP-6, IMP-10, IMP-25, and IMP-78 were significantly lower than those obtained for MEM. Molecular dynamic simulations reveal that V67F and S262G substitutions (found in IMP-78) reposition active site loop 3, ASL-3, to better accommodate the bicyclic azetidine thiazole side chain, allowing microbiological/catalytic activity to approach that of comparison MBLs used in this study. These findings suggest that modifying the R2 side chain of carbapenems can significantly impact hydrolytic stability. Furthermore, changes in conformational dynamics due to single amino acid substitutions should be used to inform drug design of novel carbapenems.

Chiral inversion induced by aromatic interactions in short peptide assembly

Although hydrophobic interactions provide the main driving force for initial peptide aggregation, their role in regulating suprastructure handedness of higher-order architectures remains largely unknown. We here interrogate the effects of hydrophobic amino acids on handedness at various assembly stages of peptide amphiphiles. Our studies reveal that relative to aliphatic side chains, aromatic side chains set the twisting directions of single β-strands due to their strong steric repulsion to the backbone, and upon packing into multi-stranded β-sheets, the side-chain aromatic interactions between strands form the aromatic ladders with a directional preference. This ordering not only leads to parallel β-sheet arrangements but also induces the chiral flipping over of single β-strands within a β-sheet. In contrast, the lack of orientational hydrophobic interactions in the assembly of aliphatic peptides implies no chiral inversion upon packing into β-sheets. This study opens an avenue to harness peptide aggregates with targeted handedness via aromatic side-chain interactions.

Structural Engineering of l-Aspartic Amphiphilic Polyesters for Enzyme-Responsive Drug Delivery and Bioimaging in Cancer Cells.

Design and development of amphiphilic polyesters based on bioresources are very important to cater to the ever-growing need for biodegradable polymers in biomedical applications. Here, we report structural engineering of enzyme-responsive amphiphilic polyesters based on l-amino acid bioresources and study their drug delivery aspects in the cancer cell line. For this purpose, an l-aspartic acid-based polyester platform is chosen, and two noncovalent forces such as hydrogen bonding and side-chain hydrophobic interactions are introduced to study their effect on the aqueous self-assembly of nanoparticles. The synthetic strategy involves the development of l-aspartic acid-based dimethyl ester monomers with acetal and stearate side chains and subjecting them to solvent-free melt polycondensation reactions to produce side-chain-functionalized polyesters in the entire composition range. Postpolymerization acid catalyst deprotection of acetal yielded hydroxyl-functionalized polyesters. Amphiphilicity of the polymer is carefully fine-tuned by varying the composition of the stearate and hydroxyl units in the polymer chains to produce self-assembly in water. Various drugs such as camptothecin (CPT), curcumin (CUR), and doxorubicin (DOX) and biomarkers like 8-hydroxypyrene-1,3,6-trisulfonic acid trisodium salt (HPTS), rose bengal (RB), and Nile red (NR) are successfully encapsulated in the polymer nanoparticles. Cytotoxicity of biodegradable polymer nanoparticles is tested in normal and breast cancer cell lines. The polymer nanoparticles are found to be highly biocompatible and delivered the anticancer drugs in the intracellular compartments of the cells.

High conductivity and dimensional stability brominated poly (aryl piperidinium) anion exchange membranes: Synergistic interaction of hydrophilic and superhydrophobic fluorinated side chains

Anion-exchange membrane fuel cells (AEMFCs) is a promising renewable energy device for hydrogen production. Anion exchange membranes (AEMs), the key component of AEMFCs, is generally needed for excellent dimensional stability and conductivity. Here, we synthesised a variety of AEMs with different number of C–F bonds to enhance the hydrophilic/hydrophobic micro-phase separation structure while maintaining dimensional stability. At 80 °C, the QPBTTP-PEG-10-F9 shows high OH− conductivity (∼179.02 mS cm−1), excellent dimensional stability (swelling ratio <19 %), and mechanical properties (elongation at break: 37 %, tensile strength: 71 MPa). In addition, QPBTTP-PEG-10-F9 shows good alkaline stability (conductivity remaining of 92.52 %) in 3 M NaOH solution for 1800 h. The peak power density (PPD) of a single cell assembled with QPBTTP-PEG-10-F9 can reach 534.79 mW cm−2 with a current density of 951 mA cm−2.

Three-Dimensional Interaction Homology: Deconstructing Residue-Residue and Residue-Lipid Interactions in Membrane Proteins.

A method is described to deconstruct the network of hydropathic interactions within and between a protein's sidechain and its environment into residue-based three-dimensional maps. These maps encode favorable and unfavorable hydrophobic and polar interactions, in terms of spatial positions for optimal interactions, relative interaction strength, as well as character. In addition, these maps are backbone angle-dependent. After map calculation and clustering, a finite number of unique residue sidechain interaction maps exist for each backbone conformation, with the number related to the residue's size and interaction complexity. Structures for soluble proteins (~749,000 residues) and membrane proteins (~387,000 residues) were analyzed, with the latter group being subdivided into three subsets related to the residue's position in the membrane protein: soluble domain, core-facing transmembrane domain, and lipid-facing transmembrane domain. This work suggests that maps representing residue types and their backbone conformation can be reassembled to optimize the medium-to-high resolution details of a protein structure. In particular, the information encoded in maps constructed from the lipid-facing transmembrane residues appears to paint a clear picture of the protein-lipid interactions that are difficult to obtain experimentally.

Alkyl chain length affecting the structure and interactions of choline-based amino acid ionic liquids: Insights from a computational study

The influences of the alkyl chain length of amino acid anions on the structure and interactions of choline-based amino acid ([Cho][AA]) ionic liquids were investigated based on density functional theory (DFT) and molecular dynamics (MD) simulations. The most stable configuration of the individual [Cho][AA] ion pair was determined, revealing negligible variations in the binding energy between choline cation ([Cho]+) and different amino acid anions ([AA]−). To further characterize the intermolecular interactions between charged entities, the atoms in molecules theory (AIM) and independent gradient model based on Hirshfeld partition (IGMH) were employed, which highlight the presence of robust hydrogen bonds between [Cho]+ and [AA] −, underscoring their significance in shaping the structure of [Cho][AA] ion pairs. The MD studies exhibited the aggregation of anions in [Cho][AA] ILs, which was attributed to side chain interactions of anions. This aggregation effect becomes more pronounced with the elongation of anionic alkyl chain. Simultaneously, radial distribution functions (RDF) analysis revealed a higher probability of hydrogen bonds occurring between [Cho]+ and [AA]− as the alkyl chain in [AA]− lengthened. Additionally, the hydrogen bond lifetimes roughly display an increasing trend as the alkyl chain of the [AA]− increases, which may be due to the reduced diffusion coefficients of ions. Hopefully, our results would contribute fundamental knowledge regarding the impact of alkyl chain length of anions on the structure and interactions in the choline-based amino acid ionic liquids.

Mechanistic Insight into the Mechanical Unfolding of the Integral Membrane Diacylglycerol Kinase.

The essential forces stabilizing membrane proteins and governing their folding and unfolding are difficult to decipher. Single-molecule atomic force spectroscopy mechanically unfolds individual membrane proteins and quantifies their dynamics and energetics. However, it remains challenging to structurally assign unfolding intermediates precisely and to deduce dominant interactions between specific residues that facilitate either the localized stabilization of these intermediates or the global assembly of membrane proteins. Here, we performed force spectroscopy experiments and multiscale molecular dynamics simulations to study the unfolding pathway of diacylglycerol kinase (DGK), a small trimeric multispan transmembrane enzyme. The remarkable agreement between experiments and simulations allowed precise structural assignment and interaction analysis of unfolding intermediates, bypassing existing limitations on structural mapping, and thus provided mechanistic explanations for the formation of these states. DGK unfolding was found to proceed with structural segments varying in size that do not correlate with its secondary structure. We identified intermolecular side-chain packing interactions as one of the major contributions to the stability of unfolding intermediates. Mutagenesis creating packing defects induced a dramatic decrease in the mechano-stability of corresponding intermediates and also in the thermo-stability of DGK trimer, in good agreement with predictions from simulations. Hence, the molecular determinants of the mechano- and thermo-stability of a membrane protein can be identified at residue resolution. The accurate structural assignment established and microscopic mechanism revealed in this work may substantially expand the scope of single-molecule studies of membrane proteins.

Single-molecule force spectroscopy shows that side chain interactions govern the mechanochemical response of polypeptide α-helices and prevent the formation of β-sheets

AFM single-molecule experiments on poly(l-glutamic acid) and poly(l-lysine) show that hydrophobic side chain interactions stabilize α-helices and inhibit the formation of a metastable β-sheet-like structure under mechanical deformation.

Electronic substituent effect on the conformation of a phenylalanine-incorporated cyclic peptide.

The Phe-incorporated cyclic peptide [cyclo(-Phe1-oxazoline2-d-Val3-thiazole4-Ile5-oxazoline6-d-Val7-thiazole8-)] is in a conformational equilibrium between square and folded forms in solution. In the folded form, a CH⋯π interaction between the Phe1 aromatic ring and the Oxz2 methyl group is observed. We endeavored to control the local conformation and thus modulate the CH⋯π interaction and flexibility of the Phe1 side chain by controlling the electronic substituent effects at the 4-position of the aromatic ring of the Phe1 residue. The effect of the 4-substituent on the global conformation was indicated by the linear relationship between the conformational free energies (ΔGo) determined through NMR-based quantification and the Hammett constants (σ). Electron-donating substituents, which had relatively strong CH⋯π interactions, promoted peptide folding by restraining the loss in entropy. Local control by the 4-substituent effects suggested that the Phe side chain exerts an entropic influence on the folding of these cyclic peptides.

Self-assembly prediction of architecture-controlled bottlebrush copolymers in solution using graph convolutional networks.

The investigation of bottlebrush copolymer self-assembly in solution involves a comprehensive approach integrating simulation and experimental research, due to their unique physical characteristics. However, the intricate architecture of bottlebrush copolymers and the diverse solvent conditions introduce a wide range of parameter spaces. In this study, we investigated the solution self-assembly behavior of bottlebrush copolymers by combining dissipative particle dynamics (DPD) simulation results and machine learning (ML) including graph convolutional networks (GCNs). The architecture of bottlebrush copolymers is encoded by graphs including connectivity, side chain length, bead types, and interaction parameters of DPD simulation. Using GCN, we accurately predicted the single chain properties of bottlebrush copolymers with over 95% accuracy. Furthermore, phase behavior was precisely predicted using these single chain properties. Shapley additive explanations (SHAP) values of single chain properties to the various self-assembly morphologies were calculated to investigate the correlation between single chain properties and morphologies. In addition, we analyzed single chain properties and phase behavior as a function of DPD interaction parameters, extracting relevant physical properties for vesicle morphology formation. This work paves the way for tailored design in solution of self-assembled nanostructures of bottlebrush copolymers, offering a GCN framework for precise prediction of self-assembly morphologies under various chain architectures and solvent conditions.

Structure-Property Relationship, Glass Transition, and Crystallization Behaviors of Conjugated Polymers.

Conjugated polymers have gained considerable interest due to their unique structures and promising applications in areas such as optoelectronics, photovoltaics, and flexible electronics. This review focuses on the structure-property relationship, glass transition, and crystallization behaviors of conjugated polymers. Understanding the relationship between the molecular structure of conjugated polymers and their properties is essential for optimizing their performance. The glass transition temperature (Tg) plays a key role in determining the processability and application of conjugated polymers. We discuss the mechanisms underlying the glass transition phenomenon and explore how side-chain interaction affects Tg. The crystallization behavior of conjugated polymers significantly impacts their mechanical and electrical properties. We investigate the nucleation and growth processes, as well as the factors that influence the crystallization process. The development of the three generations of conjugated polymers in controlling the crystalline structure and enhancing polymer ordering is also discussed. This review highlights advanced characterization techniques such as X-ray diffraction, atomic force microscopy, and thermal analysis, which provide insights into molecular ordering and polymer-crystal interfaces. This review provides an insight of the structure-property relationship, glass transition, and crystallization behaviors of conjugated polymers. It serves as a foundation for further research and development of conjugated polymer-based materials with enhanced properties and performance.

CC+ : A searchable database of validated coiled coils in PDB structures and AlphaFold2 models.

α-Helical coiled coils are common tertiary and quaternary elements of protein structure. In coiled coils, two or more α helices wrap around each other to form bundles. This apparently simple structural motif can generate many architectures and topologies. Coiled coil-forming sequences can be predicted from heptad repeats of hydrophobic and polar residues, hpphppp, although this is not always reliable. Alternatively, coiled-coil structures can be identified using the program SOCKET, which finds knobs-into-holes (KIH) packing between side chains of neighboring helices. SOCKET also classifies coiled-coil architecture and topology, thus allowing sequence-to-structure relationships to be garnered. In 2009, we used SOCKET to create a relational database of coiled-coil structures, CC+ , from the RCSB Protein Data Bank (PDB). Here, we report an update of CC+ following an update of SOCKET (to Socket2) and the recent explosion of structural data and the success of AlphaFold2 in predicting protein structures from genome sequences. With the most-stringent SOCKET parameters, CC+ contains ≈12,000 coiled-coil assemblies from experimentally determined structures, and ≈120,000 potential coiled-coil structures within single-chain models predicted by AlphaFold2 across 48 proteomes. CC+ allows these and other less-stringently defined coiled coils to be searched at various levels of structure, sequence, and side-chain interactions. The identified coiled coils can be viewed directly from CC+ using the Socket2 application, and their associated data can be downloaded for further analyses. CC+ is available freely at http://coiledcoils.chm.bris.ac.uk/CCPlus/Home.html. It will be updated automatically. We envisage that CC+ could be used to understand coiled-coil assemblies and their sequence-to-structure relationships, and to aid protein design and engineering.

Disallowed spots in protein structures

Ramachandran (ϕ, ψ) steric map was introduced in 1963 to describe available conformation space for protein structures. Subsequently, residues were observed in high-energy disallowed regions of the map. To unequivocally identify the locations of disallowed conformations of residues, we got 36 noise-free protein structures (resolution ≤1 Å, Rwork/Rfree ≤ 0.10). These stringent criteria were applied to rule out data or model errors or any crystallographic disorders. No disallowed conformation was found in the dataset. Further, we also examined disallowed conformations in a larger dataset (resolution ≤1.5 Å, devoid of any model errors, or disorders). The observed locations of disallowed residues are referred as disallowed spots. These spots include short loops of 3–5 residues, and locations where residues participate in disulfide bonding or intramolecular interactions or inter-molecular interactions with neighboring water, metals or ligands. Conformational sampling revealed that short loops in between secondary structures hardly have any opportunity to relieve from conformational strain. Residues involved in interactions, which provide energetic compensation for high-energy conformational states, were relieved from strain once the causative interaction was removed. The present study aims to identify disallowed spots in the native state of proteins, wherein residues are forced to be trapped in high-energy disallowed conformations. Moreover, it was also observed that pre-Pro, Ser, Asp, trans-Pro, Val, Asn & Gly have higher tendency to occur in disallowed conformation, which could be attributed to factors such as conformational restrictions, residue propensity of secondary structures and compensating sidechain and mainchain interactions, stabilizing turn-mimics.

Effect of sliding velocity on the nanoscale friction behaviour of articular cartilage contact interface: insights from all-atom molecular dynamics investigation

ABSTRACT This study employs molecular dynamics simulations to explore nanoscale friction behaviour as a function of varying loading and sliding speeds on a developed top-layer articular cartilage contact interface atomistic model. To investigate nanotribological behaviour, sliding speed variations on the normal force, friction force, non-bonded interaction energy and interface temperature is obtained at the inter-cartilage interface. Analysis conducted at high velocity in a simplified tissue-like hydrated environment revealed ice-like dynamic smooth sliding behaviour of protein chains when cartilage interfaces are even 3.8 Å apart. With an increase in velocity, the coefficient of friction (COF) increases significantly in a hydrated environment. Additionally, at lower loads, the effect of sliding velocity is more pronounced than at higher loads. However, results show that articular cartilage adapts to higher load and speed sliding conditions exhibiting lower friction (COF-0.03–1.17) by means of interfacial water rearrangements and protein side-chain non-bonded interactions reducing heavy shear deformation. This is attributed to an alteration in the load-bearing and friction mechanism owing to water rearrangement and adsorption at nanoconfined biointerfaces. This study provides mechanistic insights into friction mechanisms at the cartilage interface which could lead to wear-like conditions under physiological sliding contact conditions, thereby facilitating the design of better implants.

Synthesis of a 13C-methylene-labeled isoleucine precursor as a useful tool for studying protein side-chain interactions and dynamics

In this study, we present the synthesis and incorporation of a metabolic isoleucine precursor compound for selective methylene labeling. The utility of this novel α-ketoacid isotopologue is shown by incorporation into the protein Brd4-BD1, which regulates gene expression by binding to acetylated histones. High quality single quantum 13C−1 H-HSQC were obtained, as well as triple quantum HTQC spectra, which are superior in terms of significantly increased 13C-T2 times. Additionally, large chemical shift perturbations upon ligand binding were observed. Our study thus proves the great sensitivity of this precursor as a reporter for side-chain dynamic studies and for investigations of CH-π interactions in protein-ligand complexes.