Methods for Genome-Wide Chromatin Interaction Analysis.

Recent analyses revealed the essential function of chromatin structure in maintaining and regulating genomic information. Advancements in microscopy, nuclear structure observation techniques, and the development of methods utilizing next-generation sequencers (NGSs) have significantly progressed these discoveries. Methods utilizing NGS enable genome-wide analysis, which is challenging with microscopy, and have elucidated concepts of important chromatin structures such as a loop structure, a domain structure called topologically associating domains (TADs), and compartments. In this chapter, I introduce chromatin interaction techniques using NGS and outline the principles and features of each method.

HOMER for Analysis of Hi-C Data and Assessment of Composite Structure of the X Chromosome.

Hi-C reads, which represent ligation events between different regions of the genome, must be processed into matrices of interaction frequencies for downstream analysis. Here, I describe a procedure for mapping Hi-C reads to the genome and conversion of mapped reads into the HOMER tag directory format and interaction matrix format for visualization with Juicebox. The method is demonstrated for the mouse composite X chromosome in which reads from the active and inactive X chromosomes are combined after mock DMSO treatment or targeted degradation of cohesin.

Single-Cell Hi-C Analysis Workflow with Pairtools.

Single-cell Hi-C (scHi-C) is a collection of protocols for studying genomic interactions within individual cells. Although data analysis for scHi-C resembles data analysis for bulk Hi-C, the unique challenges of scHi-C, such as high noise and protocol-specific biases, require specialized data processing strategies. In this tutorial chapter, we focus on using pairtools, a suite of tools optimized for scHi-C data, demonstrating its application on a Drosophila snHi-C dataset. While centered on pairtools for snHi-C data, the principles outlined are applicable across scHi-C variants with minor adjustments. This educational chapter aims to guide researchers in using open-source tools for scHi-C analysis, emphasizing critical steps of contact pair extraction, detection of ligation junctions, filtration, and deduplication.

Synthesis and conformational analysis of pyran inter-halide analogues of á´ -talose.

In this work, we describe the synthesis of halogenated pyran analogues of ᴅ-talose using a halo-divergent strategy from known 1,6-anhydro-2,3-dideoxy-2,3-difluoro-ß-ᴅ-mannopyranose. In solution and in the solid-state, all analogues adopt standard 4 C 1-like conformations despite 1,3-diaxial repulsion between the F2 and the C4 halogen. Moreover, the solid-state conformational analysis of halogenated pyrans reveals deviation in the intra-annular torsion angles arising from repulsion between the axial fluorine at C2 and the axial halogen at C4, which increases with the size of the halogen at C4 (F < Cl < Br < I). Crystal packing arrangements of pyran inter-halides show hydrogen bond acceptor and nonbonding interactions for the halogen at C4. Finally, density functional theory (DFT) calculations corroborate the preference of talose analogues to adopt a 4 C 1-like conformation and a natural bonding orbital (NBO) analysis demonstrates the effects of hyperconjugation from C-F antibonding orbitals.

Harnessing conformational drivers in drug design.

This review article explores the pivotal role of conformational drivers in the discovery of drug-like molecules and illustrates their significance through real-life examples. Understanding molecular conformation is paramount to drug hunting as it can impact on- and off-target potency, metabolism, permeability, and solubility. Each conformational driver or effector is described and exemplified in a separate section. The final section is dedicated to NMR spectroscopy and illustrates its utility as an essential tool for conformational design.

Identification of Pathogenic Missense Mutations of NF1 Using Computational Approaches.

Neurofibromatosis type 1 (NF1) is a prevalent autosomal dominant disorder caused by mutations in the NF1 gene, leading to multisystem disorders. Given the critical role of cysteine residues in protein stability and function, we aimed to identify key NF1 mutations affecting cysteine residues that significantly contribute to neurofibromatosis pathology. To identify the most critical mutations in the NF1 gene that contribute to the pathology of neurofibromatosis, we employed a sophisticated computational pipeline specifically designed to detect significant mutations affecting the NF1 gene. Our approach involved an exhaustive search of databases such as the Human Gene Mutation Database (HGMD), UniProt, and ClinVar for information on missense mutations associated with NF1. Our search yielded a total of 204 unique cysteine missense mutations. We then employed in silico prediction tools, including PredictSNP, iStable, and Align GVGD, to assess the impact of these mutations. Among the mutations, C379R, R1000C, and C1016Y stood out due to their deleterious effects on the biophysical properties of the neurofibromin protein, significantly destabilizing its structure. These mutations were subjected to further phenotyping analysis using SNPeffect 4.0, which predicted disturbances in the protein's chaperone binding sites and overall structural stability. Furthermore, to directly visualize the impact of these mutations on protein structure, we utilized AlphaFold3 to simulate both the wild-type and mutant NF1 structures, revealing the significant effects of the R1000C mutation on the protein's conformation. In conclusion, the identification of these mutations can play a pivotal role in advancing the field of precision medicine and aid in the development of effective drugs for associated diseases.

