Edmond Y. Lau

A bimodal distribution of two distinct categories of intrinsically-disordered structures with separate functions in FG nucleoporins

Nuclear pore complexes (NPCs) gate the only conduits for nucleocytoplasmic transport in eukaryotes. Their gate is formed by nucleoporins containing large intrinsically disordered domains with multiple phenylalanine-glycine repeats (FG domains). In combination, these are hypothesized to form a structurally and chemically homogeneous network of random coils at the NPC center, which sorts macromolecules by size and hydrophobicity. Instead, we found that FG domains are structurally and chemically heterogeneous. They adopt distinct categories of intrinsically disordered structures in non-random distributions. Some adopt globular, collapsed coil configurations and are characterized by a low charge content. Others are highly charged and adopt more dynamic, extended coil conformations. Interestingly, several FG nucleoporins feature both types of structures in a bimodal distribution along their polypeptide chain. This distribution functionally correlates with the attractive or repulsive character of their interactions with collapsed coil FG domains displaying cohesion toward one another and extended coil FG domains displaying repulsion. Topologically, these bipartite FG domains may resemble sticky molten globules connected to the tip of relaxed or extended coils. Within the NPC, the crowding of FG nucleoporins and the segregation of their disordered structures based on their topology, dimensions, and cohesive character could force the FG domains to form a tubular gate structure or transporter at the NPC center featuring two separate zones of traffic with distinct physicochemical properties.

A Method to Predict Blood-Brain Barrier Permeability of Drug-Like Compounds Using Molecular Dynamics Simulations

The blood-brain barrier (BBB) is formed by specialized tight junctions between endothelial cells that line brain capillaries to create a highly selective barrier between the brain and the rest of the body. A major problem to overcome in drug design is the ability of the compound in question to cross the BBB. Neuroactive drugs are required to cross the BBB to function. Conversely, drugs that target other parts of the body ideally should not cross the BBB to avoid possible psychotropic side effects. Thus, the task of predicting the BBB permeability of new compounds is of great importance. Two gold-standard experimental measures of BBB permeability are logBB (the concentration of drug in the brain divided by concentration in the blood) and logPS (permeability surface-area product). Both methods are time-consuming and expensive, and although logPS is considered the more informative measure, it is lower throughput and more resource intensive. With continual increases in computer power and improvements in molecular simulations, in silico methods may provide viable alternatives. Computational predictions of these two parameters for a sample of 12 small molecule compounds were performed. The potential of mean force for each compound through a 1,2-dioleoyl-sn-glycero-3-phosphocholine bilayer is determined by molecular dynamics simulations. This system setup is often used as a simple BBB mimetic. Additionally, one-dimensional position-dependent diffusion coefficients are calculated from the molecular dynamics trajectories. The diffusion coefficient is combined with the free energy landscape to calculate the effective permeability (Peff) for each sample compound. The relative values of these permeabilities are compared to experimentally determined logBB and logPS values. Our computational predictions correlate remarkably well with both logBB (R(2) = 0.94) and logPS (R(2) = 0.90). Thus, we have demonstrated that this approach may have the potential to provide reliable, quantitatively predictive BBB permeability, using a relatively quick, inexpensive method.

Toward a small molecule, biomimetic carbonic anhydrase model: theoretical and experimental investigations of a panel of zinc(II) aza-macrocyclic catalysts.

