Recent studies suggest CNNM2 (cyclin M2) to be part of the long-sought basolateral Mg2+ extruder at the renal distal convoluted tubule, or its regulator. In the present study, we explore structural features and ligand-binding capacities of the Bateman module of CNNM2 (residues 429-584), an intracellular domain structurally equivalent to the region involved in Mg2+ handling by the bacterial Mg2+ transporter MgtE, and AMP binding by the Mg2+ efflux protein CorC. Additionally, we studied the structural impact of the pathogenic mutation T568I located in this region. Our crystal structures reveal that nucleotides such as AMP, ADP or ATP bind at only one of the two cavities present in CNNM2429-584. Mg2+ favours ATP binding by alleviating the otherwise negative charge repulsion existing between acidic residues and the polyphosphate group of ATP. In crystals CNNM2429-584 forms parallel dimers, commonly referred to as CBS (cystathionine β-synthase) modules. Interestingly, nucleotide binding triggers a conformational change in the CBS module from a twisted towards a flat disc-like structure that mostly affects the structural elements connecting the Bateman module with the transmembrane region. We furthermore show that the T568I mutation, which causes dominant hypomagnesaemia, mimics the structural effect induced by nucleotide binding. The results of the present study suggest that the T568I mutation exerts its pathogenic effect in humans by constraining the conformational equilibrium of the CBS module of CNNM2, which becomes 'locked' in its flat form.

RNA hairpins capped by 5'-GNRA-3' or 5'-UNCG-3' tetraloops (TLs) are prominent RNA structural motifs. Despite their small size, a wealth of experimental data, and recent progress in theoretical simulations of their structural dynamics and folding, our understanding of the folding and unfolding processes of these small RNA elements is still limited. Theoretical description of the folding and unfolding processes requires robust sampling, which can be achieved by either an exhaustive time scale in standard molecular dynamics simulations or sophisticated enhanced sampling methods, using temperature acceleration or biasing potentials. Here, we study structural dynamics of 5'-GNRA-3' and 5'-UNCG-3' TLs by 15-μs-long standard simulations and a series of well-tempered metadynamics, attempting to accelerate sampling by bias in a few chosen collective variables (CVs). Both methods provide useful insights. The unfolding and refolding mechanisms of the GNRA TL observed by well-tempered metadynamics agree with the (reverse) folding mechanism suggested by recent replica exchange molecular dynamics simulations. The orientation of the glycosidic bond of the GL4 nucleobase is critical for the UUCG TL folding pathway, and our data strongly support the hypothesis that GL4-anti forms a kinetic trap along the folding pathway. Along with giving useful insight, our study also demonstrates that using only a few CVs apparently does not capture the full folding landscape of the RNA TLs. Despite using several sophisticated selections of the CVs, formation of the loop appears to remain a hidden variable, preventing a full convergence of the metadynamics. Finally, our data suggest that the unfolded state might be overstabilized by the force fields used.

We propose an approach that combines an extraction of collective motions of a molecular system with a sampling of its free energy surface. A recently introduced method of metadynamics allows exploration of the free energy surface of a molecular system by means of coarse-grained dynamics with flooding of free energy minima. This free energy surface is defined as a function of a set of collective variables (e.g., interatomic distances, angles, torsions, and others). In this study, essential coordinates determined by essential dynamics (principle component analysis) were used as collective variables in metadynamics. First, dynamics of the model system (explicitly solvated alanine dipeptide, Ace-Ala-Nme) was simulated by a classical molecular dynamics simulation. The trajectory (1 ns) was then analyzed by essential dynamics to obtain essential coordinates. The free energy surface as a function of the first and second essential coordinates was then explored by metadynamics. The resulting free energy surface is in agreement with other studies of this system. We propose that a combination of these two methods (metadynamics and essential dynamics) has great potential in studies of conformational changes in peptides and proteins.

