Conformational dynamics of capping protein and interaction partners: Simulation studies

Structural dynamics and allostery of Rab proteins: strategies for drug discovery and design

Rab proteins represent the largest family of the Rab superfamily guanosine triphosphatase (GTPase). Aberrant human Rab proteins are associated with multiple diseases, including cancers and neurological disorders. Rab subfamily members display subtle conformational variations that render specificity in their physiological functions and can be targeted for subfamily-specific drug design. However, drug discovery efforts have not focused much on targeting Rab allosteric non-nucleotide binding sites which are subjected to less evolutionary pressures to be conserved, hence are likely to offer subfamily specificity and may be less prone to undesirable off-target interactions and side effects. To discover druggable allosteric binding sites, Rab structural dynamics need to be first incorporated using multiple experimentally and computationally obtained structures. The high-dimensional structural data may necessitate feature extraction methods to identify manageable representative structures for subsequent analyses. We have detailed state-of-the-art computational methods to (i) identify binding sites using data on sequence, shape, energy, etc., (ii) determine the allosteric nature of these binding sites based on structural ensembles, residue networks and correlated motions and (iii) identify small molecule binders through structure- and ligand-based virtual screening. To benefit future studies for targeting Rab allosteric sites, we herein detail a refined workflow comprising multiple available computational methods, which have been successfully used alone or in combinations. This workflow is also applicable for drug discovery efforts targeting other medically important proteins. Depending on the structural dynamics of proteins of interest, researchers can select suitable strategies for allosteric drug discovery and design, from the resources of computational methods and tools enlisted in the workflow.

Allosteric Coupling of CARMIL and V-1 Binding to Capping Protein Revealed by Hydrogen-Deuterium Exchange

Actin assembly is important for cell motility. The ability of actin subunits to join or leave filaments via the barbed end is critical to actin dynamics. Capping protein (CP) binds to barbed ends to prevent subunit gain and loss and is regulated by proteins that include V-1 and CARMIL. V-1 inhibits CP by sterically blocking one binding site for actin. CARMILs bind at a distal site and decrease the affinity of CP for actin, suggested to be caused by conformational changes. We used hydrogen-deuterium exchange with mass spectrometry (HDX-MS) to probe changes in structural dynamics induced by V-1 and CARMIL binding to CP. V-1 and CARMIL induce changes in both proteins' binding sites on the surface of CP, along with a set of internal residues. Both also affect the conformation of CP's ββ subunit "tentacle," a second distal actin-binding site. Concerted regulation of actin assembly by CP occurs through allosteric couplings between CP modulator and actin binding sites.

Allosteric binding sites in Rab11 for potential drug candidates

Rab11 is an important protein subfamily in the RabGTPase family. These proteins physiologically function as key regulators of intracellular membrane trafficking processes. Pathologically, Rab11 proteins are implicated in many diseases including cancers, neurodegenerative diseases and type 2 diabetes. Although they are medically important, no previous study has found Rab11 allosteric binding sites where potential drug candidates can bind to. In this study, by employing multiple clustering approaches integrating principal component analysis, independent component analysis and locally linear embedding, we performed structural analyses of Rab11 and identified eight representative structures. Using these representatives to perform binding site mapping and virtual screening, we identified two novel binding sites in Rab11 and small molecules that can preferentially bind to different conformations of these sites with high affinities. After identifying the binding sites and the residue interaction networks in the representatives, we computationally showed that these binding sites may allosterically regulate Rab11, as these sites communicate with switch 2 region that binds to GTP/GDP. These two allosteric binding sites in Rab11 are also similar to two allosteric pockets in Ras that we discovered previously.

Structural analysis of protein tyrosine phosphatase 1B reveals potentially druggable allosteric binding sites

Catalytic proteins such as human protein tyrosine phosphatase 1B (PTP1B), with conserved and highly polar active sites, warrant the discovery of druggable nonactive sites, such as allosteric sites, and potentially, therapeutic small molecules that can bind to these sites. Catalyzing the dephosphorylation of numerous substrates, PTP1B is physiologically important in intracellular signal transduction pathways in diverse cell types and tissues. Aberrant PTP1B is associated with obesity, diabetes, cancers, and neurodegenerative disorders. Utilizing clustering methods (based on root mean square deviation, principal component analysis, nonnegative matrix factorization, and independent component analysis), we have examined multiple PTP1B structures. Using the resulting representative structures in different conformational states, we determined consensus clustroids and used them to identify both known and novel binding sites, some of which are potentially allosteric. We report several lead compounds that could potentially bind to the novel PTP1B binding sites and can be further optimized. Considering the possibility for drug repurposing, we discovered homologous binding sites in other proteins, with ligands that could potentially bind to the novel PTP1B binding sites.

