Implicit model to capture electrostatic features of membrane environment.

Membrane protein structure prediction and design are challenging due to the complexity of capturing the interactions in the lipid layer, such as those arising from electrostatics. Accurately capturing electrostatic energies in the low-dielectric membrane often requires expensive Poisson-Boltzmann calculations that are not scalable for membrane protein structure prediction and design. In this work, we have developed a fast-to-compute implicit energy function that considers the realistic characteristics of different lipid bilayers, making design calculations tractable. This method captures the impact of the lipid head group using a mean-field-based approach and uses a depth-dependent dielectric constant to characterize the membrane environment. This energy function Franklin2023 (F23) is built upon Franklin2019 (F19), which is based on experimentally derived hydrophobicity scales in the membrane bilayer. We evaluated the performance of F23 on five different tests probing (1) protein orientation in the bilayer, (2) stability, and (3) sequence recovery. Relative to F19, F23 has improved the calculation of the tilt angle of membrane proteins for 90% of WALP peptides, 15% of TM-peptides, and 25% of the adsorbed peptides. The performances for stability and design tests were equivalent for F19 and F23. The speed and calibration of the implicit model will help F23 access biophysical phenomena at long time and length scales and accelerate the membrane protein design pipeline.

Fate and Efficacy of Engineered Allogeneic Stem Cells Targeting Cell Death and Proliferation Pathways in Primary and Brain Metastatic Lung Cancer.

Primary and metastatic lung cancer is a leading cause of cancer-related death and novel therapies are urgently needed. Epidermal growth factor receptor (EGFR) and death receptor (DR) 4/5 are both highly expressed in primary and metastatic non-small cell lung cancer (NSCLC); however, targeting these receptors individually has demonstrated limited therapeutic benefit in patients. In this study, we created and characterized diagnostic and therapeutic stem cells (SC), expressing EGFR-targeted nanobody (EV) fused to the extracellular domain of death DR4/5 ligand (DRL) (EVDRL) that simultaneously targets EGFR and DR4/5, in primary and metastatic NSCLC tumor models. We show that EVDRL targets both cell surface receptors, and induces caspase-mediated apoptosis in a broad spectrum of NSCLC cell lines. Utilizing real-time dual imaging and correlative immunohistochemistry, we show that allogeneic SCs home to tumors and when engineered to express EVDRL, alleviate tumor burden and significantly increase survival in primary and brain metastatic NSCLC. This study reports mechanistic insights into simultaneous targeting of EGFR- and DR4/5 in lung tumors and presents a promising approach for translation into the clinical setting.

Hallucinating structure-conditioned antibody libraries for target-specific binders.

Antibodies are widely developed and used as therapeutics to treat cancer, infectious disease, and inflammation. During development, initial leads routinely undergo additional engineering to increase their target affinity. Experimental methods for affinity maturation are expensive, laborious, and time-consuming and rarely allow the efficient exploration of the relevant design space. Deep learning (DL) models are transforming the field of protein engineering and design. While several DL-based protein design methods have shown promise, the antibody design problem is distinct, and specialized models for antibody design are desirable. Inspired by hallucination frameworks that leverage accurate structure prediction DL models, we propose the FvHallucinator for designing antibody sequences, especially the CDR loops, conditioned on an antibody structure. Such a strategy generates targeted CDR libraries that retain the conformation of the binder and thereby the mode of binding to the epitope on the antigen. On a benchmark set of 60 antibodies, FvHallucinator generates sequences resembling natural CDRs and recapitulates perplexity of canonical CDR clusters. Furthermore, the FvHallucinator designs amino acid substitutions at the VH-VL interface that are enriched in human antibody repertoires and therapeutic antibodies. We propose a pipeline that screens FvHallucinator designs to obtain a library enriched in binders for an antigen of interest. We apply this pipeline to the CDR H3 of the Trastuzumab-HER2 complex to generate in silico designs predicted to improve upon the binding affinity and interfacial properties of the original antibody. Thus, the FvHallucinator pipeline enables generation of inexpensive, diverse, and targeted antibody libraries enriched in binders for antibody affinity maturation.

Affinity maturation of TCR-like antibodies using phage display guided by structural modeling.

