Revealing elusive conformations of sucrose from hydrogen bond J-coupling in H2O: A combined NMR and quantum mechanics study.

Hydrogen bonding is a crucial feature of biomolecules, but its characterization in glycans dissolved in aqueous solutions is challenging due to rapid hydrogen exchange between hydroxyl groups and H2O. In principle, the scalar (J) coupling constant can reveal the relative orientation of the atoms in the molecule. In contrast to J-coupling through H-bonds reported in proteins and nucleic acids, research on J-coupling through H-bonds in glycans dissolved in water is lacking. Here, we use sucrose as a model system for H-bonding studies; its structure, which consists of glucose (Glc) and fructose (Frc), is well-studied, and it is readily available. We apply the in-phase, antiphase-HSQC-TOCSY and quantify previously unreported through H-bond J-values for Frc-OH1-Glc-OH2 in H2O. While earlier reports of Brown and Levy indicate this H-bond as having only a single direction, our reported findings indicate the potential presence of two involving these same atoms, namely, G2OH âž” F1O and F1OH âž” G2O (where F and G stand for Frc and Glc, respectively). The calculated density functional theory J-values for the G2OH âž” F1O agree with the experimental values. Additionally, we detected four other possible H-bonds in sucrose, which require different phi, psi (ϕ, ψ) torsion angles. The ϕ, ψ values are consistent with previous predictions of du Penhoat et al. and Venable et al. Our results will provide new insights into the molecular structure of sucrose and its interactions with proteins.

In situ serial crystallography facilitates 96-well plate structural analysis at low symmetry.

The advent of serial crystallography has rejuvenated and popularized room-temperature X-ray crystal structure determination. Structures determined at physiological temperature reveal protein flexibility and dynamics. In addition, challenging samples (e.g. large complexes, membrane proteins and viruses) form fragile crystals that are often difficult to harvest for cryo-crystallography. Moreover, a typical serial crystallography experiment requires a large number of microcrystals, mainly achievable through batch crystallization. Many medically relevant samples are expressed in mammalian cell lines, producing a meager quantity of protein that is incompatible with batch crystallization. This can limit the scope of serial crystallography approaches. Direct in situ data collection from a 96-well crystallization plate enables not only the identification of the best diffracting crystallization condition but also the possibility for structure determination under ambient conditions. Here, we describe an in situ serial crystallography (iSX) approach, facilitating direct measurement from crystallization plates mounted on a rapidly exchangeable universal plate holder deployed at a microfocus beamline, ID23-2, at the European Synchrotron Radiation Facility. We applied our iSX approach on a challenging project, autotaxin, a therapeutic target expressed in a stable human cell line, to determine the structure in the lowest-symmetry P1 space group at 3.0 Šresolution. Our in situ data collection strategy provided a complete dataset for structure determination while screening various crystallization conditions. Our data analysis reveals that the iSX approach is highly efficient at a microfocus beamline, improving throughput and demonstrating how crystallization plates can be routinely used as an alternative method of presenting samples for serial crystallography experiments at synchrotrons.

Human cytochrome C natural variants: Studying the membrane binding properties of G41S and Y48H by fluorescence energy transfer and molecular dynamics.

Cytochrome C (cyt C), the protein involved in oxidative phosphorylation, plays several other crucial roles necessary for both cell life and death. Studying natural variants of cyt C offers the possibility to better characterize the structure-to-function relationship that modulates the different activities of this protein. Naturally mutations in human cyt C (G41S and Y48H) occur in the protein central Ω-loop and cause thrombocytopenia 4. In this study, we have investigated the binding of such variants and of wild type (wt) cyt C to synthetic cardiolipin-containing vesicles. The mutants have a lower propensity in membrane binding, displaying higher dissociation constants with respect to the wt protein. Compressibility measurements reveal that both variants are more flexible than the wt, suggesting that the native central Ω-loop is important for the interaction with membranes. Such hypothesis is supported by molecular dynamics simulations. A minimal distance analysis indicates that in the presence of cardiolipin the central Ω-loop of the mutants is no more in contact with the membrane, as it happens instead in the case of wt cyt C. Such finding might provide a hint for the reduced membrane binding capacity of the variants and their enhanced peroxidase activity in vivo.

