Type 2 Diabetes Variants Disrupt Function of SLC16A11 through Two Distinct Mechanisms.

Type 2 diabetes (T2D) affects Latinos at twice the rate seen in populations of European descent. We recently identified a risk haplotype spanning SLC16A11 that explains ∼20% of the increased T2D prevalence in Mexico. Here, through genetic fine-mapping, we define a set of tightly linked variants likely to contain the causal allele(s). We show that variants on the T2D-associated haplotype have two distinct effects: (1) decreasing SLC16A11 expression in liver and (2) disrupting a key interaction with basigin, thereby reducing cell-surface localization. Both independent mechanisms reduce SLC16A11 function and suggest SLC16A11 is the causal gene at this locus. To gain insight into how SLC16A11 disruption impacts T2D risk, we demonstrate that SLC16A11 is a proton-coupled monocarboxylate transporter and that genetic perturbation of SLC16A11 induces changes in fatty acid and lipid metabolism that are associated with increased T2D risk. Our findings suggest that increasing SLC16A11 function could be therapeutically beneficial for T2D. VIDEO ABSTRACT.

A method to apply bond-angle constraints in molecular dynamics simulations.

An algorithm to apply bond-angle constraints in molecular dynamics simulations of macromolecules or molecular liquids is presented. It uses Cartesian coordinates and determines the Lagrange multipliers required for maintaining the constraints iteratively. It constitutes an alternative to the use of only distance constraints (DCs) between particles to maintain a particular geometry. DCs are unsuitable to maintain particular, for example, linear or flat, geometries of molecules. The proposed algorithm can easily handle bond-length, bond-angle, and dihedral-angle constraints simultaneously, as when calculating a potential of mean force along a dihedral-angle degree of freedom.

Rotational Profiler: A Fast, Automated, and Interactive Server to Derive Torsional Dihedral Potentials for Classical Molecular Simulations.

Rotational Profiler provides an analytical algorithm to compute sets of classical torsional dihedral parameters by fitting an empirical energy profile to a reference one that can be obtained experimentally or by quantum-mechanical methods. The resulting profiles are compatible with the functional forms in the most widely used biomolecular force fields (e.g., GROMOS, AMBER, OPLS, and CHARMM). The linear least-squares regression method is used to generate sets of parameters that best satisfy the fitting. Rotational Profiler is free to use, analytical, and force field/package independent. The formalism is herein described, and its usage, in an interactive and automated manner, is made available as a Web server at http://rotprof.lncc.br.

Genetic modifiers of EGFR dependence in non-small cell lung cancer.

Lung adenocarcinomas harboring activating mutations in the epidermal growth factor receptor (EGFR) represent a common molecular subset of non-small cell lung cancer (NSCLC) cases. EGFR mutations predict sensitivity to EGFR tyrosine kinase inhibitors (TKIs) and thus represent a dependency in NSCLCs harboring these alterations, but the genetic basis of EGFR dependence is not fully understood. Here, we applied an unbiased, ORF-based screen to identify genetic modifiers of EGFR dependence in EGFR-mutant NSCLC cells. This approach identified 18 kinase and kinase-related genes whose overexpression can substitute for EGFR in EGFR-dependent PC9 cells, and these genes include seven of nine Src family kinase genes, FGFR1, FGFR2, ITK, NTRK1, NTRK2, MOS, MST1R, and RAF1. A subset of these genes can complement loss of EGFR activity across multiple EGFR-dependent models. Unbiased gene-expression profiling of cells overexpressing EGFR bypass genes, together with targeted validation studies, reveals EGFR-independent activation of the MEK-ERK and phosphoinositide 3-kinase (PI3K)-AKT pathways. Combined inhibition of PI3K-mTOR and MEK restores EGFR dependence in cells expressing each of the 18 EGFR bypass genes. Together, these data uncover a broad spectrum of kinases capable of overcoming dependence on EGFR and underscore their convergence on the PI3K-AKT and MEK-ERK signaling axes in sustaining EGFR-independent survival.

Investigation of the structural preference and flexibility of the loop residues in amyloid fibrils of the HET-s prion.