Human cytosolic thymidine kinase (hTK1) has proven to be a suitable target for noninvasive imaging of cancer cell proliferation using radiolabeled substrates such as [ (18)F]fluorothymidine ([ (18)F]FLT). However, a thymidine tracer useful for single photon emission tomography (SPECT) based on the inexpensive radionuclide technetium-99m would be of significant interest. In this work, a series of thymidine derivatives labeled with the organometallic [M(CO) 3] (+) core (M = (99m)Tc, Re) were synthesized. Neutral, cationic, and anionic complexes were readily formed in aqueous media, and all were substrates of recombinant hTK1 when incubated with ATP. The neutral complexes were phosphorylated to a greater extent than the charged complexes. The extent of phosphorylation was further improved by increasing the spacer length separating thymidine and the organometallic core. A molecular dynamics simulation study performed with a modified hTK1 structure supported the experimental findings. In vitro cell internalization experiments performed in a human neuroblastoma cell line (SKNMC) showed low uptake of the charged complexes but significant uptake for the neutral, lipophilic complexes with a log P value >1.

Many carbohydrate-binding proteins contain aromatic amino acid residues in their binding sites. These residues interact with carbohydrates in a stacking geometry via CH/π interactions. These interactions can be found in carbohydrate-binding proteins, including lectins, enzymes and carbohydrate transporters. Besides this, many non-protein aromatic molecules (natural as well as artificial) can bind saccharides using these interactions. Recent computational and experimental studies have shown that carbohydrate–aromatic CH/π interactions are dispersion interactions, tuned by electrostatics and partially stabilized by a hydrophobic effect in solvated systems.

Galectin-3 (Gal-3) is recognized as a prognostic marker in several cancer types. Its involvement in tumor development and proliferation makes this lectin a promising target for early cancer diagnosis and anti-cancer therapies. Gal-3 recognizes poly-N-acetyllactosamine (LacNAc)-based carbohydrate motifs of glycoproteins and glycolipids with a high specificity for internal LacNAc epitopes. This study analyzes the mode and kinetics of binding of Gal-3 to a series of multivalent neo-glycoproteins presenting complex poly-LacNAc-based oligosaccharide ligands on a scaffold of bovine serum albumin. These neo-glycoproteins rank among the strongest Gal-3 ligands reported, with Kd reaching sub-nanomolar values as determined by surface plasmon resonance. Significant differences in the binding kinetics were observed within the ligand series, showing the tetrasaccharide capped with N,N′-diacetyllactosamine (LacdiNAc) as the strongest ligand of Gal-3 in this study. A molecular model of the Gal-3 carbohydrate recognition domain with docked oligosaccharide ligands is presented that shows the relations in the binding site at the molecular level. The neo-glycoproteins presented herein may be applied for selective recognition of Gal-3 both on the cell surface and in blood serum.

The free-energy surface (FES) of glycyl-phenylalanyl-alanine (GFA) tripeptide was explored by molecular dynamics (MD) simulations in combination with high-level correlated ab initio quantum chemical calculations and metadynamics. Both the MD and metadynamics employed the tight-binding DFT-D method instead of the AMBER force field, which yielded inaccurate results. We classified the minima localised in the FESs as follows: a) the backbone-conformational arrangement; and b) the existence of a COOH...OC intramolecular H-bond (families CO(2)H(free) and CO(2)H(bonded)). Comparison with experimental results showed that the most stable minima in the FES correspond to the experimentally observed structures. Remarkably, however, we did not observe experimentally the CO(2)H(bonded) family (also predicted by metadynamics), although its stability is comparable to that of the CO(2)H(free) structures. This fact was explained by the former's short excited-state lifetime. We also carried out ab initio calculations using DFT-D and the M06-2X functional. The importance of the dispersion energy in stabilising peptide conformers is well reflected by our pioneer analysis using the DFT-SAPT method to explore the nature of the backbone/side-chain interactions.

No abstract