A panel of five zinc-chelated aza-macrocycle ligands and their ability to catalyze the hydration of carbon dioxide to bicarbonate, H(2)O + CO(2) → H(+) + HCO(3)(–), was investigated using quantum-mechanical methods and stopped-flow experiments. The key intermediates in the reaction coordinate were optimized using the M06-2X density functional with aug-cc-pVTZ basis set. Activation energies for the first step in the catalytic cycle, nucleophilic CO(2) addition, were calculated from gas-phase optimized transition-state geometries. The computationally derived trend in activation energies was found to not correspond with the experimentally observed rates. However, activation energies for the second, bicarbonate release step, which were estimated using calculated bond dissociation energies, provided good agreement with the observed trend in rate constants. Thus, the joint theoretical and experimental results provide evidence that bicarbonate release, not CO(2) addition, may be the rate-limiting step in CO(2) hydration by zinc complexes of aza-macrocyclic ligands. pH-independent rate constants were found to increase with decreasing Lewis acidity of the ligand-Zn complex, and the trend in rate constants was correlated with molecular properties of the ligands. It is suggested that tuning catalytic efficiency through the first coordination shell of Zn(2+) ligands is predominantly a balance between increasing charge-donating character of the ligand and maintaining the catalytically relevant pK(a) below the operating pH.

Intramolecular Cohesion of Coils Mediated by Phenylalanine–Glycine Motifs in the Natively Unfolded Domain of a Nucleoporin

The nuclear pore complex (NPC) provides the sole aqueous conduit for macromolecular exchange between the nucleus and the cytoplasm of cells. Its diffusion conduit contains a size-selective gate formed by a family of NPC proteins that feature large, natively unfolded domains with phenylalanine–glycine repeats (FG domains). These domains of nucleoporins play key roles in establishing the NPC permeability barrier, but little is known about their dynamic structure. Here we used molecular modeling and biophysical techniques to characterize the dynamic ensemble of structures of a representative FG domain from the yeast nucleoporin Nup116. The results showed that its FG motifs function as intramolecular cohesion elements that impart order to the FG domain and compact its ensemble of structures into native premolten globular configurations. At the NPC, the FG motifs of nucleoporins may exert this cohesive effect intermolecularly as well as intramolecularly to form a malleable yet cohesive quaternary structure composed of highly flexible polypeptide chains. Dynamic shifts in the equilibrium or competition between intra- and intermolecular FG motif interactions could facilitate the rapid and reversible structural transitions at the NPC conduit needed to accommodate passing karyopherin–cargo complexes of various shapes and sizes while simultaneously maintaining a size-selective gate against protein diffusion.

The Use of One-Bead One-Compound Combinatorial Library Technology to Discover High-Affinity αvβ3 Integrin and Cancer Targeting Arginine-Glycine-Aspartic Acid Ligands with a Built-in Handle

The αvβ3 integrin, expressed on the surface of various normal and cancer cells, is involved in numerous physiologic processes such as angiogenesis, apoptosis, and bone resorption. Because this integrin plays a key role in angiogenesis and metastasis of human tumors, αvβ3 integrin ligands are of great interest to advances in targeted therapy and cancer imaging. In this report, one-bead one-compound (OBOC) combinatorial libraries containing the arginine-glycine-aspartic acid (RGD) motif were designed and screened against K562 myeloid leukemia cells that had been transfected with the human αvβ3 integrin gene. Cyclic peptide LXW7 was identified as a leading ligand with a built-in handle that binds specifically to αvβ3 and showed comparable binding affinity (IC50 = 0.68 ± 0.08 μmol/L) to some of the well-known RGD “head-to-tail” cyclic pentapeptide ligands reported in the literature. The biotinylated form of LXW7 ligand showed similar binding strength as LXW7 against αvβ3 integrin, whereas biotinylated RGD cyclopentapeptide ligands revealed a 2- to 8-fold weaker binding affinity than their free forms. LXW7 was able to bind to both U-87MG glioblastoma and A375M melanoma cell lines, both of which express high levels of αvβ3 integrin. In vivo and ex vivo optical imaging studies with the biotinylated ligand/streptavidin-Cy5.5 complex in nude mice bearing U-87MG or A375M xenografts revealed preferential uptake of biotinylated LXW7 in tumor. When compared with biotinylated RGD cyclopentapeptide ligands, biotinylated LXW7 showed higher tumor uptake but lower liver uptake. Mol Cancer Ther; 9(10); 2714–23. ©2010 AACR.

Selectively Targeting T- and B-Cell Lymphomas: A Benzothiazole Antagonist of α4β1 Integrin

