Molecular dynamics investigation of cysteine mutations: Effects on calcium ion affinity and structural stability in the RET cysteine-rich domain.

Continuous Molecular Evolution of Protein-Domain Structures by Single Amino Acid Changes

Advances in molecular dynamics simulations and free-energy calculations relevant for drug design

Serotonin transporter (SERT) regulates serotonergic signals by reuptaking serotonin from the synaptic clefts back into the presynaptic neurons. The recent resolution of the serotonin-SERT complex in multiple conformational states outlined the complete serotonin import cycle. However, a detailed functional appreciation of SERT also involves deciphering the coupling between global structural changes in the transport cycle to the bound chemicals to be transported. By employing molecular dynamics (MD) simulations and free energy calculations in different ligand binding states, here, we reveal how serotonin binding to SERT initiates the global conformational changes essential for serotonin import. Only when serotonin is bound to the central binding site, wedged between transmembrane helices (TMs) 3 and 8, can the system form an interaction network that bridges the two helical domains of the protein, thereby promoting the closure of an extracellular hydrophobic gate and sealing the bound serotonin. To test the role of this hydrophobic gate closure, we designed a series of nonequilibrium MD simulations to steer the outward-facing ↔ occluded transition with different gating configurations. The difference in nonequilibrium work required to fuel the transition indicates that the transition is more likely to happen when the extracellular gate is closed. The transition is not promoted when the gate is open or when 5-HT moves away from TM3 and TM8 toward an alternate pose. Such a local-global coupling is likely shared by other monoamine transporters considering the conservation of all involved structural elements.

RET is a tyrosine kinase receptor involved in numerous cellular mechanisms including proliferation, neuronal navigation, migration, and differentiation upon binding with glial cell derived neurotrophic factor family ligands. RET is an atypical tyrosine kinase receptor containing four cadherin domains in its extracellular part. Furthermore, it has been shown to act as a dependence receptor. Such a receptor is active in the absence of ligand, triggering apoptosis through a mechanism that requires receptor intracellular caspase cleavage. However, different data suggest that RET is not always associated with the cell death/survival balance but rather provides positional information. We demonstrate here that caspase cleavage of RET is involved in the regulation of adhesion in sympathetic neurons. The cleavage of RET generates an N-terminal truncated fragment that functions as a cadherin accessory protein, modifying cadherin environment and potentiating cadherin-mediated cell aggregation. Thus, the caspase cleavage of RET generates two RET fragments: one intracellular domain that can trigger cell death in apoptotic permissive settings, and one membrane-anchored ectodomain with cadherin accessory activity. We propose that this latter function may notably be important for the adequate development of the superior cervical ganglion.

RET Recognition of GDNF-GFRα1 Ligand by a Composite Binding Site Promotes Membrane-Proximal Self-Association

Signaling through the receptor tyrosine kinase RET is essential during normal development. Both gain- and loss-of-function mutations are involved in a variety of diseases, yet the molecular details of receptor activation have remained elusive. We have reconstituted the complete extracellular region of the RET signaling complex together with Neurturin (NRTN) and GFRα2 and determined its structure at 5.7-Å resolution by cryo-EM. The proteins form an assembly through RET-GFRα2 and RET-NRTN interfaces. Two key interaction points required for RET extracellular domain binding were observed: (i) the calcium-binding site in RET that contacts GFRα2 domain 3 and (ii) the RET cysteine-rich domain interaction with NRTN. The structure highlights the importance of the RET cysteine-rich domain and allows proposition of a model to explain how complex formation leads to RET receptor dimerization and its activation. This provides a framework for targeting RET activity and for further exploration of mechanisms underlying neurological diseases.