TCR-like antibodies represent a unique type of engineered antibodies with specificity toward pHLA, a ligand normally restricted to the sensitive recognition by T cells. Here, we report a phage display-based sequential development path of such antibodies. The strategy goes from initial lead identification through in silico informed CDR engineering in combination with framework engineering for affinity and thermostability optimization, respectively. The strategy allowed the identification of HLA-DQ2.5 gluten peptide-specific TCR-like antibodies with low picomolar affinity. Our method outlines an efficient and general method for development of this promising class of antibodies, which should facilitate their utility including translation to human therapy.

The structural basis of PTEN regulation by multi-site phosphorylation.

Phosphatase and tensin homolog (PTEN) is a phosphatidylinositol-3,4,5-triphosphate (PIP3) phospholipid phosphatase that is commonly mutated or silenced in cancer. PTEN's catalytic activity, cellular membrane localization and stability are orchestrated by a cluster of C-terminal phosphorylation (phospho-C-tail) events on Ser380, Thr382, Thr383 and Ser385, but the molecular details of this multi-faceted regulation have remained uncertain. Here we use a combination of protein semisynthesis, biochemical analysis, NMR, X-ray crystallography and computational simulations on human PTEN and its sea squirt homolog, VSP, to obtain a detailed picture of how the phospho-C-tail forms a belt around the C2 and phosphatase domains of PTEN. We also visualize a previously proposed dynamic N-terminal α-helix and show that it is key for PTEN catalysis but disordered upon phospho-C-tail interaction. This structural model provides a comprehensive framework for how C-tail phosphorylation can impact PTEN's cellular functions.

Shotgun scanning glycomutagenesis: A simple and efficient strategy for constructing and characterizing neoglycoproteins.

As a common protein modification, asparagine-linked (N-linked) glycosylation has the capacity to greatly influence the biological and biophysical properties of proteins. However, the routine use of glycosylation as a strategy for engineering proteins with advantageous properties is limited by our inability to construct and screen large collections of glycoproteins for cataloguing the consequences of glycan installation. To address this challenge, we describe a combinatorial strategy termed shotgun scanning glycomutagenesis in which DNA libraries encoding all possible glycosylation site variants of a given protein are constructed and subsequently expressed in glycosylation-competent bacteria, thereby enabling rapid determination of glycosylatable sites in the protein. The resulting neoglycoproteins can be readily subjected to available high-throughput assays, making it possible to systematically investigate the structural and functional consequences of glycan conjugation along a protein backbone. The utility of this approach was demonstrated with three different acceptor proteins, namely bacterial immunity protein Im7, bovine pancreatic ribonuclease A, and human anti-HER2 single-chain Fv antibody, all of which were found to tolerate N-glycan attachment at a large number of positions and with relatively high efficiency. The stability and activity of many glycovariants was measurably altered by N-linked glycans in a manner that critically depended on the precise location of the modification. Structural models suggested that affinity was improved by creating novel interfacial contacts with a glycan at the periphery of a protein-protein interface. Importantly, we anticipate that our glycomutagenesis workflow should provide access to unexplored regions of glycoprotein structural space and to custom-made neoglycoproteins with desirable properties.

Colicin-Mediated Transport of DNA through the Iron Transporter FepA.

Colicins are protein antibiotics deployed by Escherichia coli to eliminate competing strains. Colicins frequently exploit outer membrane (OM) nutrient transporters to penetrate the selectively permeable bacterial cell envelope. Here, by applying live-cell fluorescence imaging, we were able to monitor the entry of the pore-forming toxin colicin B (ColB) into E. coli and localize it within the periplasm. We further demonstrate that single-stranded DNA coupled to ColB can also be transported to the periplasm, emphasizing that the import routes of colicins can be exploited to carry large cargo molecules into bacteria. Moreover, we characterize the molecular mechanism of ColB association with its OM receptor FepA by applying a combination of photoactivated cross-linking, mass spectrometry, and structural modeling. We demonstrate that complex formation is coincident with large-scale conformational changes in the colicin. Thereafter, active transport of ColB through FepA involves the colicin taking the place of the N-terminal half of the plug domain that normally occludes this iron transporter. IMPORTANCE Decades of excessive use of readily available antibiotics has generated a global problem of antibiotic resistance and, hence, an urgent need for novel antibiotic solutions. Bacteriocins are protein-based antibiotics produced by bacteria to eliminate closely related competing bacterial strains. Bacteriocin toxins have evolved to bypass the complex cell envelope in order to kill bacterial cells. Here, we uncover the cellular penetration mechanism of a well-known but poorly understood bacteriocin called colicin B that is active against Escherichia coli. Moreover, we demonstrate that the colicin B-import pathway can be exploited to deliver conjugated DNA cargo into bacterial cells. Our work leads to a better understanding of the way bacteriocins, as potential alternative antibiotics, execute their mode of action as well as highlighting how they might even be exploited in the genomic manipulation of Gram-negative bacteria.

