Systematically testing human HMBS missense variants to reveal mechanism and pathogenic variation

Defects in hydroxymethylbilane synthase (HMBS) can cause acute intermittent porphyria (AIP), an acute neurological disease. Although sequencing-based diagnosis can be definitive, ∼⅓ of clinical HMBS variants are missense variants, and most clinically reported HMBS missense variants are designated as "variants of uncertain significance" (VUSs). Using saturation mutagenesis, en masse selection, and sequencing, we applied a multiplexed validated assay to both the erythroid-specific and ubiquitous isoforms of HMBS, obtaining confident functional impact scores for >84% of all possible amino acid substitutions. The resulting variant effect maps generally agreed with biochemical expectations and provide further evidence that HMBS can function as a monomer. Additionally, the maps implicated specific residues as having roles in active site dynamics, which was further supported by molecular dynamics simulations. Most importantly, these maps can help discriminate pathogenic from benign HMBS variants, proactively providing evidence even for yet-to-be-observed clinical missense variants.

Identification of Catechins' Binding Sites in Monomeric Aβ42 through Ensemble Docking and MD Simulations

The assembly of the amyloid-β peptide (Aβ) into toxic oligomers and fibrils is associated with Alzheimer's disease and dementia. Therefore, disrupting amyloid assembly by direct targeting of the Aβ monomeric form with small molecules or antibodies is a promising therapeutic strategy. However, given the dynamic nature of Aβ, standard computational tools cannot be easily applied for high-throughput structure-based virtual screening in drug discovery projects. In the current study, we propose a computational pipeline-in the framework of the ensemble docking strategy-to identify catechins' binding sites in monomeric Aβ₄₂. It is shown that both hydrophobic aromatic interactions and hydrogen bonding are crucial for the binding of catechins to Aβ₄₂. Additionally, it has been found that all the studied ligands, especially EGCG, can act as potent inhibitors against amyloid aggregation by blocking the central hydrophobic region of Aβ. Our findings are evaluated and confirmed with multi-microsecond MD simulations. Finally, it is suggested that our proposed pipeline, with low computational cost in comparison with MD simulations, is a suitable approach for the virtual screening of ligand libraries against Aβ.

Systematically testing human HMBS missense variants to reveal mechanism and pathogenic variation

Defects in hydroxymethylbilane synthase (HMBS) can cause Acute Intermittent Porphyria (AIP), an acute neurological disease. Although sequencing-based diagnosis can be definitive, ~⅓ of clinical HMBS variants are missense variants, and most clinically-reported HMBS missense variants are designated as "variants of uncertain significance" (VUS). Using saturation mutagenesis, en masse selection, and sequencing, we applied a multiplexed validated assay to both the erythroid-specific and ubiquitous isoforms of HMBS, obtaining confident functional impact scores for >84% of all possible amino-acid substitutions. The resulting variant effect maps generally agreed with biochemical expectation. However, the maps showed variants at the dimerization interface to be unexpectedly well tolerated, and suggested residue roles in active site dynamics that were supported by molecular dynamics simulations. Most importantly, these HMBS variant effect maps can help discriminate pathogenic from benign variants, proactively providing evidence even for yet-to-be-observed clinical missense variants.

Dawn of a New Era for Membrane Protein Design

A major advancement has recently occurred in the ability to predict protein secondary structure from sequence using artificial neural networks. This new accessibility to high-quality predicted structures provides a big opportunity for the protein design community. It is particularly welcome for membrane protein design, where the scarcity of solved structures has been a major limitation of the field for decades. Here, we review the work done to date on the membrane protein design and set out established and emerging tools that can be used to most effectively exploit this new access to structures.

Free Energy and Stacking of Eumelanin Nanoaggregates

Eumelanin, a member of the melanin family, is a black-brown insoluble pigment. It possesses a broad range of properties such as antioxidation, free radical scavenging, photoprotection, and charge carrier transportation. Surprisingly, the exact molecular structure of eumelanin remains undefined. It is, however, generally considered to consist of two main building blocks, 5,6-dihydroxyindole (DHI) and 5,6- dihydroxyindole carboxylic acid (DHICA). We focus on DHI and report, for the first time, a computational investigation of the structural properties of DHI-eumelanin aggregates in aqueous solutions. First, multimicrosecond molecular dynamics (MD) simulations at different concentrations were performed to investigate the aggregation and ordering of tetrameric DHI-eumelanin protomolecules. This was followed by umbrella sampling (US) and density functional theory (DFT) calculations to study the physical mechanisms of stacking. Aggregation occurs through formation of nanoscale stacks and was observed in all systems. Further analyses showed that aggregation and coarsening of the domains is due to a decrease in hydrogen bonds between the eumelanins and water; while domains exist, there is no long-range order. The results show noncovalent stacks with the interlayer distance between eumelanin protomolecules being less than 3.5 Å. This is in good agreement with transmission electron microscopy data. Both free energy calculations and DFT revealed strong stacking interactions. The electrostatic potential map provides an explanation and a rationale for the slightly sheared relative orientations and, consequently, for the curved shapes of the nanoscale domains.