Dynamic conformation: Marching toward circular RNA function and application.

Circular RNA is a group of covalently closed, single-stranded transcripts with unique biogenesis, stability, and conformation that play distinct roles in modulating cellular functions and also possess a great potential for developing circular RNA-based therapies. Importantly, due to its circular conformation, circular RNA generates distinct intramolecular base pairing that is different from the linear transcript. In this perspective, we review how circular RNA conformation can affect its turnover and modes of action, as well as what factors can modulate circular RNA conformation. We also discuss how understanding circular RNA conformation can facilitate learning about their functions as well as the remaining technological issues to further address their conformation. These efforts will ultimately inform the design of circular RNA-based platforms for biomedical applications.

Deciphering the Conformations of Glutathione Oxidized Peptide: A Comparative NMR Study in Solution and Solid-State Environments.

Glutathione (GSH) and its oxidized dimer (GSSG) play an important role in living systems as an antioxidant, balancing the presence of reactive oxygen species (ROS). The central thiol (-S-S-) bond in GSSG can undergo free rotation, providing multiple conformations with respect to the S-S bridge. The six titratable sites of GSSG, which are influenced by pH variations, affect these conformations in solution, whereas in solids, additionally crystal packing effects come into play. In view of differing reports about the structure of GSSG in literature, we have here conducted an extensive reexamination of its conformations using NMR, and contrasting results have been obtained for solution and solid state. In solution, the existence of more than one antiparallel orientation of the monomer unit with different hydrogen bonding schemes has been indicated by NOE and amide temperature coefficient results. On the other hand, in the solid-state, a 1H-1H double-quantum (DQ) to 13C single-quantum (SQ) correlation study has confirmed a parallel orientation, consistent with the reported X-ray crystal structure. Experimentally assigned solid-state NMR resonances have been validated using GIPAW calculations incorporated in the Quantum ESPRESSO package.

NMR Coupling Constants, Karplus Equations, and Adjusted MD Statistics: Detecting Diagnostic Torsion Angles for the Solution Geometry of 6-[α-d-Mannopyranosyl]-d-Mannopyranose (Mannobiose).

The quantitative solution conformations of 2-(hydroxymethyl)-tetrahydropyran, α-methyl-d-mannopyranoside, and 6-[α-d-mannopyranosyl]-d-mannopyranose (mannobiose) are described. Parametrized Karplus equations for redundant spin pairs across the terminal ω-torsion and the glycosidic ω-torsion for mannobiose are developed, including ω/θ-hypersurfaces for the terminal hydroxymethylene group. Experimental NMR data, algorithmic spectral simulation (clustered Hamiltonian method), molecular dynamics (MD) simulations (GLYCAM06), energy minimizations by DFT, and adjusted torsion angle populations weighted over the Karplus-type equations are used. We demonstrate that spectral simulation is a powerful tool in the refinement of initial J values obtained from static GAIO DFT calculations. We also show that only as few as one of multiple redundant torsions can be diagnostic for conformational analysis of the disaccharide.

Genome-Wide Association Study on Body Conformation Traits in Xinjiang Brown Cattle.

Body conformation traits are linked to the health, longevity, reproductivity, and production performance of cattle. These traits are also crucial for herd selection and developing new breeds. This study utilized pedigree information and phenotypic (1185 records) and genomic (The resequencing of 496 Xinjiang Brown cattle generated approximately 74.9 billion reads.) data of Xinjiang Brown cattle to estimate the genetic parameters, perform factor analysis, and conduct a genome-wide association study (GWAS) for these traits. Our results indicated that most traits exhibit moderate to high heritability. The principal factors, which explained 59.12% of the total variance, effectively represented body frame, muscularity, rump, feet and legs, and mammary system traits. Their heritability estimates range from 0.17 to 0.73, with genetic correlations ranging from -0.53 to 0.33. The GWAS identified 102 significant SNPs associated with 12 body conformation traits. A few of the SNPs were located near previously reported genes and quantitative trait loci (QTLs), while others were novel. The key candidate genes such as LCORL, NCAPG, and FAM184B were annotated within 500 Kb upstream and downstream of the significant SNPs. Therefore, factor analysis can be used to simplify multidimensional conformation traits into new variables, thus reducing the computational burden. The identified candidate genes from GWAS can be incorporated into the genomic selection of Xinjiang Brown cattle, enhancing the reliability of breeding programs.