Calcium-dependent lectin I from Pseudomonas aeruginosa (PA-IL) binds specifically to oligosaccharides presenting an alpha-galactose residue at their nonreducing end, such as the disaccharides alphaGal1-2betaGalOMe, alphaGal1-3betaGalOMe, and alphaGal1-4betaGalOMe. This provides a unique model for studying the effect of the glycosidic linkage of the ligands on structure and thermodynamics of the complexes by means of experimental and theoretical tools. The structural features of PA-IL in complex with the three disaccharides were established by docking and molecular dynamics simulations and compared with those observed in available crystal structures, including PA-IL.alphaGal1-2betaGalOMe complex, which was solved at 2.4 A resolution and reported herein. The role of a structural bridge water molecule in the binding site of PA-IL was also elucidated through molecular dynamics simulations and free energy calculations. This water molecule establishes three very stable hydrogen bonds with O6 of nonreducing galactose, oxygen from Pro-51 main chain, and nitrogen from Gln-53 main chain of the lectin binding site. Binding free energies for PA-IL in complex with the three disaccharides were investigated, and the results were compared with the experimental data determined by titration microcalorimetry. When the bridge water molecule was included in the free energy calculations, the simulations predicted the correct binding affinity trends with the 1-2-linked disaccharide presenting three times stronger affinity ligand than the other two. These results highlight the role of the water molecule in the binding site of PA-IL and indicate that it should be taken into account when designing glycoderivatives active against P. aeruginosa adhesion.

We have studied the binding of two organic cations, an iminium (IM) and a guanidinium (GU), to a cyclophane host P4--4Na+, using molecular dynamics simulations and free energy calculations. A proper treatment of the long-range electrostatic forces is essential for the stability of these highly charged complexes, and a simple cutoff at 12 A results in an artifactual dissociation of the IM−P4--4Na+ complex. Since the host is highly aromatic and the guests cationic, cation−π interactions play an important role in the complex stability. In free energy calculations, using a simple additive force field, we calculate that the relative free energy of association of IM and GU binding to the host is 2.3 kcal/mol favoring IM, which is of the correct sign but 1.4 kcal/mol too small in magnitude. Differences in van der Waals interaction energies are mainly responsible for the different binding strengths, and the host adopts different shapes when accommodating IM compared to GU. To approximately estimate the contributi...

4-(Phenylamino)-pyrrolo[2,1-f][1,2,4]triazines have been discovered as inhibitors of p38alpha. Experimental assays have proven that the configuration of alpha-Me-benzyl connected with amide at C6 is essential for the binding affinity. The S-configured inhibitor (11j) displays 80 times more potency than the R-configured one (11k). Here we investigated the mechanism how different configurations influence the binding affinity using molecular dynamics simulations, free energy calculations and free energy decomposition analysis. We found that the van der Waals interactions play the most important role in differentiating the activities between 11j and 11k with p38alpha. The difference of the van der Waals interactions is primarily determined by two residues, LEU108 and LEU167. Consequently stabilization of pyrrolo[2,1-f][1,2,4]triazine ring is important for the activities of inhibitors. Meanwhile we observed that the different configuration of the alpha-Me-benzyl group leads to the difference of binding between 11j and 11k. In conclusion, our work shows that it is feasible to analyze the chirality effect of inhibitors with different configurations by molecular dynamics simulations and free energy calculations, and provides useful information for drug design.

Free energy calculations of RNA interactions.

Cannulae are tubular protein filaments that accumulate on the extracellular surface of the hyperthermophilic archaeon Pyrodictium abyssi during cell division. Cannulae have been postulated to act as a primitive extracellular matrix through which cells could communicate or exchange material, although their native biological function remains obscure. Here, we report cryoEM structural analyses of ex vivo cannulae and of in vitro protein assemblies derived from recombinant cannula-like proteins. Three-dimensional reconstructions of P. abyssi cannulae revealed that the structural interactions between protomers in the native and recombinant filaments were based on donor strand complementation, a form of non-covalent polymerization in which a donor β-strand from one subunit is inserted into an acceptor groove in a β-sheet of a neighboring subunit. Donor strand complementation in cannulae is reinforced through calcium ion coordination at the interfaces between structural subunits in the respective assemblies. While donor strand complementation occurs during the assembly of chaperone-usher pili, this process requires the participation of accessory proteins that are localized in the outer membrane. In contrast, we demonstrate that calcium ions can induce assembly of cannulae in the absence of other co-factors. Crystallographic analysis of a recombinant cannula-like protein monomer provided evidence that calcium ion binding primes the precursor for donor strand invasion through unblocking of the acceptor groove. Bioinformatic analysis suggested that structurally homologous cannula-like proteins occurred within the genomes of other hyperthermophilic archaea and were encompassed within the TasA superfamily of biomatrix proteins. CryoEM structural analyses of tubular filaments derived from in vitro assembly of a recombinant cannula-like protein from an uncultured Hyperthermus species revealed a common mode of assembly to the Pyrodictium cannulae, in which donor strand complementation and calcium ion binding stabilized longitudinal and lateral assembly in tubular 2D sheets.