A high-affinity human TCR-like antibody detects celiac disease gluten peptide-MHC complexes and inhibits T cell activation.

Antibodies specific for peptides bound to human leukocyte antigen (HLA) molecules are valuable tools for studies of antigen presentation and may have therapeutic potential. Here, we generated human T cell receptor (TCR)-like antibodies toward the immunodominant signature gluten epitope DQ2.5-glia-α2 in celiac disease (CeD). Phage display selection combined with secondary targeted engineering was used to obtain highly specific antibodies with picomolar affinity. The crystal structure of a Fab fragment of the lead antibody 3.C11 in complex with HLA-DQ2.5:DQ2.5-glia-α2 revealed a binding geometry and interaction mode highly similar to prototypic TCRs specific for the same complex. Assessment of CeD biopsy material confirmed disease specificity and reinforced the notion that abundant plasma cells present antigen in the inflamed CeD gut. Furthermore, 3.C11 specifically inhibited activation and proliferation of gluten-specific CD4+ T cells in vitro and in HLA-DQ2.5 humanized mice, suggesting a potential for targeted intervention without compromising systemic immunity.

Structural basis for peptide substrate specificities of glycosyltransferase GalNAc-T2.

The polypeptide N-acetylgalactosaminyl transferase (GalNAc-T) enzyme family initiates O-linked mucin-type glycosylation. The family constitutes 20 isoenzymes in humans. GalNAc-Ts exhibit both redundancy and finely tuned specificity for a wide range of peptide substrates. In this work, we deciphered the sequence and structural motifs that determine the peptide substrate preferences for the GalNAc-T2 isoform. Our approach involved sampling and characterization of peptide-enzyme conformations obtained from Rosetta Monte Carlo-minimization-based flexible docking. We computationally scanned 19 amino acid residues at positions -1 and +1 of an eight-residue peptide substrate, which comprised a dataset of 361 (19x19) peptides with previously characterized experimental GalNAc-T2 glycosylation efficiencies. The calculations recapitulated experimental specificity data, successfully discriminating between glycosylatable and non-glycosylatable peptides with a probability of 96.5% (ROC-AUC score), a balanced accuracy of 85.5% and a false positive rate of 7.3%. The glycosylatable peptide substrates viz. peptides with proline, serine, threonine, and alanine at the -1 position of the peptide preferentially exhibited cognate sequon-like conformations. The preference for specific residues at the -1 position of the peptide was regulated by enzyme residues R362, K363, Q364, H365 and W331, which modulate the pocket size and specific enzyme-peptide interactions. For the +1 position of the peptide, enzyme residues K281 and K363 formed gating interactions with aromatics and glutamines at the +1 position of the peptide, leading to modes of peptide-binding sub-optimal for catalysis. Overall, our work revealed enzyme features that lead to the finely tuned specificity observed for a broad range of peptide substrates for the GalNAc-T2 enzyme. We anticipate that the key sequence and structural motifs can be extended to analyze specificities of other isoforms of the GalNAc-T family and can be used to guide design of variants with tailored specificity.

Novel sampling strategies and a coarse-grained score function for docking homomers, flexible heteromers, and oligosaccharides using Rosetta in CAPRI rounds 37-45.

Critical Assessment of PRediction of Interactions (CAPRI) rounds 37 through 45 introduced larger complexes, new macromolecules, and multistage assemblies. For these rounds, we used and expanded docking methods in Rosetta to model 23 target complexes. We successfully predicted 14 target complexes and recognized and refined near-native models generated by other groups for two further targets. Notably, for targets T110 and T136, we achieved the closest prediction of any CAPRI participant. We created several innovative approaches during these rounds. Since round 39 (target 122), we have used the new RosettaDock 4.0, which has a revamped coarse-grained energy function and the ability to perform conformer selection during docking with hundreds of pregenerated protein backbones. Ten of the complexes had some degree of symmetry in their interactions, so we tested Rosetta SymDock, realized its shortcomings, and developed the next-generation symmetric docking protocol, SymDock2, which includes docking of multiple backbones and induced-fit refinement. Since the last CAPRI assessment, we also developed methods for modeling and designing carbohydrates in Rosetta, and we used them to successfully model oligosaccharide-protein complexes in round 41. Although the results were broadly encouraging, they also highlighted the pressing need to invest in (a) flexible docking algorithms with the ability to model loop and linker motions and in (b) new sampling and scoring methods for oligosaccharide-protein interactions.