Structural Insights and Catalytic Mechanism of 3-Hydroxybutyryl-CoA Dehydrogenase from Faecalibacterium Prausnitzii A2-165.

Atopic dermatitis (AD) is characterized by a T-helper cell type 2 (Th2) inflammatory response leading to skin damage with erythema and edema. Comparative fecal sample analysis has uncovered a strong correlation between AD and Faecalibacterium prausnitzii strain A2-165, specifically associated with butyrate production. Therefore, understanding the functional mechanisms of crucial enzymes in the butyrate pathway, such as 3-hydroxybutyryl-CoA dehydrogenase of A2-165 (A2HBD), is imperative. Here, we have successfully elucidated the three-dimensional structure of A2HBD in complex with acetoacetyl-CoA and NAD+ at a resolution of 2.2Å using the PAL-11C beamline (third generation). Additionally, X-ray data of A2HBD in complex with acetoacetyl-CoA at a resolution of 1.9 Å were collected at PAL-XFEL (fourth generation) utilizing Serial Femtosecond Crystallography (SFX). The monomeric structure of A2HBD consists of two domains, N-terminal and C-terminal, with cofactor binding occurring at the N-terminal domain, while the C-terminal domain facilitates dimerization. Our findings elucidate the binding mode of NAD+ to A2HBD. Upon acetoacetyl-CoA binding, the crystal structure revealed a significant conformational change in the Clamp-roof domain (root-mean-square deviation of 2.202 Å). Notably, residue R143 plays a critical role in capturing the adenine phosphate ring, underlining its significance in substrate recognition and catalytic activity. The binding mode of acetoacetyl-CoA was also clarified, indicating its lower stability compared to NAD+. Furthermore, the conformational change of hydrophobic residues near the catalytic cavity upon substrate binding resulted in cavity shrinkage from an open to closed conformation. This study confirms the conformational changes of catalytic triads involved in the catalytic reaction and presents a proposed mechanism for substrate reduction based on structural observations.

Structural dissection of the CMP-pseudaminic acid synthetase, PseF.

Pseudaminic acid is a non-mammalian sugar found in the surface glycoconjugates of many bacteria, including several human pathogens, and is a virulence factor thought to facilitate immune evasion. The final step in the biosynthesis of the nucleotide activated form of the sugar, CMP-Pse5Ac7Ac is performed by a CMP-Pse5Ac7Ac synthetase (PseF). Here we present the biochemical and structural characterization of PseF from Aeromonas caviae (AcPseF), with AcPseF displaying metal-dependent activity over a broad pH and temperature range. Upon binding to CMP-Pse5Ac7Ac, AcPseF undergoes dynamic movements akin to other CMP-ulosonic acid synthetases. The enzyme clearly discriminates Pse5Ac7Ac from other ulosonic acids, through active site interactions with side-chain functional groups and by positioning the molecule in a hydrophobic pocket. Finally, we show that AcPseF binds the CMP-Pse5Ac7Ac side chain in the lowest energy conformation, a trend that we observed in the structures of other enzymes of this class.

Potential allosteric pockets identification of glucagon receptor based on molecular dynamics simulations.

Glucagon receptor (GCGR) is an important target for the treatment of type 2 diabetes mellitus. Although several small molecules with antagonistic activity have been discovered, so far, only one small molecule binding site has been resolved. To discover more novel allosteric pockets and allosteric molecules, we started with the unique full-length inactive conformation of GCGR and applied all-atom molecular dynamics (MD) simulations to obtain extensive dynamic conformations of the GCGR/glucagon complex. For the first time, MDpocket, FTMove and FTMap were used to detect allosteric pockets in simulation trajectories, selecting 4 stable pockets with a total of 14 structures as templates for virtual screening. From the results of virtual screening, 14 compounds were ultimately selected after a series of filtering steps. The cAMP accumulation assay indicated that compound gs6 has antagonistic activity, and MD simulations further revealed the allosteric mechanism of gs6. We are the first to identify new allosteric pockets and allosteric molecules in simulation trajectories of the GCGR/glucagon complex, providing a reference for research on other G-protein-coupled receptors (GPCR). However, there is still considerable room for improvement, such as using more simulation methods to obtain a richer set of dynamic conformations.

