Thermal isomerization rates in retinal analogues using Ab-Initio molecular dynamics

For a detailed understanding of chemical processes in nature and industry, we need accurate models of chemical reactions in complex environments. While Eyring transition state theory is commonly used for modeling chemical reactions, it is most accurate for small molecules in the gas phase. A wide range of alternative rate theories exist that can better capture reactions involving complex molecules and environmental effects. However, they require that the chemical reaction is sampled by molecular dynamics simulations. This is a formidable challenge since the accessible simulation timescales are many orders of magnitude smaller than typical timescales of chemical reactions. To overcome these limitations, rare event methods involving enhanced molecular dynamics sampling are employed. In this work, thermal isomerization of retinal is studied using tight-binding density functional theory. Results from transition state theory are compared to those obtained from enhanced sampling. Rates obtained from dynamical reweighting using infrequent metadynamics simulations were in close agreement with those from transition state theory. Meanwhile, rates obtained from application of Kramers' rate equation to a sampled free energy profile along a torsional dihedral reaction coordinate were found to be up to three orders of magnitude higher. This discrepancy raises concerns about applying rate methods to one-dimensional reaction coordinates in chemical reactions.

Impact of glycan nature on structure and viscoelastic properties of glycopeptide hydrogels

Mucus is a complex biological hydrogel that acts as a barrier for almost everything entering or exiting the body. It is therefore of emerging interest for biomedical and pharmaceutical applications. Besides water, the most abundant components are the large and densely glycosylated mucins, glycoproteins of up to 20 MDa and carbohydrate content of up to 80 wt%. Here, we designed and explored a library of glycosylated peptides to deconstruct the complexity of mucus. Using the well-characterized hFF03 coiled-coil system as a hydrogel-forming peptide scaffold, we systematically probed the contribution of single glycans to the secondary structure as well as the formation and viscoelastic properties of the resulting hydrogels. We show that glycan-decoration does not affect α-helix and coiled-coil formation while it alters gel stiffness. By using oscillatory macrorheology, dynamic light scattering microrheology, and fluorescence lifetime-based nanorheology, we characterized the glycopeptide materials over several length scales. Molecular simulations revealed that the glycosylated linker may extend into the solvent, but more frequently interacts with the peptide, thereby likely modifying the stability of the self-assembled fibers. This systematic study highlights the interplay between glycan structure and hydrogel properties and may guide the development of synthetic mucus mimetics.

How Chromophore Labels Shape the Structure and Dynamics of a Peptide Hydrogel

Biocompatible and functionalizable hydrogels have a wide range of (potential) medicinal applications. The hydrogelation process, particularly for systems with very low polymer weight percentages (<1 wt %), remains poorly understood, making it challenging to predict the self-assembly of a given molecular building block into a hydrogel. This severely hinders the rational design of self-assembled hydrogels. In this study, we demonstrate the impact of an N-terminal group on the self-assembly and rheology of the peptide hydrogel hFF03 (hydrogelating, fibril forming peptide 03) using molecular dynamics simulations, oscillatory shear rheology, and circular dichroism spectroscopy. We find that the chromophore and even its specific regioisomers have a significant influence on the microscopic structure and dynamics of the self-assembled fibril, and on the macroscopic mechanical properties. This is because the chromophore influences the possible salt bridges, which form and stabilize the fibril formation. Furthermore, we find that the solvation shell fibrils by itself cannot explain the viscoelasticity of hFF03 hydrogels. Our atomistic model of the hFF03 fibril formation enables a more rational design of these hydrogels. In particular, altering the N-terminal chromophore emerges as a design strategy to tune the mechanic properties of these self-assembled peptide hydrogels.

Biosynthetic incorporation of fluorinated amino acids into the nonribosomal peptide gramicidin S

Fluorine is a key element in medicinal chemistry, as it can significantly enhance the pharmacological properties of drugs. In this study, we aimed to biosynthetically produce fluorinated analogues of the antimicrobial cyclic decapeptide gramicidin S (GS). However, our results show that the A-domain of the NRPS module GrsA rejects 4-fluorinated analogues of its native substrate Phe due to an interrupted T-shaped aromatic interaction in the binding pocket. We demonstrate that GrsA mutant W239S improves the incorporation of 4-fluorinated Phe into GS both in vitro and in vivo. Our findings provide new insights into the behavior of NRPSs towards fluorinated amino acids and strategies for the engineered biosynthesis of fluorinated peptides.