Biphasic Proton Transport Mechanism for Uncoupling Proteins

It has been suggested that uncoupling proteins (UCPs) transport protons via interconversion between two conformational states: one in the "cytoplasmic state" and the other in the "matrix state". Matrix and cytoplasmic salt-bridge networks are key controllers of these states. This study proposes a mechanism for proton transport in tetrameric UCP2, with focus on the role of the matrix network. Eleven mutants were prepared to disrupt (K → Q or D → N mutations) or alter (K → D and D → K mutations) the salt-bridges in the matrix network. Proteins were recombinantly expressed in Escherichia coli membrane, reconstituted in model lipid membranes, and their structures and functions were analyzed by gel electrophoresis, circular dichroism spectroscopy, fluorescence assays, as well as molecular dynamics simulations. It is shown that the UCP2 matrix network contains five salt-bridges (rather than the previously reported three), and the matrix network can regulate the proton transport by holding the protein's transmembrane helices in close proximity, limiting the movement of the activator fatty acid(s). A biphasic two-state molecular model is proposed for proton transport in tetrameric (a dimer of stable dimers) UCP2, in which all the monomers are functional, and monomers in each dimer are in the same transport mode. Purine nucleotide (e.g., ATP) can occlude the internal pore of the monomeric units of UCP tetramers via interacting with positive residues at or in the proximity of the matrix network (K38, K141, K239, R88, R185, and R279) and prevent switching between cytoplasmic and matrix states, thus inhibiting the proton transport. This study provides new insights into the mechanism of proton transport and regulation in UCPs.

Functional Oligomeric Forms of Uncoupling Protein 2: Strong Evidence for Asymmetry in Protein and Lipid Bilayer Systems

Stoichiometry of uncoupling proteins (UCPs) and their coexistence as functional monomeric and associated forms in lipid membranes remain intriguing open questions. In this study, tertiary and quaternary structures of UCP2 were analyzed experimentally and through molecular dynamics (MD) simulations. UCP2 was overexpressed in the inner membrane of Escherichia coli, then purified and reconstituted in lipid vesicles. Structure and proton transport function of UCP2 were characterized by circular dichroism (CD) spectroscopy and fluorescence methods. Findings suggest a tetrameric functional form for UCP2. MD simulations conclude that tetrameric UCP2 is a dimer of dimers, is more stable than its monomeric and dimeric forms, is asymmetrical and induces asymmetry in the membrane's lipid structure, and a biphasic on-off switch between the dimeric units is its possible mode of transport. MD simulations also show that the water density inside the UCP2 monomer is asymmetric, with the cytoplasmic side having a higher water density and a wider radius. In contrast, the structurally comparable adenosine 5'-diphosphate (ADP)/adenosine 5'-triphosphate (ATP) carrier (AAC1) did not form tetramers, implying that tetramerization cannot be generalized to all mitochondrial carriers.

Manipulation of Diatomic Molecules with Oriented External Electric Fields: Linear Correlations in Atomic Properties Lead to Nonlinear Molecular Responses

Oriented external electric fields (OEEFs) have been shown to have great potential in being able to provide unprecedented control of chemical reactions, catalysis, and selectivity with applications ranging from H₂ storage to molecular machines. We report a theoretical study of the atomic origins of molecular changes because of OEEFs since understanding the characteristics of OEEF-induced couplings between atomic and molecular properties is an important step toward comprehensive understanding of the effects of strong external fields on the molecular structure, stability, and reactivity. We focus on the atomic and molecular (bond) properties of a set of homo- (H₂, N₂, O₂, F₂, and Cl₂) and heterodiatomic (HF, HCl, CO, and NO) molecules under intense external electric fields in the context of quantum theory of atoms in molecules (QTAIM). It is shown that the atomic properties (atomic charges, energies, and localization indices) correlate linearly with the field strengths, but molecular properties (bond length, electron density at the bond critical point, and electron delocalization index) exhibit nonlinear responses to the imposed fields. In particular, the changes in the electron density distribution alter the shapes and locations of the zero-flux surfaces, atomic volumes, atomic electron population, and localization/delocalization indices. The topography and topology of the molecular electrostatic potential undergo dramatic changes. External fields also perturb the covalent-polar-ionic characteristic of the studied chemical bonds, hallmarking the impact of electric fields on the stability and reactivity of chemical compounds. The findings are well-rationalized within the framework of the QTAIM and form a coherent conceptual understanding of these effects in prototypical diatomic molecules.