Synthesis of 6-deoxy-d-ido-heptopyranose-containing fragments of the Campylobacter jejuni strain CG8486 capsular polysaccharide.

Campylobacters are important causes of gastrointestinal illness and the capsular polysaccharides (CPS) they produce are key virulence factors and targets for vaccine development. We report here the synthesis of two fragments of the Campylobacter jejuni CG8486 strain CPS that contain a rare 6-deoxy-d-ido-heptopyranose residue and, in one target, two O-methyl phosphoramidate (MeOPN) motifs. The synthetic approach features the stereoselective construction of the ß-d-ido-heptopyranoside linkage via glycosylation with a ß-d-galacto-heptopyranoside donor followed by a one-pot sequential C-2 and C-3 inversion. During the syntheses, we uncovered a number of interesting conformational effects with regard to the 6-deoxy-ido-heptopyranose ring, the glycosidic linkage connecting the two monosaccharides, and the MeOPN groups.

Conformational Dynamics of Lipoxygenases and Their Interaction with Biological Membranes.

Lipoxygenases (LOXs) are a family of enzymes that includes different fatty acid oxygenases with a common tridimensional structure. The main functions of LOXs are the production of signaling compounds and the structural modifications of biological membranes. These features of LOXs, their widespread presence in all living organisms, and their involvement in human diseases have attracted the attention of the scientific community over the last decades, leading to several studies mainly focused on understanding their catalytic mechanism and designing effective inhibitors. The aim of this review is to discuss the state-of-the-art of a different, much less explored aspect of LOXs, that is, their interaction with lipid bilayers. To this end, the general architecture of six relevant LOXs (namely human 5-, 12-, and 15-LOX, rabbit 12/15-LOX, coral 8-LOX, and soybean 15-LOX), with different specificity towards the fatty acid substrates, is analyzed through the available crystallographic models. Then, their putative interface with a model membrane is examined in the frame of the conformational flexibility of LOXs, that is due to their peculiar tertiary structure. Finally, the possible future developments that emerge from the available data are discussed.

SAXS/WAXS data of conformationally flexible ribose binding protein.

Modern artificial intelligence-based protein structure prediction methods, such as Alphafold2, can predict structures of folded proteins with reasonable accuracy. However, Alphafold2 provides a static view of a protein, which does not show the conformational variability of the protein, domain movement in a multi-domain protein, or ligand-induced conformational changes it might undergo in solution. Small-angle X-ay scattering (SAXS) and wide-angle X-ray scattering (WAXS) are solution techniques that can aid in integrative modeling of conformationally flexible proteins, or in validating their predicted ensemble structures. While SAXS is sensitive to global structural features, WAXS can expand the scope of structural modeling by including information about local structural changes. We present SAXS and WAXS datasets obtained from conformationally flexible d-ribose binding protein (RBP) from Escherichia coli in the ribose bound and unbound forms. SAXS/WAXS datasets of RBP provided here may aid in method development efforts for more accurate prediction of structural ensembles of conformationally flexible proteins, and their conformational changes.

Isolated and solvated chiroptical behavior in conformationally flexible butanamines.

The nonresonant optical activity of two highly flexible aliphatic amines, (2R)-3-methyl-2-butanamine (R-MBA) and (2R)-(3,3)-dimethyl-2-butanamine (R-DMBA), has been probed under isolated and solvated conditions to examine the roles of conformational isomerism and to explore the influence of extrinsic perturbations. The optical rotatory dispersion (ORD) measured in six solvents presented uniformly negative rotatory powers over the 320-590 nm region, with the long-wavelength magnitude of chiroptical response growing nearly monotonically as the dielectric constant of the surroundings diminished. The intrinsic specific optical rotation, α λ T (in deg dm-1 [g/mL]-1 ), extracted for ambient vapor-phase samples of R-MBA [-11.031(98) and -2.29 (11)] and R-DMBA [-9.434 (72) and -1.350 (48)] at 355 and 633 nm were best reproduced by counterintuitive solvents of high polarity (yet low polarizability) like acetonitrile and methanol. Attempts to interpret observed spectral signatures quantitatively relied on the linear-response frameworks of density-functional theory (B3LYP, cam-B3LYP, and dispersion-corrected analogs) and coupled-cluster theory (CCSD), with variants of the polarizable continuum model (PCM) deployed to account for the effects of implicit solvation. Building on the identification of several low-lying equilibrium geometries (nine for R-MBA and three for R-DMBA), ensemble-averaged ORD profiles were calculated at T = 300 K by means of the independent-conformer ansatz, which enabled response properties predicted for the optimized structure of each isomer to be combined through Boltzmann-weighted population fractions derived from corresponding relative internal-energy or free-energy values, the latter of which stemmed from composite CBS-APNO and G4 analyses. Although reasonable accord between theory and experiment was realized for the isolated (vapor-phase) species, the solution-phase results were less satisfactory and tended to degrade progressively as the solvent polarity increased. These trends were attributed to solvent-mediated changes in structural parameters and energy metrics for the transition states that separate and putatively isolate the equilibrium conformations supported by the ground electronic potential-energy surface, with the resulting displacement of barrier locations and/or decrease of barrier heights compromising the underlying premise of the independent-conformer ansatz.