Current cancer chemotherapeutic agents clinically deployed today are designed to be indiscriminately cytotoxic, however, achieving selective targeting of cancer malignancies would allow for improved diagnostic and chemotherapeutic tools. Integrin alpha(4)beta(1), a heterodimeric cell surface receptor, is believed to have a low-affinity conformation in resting normal lymphocytes and an activated high-affinity conformation in cancerous cells, specifically T- and B-cell lymphomas. This highly attractive yet poorly understood receptor has been selectively targeted with the bisaryl urea peptidomimetic antagonist 1. However, concerns regarding its preliminary pharmacokinetic (PK) profile provided an impetus to change the pharmacophore from a bisaryl urea to a 2-arylaminobenzothiazole moiety, resulting in an analogue with improved physicochemical properties, solubility, and kidney:tumor ratio while maintaining potency (6; IC(50) = 53 pM). The results presented herein utilized heterocyclic and solid-phase chemistry, cell adhesion assay, and in vivo optical imaging using the cyanine dye Cy5.5 conjugate.

Effects of fluorine substitution on the structure and dynamics of complexes of dihydrofolate reductase (Escherichia coli).

Fluorine NMR experiments with a protein containing fluorinated amino acid analogs can often be used to probe structure and dynamics of the protein as well as conformational changes produced by binding of small molecules. The relevance of NMR experiments with fluorine-containing materials to characteristics of the corresponding native (nonfluorinated) proteins depends upon the extent to which these characteristics are altered by the presence of fluorine. The present work uses molecular dynamics simulations to explore the effects of replacement of tryptophan by 6-fluorotryptophan in folate and methotrexate complexes of the enzyme dihydrofolate reductase (DHFR) (Escherichia coli). Simulations of the folate-native enzyme complex produce local correlation times and order parameters that are generally in good agreement with experimental values. Simulations of the corresponding fluorotryptophan-containing system indicate that the structure and dynamics of this complex are scarcely changed by the presence of fluorinated amino acids. Calculations with the pharmacologically important methotrexate-enzyme complex predict dynamical behavior of the protein similar to that of the folate complex for both the fluorinated and native enzyme. It thus appears that, on the time scale sampled by these computer simulations, substitution of 6-fluorotryptophan for tryptophan has little effect on either the structures or dynamics of DHFR in these complexes.

Consequences of breaking the Asp-His hydrogen bond of the catalytic triad: effects on the structure and dynamics of the serine esterase cutinase.

The objective of this study has been to investigate the effects on the structure and dynamics that take place with the breaking of the Asp-His hydrogen bond in the catalytic triad Asp175-His188-Ser120 of the serine esterase cutinase in the ground state. Four molecular dynamics simulations were performed on this enzyme in solution. The starting structures in two simulations had the Asp175-His188 hydrogen bond intact, and in two simulations the Asp175-His188 hydrogen bond was broken. Conformations of the residues comprising the catalytic triad are well behaved during both simulations containing the intact Asp175-His188 hydrogen bond. Short contacts of less than 2.6 A were observed in 1.2% of the sampled distances between the carboxylate oxygens of Asp175 and the NE2 of His188. The simulations showed that the active site residues exhibit a great deal of mobility when the Asp175-His188 hydrogen bond is broken. In the two simulations in which the Asp175-His188 hydrogen bond is not present, the final geometries for the residues in the catalytic triad are not in catalytically productive conformations. In both simulations, Asp175 and His188 are more than 6 A apart in the final structure from dynamics, and the side chains of Ser120 and Asp175 are in closer proximity to the NE2 of His188 than to ND1. Nonlocal effects on the structure of cutinase were observed. A loop formed by residues 26-31, which is on the opposite end of the protein relative to the active site, was greatly affected. Further changes in the dynamics of cutinase were determined from quasiharmonic mode analysis. The frequency of the second lowest mode was greatly reduced when the Asp175-His188 hydrogen bond was broken, and several higher modes showed lower frequencies. All four simulations showed that the oxyanion hole, composed of residues Ser42 and Gln121, is stable. Only one of the hydrogen bonds (Ser42 OG to Gln121 NE2) observed in the crystal structure that stabilize the conformation of Ser42 OG persisted throughout the simulations. This hydrogen bond appears to be enough for the oxyanion hole to retain its structural integrity.

