Dead-End Elimination with a Polarizable Force Field Repacks PCNA Structures

GJA3 Genetic Variation and Autosomal Dominant Congenital Cataracts and Glaucoma Following Cataract Surgery

Importance

The p.Asp67Tyr genetic variant in the GJA3 gene is responsible for congenital cataracts in a family with a high incidence of glaucoma following cataract surgery.

Objective

To describe the clinical features of a family with a strong association between congenital cataracts and glaucoma following cataract surgery secondary to a genetic variant in the GJA3 gene (NM_021954.4:c.199G>T, p.Asp67Tyr).

Design, Setting, and Participants

This was a retrospective, observational, case series, genetic association study from the University of Iowa spanning 61 years. Examined were the ophthalmic records from 1961 through 2022 of the family members of a 4-generation pedigree with autosomal dominant congenital cataracts.

Main Outcomes and Measures

Frequency of glaucoma following cataract surgery and postoperative complications among family members with congenital cataract due to the p.Asp67Tyr GJA3 genetic variant.

Results

Medical records were available from 11 of 12 family members (7 male [63.6%]) with congenital cataract with a mean (SD) follow-up of 30 (21.7) years (range, 0.2-61 years). Eight of 9 patients with congenital cataracts developed glaucoma, and 8 of 8 patients who had cataract surgery at age 2 years or younger developed glaucoma following cataract surgery. The only family member with congenital cataracts who did not develop glaucoma had delayed cataract surgery until 12 and 21 years of age. Five of 11 family members (45.5%) had retinal detachments after cataract extraction and vitrectomy. No patients developed retinal detachments after prophylactic 360-degree endolaser.

Conclusions and Relevance

The GJA3 genetic variant, p.Asp67Tyr, was identified in a 4-generation congenital cataract pedigree from Iowa. This report suggests that patients with congenital cataract due to some GJA3 genetic variants may be at especially high risk for glaucoma following cataract surgery. Retinal detachments after cataract extraction in the first 2 years of life were also common in this family, and prophylactic retinal endolaser may be indicated at the time of surgery.

Implicit Solvents for the Polarizable Atomic Multipole AMOEBA Force Field

Computational protein design, ab initio protein/RNA folding, and protein-ligand screening can be too computationally demanding for explicit treatment of solvent. For these applications, implicit solvent offers a compelling alternative, which we describe here for the polarizable atomic multipole AMOEBA force field based on three treatments of continuum electrostatics: numerical solutions to the nonlinear and linearized versions of the Poisson-Boltzmann equation (PBE), the domain-decomposition conductor-like screening model (ddCOSMO) approximation to the PBE, and the analytic generalized Kirkwood (GK) approximation. The continuum electrostatics models are combined with a nonpolar estimator based on novel cavitation and dispersion terms. Electrostatic model parameters are numerically optimized using a least-squares style target function based on a library of 103 small-molecule solvation free energy differences. Mean signed errors for the adaptive Poisson-Boltzmann solver (APBS), ddCOSMO, and GK models are 0.05, 0.00, and 0.00 kcal/mol, respectively, while the mean unsigned errors are 0.70, 0.63, and 0.58 kcal/mol, respectively. Validation of the electrostatic response of the resulting implicit solvents, which are available in the Tinker (or Tinker-HP), OpenMM, and Force Field X software packages, is based on comparisons to explicit solvent simulations for a series of proteins and nucleic acids. Overall, the emergence of performative implicit solvent models for polarizable force fields opens the door to their use for folding and design applications.

Syngeneic homograft of framework regions enhances the affinity of the mouse anti-human epidermal receptor 2 single-chain antibody e23sFv