FLEXR: automated multi-conformer model building using electron-density map sampling.

Protein conformational dynamics that may inform biology often lie dormant in high-resolution electron-density maps. While an estimated ∼18% of side chains in high-resolution models contain alternative conformations, these are underrepresented in current PDB models due to difficulties in manually detecting, building and inspecting alternative conformers. To overcome this challenge, we developed an automated multi-conformer modeling program, FLEXR. Using Ringer-based electron-density sampling, FLEXR builds explicit multi-conformer models for refinement. Thereby, it bridges the gap of detecting hidden alternate states in electron-density maps and including them in structural models for refinement, inspection and deposition. Using a series of high-quality crystal structures (0.8-1.85 Šresolution), we show that the multi-conformer models produced by FLEXR uncover new insights that are missing in models built either manually or using current tools. Specifically, FLEXR models revealed hidden side chains and backbone conformations in ligand-binding sites that may redefine protein-ligand binding mechanisms. Ultimately, the tool facilitates crystallographers with opportunities to include explicit multi-conformer states in their high-resolution crystallographic models. One key advantage is that such models may better reflect interesting higher energy features in electron-density maps that are rarely consulted by the community at large, which can then be productively used for ligand discovery downstream. FLEXR is open source and publicly available on GitHub at https://github.com/TheFischerLab/FLEXR.

Pan-HSP90 ligand binding reveals isoform-specific differences in plasticity and water networks.

Isoforms of heat shock protein 90 (HSP90) fold oncoproteins that facilitate all 10 hallmarks of cancer. However, its promise as a therapeutic target remains unfulfilled as there is still no FDA-approved drug targeting HSP90 in disease. Among the reasons hindering progress are side effects caused by pan-HSP90 inhibition. Selective targeting of the four isoforms is challenging due to high sequence and structural similarity. Surprisingly, while decades of drug discovery efforts have produced almost 400 human HSP90 structures, no single ligand has been structurally characterized across all four human isoforms to date, which could reveal structural differences to achieve selectivity. To better understand the HSP90 landscape relevant for ligand binding and design we take a three-pronged approach. First, we solved the first complete set of structures of a single ligand bound to all four human isoforms. This enabled a systematic comparison of how side-chains and water networks respond to ligand binding across isoforms. Second, we expanded our analysis to publicly available, incomplete isoform-ligand series with distinct ligand chemistry. This highlighted general trends of protein and water mobility that differ among isoforms and impact ligand binding. Third, we further probed the Hsp90α conformational landscape for accommodating a congeneric series containing the purine scaffold common to HSP90 inhibitors. This revealed how minor ligand modifications flip ligand poses and perturb water and protein conformations. Taken together, this work illustrates how a systematic approach can shed new light on an "old" target and reveal hidden isoform-specific accommodations of congeneric ligands that may be exploited in ligand discovery and design.

Diverse strategies adopted by nature for regulating purine biosynthesis via fine-tuning of purine metabolic enzymes.