The structural variability of a 16-residue loop (residues 246-261) which is in part disordered and connects two layers of the ß-solenoid formed by the prion-form of HET-s and its prion domain HET-s(218-289) is investigated using molecular dynamics computer simulation. A system of three HET-s(218-289) molecules in a ß-sheet structure as in the fibril is simulated in aqueous solution. The trajectory structures appear to be consistent with the Cα chemical shift data obtained. In order to delineate the influence of the ß-sheet core of the fibril upon the structural variability of the loop, the latter is also simulated without the ß-sheet core, but with its N- and C-terminal residues restrained at their positions in the fibril. The analysis of the trajectories shows that the structural variability of the loop is restricted by the ß-sheet core, least at its N-terminal end and most in the middle of the trimer.

Deriving Structural Information from Experimentally Measured Data on Biomolecules.

During the past half century, the number and accuracy of experimental techniques that can deliver values of observables for biomolecular systems have been steadily increasing. The conversion of a measured value Qexp of an observable quantity Q into structural information is, however, a task beset with theoretical and practical problems: 1) insufficient or inaccurate values of Qexp , 2) inaccuracies in the function Q(r→) used to relate the quantity Q to structure r→ , 3) how to account for the averaging inherent in the measurement of Qexp , 4) how to handle the possible multiple-valuedness of the inverse r→(Q) of the function Q(r→) , to mention a few. These apply to a variety of observable quantities Q and measurement techniques such as X-ray and neutron diffraction, small-angle and wide-angle X-ray scattering, free-electron laser imaging, cryo-electron microscopy, nuclear magnetic resonance, electron paramagnetic resonance, infrared and Raman spectroscopy, circular dichroism, Förster resonance energy transfer, atomic force microscopy and ion-mobility mass spectrometry. The process of deriving structural information from measured data is reviewed with an eye to non-experts and newcomers in the field using examples from the literature of the effect of the various choices and approximations involved in the process. A list of choices to be avoided is provided.

Assessing protein conformational sampling and structural stability via de novo design and molecular dynamics simulations.

Molecular dynamics and de novo techniques, associated to quality parameter sets, have excelled at determining the structure of small proteins with high accuracy. To achieve a detailed description of protein conformations, these methods must critically assess the thermodynamic features of the molecular ensembles. Here, a comparison of the conformational ensemble generated by molecular dynamics and de novo techniques were carried out for six Top7-based proteins carrying gp41 HIV-1 epitopes. The native Top7, a highly stable computationally designed protein, was used as benchmark. Structural stability, flexibility, and secondary structure content were assessed. The consistency of the latter was compared to experimental circular dichroism spectra for all proteins. While both methods are capable to identify the stable from unstable chimeric proteins, the sampled conformational space and flexibility differ significantly in both methods. Molecular dynamics simulations seem to better describe secondary structure content and identify regions responsible for conformational instability. The de novo method, as implemented in Rosetta-a prime tool for protein design, overestimates secondary structure content. On the other hand, its empirical energy function is capable to predict the threshold for protein stability.

Alpha-hemolysin nanopore allows discrimination of the microcystins variants.

Microcystins (MCs) are a class of cyclic heptapeptides with more than 100 variants produced by cyanobacteria present in surface waters. MCs are potent hepatotoxic agents responsible for fatal poisoning in animals and humans. Several techniques are employed in the detection of MCs, however, there is a shortage of methods capable of discriminating variants of MCs. In this work we demonstrate that the α-hemolysin (αHL) nanopore can detect and discriminate the variants (LR, YR and RR) of MCs in aqueous solution. The discrimination process is based on the analysis of the residence times of each variant of MCs within the unitary nanopore, as well as, on the amplitudes of the blockages in the ionic current flowing through it. Simulations of molecular dynamics and calculation of the electrostatic potential revealed that the variants of MCs present different charge distribution and correlated with the three patterns on the amplitudes of the blockages in the ionic current. Additionally, molecular docking analysis indicates different patterns of interaction of the variants of MCs with two specific regions of the nanopore. We conclude that αHL nanopore can discriminate variants of microcystins by a mechanism based mainly on electrostatic interaction. Finally, we propose the use of nanopore-based technology as a promising method for analyzing microcystins in aqueous solutions.

The patterns of base sequences in the nucleic acids of prokaryotes and eukaryotes reflect features of their abiotic past.