e23sFv is a HER2-targeted single-chain variable fragment (scFV) that was characterized as the targeting portion of a HER2-targeted tumour proapoptotic molecule in our previous study. In vitro antibody affinity maturation is a method to enhance antibody affinity either by complementarity-determining region (CDR) mutagenesis or by framework region (FR) engraftment. In the present study, the affinity of e23sFv was enhanced using two strategies. In one approach, site-directed mutations were introduced into the FRs of e23sFv (designated EMEY), and in the other approach e23sFv FRs were substituted with FRs from the most homologous screened antibodies (designated EX1 and EX2). Notably, EX1 derived from the FR engraftment strategy demonstrated a 4-fold higher affinity for HER2 compared with e23sFv and was internalized into HER2-overexpressing cells; however, EMEY and EX2 exhibited reduced affinity for HER2 and decreased internalization potential compared with EX1. The 3D structure of EX1 and the HER2-EX1 complex was acquired using molecular homology modelling and docking and the HER2 epitopes of EX1 and the molecular interaction energy of the EX1-HER2 complex were predicted. In the present study, it was demonstrated that scFv affinity improvement based on sequence alignment was feasible and effective. Moreover, the FR grafting strategy was indicated to be more effective and simple compared with site-directed mutagenesis to improve e23sFv affinity. In conclusion, it was indicated that the affinity-improved candidate EX1 may present a great potential for the diagnosis and treatment of HER2-overexpressing tumours.

Effective mismatch repair depends on timely control of PCNA retention on DNA by the Elg1 complex

Proliferating cell nuclear antigen (PCNA) is a sliding clamp that acts as a central co-ordinator for mismatch repair (MMR) as well as DNA replication. Loss of Elg1, the major subunit of the PCNA unloader complex, causes over-accumulation of PCNA on DNA and also increases mutation rate, but it has been unclear if the two effects are linked. Here we show that timely removal of PCNA from DNA by the Elg1 complex is important to prevent mutations. Although premature unloading of PCNA generally increases mutation rate, the mutator phenotype of elg1Δ is attenuated by PCNA mutants PCNA-R14E and PCNA-D150E that spontaneously fall off DNA. In contrast, the elg1Δ mutator phenotype is exacerbated by PCNA mutants that accumulate on DNA due to enhanced electrostatic PCNA–DNA interactions. Epistasis analysis suggests that PCNA over-accumulation on DNA interferes with both MMR and MMR-independent process(es). In elg1Δ, over-retained PCNA hyper-recruits the Msh2–Msh6 mismatch recognition complex through its PCNA-interacting peptide motif, causing accumulation of MMR intermediates. Our results suggest that PCNA retention controlled by the Elg1 complex is critical for efficient MMR: PCNA needs to be on DNA long enough to enable MMR, but if it is retained too long it interferes with downstream repair steps.

Protein Design by Provable Algorithms

Protein design algorithms can leverage provable guarantees of accuracy to provide new insights and unique optimized molecules.

Structural Insights into Hearing Loss Genetics from Polarizable Protein Repacking

Hearing loss is associated with ∼8100 mutations in 152 genes, and within the coding regions of these genes are over 60,000 missense variants. The majority of these variants are classified as "variants of uncertain significance" to reflect our inability to ascribe a phenotypic effect to the observed amino acid change. A promising source of pathogenicity information is biophysical simulation, although input protein structures often contain defects because of limitations in experimental data and/or only distant homology to a template. Here, we combine the polarizable atomic multipole optimized energetics for biomolecular applications force field, many-body optimization theory, and graphical processing unit acceleration to repack all deafness-associated proteins and thereby improve average structure MolProbity score from 2.2 to 1.0. We then used these optimized wild-type models to create over 60,000 structures for missense variants in the Deafness Variation Database, which are being incorporated into the Deafness Variation Database to inform deafness pathogenicity prediction. Finally, this work demonstrates that advanced polarizable atomic multipole force fields are efficient enough to repack the entire human proteome.

Conformational Flexibility of Ubiquitin-Modified and SUMO-Modified PCNA Shown by Full-Ensemble Hybrid Methods