Optimizing molecular potential models by imposing kinetic constraints with path reweighting

Empirical force fields employed in molecular dynamics simulations of complex systems are often optimized to reproduce experimentally determined structural and thermodynamic properties. In contrast, experimental knowledge about the interconversion rates between metastable states in such systems is hardly ever incorporated in a force field due to a lack of an efficient approach. Here, we introduce such a framework based on the relationship between dynamical observables, such as rate constants, and the underlying molecular model parameters using the statistical mechanics of trajectories. Given a prior ensemble of molecular dynamics trajectories produced with imperfect force field parameters, the approach allows for the optimal adaption of these parameters such that the imposed constraint of equally predicted and experimental rate constant is obeyed. To do so, the method combines the continuum path ensemble maximum caliber approach with path reweighting methods for stochastic dynamics. When multiple solutions are found, the method selects automatically the combination that corresponds to the smallest perturbation of the entire path ensemble, as required by the maximum entropy principle. To show the validity of the approach, we illustrate the method on simple test systems undergoing rare event dynamics. Next to simple 2D potentials, we explore particle models representing molecular isomerization reactions and protein-ligand unbinding. Besides optimal interaction parameters, the methodology gives physical insights into what parts of the model are most sensitive to the kinetics. We discuss the generality and broad implications of the methodology.

Grid-based state space exploration for molecular binding

Binding processes are difficult to sample with molecular-dynamics (MD) simulations. In particular, the state space exploration is often incomplete. Evaluating the molecular interaction energy on a grid circumvents this problem but is heavily limited by state space dimensionality. Here, we make the first steps towards a low-dimensional grid-based model of molecular binding. We discretise the state space of relative positions and orientations of the two molecules under the rigid body assumption. The corresponding program is published as the Python package molgri. For the rotational component of the grids, we test algorithms based on Euler angles, polyhedra and quaternions, of which the polyhedra-based are the most uniform. The program outputs a sequence of molecular structures that can be easily processed by standard MD programs to calculate grid point energies. We demonstrate the grid-based approach on two molecular systems: a water dimer and a coiled-coil protein interacting with a chloride anion. For the second system we relax the rigid-body assumption and improve the accuracy of the grid point energies by an energy minimisation. In both cases, oriented bonding patterns and energies confirm expectations from chemical intuition and MD simulations. We also demonstrate how analysis of energy contributions on a grid can be performed and demonstrate that electrostatically-driven association is sufficiently resolved by point-energy calculations. Overall, grid-based models of molecular binding are potentially a powerful complement to molecular sampling approaches, and we see the potential to expand the method to quantum chemistry and flexible docking applications.

CommonNNClustering─A Python Package for Generic Common-Nearest-Neighbor Clustering

Density-based clustering procedures are widely used in a variety of data science applications. Their advantage lies in the capability to find arbitrarily shaped and sized clusters and robustness against outliers. In particular, they proved effective in the analysis of molecular dynamics simulations, where they serve to identify relevant, low-energetic molecular conformations. As such, they can provide a convenient basis for the construction of kinetic (core-set) Markov-state models. Here we present the open-source Python project CommonNNClustering, which provides an easy-to-use and efficient reimplementation of the common-nearest-neighbor (CommonNN) method. The package provides functionalities for hierarchical clustering and an evaluation of the results. We put our emphasis on a generic API design to keep the implementation flexible and open for customization.

Target Recognition in Tandem WW Domains: Complex Structures for Parallel and Antiparallel Ligand Orientation in h-FBP21 Tandem WW

Protein-protein interactions often rely on specialized recognition domains, such as WW domains, which bind to specific proline-rich sequences. The specificity of these protein-protein interactions can be increased by tandem repeats, i.e., two WW domains connected by a linker. With a flexible linker, the WW domains can move freely with respect to each other. Additionally, the tandem WW domains can bind in two different orientations to their target sequences. This makes the elucidation of complex structures of tandem WW domains extremely challenging. Here, we identify and characterize two complex structures of the tandem WW domain of human formin-binding protein 21 and a peptide sequence from its natural binding partner, the core-splicing protein SmB/B'. The two structures differ in the ligand orientation and, consequently, also in the relative orientation of the two WW domains. We analyze and probe the interactions in the complexes by molecular simulations and NMR experiments. The workflow to identify the complex structures uses molecular simulations, density-based clustering, and peptide docking. It is designed to systematically generate possible complex structures for repeats of recognition domains. These structures will help us to understand the synergistic and multivalency effects that generate the astonishing versatility and specificity of protein-protein interactions.