The structure of the colorectal cancer-associated enzyme GalNAc-T12 reveals how nonconserved residues dictate its function.

Polypeptide N-acetylgalactosaminyl transferases (GalNAc-Ts) initiate mucin type O-glycosylation by catalyzing the transfer of N-acetylgalactosamine (GalNAc) to Ser or Thr on a protein substrate. Inactive and partially active variants of the isoenzyme GalNAc-T12 are present in subsets of patients with colorectal cancer, and several of these variants alter nonconserved residues with unknown functions. While previous biochemical studies have demonstrated that GalNAc-T12 selects for peptide and glycopeptide substrates through unique interactions with its catalytic and lectin domains, the molecular basis for this distinct substrate selectivity remains elusive. Here we examine the molecular basis of the activity and substrate selectivity of GalNAc-T12. The X-ray crystal structure of GalNAc-T12 in complex with a di-glycosylated peptide substrate reveals how a nonconserved GalNAc binding pocket in the GalNAc-T12 catalytic domain dictates its unique substrate selectivity. In addition, the structure provides insight into how colorectal cancer mutations disrupt the activity of GalNAc-T12 and illustrates how the rules dictating GalNAc-T12 function are distinct from those for other GalNAc-Ts.

Targeting the CoREST complex with dual histone deacetylase and demethylase inhibitors.

Here we report corin, a synthetic hybrid agent derived from the class I HDAC inhibitor (entinostat) and an LSD1 inhibitor (tranylcypromine analog). Enzymologic analysis reveals that corin potently targets the CoREST complex and shows more sustained inhibition of CoREST complex HDAC activity compared with entinostat. Cell-based experiments demonstrate that corin exhibits a superior anti-proliferative profile against several melanoma lines and cutaneous squamous cell carcinoma lines compared to its parent monofunctional inhibitors but is less toxic to melanocytes and keratinocytes. CoREST knockdown, gene expression, and ChIP studies suggest that corin's favorable pharmacologic effects may rely on an intact CoREST complex. Corin was also effective in slowing tumor growth in a melanoma mouse xenograft model. These studies highlight the promise of a new class of two-pronged hybrid agents that may show preferential targeting of particular epigenetic regulatory complexes and offer unique therapeutic opportunities.

Web-accessible molecular modeling with Rosetta: The Rosetta Online Server that Includes Everyone (ROSIE).

The Rosetta molecular modeling software package provides a large number of experimentally validated tools for modeling and designing proteins, nucleic acids, and other biopolymers, with new protocols being added continually. While freely available to academic users, external usage is limited by the need for expertise in the Unix command line environment. To make Rosetta protocols available to a wider audience, we previously created a web server called Rosetta Online Server that Includes Everyone (ROSIE), which provides a common environment for hosting web-accessible Rosetta protocols. Here we describe a simplification of the ROSIE protocol specification format, one that permits easier implementation of Rosetta protocols. Whereas the previous format required creating multiple separate files in different locations, the new format allows specification of the protocol in a single file. This new, simplified protocol specification has more than doubled the number of Rosetta protocols available under ROSIE. These new applications include pKa determination, lipid accessibility calculation, ribonucleic acid redesign, protein-protein docking, protein-small molecule docking, symmetric docking, antibody docking, cyclic toxin docking, critical binding peptide determination, and mapping small molecule binding sites. ROSIE is freely available to academic users at http://rosie.rosettacommons.org.

Glycoengineering of Esterase Activity through Metabolic Flux-Based Modulation of Sialic Acid.

This report describes the metabolic glycoengineering (MGE) of intracellular esterase activity in human colon cancer (LS174T) and Chinese hamster ovary (CHO) cells. In silico analysis of carboxylesterases CES1 and CES2 suggested that these enzymes are modified with sialylated N-glycans, which are proposed to stabilize the active multimeric forms of these enzymes. This premise was supported by treating cells with butanolylated ManNAc to increase sialylation, which in turn increased esterase activity. By contrast, hexosamine analogues not targeted to sialic acid biosynthesis (e.g., butanoylated GlcNAc or GalNAc) had minimal impact. Measurement of mRNA and protein confirmed that esterase activity was controlled through glycosylation and not through transcription or translation. Azide-modified ManNAc analogues widely used in MGE also enhanced esterase activity and provided a way to enrich targeted glycoengineered proteins (such as CES2), thereby providing unambiguous evidence that the compounds were converted to sialosides and installed into the glycan structures of esterases as intended. Overall, this study provides a pioneering example of the modulation of intracellular enzyme activity through MGE, which expands the value of this technology from its current status as a labeling strategy and modulator of cell surface biological events.