Ubiquitin-modified proliferating cell nuclear antigen (PCNA) and small ubiquitin-like modifier (SUMO)-modified PCNA regulate DNA damage tolerance pathways. X-ray crystal structures of these proteins suggested that they do not have much conformational flexibility because the modifiers have preferred binding sites on the surface of PCNA. By contrast, small-angle X-ray scattering analyses of these proteins suggested that they have different degrees of conformational flexibility, with SUMO-modified PCNA being more flexible. These conclusions were based on minimal-ensemble hybrid approaches, which produce unrealistic models by representing flexible proteins with only a few static structures. To overcome the limitations of minimal-ensemble hybrid approaches and to determine the degree of conformational flexibility of ubiquitin-modified PCNA and SUMO-modified PCNA, we utilized a novel full-ensemble hybrid approach. We carried out molecular simulations and small-angle X-ray scattering analyses of both proteins and obtained outstanding agreement between the full ensembles generated by the simulations and the experimental data. We found that both proteins have a high degree of conformational flexibility. The modifiers occupy many positions around the back and side of the PCNA ring. Moreover, we found no preferred ubiquitin-binding or SUMO-binding sites on PCNA. This conformational flexibility likely facilitates the recognition of downstream effector proteins and the formation of PCNA tool belts.

Tinker 8: Software Tools for Molecular Design

The Tinker software, currently released as version 8, is a modular molecular mechanics and dynamics package written primarily in a standard, easily portable dialect of Fortran 95 with OpenMP extensions. It supports a wide variety of force fields, including polarizable models such as the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field. The package runs on Linux, macOS, and Windows systems. In addition to canonical Tinker, there are branches, Tinker-HP and Tinker-OpenMM, designed for use on message passing interface (MPI) parallel distributed memory supercomputers and state-of-the-art graphical processing units (GPUs), respectively. The Tinker suite also includes a tightly integrated Java-based graphical user interface called Force Field Explorer (FFE), which provides molecular visualization capabilities as well as the ability to launch and control Tinker calculations.

Computer-aided design of amino acid-based therapeutics: a review

During the last two decades, the pharmaceutical industry has progressed from detecting small molecules to designing biologic-based therapeutics. Amino acid-based drugs are a group of biologic-based therapeutics that can effectively combat the diseases caused by drug resistance or molecular deficiency. Computational techniques play a key role to design and develop the amino acid-based therapeutics such as proteins, peptides and peptidomimetics. In this study, it was attempted to discuss the various elements for computational design of amino acid-based therapeutics. Protein design seeks to identify the properties of amino acid sequences that fold to predetermined structures with desirable structural and functional characteristics. Peptide drugs occupy a middle space between proteins and small molecules and it is hoped that they can target "undruggable" intracellular protein-protein interactions. Peptidomimetics, the compounds that mimic the biologic characteristics of peptides, present refined pharmacokinetic properties compared to the original peptides. Here, the elaborated techniques that are developed to characterize the amino acid sequences consistent with a specific structure and allow protein design are discussed. Moreover, the key principles and recent advances in currently introduced computational techniques for rational peptide design are spotlighted. The most advanced computational techniques developed to design novel peptidomimetics are also summarized.

LADD syndrome with glaucoma is caused by a novel gene

PURPOSE

Lacrimo-auriculo-dento-digital (LADD) syndrome is an autosomal dominant disorder displaying variable expression of multiple congenital anomalies including hypoplasia or aplasia of the lacrimal and salivary systems causing abnormal tearing and dry mouth. Mutations in the FGF10, FGFR2, and FGFR3 genes were found to cause some cases of LADD syndrome in prior genetic studies. The goal of this study is to identify the genetic basis of a case of LADD syndrome with glaucoma and thin central corneal thickness (CCT).

METHODS

Whole exome sequencing was performed, and previously described disease-causing genes (FGF10, FGFR2, and FGFR3) were first evaluated for mutations. Fifty-eight additional prioritized candidate genes were identified by searching gene annotations for features of LADD syndrome. The potential pathogenicity of the identified mutations was assessed by determining their frequency in large public exome databases; through sequence analysis using the Blosum62 matrix, PolyPhen2, and SIFT algorithms; and through homology analyses. A structural analysis of the effects of the top candidate mutation in tumor protein 63 (TP63) was also conducted by superimposing the mutation over the solved crystal structure.