Covalent modification of protein by chemical probe in living cells for structural and interaction studies.

Cellular activities are predominantly carried out by proteins that can dynamically adopt different structural conformations and differentially interact with other biomolecules according to cellular needs. Chemical probes are small molecules used to selectively interact and modulate the activities of specific proteins to study their functions such as the validation of potential drug targets. The remarkable performance of AlphaFold algorithms in the prediction of protein structures has pivoted interest toward elucidating the intracellular dynamics of protein structural conformation where covalent modification of proteins by chemical probes could be used to shed light upon. However, due to the barrier to entry by cell membrane and the general unfavorable reactive conditions of the intracellular environment, most studies using reactive chemical probes are still conducted on purified proteins and cell lysates. Nevertheless, recent progresses have been made in designing chemical probes with improved membrane permeability, stability and reactivity. This paper surveys the literature on recent advancements in membrane-permeable chemical probes and their applications with protein mass spectrometry for the intracellular studies of protein structural conformations and biomolecular interactions.

Bis(amino acidato)copper(II) compounds in blood plasma: a review of computed structural properties and amino acid affinities for Cu2+ informing further pharmacological research.

Neutral bis(amino acidato)copper(II) [Cu(aa)2] coordination compounds are the physiological species of copper(II) amino acid compounds in blood plasma taking the form of bis(l-histidinato)copper(II) and mixed ternary copper(II)-l-histidine complexes, preferably with l-glutamine, l-threonine, l-asparagine, and l-cysteine. These amino acids have three functional groups that can bind metal ions: the common α-amino and carboxylate groups and a side-chain polar group. In Cu(aa)2, two coordinating groups per amino acid bind to copper(II) in-plane, while the third group can bind apically, which yields many possibilities for axial and planar bonds, that is, for bidentate and tridentate binding. So far, the experimental studies of physiological Cu(aa)2 compounds in solutions have not specified their complete geometries. This paper provides a brief review of my group's research on structural properties of physiological Cu(aa)2 calculated using the density functional theory (DFT) to locate low-energy conformers that can coexist in aqueous solutions. These DFT investigations have revealed high conformational flexibility of ternary Cu(aa)2 compounds for tridentate or bidentate chelation, which may explain copper(II) exchange reactions in the plasma and inform the development of small multifunctional copper(II)-binding drugs with several possible copper(II)-binding groups. Furthermore, our prediction of metal ion affinities for Cu2+ binding with amino-acid ligands in low-energy conformers with different coordination modes of five physiological Cu(aa)2 in aqueous solution supports the findings of their abundance in human plasma obtained with chemical speciation modelling.

Injury to the palmar supporting structures of the fetlock alters limb stiffness and fetlock angle.

BACKGROUND: In vivo measurement of limb stiffness and conformation provides a non-invasive proxy assessment of superficial digital flexor tendon (SDFT) and suspensory ligament (SL) function. Here, we compared it in fore and hindlimbs and after injury. OBJECTIVES: To compare the limb stiffness and conformation in forelimbs and hindlimbs, changes with age, and following injury to the SDFT and SL. STUDY DESIGN: Retrospective cohort study. METHODS: Limb stiffness was calculated using floor scales and an electrogoniometer taped to the dorsal fetlock. The fetlock angle and weight were simultaneously recorded five times with the limb weight-bearing and when the opposite limb was picked up (increased load). Limb stiffness of both limbs was calculated from the gradient of the regression line of angle versus load. Fetlock angle when the weight was zero was extrapolated from the graph and used as a measure of conformation. Limb stiffness was measured in uninjured forelimbs (n = 42 limbs), hindlimbs (n = 19 limbs), forelimbs with SDFT injury (n = 18) and hindlimbs with SL injury (n = 5). RESULTS: Limb stiffness correlated with weight in forelimbs as shown previously (p < 0.001) but also in hindlimbs (p = 0.006). When normalised to the horse's weight (503 kg, IQR 471.5-560), forelimb stiffness was significantly higher (22.3 [±4.5] × 10-3 degree-1) than for the hindlimb (16.4 [±4.0] × 10-3 degree-1; p < 0.001). While there were no significant differences between forelimb and hindlimb conformation in unaffected or SDFT injury, both limb stiffness and conformation was significantly greater in limbs with SL injury (p = 0.009 and p = 0.002, respectively). MAIN LIMITATIONS: Small sample size, lack of clinical data including lameness and quantification of injuries. CONCLUSIONS: Injury to the forelimb SDFT does not alter limb stiffness or conformation in the long-term, while hindlimb SL injury simultaneously increases limb stiffness and fetlock angle, suggesting an increase in SL length following injury.