Functional loss of semaphorin 3C and/or semaphorin 3D and their epistatic interaction with ret are critical to Hirschsprung disease liability.

Innervation of the gut is segmentally lost in Hirschsprung disease (HSCR), a consequence of cell-autonomous and non-autonomous defects in enteric neuronal cell differentiation, proliferation, migration, or survival. Rare, high-penetrance coding variants and common, low-penetrance non-coding variants in 13 genes are known to underlie HSCR risk, with the most frequent variants in the ret proto-oncogene (RET). We used a genome-wide association (220 trios) and replication (429 trios) study to reveal a second non-coding variant distal to RET and a non-coding allele on chromosome 7 within the class 3 Semaphorin gene cluster. Analysis in Ret wild-type and Ret-null mice demonstrates specific expression of Sema3a, Sema3c, and Sema3d in the enteric nervous system (ENS). In zebrafish embryos, sema3 knockdowns show reduction of migratory ENS precursors with complete ablation under conjoint ret loss of function. Seven candidate receptors of Sema3 proteins are also expressed within the mouse ENS and their expression is also lost in the ENS of Ret-null embryos. Sequencing of SEMA3A, SEMA3C, and SEMA3D in 254 HSCR-affected subjects followed by in silico protein structure modeling and functional analyses identified five disease-associated alleles with loss-of-function defects in semaphorin dimerization and binding to their cognate neuropilin and plexin receptors. Thus, semaphorin 3C/3D signaling is an evolutionarily conserved regulator of ENS development whose dys-regulation is a cause of enteric aganglionosis.

Adding diverse noncanonical backbones to rosetta: enabling peptidomimetic design.