Water Network in the Binding Pocket of Fluorinated BPTI-Trypsin Complexes─Insights from Simulation and Experiment

Structural waters in the S1 binding pocket of β-trypsin are critical for the stabilization of the complex of β-trypsin with its inhibitor bovine pancreatic trypsin inhibitor (BPTI). The inhibitor strength of BPTI can be modulated by replacing the critical lysine residue at the P1 position by non-natural amino acids. We study BPTI variants in which the critical Lys15 in BPTI has been replaced by α-aminobutyric acid (Abu) and its fluorinated derivatives monofluoroethylglycine (MfeGly), difluoroethylglycine (DfeGly), and trifluoroethylglycine (TfeGly). We investigate the hypothesis that additional water molecules in the binding pocket can form specific noncovalent interactions with the fluorinated side chains and thereby act as an extension of the inhibitors. We report potentials of mean force (PMF) of the unbinding process for all four complexes and enzyme activity inhibition assays. Additionally, we report the protein crystal structure of the Lys15MfeGly-BPTI-β-trypsin complex (pdb: 7PH1). Both experimental and computational data show a stepwise increase in inhibitor strength with increasing fluorination of the Abu side chain. The PMF additionally shows a minimum for the encounter complex and an intermediate state just before the bound state. In the bound state, the computational analysis of the structure and dynamics of the water molecules in the S1 pocket shows a highly dynamic network of water molecules that does not indicate a rigidification or stabilizing trend in regard to energetic properties that could explain the increase in inhibitor strength. The analysis of the energy and the entropy of the water molecules in the S1 binding pocket using grid inhomogeneous solvation theory confirms this result. Overall, fluorination systematically changes the binding affinity, but the effect cannot be explained by a persistent water network in the binding pocket. Other effects, such as the hydrophobicity of fluorinated amino acids and the stability of the encounter complex as well as the additional minimum in the potential of mean force in the bound state, likely influence the affinity more directly.

The occurrence of ansamers in the synthesis of cyclic peptides

α-Amanitin is a bicyclic octapeptide composed of a macrolactam with a tryptathionine cross-link forming a handle. Previously, the occurrence of isomers of amanitin, termed atropisomers has been postulated. Although the total synthesis of α-amanitin has been accomplished this aspect still remains unsolved. We perform the synthesis of amanitin analogs, accompanied by in-depth spectroscopic, crystallographic and molecular dynamics studies. The data unambiguously confirms the synthesis of two amatoxin-type isomers, for which we propose the term ansamers. The natural structure of the P-ansamer can be ansa-selectively synthesized using an optimized synthetic strategy. We believe that the here described terminology does also have implications for many other peptide structures, e.g. norbornapeptides, lasso peptides, tryptorubins and others, and helps to unambiguously describe conformational isomerism of cyclic peptides.

A Formylglycine-Peptide for the Site-Directed Identification of Phosphotyrosine-Mimetic Fragments

Discovery of protein-binding fragments for precisely defined binding sites is an unmet challenge to date. Herein, formylglycine is investigated as a molecular probe for the sensitive detection of fragments binding to a spatially defined protein site . Formylglycine peptide 3 was derived from a phosphotyrosine-containing peptide substrate of protein tyrosine phosphatase PTP1B by replacing the phosphorylated amino acid with the reactive electrophile. Fragment ligation with formylglycine occurred in situ in aqueous physiological buffer. Structures and kinetics were validated by NMR spectroscopy. Screening and hit validation revealed fluorinated and non-fluorinated hit fragments being able to replace the native phosphotyrosine residue. The formylglycine probe identified low-affinity fragments with high spatial resolution as substantiated by molecular modelling. The best fragment hit, 4-amino-phenyl-acetic acid, was converted into a cellularly active, nanomolar inhibitor of the protein tyrosine phosphatase SHP2.

Pentafluorophosphato‐Phenylalanines: Amphiphilic Phosphotyrosine Mimetics Displaying Fluorine‐Specific Protein Interactions