RESULTS

No mutations were detected in FGF10, FGFR2, or FGFR3. The LADD syndrome patient's exome data was searched for mutations in the 58 candidate genes and only one mutation was detected, an Arg343Trp mutation in the tumor protein 63 (TP63) gene. This TP63 mutation is absent from the gnomAD sequence database. Analysis of the Arg343Trp mutation with Blosum62, PolyPhen2, and SIFT all suggest it is pathogenic. This arginine residue is highly conserved in orthologous genes. Finally, crystal structure analysis showed that the Arg343Trp mutation causes a significant alteration in the structure of TP63's DNA binding domain.

CONCLUSIONS

We report a patient with no mutations in known LADD syndrome genes (FGF10, FGFR2, and FGFR3). Our analysis provides strong evidence that the Arg343Trp mutation in TP63 caused LADD syndrome in our patient and that TP63 is a fourth gene contributing to this condition. TP63 encodes a transcription factor involved in the development and differentiation of tissues affected by LADD syndrome. These data suggest that TP63 is a novel LADD syndrome gene and may also influence corneal thickness and risk for open-angle glaucoma.

Toward polarizable AMOEBA thermodynamics at fixed charge efficiency using a dual force field approach: application to organic crystals

First principles prediction of the structure, thermodynamics and solubility of organic molecular crystals, which play a central role in chemical, material, pharmaceutical and engineering sciences, challenges both potential energy functions and sampling methodologies. Here we calculate absolute crystal deposition thermodynamics using a novel dual force field approach whose goal is to maintain the accuracy of advanced multipole force fields (e.g. the polarizable AMOEBA model) while performing more than 95% of the sampling in an inexpensive fixed charge (FC) force field (e.g. OPLS-AA). Absolute crystal sublimation/deposition phase transition free energies were determined using an alchemical path that grows the crystalline state from a vapor reference state based on sampling with the OPLS-AA force field, followed by dual force field thermodynamic corrections to change between FC and AMOEBA resolutions at both end states (we denote the three step path as AMOEBA/FC). Importantly, whereas the phase transition requires on the order of 200 ns of sampling per compound, only 5 ns of sampling was needed for the dual force field thermodynamic corrections to reach a mean statistical uncertainty of 0.05 kcal mol-1. For five organic compounds, the mean unsigned error between direct use of AMOEBA and the AMOEBA/FC dual force field path was only 0.2 kcal mol-1 and not statistically significant. Compared to experimental deposition thermodynamics, the mean unsigned error for AMOEBA/FC (1.4 kcal mol-1) was more than a factor of two smaller than uncorrected OPLS-AA (3.2 kcal mol-1). Overall, the dual force field thermodynamic corrections reduced condensed phase sampling in the expensive force field by a factor of 40, and may prove useful for protein stability or binding thermodynamics in the future.

LUTE (Local Unpruned Tuple Expansion): Accurate Continuously Flexible Protein Design with General Energy Functions and Rigid Rotamer-Like Efficiency

Most protein design algorithms search over discrete conformations and an energy function that is residue-pairwise, that is, a sum of terms that depend on the sequence and conformation of at most two residues. Although modeling of continuous flexibility and of non-residue-pairwise energies significantly increases the accuracy of protein design, previous methods to model these phenomena add a significant asymptotic cost to design calculations. We now remove this cost by modeling continuous flexibility and non-residue-pairwise energies in a form suitable for direct input to highly efficient, discrete combinatorial optimization algorithms such as DEE/A* or branch-width minimization. Our novel algorithm performs a local unpruned tuple expansion (LUTE), which can efficiently represent both continuous flexibility and general, possibly nonpairwise energy functions to an arbitrary level of accuracy using a discrete energy matrix. We show using 47 design calculation test cases that LUTE provides a dramatic speedup in both single-state and multistate continuously flexible designs.