Peptidomimetics are classes of molecules that mimic structural and functional attributes of polypeptides. Peptidomimetic oligomers can frequently be synthesized using efficient solid phase synthesis procedures similar to peptide synthesis. Conformationally ordered peptidomimetic oligomers are finding broad applications for molecular recognition and for inhibiting protein-protein interactions. One critical limitation is the limited set of design tools for identifying oligomer sequences that can adopt desired conformations. Here, we present expansions to the ROSETTA platform that enable structure prediction and design of five non-peptidic oligomer scaffolds (noncanonical backbones), oligooxopiperazines, oligo-peptoids, [Formula: see text]-peptides, hydrogen bond surrogate helices and oligosaccharides. This work is complementary to prior additions to model noncanonical protein side chains in ROSETTA. The main purpose of our manuscript is to give a detailed description to current and future developers of how each of these noncanonical backbones was implemented. Furthermore, we provide a general outline for implementation of new backbone types not discussed here. To illustrate the utility of this approach, we describe the first tests of the ROSETTA molecular mechanics energy function in the context of oligooxopiperazines, using quantum mechanical calculations as comparison points, scanning through backbone and side chain torsion angles for a model peptidomimetic. Finally, as an example of a novel design application, we describe the automated design of an oligooxopiperazine that inhibits the p53-MDM2 protein-protein interaction. For the general biological and bioengineering community, several noncanonical backbones have been incorporated into web applications that allow users to freely and rapidly test the presented protocols (http://rosie.rosettacommons.org). This work helps address the peptidomimetic community's need for an automated and expandable modeling tool for noncanonical backbones.

Computed structures of point deletion mutants and their enzymatic activities.

Point deletions in enzymes can vary in effect from negligible to complete loss of activity; however, these effects are not generally predictable. Deletions are widely observed in nature and often result in diseases such as cancer, cystic fibrosis, or osteogenesis imperfecta. Here, we have developed an algorithm to model the perturbed structures of deletion mutants with the ultimate goal of predicting their activities. The algorithm works by deleting the specified residue from the wild-type structure, creating a gap that is closed using a combination of local and global moves that change the backbone torsion angles of the protein structure. On a set of five proteins for which both wild-type and deletion mutant x-ray crystal structures are available, the algorithm produces deep, narrow energy funnels within 1.5 Å of the crystal structure for the deletion mutants. To assess the ability of our algorithm to predict activity from the predicted structures, we tested the correlation of experimental activity with several measures of the predicted structure ensemble using a set of 45 point deletions from ricin. Estimates incorporating likely prevalence of active and inactive deletion sites suggest that activity can be predicted correctly over 60% of the time from the active site root-mean squared deviation of the lowest energy predicted structures. The predictions are stronger than simple sequence organization measures, but more fundamental work is required in structure prediction and enzyme activity determination to allow consistent prediction of activity.

Prediction of calcite morphology from computational and experimental studies of mutations of a de novo-designed peptide.

Many organisms use macromolecules, often proteins or peptides, to control the growth of inorganic crystals into complex materials. The ability to model peptide-mineral interactions accurately could allow for the design of novel peptides to produce materials with desired properties. Here, we tested a computational algorithm developed to predict the structure of peptides on mineral surfaces. Using this algorithm, we analyzed energetic and structural differences between a 16-residue peptide (bap4) designed to interact with a calcite growth plane and single- and double-point mutations of the charged residues. Currently, no experimental method is available to resolve the structures of proteins on solid surfaces, which precludes benchmarking for computational models. Therefore, to test the models, we chemically synthesized each peptide and analyzed its effects on calcite crystal growth. Whereas bap4 affected the crystal growth by producing heavily stepped corners and edges, point mutants had variable influences on morphology. Calculated residue-specific binding energies correlated with experimental observations; point mutations of residues predicted to be crucial to surface interactions produced morphologies most similar to unmodified calcite. These results suggest that peptide conformation plays a role in mineral interactions and that the computational model supplies valid energetic and structural data that can provide information about expected crystal morphology.

Structurally distinct toxicity inhibitors bind at common loci on ß-amyloid fibril.

The accumulation of aggregated ß-Amyloid (Aß) in the brain is a hallmark of Alzheimer's disease and is thought to play a role in the neurotoxicity associated with the disease. The mechanism by which Aß aggregates induce toxicity is uncertain. Nonetheless, several small molecules have been found to interact with Aß fibrils and to prevent their toxicity. In this paper we studied the binding of these known toxicity inhibitors to Aß fibrils, as a means to explore surfaces or loci on Aß aggregates that may be significant in the mechanism of action of these inhibitors. We believe knowledge of these binding loci will provide insight into surfaces on the Aß fibrils important in Aß biological activity. The program DOCK was used to computationally dock the inhibitors to an Aß fibril. The inhibitors docked at two shared binding loci, near Lys28 and at the C-termini near Asn27 and Val39. The docking predictions were experimentally verified using lysine specific chemical modifications and Aß fibrils mutated at Asn27. We found that both Congo red and Myricetin, despite being structurally different, bound at the same two sites. Additionally, our data suggests that three additional Aß toxicity inhibitors may also bind in one of the sites. Identification of these common binding loci provides targets on the Aß fibril surface that can be tested in the future for their role in Aß biological activity.

De novo design of peptide-calcite biomineralization systems.

Many organisms produce complex, hierarchically structured, inorganic materials via protein-influenced crystal growth--a process known as biomineralization. Understanding this process would shed light on hard-tissue formation and guide efforts to develop biomaterials. We created and tested a computational method to design protein-biomineralization systems. The algorithm folds a protein from a fully extended structure and simultaneously optimizes the fold, orientation, and sequence of the protein adsorbed to a crystal surface. We used the algorithm to design peptides (16 residues) to modify calcite (CaCO(3)) crystallization. We chemically synthesized six peptides that were predicted to bind different states of a calcite growth plane. All six peptides dramatically affected calcite crystal growth (as observed by scanning electron microscopy), and the effects were dependent on the targeted state of the {001} growth plane. Additionally, we synthesized and assayed scrambled variants of all six designed peptides to distinguish cases where sequence composition determines the interactions versus cases where sequence order (and presumably structure) plays a role. Scrambled variants of negatively charged peptides also had dramatic effects on calcite crystallization; in contrast, scrambled variants of positively charged peptides had a variable effect on crystallization, ranging from dramatic to mild. Special emphasis is often placed on acidic protein residues in calcified tissue mineralization; the work presented here suggests an important role for basic residues as well. In particular, this work implicates a potential role for basic residues in sequence-order specificity for peptide-mineral interactions.