Association mechanism of bicalutamide and human serum albumin for potential clinical implications.

The binding mechanism of molecular interaction between bicalutamide and human serum albumin (HSA) in a pH 7.4 phosphate buffer was studied using various spectroscopic techniques in combination with molecular modeling. Fluorescence data revealed that the fluorescence quenching of HSA by bicalutamide was a static quenching procedure. The binding constants and number of binding sites were evaluated at different temperatures. The thermodynamic parameters, ΔH and ΔS, were calculated to be 4.30 × 104 J·mol-1 and 245 J·mol-1·K-1, respectively, suggesting that the binding of bicalutamide to HSA was driven mainly by hydrophobic interactions and hydrogen bonds. The displacement studies indicated neither Sudlow's site I nor II but subdomain IB as the main binding site for bicalutamide on HSA. The binding distance between bicalutamide and HSA was determined to be 3.54 nm based on the Förster theory. Analysis of circular dichroism, synchronous, and 3D fluorescence spectra demonstrated that HSA conformation was slightly altered in the presence of bicalutamide.

Crystal structure, Hirshfeld surface analysis, DFT and the mol-ecular docking studies of 3-(2-chloro-acet-yl)-2,4,6,8-tetra-phenyl-3,7-di-azabicyclo-[3.3.1]nonan-9-one.

In the title compound, C33H29ClN2O2, the two piperidine rings of the di-aza-bicyclo moiety adopt distorted-chair conformations. Inter-molecular C-H⋯π inter-actions are mainly responsible for the crystal packing. The inter-molecular inter-actions were qu-anti-fied and analysed using Hirshfeld surface analysis, revealing that H⋯H inter-actions contribute most to the crystal packing (52.3%). The mol-ecular structure was further optimized by density functional theory (DFT) at the B3LYP/6-31 G(d,p) level and is compared with the experimentally determined mol-ecular structure in the solid state.

Shape and joint angle data for seven European horse breeds and their repeatability.

Conformation traits are important in the selection and distinction between horse breeds, but tend to be evaluated subjectively within a breed and cannot be compared between them. The horse shape space model, using a combination of 253 landmarks and semi-landmarks, provides objective information on the shape of a horse photographed from the side that can be compared between breeds. In this dataset, we are providing the full set of 253 landmarks for 1241 horses from seven breeds, including an R code file to extract joint angle information and transform the raw data into csv files for further analysis, such as breed comparisons, heritability or genome-wide association studies (single- or multibreed). The repeatability of the joint angles are also reported.

Genome-Wide Association Study of Conformation Traits in Brazilian Holstein Cattle.

The linear conformation of animals exerts an influence on health, reproduction, production, and welfare, in addition to longevity, which directly affects the profitability of milk-producing farms. The objectives of this study were (1) to perform genome-wide association studies (GWASs) of conformation traits, namely the Rump, Feet and Legs, Mammary System, Dairy Strength, and Final Classification traits, and (2) to identify genes and related pathways involved in physiological processes associated with conformation traits in Brazilian Holstein cattle. Phenotypic and genotypic data from 2339 Holstein animals distributed across the states of Rio Grande do Sul, Paraná, São Paulo, and Minas Gerais were used. The genotypic data were obtained with a 100 K SNP marker panel. The single-step genome-wide association study (ssGWAS) method was employed in the analyses. Genes close to a significant SNP were identified in an interval of 100 kb up- and downstream using the Ensembl database available in the BioMart tool. The DAVID database was used to identify the main metabolic pathways and the STRING program was employed to create the gene regulatory network. In total, 36 significant SNPs were found on 15 chromosomes; 27 of these SNPs were linked to genes that may influence the traits studied. Fourteen genes most closely related to the studied traits were identified, as well as four genes that showed interactions in important metabolic pathways such as myogenesis, adipogenesis, and angiogenesis. Among the total genes, four were associated with myogenesis (TMOD2, TMOD3, CCND2, and CTBP2), three with angiogenesis (FGF23, FGF1, and SCG3), and four with adipogenesis and body size and development (C5H12orf4, CCND2, EMILIN1, and FGF6). These results contribute to a better understanding of the biological mechanisms underlying phenotypic variability in conformation traits in Brazilian Holstein cattle.