Capping protein regulatory cycle driven by CARMIL and V-1 may promote actin network assembly at protruding edges

Although capping protein (CP) terminates actin filament elongation, it promotes Arp2/3-dependent actin network assembly and accelerates actin-based motility both in vitro and in vivo. In vitro, capping protein Arp2/3 myosin I linker (CARMIL) antagonizes CP by reducing its affinity for the barbed end and by uncapping CP-capped filaments, whereas the protein V-1/myotrophin sequesters CP in an inactive complex. Previous work showed that CARMIL can readily retrieve CP from the CP:V-1 complex, thereby converting inactive CP into a version with moderate affinity for the barbed end. Here we further clarify the mechanism of this exchange reaction, and we demonstrate that the CP:CARMIL complex created by complex exchange slows the rate of barbed-end elongation by rapidly associating with, and dissociating from, the barbed end. Importantly, the cellular concentrations of V-1 and CP determined here argue that most CP is sequestered by V-1 at steady state in vivo. Finally, we show that CARMIL is recruited to the plasma membrane and only at cell edges undergoing active protrusion. Assuming that CARMIL is active only at this location, our data argue that a large pool of freely diffusing, inactive CP (CP:V-1) feeds, via CARMIL-driven complex exchange, the formation of weak-capping complexes (CP:CARMIL) at the plasma membrane of protruding edges. In vivo, therefore, CARMIL should promote Arp2/3-dependent actin network assembly at the leading edge by promoting barbed-end capping there.

Mapping the structural and dynamical features of multiple p53 DNA binding domains: insights into loop 1 intrinsic dynamics

The transcription factor p53 regulates cellular integrity in response to stress. p53 is mutated in more than half of cancerous cells, with a majority of the mutations localized to the DNA binding domain (DBD). In order to map the structural and dynamical features of the DBD, we carried out multiple copy molecular dynamics simulations (totaling 0.8 μs). Simulations show the loop 1 to be the most dynamic element among the DNA-contacting loops (loops 1-3). Loop 1 occupies two major conformational states: extended and recessed; the former but not the latter displays correlations in atomic fluctuations with those of loop 2 (~24 Å apart). Since loop 1 binds to the major groove whereas loop 2 binds to the minor groove of DNA, our results begin to provide some insight into the possible mechanism underpinning the cooperative nature of DBD binding to DNA. We propose (1) a novel mechanism underlying the dynamics of loop 1 and the possible tread-milling of p53 on DNA and (2) possible mutations on loop 1 residues to restore the transcriptional activity of an oncogenic mutation at a distant site.

Molecular dynamics simulations of a new branched antimicrobial peptide: a comparison of force fields

Branched antimicrobial peptides are promising as a new class of antibiotics displaying high activity and low toxicity and appear to work through a unique mechanism of action. We explore the structural dynamics of a covalently branched 18 amino acid peptide (referred to as B2088) in aqueous and membrane mimicking environments through molecular dynamics (MD) simulations. Towards this, we carry out conventional MD simulations and supplement these with replica exchange simulations. The simulations are carried out using four different force fields that are commonly employed for simulating biomolecular systems. These force fields are GROMOS53a6, CHARMM27 with cMAP, CHARMM27 without cMAP and AMBER99sb. The force fields are benchmarked against experimental data available from circular dichroism and nuclear magnetic resonance spectroscopies, and show that CHARMM27 without cMAP correction is the most successful in reproducing the structural dynamics of B2088 both in water and in the presence of micelles. Although the four force fields predict different structures of B2088, they all show that B2088 stabilizes against the head group of the lipid through hydrogen bonding of its Lys and Arg side chains. This leads us to hypothesize that B2088 is unlikely to penetrate into the hydrophobic region of the membrane owing to the high free energy costs of transfer from water, and possibly acts by carpeting and thus disrupting the membrane.