Purine nucleotides, generated by de novo synthesis and salvage pathways, are essential for metabolism and act as building blocks of genetic material. To avoid an imbalance in the nucleotide pool, nature has devised several strategies to regulate/tune the catalytic performance of key purine metabolic enzymes. Here, we discuss some recent examples, such as stress-regulating alarmones that bind to select pathway enzymes, huge ensembles like dynamic metabolons and self-assembled filaments that highlight the layered fine-control prevalent in the purine metabolic pathway to fulfill requisite purine demands. Examples of enzymes that turn-on only under allosteric control, are regulated via long-distance communication that facilitates transient conduits have additionally been explored.

Dissecting the general mechanisms of protein cage self-assembly by coarse-grained simulations.

The development of artificial protein cages has recently gained massive attention due to their promising application prospect as novel delivery vehicles for therapeutics. These nanoparticles are formed through a process called self-assembly, in which individual subunits spontaneously arrange into highly ordered patterns via non-covalent but specific interactions. Therefore, the first step toward the design of novel engineered protein cages is to understand the general mechanisms of their self-assembling dynamics. Here we have developed a new computational method to tackle this problem. Our method is based on a coarse-grained model and a diffusion-reaction simulation algorithm. Using a tetrahedral cage as test model, we showed that self-assembly of protein cage requires of a seeding process in which specific configurations of kinetic intermediate states are identified. We further found that there is a critical concentration to trigger self-assembly of protein cages. This critical concentration allows that cages can only be successfully assembled under a persistently high concentration. Additionally, phase diagram of self-assembly has been constructed by systematically testing the model across a wide range of binding parameters. Finally, our simulations demonstrated the importance of protein's structural flexibility in regulating the dynamics of cage assembly. In summary, this study throws lights on the general principles underlying self-assembly of large cage-like protein complexes and thus provides insights to design new nanomaterials.

Monomer and Excimer Emission in a Conformational and Stacking-Adaptable Molecular System.

A design strategy that combines molecular conformation, alkyl chain length, and charge-transfer effects has been developed to obtain conformational and stacking-adaptable donor-acceptor-π type molecules for precisely regulating the monomer and excimer emission in a single luminous platform under different environments. These fluorophores can exhibit bright monomer emissions when they are in the dispersed state based on their planar conformation. However, when the luminous molecules with short alkyl side chains are in the crystalline state, their molecular conformation can become distorted, further inducing strong intermolecular interactions and staggered π-π stacking for bright excimer emission. More importantly, their dispersed and aggregated states can be reversibly regulated in a phase-change fatty acid matrix, to achieve temperature-responsive fluorescence for temperature monitoring and advanced information encryption.

The molecular docking and molecular dynamics study of flavonol synthase and flavonoid 3'-monooxygenase enzymes involved for the enrichment of kaempferol.

Kaempferol is a natural flavonol that shows many pharmacological properties including anti-inflammatory, antioxidant, anticancer, antidiabetic activities etc. It has been reported in many vegetables, fruits, herbs and medicinal plants. The enzyme flavonol synthase (FLS, EC 1.14.20.6) catalyses the conversion of dihydroflavonols to flavonols. Whereas flavonoid 3'-monooxygenase (F3'H, EC 1.14.14.82) catalyses the hydroxylation of dihydroflavonol, and flavonol. FLS is involved in the synthesis of the kaempferol whereas F3'H causes degradation of kaempferol. The present study aimed to analyse the binding affinity, stability and activating activity of enzyme FLS as well as inhibitory activity of enzyme F3'H involved in the enrichment of the kaempferol using the in-silico approaches. Computational study for physico-chemical properties, conserved domain identification, 3-D structure prediction and its validation, conservation analysis, molecular docking followed by molecular dynamics analysis of FLS and F3'H, protein-activator (FLS-LIG Complex) and protein-inhibitor (F3'H-LIG Complex) complexes have been performed. Other structural analyses like root mean square fluctuation (RMSF), root mean square deviation (RMSD), surface area solvent accessibility (SASA), radius of gyration (Rg), hydrogen bond analysis, principal component analysis (PCA), Poisson-Boltzmann analysis (MM_PBSA) and the dynamic cross correlation map (DCCM) analysis to explore the structural, functional and thermodynamic stability of the proteins and the complexes were also studied. The molecular docking result showed that FLS binds strongly with the activator ascorbate (CID _54670067) while F3'H binds with the inhibitor ketoconazole (CID_456201). The most powerful inhibitor (ketoconazole for F3'H) and activator (ascorbate for FLS) is determined by computing the thermodynamic binding free energy through MM_PBSA analysis. The current work provides wide-ranging structural and functional information about FLS and F3'H enzymes showing detailed molecular mechanism of kaempferol biosynthesis and its degradation and hence kaempferol enrichment. Finding of the present work opens up new possibilities for future research towards enrichment of kaempferol by using activator (ascorbate) for FLS and inhibitor (ketoconazole) for F3'H as well as for its large-scale production using in vitro approaches.Communicated by Ramaswamy H. Sarma.