Identification of a Possible Secondary Picrotoxin-Binding Site on the GABAA Receptor

The type A GABA receptors (GABARs) are ligand-gated ion channels (LGICs) found in the brain and are the major inhibitory neurotransmitter receptors. Upon binding of an agonist, the GABAR opens and increases the intraneuronal concentration of chloride ions, thus hyperpolarizing the cell and inhibiting the transmission of the nerve action potential. GABARs also contain many other modulatory binding pockets that differ from the agonist-binding site. The composition of the GABAR subunits can alter the properties of these modulatory sites. Picrotoxin is a noncompetitive antagonist for LGICs, and by inhibiting GABAR, picrotoxin can cause overstimulation and induce convulsions. We use addition of picrotoxin to probe the characteristics and possible mechanism of an additional modulatory pocket located at the interface between the ligand-binding domain and the transmembrane domain of the GABAR. Picrotoxin is widely regarded as a pore-blocking agent that acts at the cytoplasmic end of the channel. However, there are also data to suggest that there may be an additional, secondary binding site for picrotoxin. Through homology modeling, molecular docking, and molecular dynamics simulations, we show that binding of picrotoxin to this interface pocket correlates with these data, and negative modulation occurs at the pocket via a kinking of the pore-lining helices into a more closed orientation.

A theoretical examination of the factors controlling the catalytic efficiency of the DNA-(adenine-N6)-methyltransferase from Thermus aquaticus.

Ab initio and density functional calculations have been carried out to more fully understand the factors controlling the catalytic activity of the Thermus aquaticus DNA methyltransferase (MTaqI) in the N-methylation at the N6 of an adenine nucleobase. The noncatalyzed reaction was modeled as a methyl transfer from trimethylsulfonium to the N6 of adenine. Activation barriers of 32.0 kcal/mol and 24.0 kcal/mol were predicted for the noncatalyzed reaction in the gas phase by MP2/6–31+G(d,p)//HF/6–31+G(d,p) and B3LYP/6–31+G(d,p) calculations, respectively. Calculations performed to evaluate the effect of substrate positioning in the active site of MTaqI on the reaction determine the barrier to be 23.4 kcal/mol and 17.3 kcal/mol for the MP2/6–31+G(d,p)//HF/6–31+G(d,p) and B3LYP/6–31+G(d,p) gas phase calculations, respectively. The effect of hydrogen bonding between the N6 of adenine and the terminal oxygen of Asn-105 on the activation barrier was also studied. A formamide molecule was modeled into the system to mimic the function of active site residue Asn-105. The activation barrier for this reaction was found to be 21.8 kcal/mol and 15.8 kcal/mol as determined from the MP2/6–31+G(d,p)//HF/6–31+G(d,p) and B3LYP/6–31+G(d,p) calculations, respectively. This result predicts a contribution of less than 2 kcal/mol to the lowering of the activation barrier from amide hydrogen bonding between formamide and N6 of adenine. Comparison of the reaction coordinates suggest that it is not the hydrogen bonding of the Asn-105 that lends to the catalytic prowess of the enzyme since the organization of the substrates in the active site of the enzyme has a far greater effect on reducing the activation barrier. The results also suggest a stepwise mechanism for the removal of the hydrogen from the N6 of adenine as opposed to a concerted reaction in which a proton is abstracted simultaneously with the transfer of the methyl group. The hydrogen on the N6 of the intermediate methyl adenine product is far more acidic than in the reactant complex and may be subsequently abstracted by basic groups in the active site that are too weak to abstract the proton before the full sp3 hybridization of the attacking nitrogen.