Protein engineering and design principles employing the 20 standard amino acids have been extensively used to achieve stable protein scaffolds and deliver their specific activities. Although this confers some advantages, it often restricts the sequence, chemical space, and ultimately the functional diversity of proteins. Moreover, although site-specific incorporation of non-natural amino acids (nnAAs) has been proven to be a valuable strategy in protein engineering and therapeutics development, its utility in the affinity-maturation of nanobodies is not fully explored. Besides, current experimental methods do not routinely employ nnAAs due to their enormous library size and infinite combinations. To address this, we have developed an integrated computational pipeline employing structure-based protein design methodologies, molecular dynamics simulations and free energy calculations, for the binding affinity prediction of an nnAA-incorporated nanobody toward its target and selection of potent binders. We show that by incorporating halogenated tyrosines, the affinity of 9G8 nanobody can be improved toward epidermal growth factor receptor (EGFR), a crucial cancer target. Surface plasmon resonance (SPR) assays showed that the binding of several 3-chloro-l-tyrosine (3MY)-incorporated nanobodies were improved up to 6-fold into a picomolar range, and the computationally estimated binding affinities shared a Pearson's r of 0.87 with SPR results. The improved affinity was found to be due to enhanced van der Waals interactions of key 3MY-proximate nanobody residues with EGFR, and an overall increase in the nanobody's structural stability. In conclusion, we show that our method can facilitate screening large libraries and predict potent site-specific nnAA-incorporated nanobody binders against crucial disease-targets.

G-quadruplexes are higher-order DNA structures formed from guanine-rich sequences, and have been identified as attractive anticancer drug targets. Elucidating the three-dimensional structure of G-quadruplex with 9-amino acridines and the specific interactions involved in binding selectivity are the key to understanding their mechanism of action. Fluorescence titration assays, competitive dialysis and NMR studies have been used to study the binding specificity of 9-amino acridines to DNA. Structural models of the complexes with the telomeric DNA G-quadruplex based on NMR measurements were developed and further examined by molecular dynamics simulations and free energy calculations. Selective binding of 9-amino acridines for G-quadruplex sequences were observed. These compounds bind between A and G-tetrads, involving significant π-π interactions and several strong hydrogen bonds. The specific interactions between different moieties of the 9-amino acridines to the DNA were examined and shown to play a significant role in governing the overall stabilities of DNA G-quadruplex complexes. Both 9-amino acridines, with similar binding affinities to the G-quadruplex, were shown to induce different level of structural stabilization through intercalation. This unique property of altering structural stability is likely a contributing factor for affecting telomerase function and, subsequently, the observed differences in the anticancer activities between the two 9-amino acridines.

Molecular mechanism of Gd@C82(OH)22 increasing collagen expression: Implication for encaging tumor

ABTRACT The epidermal growth factor receptor (EGFR) kinase inhibitors Gefitinib, Erlotinib, Afatinib and Osimertinib have been approved for the treatments of non-small cell lung cancer patients harboring sensitive EGFR mutations, but resistance arises rapidly. To date all approved EGFR inhibitors are ATP-competitive inhibitors, highlighting the need for therapeutic agents with alternative mechanisms of action. Allosteric kinase inhibitors offer a promising new therapeutic strategy to ATP-competitive inhibitors. The mutant-selective allosteric EGFR inhibitors EAI045 exhibited higher potency for EGFRL858R&T790M compared to WT, which was also effective in EGFR-mutant models including those harboring the C797S mutation. However, it was not effective as a single-agent inhibitor, and require the co-administration of the anti-EGFR antibody Cetuximab. Further efforts produced a more potent analog JBJ-04-125-02, which can inhibit cell proliferation as a single-agent inhibitor. In the present study, molecular dynamics simulations and free energy calculations were performed and revealed the detailed inhibitory mechanism of JBJ-04-125-02 as more potent EGFR inhibitor. Moreover, the energy difference between HOMO and LUMO calculated by DFT implied the higher interaction of JBJ-04-125-02 than EAI045 in the active site of the EGFR. The identified key features obtained from the molecular modeling enabled us to design novel EGFR allosteric inhibitors. Communicated by Ramaswamy H. Sarma