The base sequences of the nucleic acids corresponding to ten proteins (aconitase, alcohol dehydrogenase, enolase, fumarase, isocitrate dehydrogenase, lactate dehydrogenase, phosphofructokinase, phosphoglycerate mutase, pyruvate kinase and succinate dehydrogenase) belonging to a total of 154 species, ranging from prokaryotes to vertebrates, were compared with the base sequences of oligoribotides whose growth rates were calculated by a chemical kinetics model. It was shown that oligoribotides grown according to the kinetics model have a fraction of repetitive bases larger than expected from random processes. The base sequences of nucleic acids of prokaryotes and eukaryotes retain, in decreasing proportions, this feature of their abiotic past. Chemically synthesized pentameric stretches with repetitive bases are slightly more abundant than those present in prokaryotes. Genetic drift and natural selection, operating as fundamental laws even for the most primitive living systems, reduced the original, chemically controlled, repetitive base frequency in prokaryotes, which was further reduced for eukaryotes.

Aromatic Rings Commonly Used in Medicinal Chemistry: Force Fields Comparison and Interactions With Water Toward the Design of New Chemical Entities.

The identification of lead compounds usually includes a step of chemical diversity generation. Its rationale may be supported by both qualitative (SAR) and quantitative (QSAR) approaches, offering models of the putative ligand-receptor interactions. In both scenarios, our understanding of which interactions functional groups can perform is mostly based on their chemical nature (such as electronegativity, volume, melting point, lipophilicity etc.) instead of their dynamics in aqueous, biological solutions (solvent accessibility, lifetime of hydrogen bonds, solvent structure etc.). As a consequence, it is challenging to predict from 2D structures which functional groups will be able to perform interactions with the target receptor, at which intensity and relative abundance in the biological environment, all of which will contribute to ligand potency and intrinsic activity. With this in mind, the aim of this work is to assess properties of aromatic rings, commonly used for drug design, in aqueous solution through molecular dynamics simulations in order to characterize their chemical features and infer their impact in complexation dynamics. For this, common aromatic and heteroaromatic rings were selected and received new atomic charge set based on the direction and module of the dipole moment from MP2/6-31G* calculations, while other topological terms were taken from GROMOS53A6 force field. Afterwards, liquid physicochemical properties were simulated for a calibration set composed by nearly 40 molecules and compared to their respective experimental data, in order to validate each topology. Based on the reliance of the employed strategy, we expanded the dataset to more than 100 aromatic rings. Properties in aqueous solution such as solvent accessible surface area, H-bonds availability, H-bonds residence time, and water structure around heteroatoms were calculated for each ring, creating a database of potential interactions, shedding light on features of drugs in biological solutions, on the structural basis for bioisosterism and on the enthalpic/entropic costs for ligand-receptor complexation dynamics.

Morphological and functional renovascular changes as cause of resistant arterial hypertension - case report and literature review.

Resistant hypertension is defined by the inability to maintain within normal limits the blood pressure values of an individual, while he is under treatment with maximal tolerated doses of three antihypertensive agents. One of the most common types of resistant hypertension is renovascular hypertension (RVH), which is caused by the narrowing of the renal arteries, in the context of existing atherosclerotic plaques at that level. We are presenting the case of a hypertensive 56-year-old man admitted in the Clinic of Cardiology for a sudden rise of his blood pressure values, despite undergoing the scheduled treatment. The abdominal bruit discovered at the clinical examination and the hypokalemia, together with the mild impairment of the renal function raised the suspicion of an existing stenosis of the main renal blood vessels. Simple grey scale kidney ultrasound, Doppler ultrasound of the renal arteries, abdominal computed tomography and magnetic resonance angiography of the renal arteries, along with invasive renal angiography demonstrated a smaller right kidney, adrenal incidentalomas, reduced vascular diameter of renal arteries due to atheromatous lesions, thrombosis of the infrarenal segment of the abdominal aorta, and reduced vascular hemodynamics in the same territories. After the renal arteries revascularization and with minimal antihypertensive treatment, the patient had a favorable outcome, with normalization of blood pressure and renal function. Atherosclerotic disease causing renal artery stenosis is essential to be taken into consideration in the etiopathogenesis of resistant hypertension especially because RVH is a potentially curable disease.

PITOMBA: Parameter Interface for Oligosaccharide Molecules Based on Atoms.