Crystal structure of ethyl 4-[(4-methyl-benz-yl)-oxy]benzoate.

The title compound, C17H18O3, crystallizes with three mol-ecules in the asymmetric unit. The mol-ecules differ in the conformation related to the eth-oxy group and in the orientation of the two phenyl rings, one of which has the eth-oxy group disordered over two positions with refined occupancies of 0.735:0.265 (9). In the crystal packing, the mol-ecules are connected by weak C-H⋯π inter-actions.

Comparison of the conformational dynamics of an N-glycan in implicit and explicit solvents.

All-atom MD simulations provide the atomic details of glycan structure in solution, but extensive sampling is required for simulating the transition between rotameric states. For carbohydrate modeling, the motions of the glycosidic linkages in a glycan are often sampled well on the sub-µs time scale. Consequently, simulations of large glycans or glycoconjugates in aqueous environments are very challenging at the atomistic level because of a huge number of water molecules, which eventually slows down the conformational sampling. An alternative to the all-atom approach is the multiscale simulation, where the glycan is in all-atom form while the solvent is treated implicitly. Here we present a comparative analysis of the implicit model with the currently used explicit model using Gaussian accelerated molecular dynamics (GaMD). Here we used a hybrid N-glycan to investigate the accuracy of the implicit solvent model. Estimating several structural parameters, namely dihedral torsional angles, puckering angles, and end-to-end distances, yields a complete classification of structural space in both solvent models. Both solvent models displayed similar conformations at the monosaccharide level, whereas minor differences were observed in the global conformation. The only major difference was observed in the inter-residue hydrogen bonds, which increased by 2-fold in the explicit solvent model. Detailed analysis of structural parameters established a balanced trade-off between the computational speed and accuracy for studying long-chain glycan molecules in an implicit solvent.

Crystal Structures of the Clostridium botulinum Neurotoxin A6 Cell Binding Domain Alone and in Complex with GD1a Reveal Significant Conformational Flexibility.

Clostridium botulinum neurotoxin A (BoNT/A) targets the soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) complex, by cleaving synaptosomal-associated protein of 25 kDa size (SNAP-25). Cleavage of SNAP-25 results in flaccid paralysis due to repression of synaptic transmission at the neuromuscular junction. This activity has been exploited to treat a range of diseases associated with hypersecretion of neurotransmitters, with formulations of BoNT/A commercially available as therapeutics. Generally, BoNT activity is facilitated by three essential domains within the molecule, the cell binding domain (HC), the translocation domain (HN), and the catalytic domain (LC). The HC, which consists of an N-terminal (HCN) and a C-terminal (HCC) subdomain, is responsible for BoNT's high target specificity where it forms a dual-receptor complex with synaptic vesicle protein 2 (SV2) and a ganglioside receptor on the surface of motor neurons. In this study, we have determined the crystal structure of botulinum neurotoxin A6 cell binding domain (HC/A6) in complex with GD1a and describe the interactions involved in ganglioside binding. We also present a new crystal form of wild type HC/A6 (crystal form II) where a large 'hinge motion' between the HCN and HCC subdomains is observed. These structures, along with a comparison to the previously determined wild type crystal structure of HC/A6 (crystal form I), reveals the degree of conformational flexibility exhibited by HC/A6.

Molecular Events behind the Selectivity and Inactivation Properties of Model NaK-Derived Ion Channels.