Phosphotyrosine residues are essential functional switches in health and disease. Thus, phosphotyrosine biomimetics are crucial for the development of chemical tools and drug molecules. We report here the discovery and investigation of pentafluorophosphato amino acids as novel phosphotyrosine biomimetics. A mild acidic pentafluorination protocol was developed and two PF₅‐amino acids were prepared and employed in peptide synthesis. Their structures, reactivities, and fluorine‐specific interactions were studied by NMR and IR spectroscopy, X‐ray diffraction, and in bioactivity assays. The mono‐anionic PF₅ motif displayed an amphiphilic character binding to hydrophobic surfaces, to water molecules, and to protein‐binding sites, exploiting charge and H−F‐bonding interactions. The novel motifs bind 25‐ to 30‐fold stronger to the phosphotyrosine binding site of the protein tyrosine phosphatase PTP1B than the best current biomimetics, as rationalized by computational methods, including molecular dynamics simulations.

Fluorine-induced polarity increases inhibitory activity of BPTI towards chymotrypsin

Substituting the P₁ position in bovine pancreatic trypsin inhibitor (BPTI) is known to heavily influence its inhibitory activity towards serine proteases. Side-chain fluorinated aliphatic amino acids have been shown to alter numerous properties of peptides and proteins and thus are of interest in the context of BPTI. In our study, we systematically investigated the site-specific incorporation of non-canonical amino acids into BPTI by microwave-assisted solid-phase peptide synthesis (SPPS). Inhibitor activity of the variants was tested towards the serine protease α-chymotrypsin. We observed enhanced inhibition of two fluorinated BPTIs compared to wild type and hydrocarbon variants. To further investigate the complexes, we performed X-ray structure analysis. Our findings underline the power fluorine offers as a tool in protein engineering to beneficially alter the effects on phenomena as protein–protein interactions.

Site-specific incorporation of fluorine into bovine pancreatic trypsin inhibitor (BPTI) through chemical protein synthesis, followed by investigations of influence of fluorine on inhibition properties towards serine protease chymotrypsin.

Iodine-Mediated Tryptathionine Formation Facilitates the Synthesis of Amanitins

Synthetic methods on the macrocyclization of peptides are of high interest since they facilitate the synthesis of various types of potentially bioactive compounds, e.g. addressing targets like protein-protein-interactions. Herein, we report on an efficient method to construct tryptathionine-cross-links in peptides between the amino acids Trp and Cys. This reaction not only is the basis for the total synthesis of the death cap toxin α-amanitin but also provides rapid access to various new amanitin analogues. This study for the first time presents a systematic compilation of structure-activity relations (SAR) of amatoxins with regard to RNA polymerase II inhibition and cytotoxicity with one amanitin derivative of superior RNAP II inhibition. The present approach paves the way for the synthesis of structurally diverse amatoxins as future payloads for antibody-toxin conjugates in cancer therapy.

Markov models from the square root approximation of the Fokker-Planck equation: calculating the grid-dependent flux

Molecular dynamics (MD) are extremely complex, yet understanding the slow components of their dynamics is essential to understanding their macroscopic properties. To achieve this, one models the MD as a stochastic process and analyses the dominant eigenfunctions of the associated Fokker-Planck operator, or of closely related transfer operators. So far, the calculation of the discretized operators requires extensive MD simulations. The square-root approximation of the Fokker-Planck equation is a method to calculate transition rates as a ratio of the Boltzmann densities of neighboring grid cells times a flux, and can in principle be calculated without a simulation. In a previous work we still used MD simulations to determine the flux. Here, we propose several methods to calculate the exact or approximate flux for various grid types, and thus estimate the rate matrix without a simulation. Using model potentials we test computational efficiency of the methods, and the accuracy with which they reproduce the dominant eigenfunctions and eigenvalues. For these model potentials, rate matrices with up to [Formula: see text] states can be obtained within seconds on a single high-performance compute server if regular grids are used.

Path probability ratios for Langevin dynamics-Exact and approximate

Path reweighting is a principally exact method to estimate dynamic properties from biased simulations-provided that the path probability ratio matches the stochastic integrator used in the simulation. Previously reported path probability ratios match the Euler-Maruyama scheme for overdamped Langevin dynamics. Since molecular dynamics simulations use Langevin dynamics rather than overdamped Langevin dynamics, this severely impedes the application of path reweighting methods. Here, we derive the path probability ratio M_L for Langevin dynamics propagated by a variant of the Langevin Leapfrog integrator. This new path probability ratio allows for exact reweighting of Langevin dynamics propagated by this integrator. We also show that a previously derived approximate path probability ratio M_approx differs from the exact M_L only by O(ξ⁴Δt⁴) and thus yields highly accurate dynamic reweighting results. (Δt is the integration time step, and ξ is the collision rate.) The results are tested, and the efficiency of path reweighting is explored using butane as an example.