A novel four-bead coarse-grained (CG) model for carbohydrates denoted PITOMBA was devised using a bottom-up approach based on the atomistic GROMOS 53A6GLYC force field and on experimental thermodynamical data. The model was developed to be used in conjunction with the SPC CG water model (J. Chem. Phys. 2011, 134, 084110) and the GROMOS force field functional form. Explicit electrostatic interactions are considered by assigning point charges to each CG bead. Validation of the model is presented to a variety of structural and thermodynamic properties for mono- and oligosaccharides in solution. In addition, the model development philosophy allows for prompt extensions to include hexopyranose chains with diverse glycosidic linkages and branches.

MDWiZ: a platform for the automated translation of molecular dynamics simulations.

A variety of popular molecular dynamics (MD) simulation packages were independently developed in the last decades to reach diverse scientific goals. However, such non-coordinated development of software, force fields, and analysis tools for molecular simulations gave rise to an array of software formats and arbitrary conventions for routine preparation and analysis of simulation input and output data. Different formats and/or parameter definitions are used at each stage of the modeling process despite largely contain redundant information between alternative software tools. Such Babel of languages that cannot be easily and univocally translated one into another poses one of the major technical obstacles to the preparation, translation, and comparison of molecular simulation data that users face on a daily basis. Here, we present the MDWiZ platform, a freely accessed online portal designed to aid the fast and reliable preparation and conversion of file formats that allows researchers to reproduce or generate data from MD simulations using different setups, including force fields and models with different underlying potential forms. The general structure of MDWiZ is presented, the features of version 1.0 are detailed, and an extensive validation based on GROMACS to LAMMPS conversion is presented. We believe that MDWiZ will be largely useful to the molecular dynamics community. Such fast format and force field exchange for a given system allows tailoring the chosen system to a given computer platform and/or taking advantage of a specific capabilities offered by different software engines.

The Effect of Temperature, Cations, and Number of Acyl Chains on the Lamellar to Non-Lamellar Transition in Lipid-A Membranes: A Microscopic View.

Lipopolysaccharides (LPS) are the main constituent of the outer bacterial membrane of Gram-negative bacteria. Lipid-A is the structural region of LPS that interacts with the innate immune system and induces inflammatory responses. It is formed by a phosphorylated ß-d-glucosaminyl-(1→6)-α-N-glucosamine disaccharide backbone containing ester-linked and amide-linked long-chain fatty acids, which may vary in length and number depending on the bacterial strains and the environment. Phenotypical variation (i.e., number of acyl chains), cation type, and temperature influence the phase transition, aggregate structure, and endotoxic activity of Lipid-A. We have applied an extension of the GROMOS force field 45a4 carbohydrate parameter set to investigate the behavior of hexa- and pentaacylated Lipid-A of Pseudomonas aeruginosa at two temperatures (300 and 328 K) and in the presence of mono- and divalent cations (represented by Ca(2+) and Na(+), respectively) through molecular dynamics simulations. The distinct phase of Lipid-A aggregates was characterized by structural properties, deuterium order parameters, the molecular shape of the lipid units (conical versus cylindrical), and molecular packing. Our results show that Na(+) ions induce a transition from the lamellar to nonlamellar phase. In contrast, the bilayer integrity is maintained in the presence of Ca(2+) ions. Through these findings, we present microscopic insights on the influence of different cations on the molecular behavior of Lipid-A associated with the lamellar to nonlamellar transition.

GROMOS 53A6GLYC, an Improved GROMOS Force Field for Hexopyranose-Based Carbohydrates.

An improved parameter set for explicit-solvent simulations of carbohydrates (referred to as GROMOS 53A6GLYC) is presented, allowing proper description of the most stable conformation of all 16 possible aldohexopyranose-based monosaccharides. This set includes refinement of torsional potential parameters associated with the determination of hexopyranose rings conformation by fitting to their corresponding quantum-mechanical profiles. Other parameters, as the rules for third and excluded neighbors, are taken directly from the GROMOS 53A6 force field. Comparisons of the herein presented parameter set to our previous version (GROMOS 45A4), the GLYCAM06 force field, and available NMR data are presented in terms of ring puckering free energies, conformational distribution of the hydroxymethyl group, and glycosidic linkage geometries for 16 selected monosaccharides and eight disaccharides. The proposed parameter modifications have shown a significant improvement for the above-mentioned quantities over the two tested force fields, while retaining full compatibility with the GROMOS 53A6 and 54A7 parameter sets for other classes of biomolecules.