Y55W mutants of non-selective NaK and partly K+-selective NaK2K channels have been used to explore the conformational dynamics at the pore region of these channels as they interact with either Na+ or K+. A major conclusion is that these channels exhibit a remarkable pore conformational flexibility. Homo-FRET measurements reveal a large change in W55-W55 intersubunit distances, enabling the selectivity filter (SF) to admit different species, thus, favoring poor or no selectivity. Depending on the cation, these channels exhibit wide-open conformations of the SF in Na+, or tight induced-fit conformations in K+, most favored in the four binding sites containing NaK2K channels. Such conformational flexibility seems to arise from an altered pattern of restricting interactions between the SF and the protein scaffold behind it. Additionally, binding experiments provide clues to explain such poor selectivity. Compared to the K+-selective KcsA channel, these channels lack a high affinity K+ binding component and do not collapse in Na+. Thus, they cannot properly select K+ over competing cations, nor reject Na+ by collapsing, as K+-selective channels do. Finally, these channels do not show C-type inactivation, likely because their submillimolar K+ binding affinities prevent an efficient K+ loss from their SF, thus favoring permanently open channel states.

Water Networks Repopulate Protein-Ligand Interfaces with Temperature.

High-resolution crystal structures highlight the importance of water networks in protein-ligand interactions. However, as these are typically determined at cryogenic temperature, resulting insights may be structurally precise but not biologically accurate. By collecting 10 matched room-temperature and cryogenic datasets of the biomedical target Hsp90α, we identified changes in water networks that impact protein conformations at the ligand binding interface. Water repositioning with temperature repopulates protein ensembles and ligand interactions. We introduce Flipper conformational barcodes to identify temperature-sensitive regions in electron density maps. This revealed that temperature-responsive states coincide with ligand-responsive regions and capture unique binding signatures that disappear upon cryo-cooling. Our results have implications for discovering Hsp90 selective ligands, and, more generally, for the utility of hidden protein and water conformations in drug discovery.

Crystal structure of human acetylcholinesterase in complex with tacrine: Implications for drug discovery.

Alzheimer's disease (AD) is one of the most common, progressive neurodegenerative disorders affecting the aged populations. Though various disease pathologies have been suggested for AD, the impairment of the cholinergic system is one of the critical factors for the disease progression. Restoration of the cholinergic transmission through acetylcholinesterase (AChE) inhibitors is a promising disease modifying therapy. Being the first marketed drug for AD, tacrine reversibly inhibits AChE and thereby slows the breakdown of the chemical messenger acetylcholine (ACh) in the brain. However, the atomic level of interactions of tacrine towards human AChE (hAChE) is unknown for years. Hence, in the current study, we report the X-ray structure of hAChE-tacrine complex at 2.85 Å resolution. The conformational heterogeneity of tacrine within the electron density was addressed with the help of molecular mechanics assisted methods and the low-energy ligand configuration is reported, which provides a mechanistic explanation for the high binding affinity of tacrine towards AChE. Additionally, structural comparison of reported hAChE structures sheds light on the conformational selection and induced fit effects of various active site residues upon binding to different ligands and provides insight for future drug design campaigns against AD where AChE is a drug target.

Raman Optical Activity of N-Acetyl-L-Cysteine in Water and in Methanol: The "Clusters-in-a-Liquid" Model and ab Initio Molecular Dynamics Simulations.

Raman and Raman Optical Activity (ROA) spectra of N-acetyl-L-cysteine (NALC), a flexible chiral molecule, were measured in water and in methanol to evaluate the solvent effects. Two different solvation approaches, that is, the DFT based "clusters-in-a-liquid" solvent model and the ab initio molecular dynamics (AIMD) simulations, were applied to simulate the Raman and ROA spectra. Systematic conformational searches were carried out using a recently developed conformational searching tool, CREST, with the inclusion of polarizable continuum model of water and of methanol. The CREST candidates of NALC and the NALC-solvent complexes were re-optimized and their Raman and ROA simulations were done at the B3LYP-D3BJ/def2-TZVP and the B3LYP-aug-cc-pVDZ//cc-pVTZ levels. Also, AIMD simulations, which includes some anharmonic effects and all intermolecular interactions in solution, were performed. By empirically weighting the computed Raman and ROA spectra of each conformer, good agreements with the experimental data were achieved with both approaches, while AIMD offered some improvements in the carbonyl and in the low wavenumber regions over the static DFT approach. The pros and cons of these two different approaches for accounting the solvent effects on Raman and ROA of this flexible chiral system will also be discussed.