Dynamical reweighting methods for Markov models

Conformational dynamics is essential to biomolecular processes. Markov State Models (MSMs) are widely used to elucidate dynamic properties of molecular systems from unbiased Molecular Dynamics (MD). However, the implementation of reweighting schemes for MSMs to analyze biased simulations is still at an early stage of development. Several dynamical reweighing approaches have been proposed, which can be classified as approaches based on (i) Kramers rate theory, (ii) rescaling of the probability density flux, (iii) reweighting by formulating a likelihood function, (iv) path reweighting. We present the state-of-the-art and discuss the methodological differences of these methods, their limitations and recent applications.

Allostery in C-type lectins

C-type lectins are the largest and most diverse family of mammalian carbohydrate-binding proteins. They share a common protein fold, which provides the unifying basis for calcium-mediated carbohydrate recognition. Their involvement in a multitude of biological functions is remarkable. Here, we review the variety of tasks these lectins are involved in alongside with the structural demands on the overall protein architecture. Subtle changes of the protein structure are implemented to cope with such diverse functional requirements. The presence of a high level of structural dynamics over a broad palette of time scales is paired with the presence of secondary binding sites and allosteric coordination of remote sites and renders this lectin fold a highly adaptable scaffold.

Exploring Rigid and Flexible Core Trivalent Sialosides for Influenza Virus Inhibition

Herein, the chemical synthesis and binding analysis of functionalizable rigid and flexible core trivalent sialosides bearing oligoethylene glycol (OEG) spacers interacting with spike proteins of influenza A virus (IAV) X31 is described. Although the flexible Tris-based trivalent sialosides achieved micromolar binding constants, a trivalent binder based on a rigid adamantane core dominated flexible tripodal compounds with micromolar binding and hemagglutination inhibition constants. Simulation studies indicated increased conformational penalties for long OEG spacers. Using a systematic approach with molecular modeling and simulations as well as biophysical analysis, these findings emphasize on the importance of the scaffold rigidity and the challenges associated with the spacer length optimization.

Estimation of the infinitesimal generator by square-root approximation

In recent years, for the analysis of molecular processes, the estimation of time-scales and transition rates has become fundamental. Estimating the transition rates between molecular conformations is-from a mathematical point of view-an invariant subspace projection problem. We present a method to project the infinitesimal generator acting on function space to a low-dimensional rate matrix. This projection can be performed in two steps. First, we discretize the conformational space in a Voronoi tessellation, then the transition rates between adjacent cells is approximated by the geometric average of the Boltzmann weights of the Voronoi cells. This method demonstrates that there is a direct relation between the potential energy surface of molecular structures and the transition rates of conformational changes. We will show also that this approximation is correct and converges to the generator of the Smoluchowski equation in the limit of infinitely small Voronoi cells. We present results for a two dimensional diffusion process and alanine dipeptide as a high-dimensional system.

Broad substrate tolerance of tubulin tyrosine ligase enables one-step site-specific enzymatic protein labeling

The broad substrate tolerance of tubulin tyrosine ligase enables its wide applicability for protein functionalization.

The broad substrate tolerance of tubulin tyrosine ligase is the basic rationale behind its wide applicability for chemoenzymatic protein functionalization. In this context, we report that the wild-type enzyme enables ligation of various unnatural amino acids that are substantially bigger than and structurally unrelated to the natural substrate, tyrosine, without the need for extensive protein engineering. This unusual substrate flexibility is due to the fact that the enzyme's catalytic pocket forms an extended cavity during ligation, as confirmed by docking experiments and all-atom molecular dynamics simulations. This feature enabled one-step C-terminal biotinylation and fluorescent coumarin labeling of various functional proteins as demonstrated with ubiquitin, an antigen binding nanobody, and the apoptosis marker Annexin V. Its broad substrate tolerance establishes tubulin tyrosine ligase as a powerful tool for in vitro enzyme-mediated protein modification with single functional amino acids in a specific structural context.

Molecular dynamics simulations data of the twenty encoded amino acids in different force fields

