Estimation of ligand efficacies of metabotropic glutamate receptors from conformational forces obtained from molecular dynamics simulations

Molecular regulation of calcium-sensing receptor (CaSR)-mediated signaling

Calcium-sensing receptor (CaSR), a family C G-protein-coupled receptor, plays a crucial role in regulating calcium homeostasis by sensing small concentration changes of extracellular Ca2+, Mg2+, amino acids (e.g., L-Trp and L-Phe), small peptides, anions (e.g., HCO3 - and PO4 3-), and pH. CaSR-mediated intracellular Ca2+ signaling regulates a diverse set of cellular processes including gene transcription, cell proliferation, differentiation, apoptosis, muscle contraction, and neuronal transmission. Dysfunction of CaSR with mutations results in diseases such as autosomal dominant hypocalcemia, familial hypocalciuric hypercalcemia, and neonatal severe hyperparathyroidism. CaSR also influences calciotropic disorders, such as osteoporosis, and noncalciotropic disorders, such as cancer, Alzheimer's disease, and pulmonary arterial hypertension. This study first reviews recent advances in biochemical and structural determination of the framework of CaSR and its interaction sites with natural ligands, as well as exogenous positive allosteric modulators and negative allosteric modulators. The establishment of the first CaSR protein-protein interactome network revealed 94 novel players involved in protein processing in endoplasmic reticulum, trafficking, cell surface expression, endocytosis, degradation, and signaling pathways. The roles of these proteins in Ca2+-dependent cellular physiological processes and in CaSR-dependent cellular signaling provide new insights into the molecular basis of diseases caused by CaSR mutations and dysregulated CaSR activity caused by its protein interactors and facilitate the design of therapeutic agents that target CaSR and other family C G-protein-coupled receptors.

DIRECT-ID: An automated method to identify and quantify conformational variations--application to β2 -adrenergic GPCR

The conformational dynamics of a macromolecule can be modulated by a number of factors, including changes in environment, ligand binding, and interactions with other macromolecules, among others. We present a method that quantifies the differences in macromolecular conformational dynamics and automatically extracts the structural features responsible for these changes. Given a set of molecular dynamics (MD) simulations of a macromolecule, the norms of the differences in covariance matrices are calculated for each pair of trajectories. A matrix of these norms thus quantifies the differences in conformational dynamics across the set of simulations. For each pair of trajectories, covariance difference matrices are parsed to extract structural elements that undergo changes in conformational properties. As a demonstration of its applicability to biomacromolecular systems, the method, referred to as DIRECT-ID, was used to identify relevant ligand-modulated structural variations in the β2 -adrenergic (β2 AR) G-protein coupled receptor. Micro-second MD simulations of the β2 AR in an explicit lipid bilayer were run in the apo state and complexed with the ligands: BI-167107 (agonist), epinephrine (agonist), salbutamol (long-acting partial agonist), or carazolol (inverse agonist). Each ligand modulated the conformational dynamics of β2 AR differently and DIRECT-ID analysis of the inverse-agonist vs. agonist-modulated β2 AR identified residues known through previous studies to selectively propagate deactivation/activation information, along with some previously unidentified ligand-specific microswitches across the GPCR. This study demonstrates the utility of DIRECT-ID to rapidly extract functionally relevant conformational dynamics information from extended MD simulations of large and complex macromolecular systems.

Cyclopropyl-containing positive allosteric modulators of metabotropic glutamate receptor subtype 5

Positive allosteric modulators (PAMs) binding to the transmembrane (TM) domain of metabotropic glutamate receptor 5 (mGluR5) are promising therapeutic agents for psychiatric disorders and traumatic brain injury (TBI). Novel PAMs based on a trans-2-phenylcyclopropane amide scaffold have been designed and synthesized. Facilitating ligand design and allowing estimation of binding affinities to the mGluR5 TM domain was the novel computational strategy, site identification by ligand competitive saturation (SILCS). The potential protective activity of the new compounds was evaluated using nitric oxide (NO) production in BV2 microglial cell cultures treated with lipopolysaccharide (LPS), and the toxicity of the new compounds tested using a cell viability assay. One of the new compounds, 3a, indicated promising activity with potency of 30 μM, which is 4.5-fold more potent than its lead compound 3,3'-difluorobenzaldazine (DFB), and showed no detectable toxicity with concentrations as high as 1000 μM. Thus this compound represents a new lead for possible development as treatment for TBI and related neurodegenerative disorders.

A consistent force field parameter set for zwitterionic amino acid residues

Isolated amino acids play an important role in biochemistry and are therefore an interesting object of study. Atomistic molecular dynamics (MD) simulations can provide a high-resolution picture of the dynamic features of these species, especially in their biological environment. Unfortunately, most standard force field packages lack libraries for isolated amino acids in their zwitterionic form. Although several studies have used ad-hoc parameterizations for single amino acids, a consistent force-field parameter set for these molecules is still missing. Here, we present such a parameter library derived from the widely used parm99SB set from the AMBER program package. The parameter derivation for all 20 proteinogenic amino acids transparently followed established procedures with histidine treated in three different protonation states. All amino acids were subjected to MD simulations in four different forms for comparison: zwitterionic, N-teminally capped with acetyl, C-terminally capped with N-methyl, and capped at both termini. Simulation results show similarities between the different forms. Five zwitterionic amino acids-arginine, glutamate, glycine, phenylalanine, leucine-were simulated in a protein environment. Proteins and ligands generally retained their initial structure. The new parameter set will thus facilitate future atomistic simulations of these species.