Re-evaluation of low-resolution crystal structures via interactive molecular-dynamics flexible fitting (iMDFF): a case study in complement C4

While the rapid proliferation of high-resolution structures in the Protein Data Bank provides a rich set of templates for starting models, it remains the case that a great many structures both past and present are built at least in part by hand-threading through low-resolution and/or weak electron density. With current model-building tools this task can be challenging, and the de facto standard for acceptable error rates (in the form of atomic clashes and unfavourable backbone and side-chain conformations) in structures based on data with dmax not exceeding 3.5 Å reflects this. When combined with other factors such as model bias, these residual errors can conspire to make more serious errors in the protein fold difficult or impossible to detect. The three recently published 3.6-4.2 Å resolution structures of complement C4 (PDB entries 4fxg, 4fxk and 4xam) rank in the top quartile of structures of comparable resolution both in terms of Rfree and MolProbity score, yet, as shown here, contain register errors in six β-strands. By applying a molecular-dynamics force field that explicitly models interatomic forces and hence excludes most physically impossible conformations, the recently developed interactive molecular-dynamics flexible fitting (iMDFF) approach significantly reduces the complexity of the conformational space to be searched during manual rebuilding. This substantially improves the rate of detection and correction of register errors, and allows user-guided model building in maps with a resolution lower than 3.5 Å to converge to solutions with a stereochemical quality comparable to atomic resolution structures. Here, iMDFF has been used to individually correct and re-refine these three structures to MolProbity scores of <1.7, and strategies for working with such challenging data sets are suggested. Notably, the improved model allowed the resolution for complement C4b to be extended from 4.2 to 3.5 Å as demonstrated by paired refinement.

Identification of New Mutations at the PCNA Subunit Interface that Block Translesion Synthesis

Proliferating cell nuclear antigen (PCNA) plays an essential role in DNA replication and repair by interacting with a large number of proteins involved in these processes. Two amino acid substitutions in PCNA, both located at the subunit interface, have previously been shown to block translesion synthesis (TLS), a pathway for bypassing DNA damage during replication. To better understand the role of the subunit interface in TLS, we used random mutagenesis to generate a set of 33 PCNA mutants with substitutions at the subunit interface. We assayed the full set of mutants for viability and sensitivity to ultraviolet (UV) radiation. We then selected a subset of 17 mutants and measured their rates of cell growth, spontaneous mutagenesis, and UV-induced mutagenesis. All except three of these 17 mutants were partially or completely defective in induced mutagenesis, which indicates a partial or complete loss of TLS. These results demonstrate that the integrity of the subunit interface of PCNA is essential for efficient TLS and that even conservative substitutions have the potential to disrupt this process.

Comparing three stochastic search algorithms for computational protein design: Monte Carlo, replica exchange Monte Carlo, and a multistart, steepest-descent heuristic

Computational protein design depends on an energy function and an algorithm to search the sequence/conformation space. We compare three stochastic search algorithms: a heuristic, Monte Carlo (MC), and a Replica Exchange Monte Carlo method (REMC). The heuristic performs a steepest-descent minimization starting from thousands of random starting points. The methods are applied to nine test proteins from three structural families, with a fixed backbone structure, a molecular mechanics energy function, and with 1, 5, 10, 20, 30, or all amino acids allowed to mutate. Results are compared to an exact, "Cost Function Network" method that identifies the global minimum energy conformation (GMEC) in favorable cases. The designed sequences accurately reproduce experimental sequences in the hydrophobic core. The heuristic and REMC agree closely and reproduce the GMEC when it is known, with a few exceptions. Plain MC performs well for most cases, occasionally departing from the GMEC by 3-4 kcal/mol. With REMC, the diversity of the sequences sampled agrees with exact enumeration where the latter is possible: up to 2 kcal/mol above the GMEC. Beyond, room temperature replicas sample sequences up to 10 kcal/mol above the GMEC, providing thermal averages and a solution to the inverse protein folding problem. © 2016 Wiley Periodicals, Inc.