We present extensive all-atom Molecular Dynamics (MD) simulation data of the twenty encoded amino acids in explicit water, simulated with different force fields. The termini of the amino acids have been capped to ensure that the dynamics of the Φ and ψ torsion angles are analogues to the dynamics within a peptide chain. We use representatives of each of the four major force field families: AMBER ff-99SBILDN [1], AMBER ff-03 [2], OPLS-AA/L [3], CHARMM27 [4] and GROMOS43a1 [5], [6]. Our data represents a library and test bed for method development for MD simulations and for force fields development. Part of the data set has been previously used for comparison of the dynamic properties of force fields (Vitalini et al., 2015) [7] and for the construction of peptide basis functions for the variational approach to molecular kinetics [8].

Complex RNA folding kinetics revealed by single-molecule FRET and hidden Markov models

We have developed a hidden Markov model and optimization procedure for photon-based single-molecule FRET data, which takes into account the trace-dependent background intensities. This analysis technique reveals an unprecedented amount of detail in the folding kinetics of the Diels–Alderase ribozyme. We find a multitude of extended (low-FRET) and compact (high-FRET) states. Five states were consistently and independently identified in two FRET constructs and at three Mg²⁺ concentrations. Structures generally tend to become more compact upon addition of Mg²⁺. Some compact structures are observed to significantly depend on Mg²⁺ concentration, suggesting a tertiary fold stabilized by Mg²⁺ ions. One compact structure was observed to be Mg²⁺-independent, consistent with stabilization by tertiary Watson–Crick base pairing found in the folded Diels–Alderase structure. A hierarchy of time scales was discovered, including dynamics of 10 ms or faster, likely due to tertiary structure fluctuations, and slow dynamics on the seconds time scale, presumably associated with significant changes in secondary structure. The folding pathways proceed through a series of intermediate secondary structures. There exist both compact pathways and more complex ones, which display tertiary unfolding, then secondary refolding, and, subsequently, again tertiary refolding.

Modulation of a ligand's energy landscape and kinetics by the chemical environment

Understanding how the chemical environment modulates the predominant conformations and kinetics of flexible molecules is a core interest of biochemistry and a prerequisite for the rational design of synthetic catalysts. This study combines molecular dynamics simulation and Markov state models (MSMs) to a systematic computational strategy for investigating the effect of the chemical environment of a molecule on its conformations and kinetics. MSMs allow quantities to be computed that are otherwise difficult to access, such as the metastable sets, their free energies, and the relaxation time scales related to the rare transitions between metastable states. Additionally, MSMs are useful to identify observables that may act as sensors for the conformational or binding state of the molecule, thus guiding the design of experiments. In the present study, the conformation dynamics of UDP-GlcNAc are studied in vacuum, water, water + Mg(2+), and in the protein UDP-GlcNAc 2-epimerase. It is found that addition of Mg(2+) significantly affects the conformational stability, thermodynamics, and kinetics of UDP-GlcNAc. In particular, the slowest structural process, puckering of the GlcNAc sugar, depends on the overall conformation of UDP-GlcNAc and may thus act as a sensor of whether Mg(2+) is bound or not. Interestingly, transferring the molecule from vacuum to water makes the protein-binding conformations UDP-GlcNAc first accessible, while adding Mg(2+) further stabilizes them by specifically associating to binding-competent conformations. While Mg(2+) is not cocrystallized in the UDP-GlcNAc 2-epimerase complex, the selectively stabilized Mg(2+)/UDP-GlcNAc complex may be a template for the bound state, and Mg(2+) may accompany the binding-competent ligand conformation to the binding pocket. This serves as a possible explanation of the enhanced epimerization rate in the presence of Mg(2+). This role of Mg(2+) has previously not been described and opens the question whether "binding co-factors" may be a concept of general relevance for protein-ligand binding.

Maxillary dental arch dimensions following pharyngeal-flap surgery

Pharyngeal-flap surgery (PFS) is a surgical procedure utilized in the treatment of hypernasal speech. The purpose of this study was to evaluate the effect of PFS on the maxillary dental arch. Sixteen patients who had undergone PFS between the ages of 5 and 7 years were selected from the longitudinal growth study of the Lancaster Cleft Palate Clinic. Eight points were marked on the occlusal surfaces of serial maxillary dental models to determine arch width and length. These data were compared to a matched sample of patients from the same growth study who did not undergo PFS. The flap group demonstrated significant reduction in both arch width (at cuspids and molars) and arch length development compared to the control sample, following PFS. These results are not inconsistent with an hypothesis of dental-arch form changes related to possible airway obstruction following PFS, although data on the latter were not available on this sample.