Membrane Disruption by Very Long Chain Fatty Acids during Necroptosis

Necroptosis is a form of regulated cell death which results in loss of plasma membrane integrity, release of intracellular contents, and an associated inflammatory response. We previously found that saturated very long chain fatty acids (VLCFAs), which contain ≥20 carbons, accumulate during necroptosis. Here, we show that genetic knockdown of Fatty Acid (FA) Elongase 7 (ELOVL7) reduces accumulation of specific very long chain FAs during necroptosis, resulting in reduced necroptotic cell death and membrane permeabilization. Conversely, increasing the expression of ELOVL7 increases very long chain fatty acids and membrane permeabilization. In vitro, introduction of the VLCFA C24 FA disrupts bilayer integrity in liposomes to a greater extent than a conventional C16 FA. To investigate the microscopic origin of these observations, atomistic Molecular Dynamics (MD) simulations were performed. MD simulations suggest that fatty acids cause clear differences in bilayers based on length and that it is the interdigitation of C24 FA between the individual leaflets that results in disorder in the region and, consequently, membrane disruption. We synthesized clickable VLCFA analogs and observed that many proteins were acylated by VLCFAs during necroptosis. Taken together, these results confirm the active role of VLCFAs during necroptosis and point to multiple potential mechanisms of membrane disruption including direct permeabilization via bilayer disruption and permeabilization by targeting of proteins to cellular membranes by fatty acylation.

Manipulation of a DNA aptamer-protein binding site through arylation of internal guanine residues

Chemically modified aptamers have the opportunity to increase aptamer target binding affinity and provide structure-activity relationships to enhance our understanding of molecular target recognition by the aptamer fold. In the current study, 8-aryl-2'-deoxyguanosine nucleobases have been inserted into the G-tetrad and central TGT loop of the thrombin binding aptamer (TBA) to determine their impact on antiparallel G-quadruplex (GQ) folding and thrombin binding affinity. The aryl groups attached to the dG nucleobase vary greatly in aryl ring size and impact on GQ stability (∼20 °C change in GQ thermal melting (Tm) values) and thrombin binding affinity (17-fold variation in dissociation constant (Kd)). At G8 of the central TGT loop that is distal from the aptamer recognition site, the probes producing the most stable GQ structure exhibited the strongest thrombin binding affinity. However, within the G-tetrad, changes to the electron density of the dG component within the modified nucleobase can diminish thrombin binding affinity. Detailed molecular dynamics (MD) simulations on the modified TBA (mTBA) and mTBA-protein complexes demonstrate how the internal 8-aryl-dG modification can manipulate the interactions between the DNA nucleobases and the amino acid residues of thrombin. These results highlight the potential of internal fluorescent nuclobase analogs (FBAs) to broaden design options for aptasensor development.

Structural Insight into the Discrimination between 8-Oxoguanine Glycosidic Conformers by DNA Repair Enzymes: A Molecular Dynamics Study of Human Oxoguanine Glycosylase 1 and Formamidopyrimidine-DNA Glycosylase

hOgg1 and FPG are the primary DNA repair enzymes responsible for removing the major guanine (G) oxidative product, namely, 7,8-dihydro-8-oxoguanine (OG), in humans and bacteria, respectively. While natural G adopts the anti conformation and forms a Watson-Crick pair with cytosine (C), OG can also adopt the syn conformation and form a Hoogsteen pair with adenine (A). hOgg1 removes OG paired with C but is inactive toward the OG:A pair. In contrast, FPG removes OG from OG:C pairs and also exhibits appreciable (although diminished) activity toward OG:A pairs. As a first step toward understanding this difference in activity, we have employed molecular dynamics simulations to examine how the anti and syn conformers of OG are accommodated in the hOgg1 and FPG active sites. When anti-OG is bound, hOgg1 active site residues are properly aligned to initiate catalytic base departure, while geometrical parameters required for the catalytic reaction are not conserved for syn-OG. On the other hand, the FPG catalytic residues are suitably aligned for both OG conformers, with anti-OG being more favorably bound. Thus, our data suggests that the differential ability of hOgg1 and FPG to accommodate the anti- and syn-OG glycosidic conformations is an important factor that contributes to the relative experimental excision rates. Nevertheless, the positions of the nucleophiles with respect to the lesion in the active sites suggest that the reactant complex is poised to initiate catalysis through a similar mechanism for both repair enzymes and supports a recently proposed mechanism in which sugar-ring opening precedes nucleoside deglycosylation.

Quantum mechanical study of the β- and δ-lyase reactions during the base excision repair process: application to FPG

The β- and δ-elimination reactions catalyzed by FPG during the base excision repair of 8-oxoguanine are intrinsically different.