LoCoHD: a metric for comparing local environments of proteins

Protein folds and the local environments they create can be compared using a variety of differently designed measures, such as the root mean squared deviation, the global distance test, the template modeling score or the local distance difference test. Although these measures have proven to be useful for a variety of tasks, each fails to fully incorporate the valuable chemical information inherent to atoms and residues, and considers these only partially and indirectly. Here, we develop the highly flexible local composition Hellinger distance (LoCoHD) metric, which is based on the chemical composition of local residue environments. Using LoCoHD, we analyze the chemical heterogeneity of amino acid environments and identify valines having the most conserved-, and arginines having the most variable chemical environments. We use LoCoHD to investigate structural ensembles, to evaluate critical assessment of structure prediction (CASP) competitors, to compare the results with the local distance difference test (lDDT) scoring system, and to evaluate a molecular dynamics simulation. We show that LoCoHD measurements provide unique information about protein structures that is distinct from, for example, those derived using the alignment-based RMSD metric, or the similarly distance matrix-based but alignment-free lDDT metric.

Polymorphic amyloid nanostructures of hormone peptides involved in glucose homeostasis display reversible amyloid formation

A large group of hormones are stored as amyloid fibrils in acidic secretion vesicles before they are released into the bloodstream and readopt their functional state. Here, we identify an evolutionarily conserved hexapeptide sequence as the major aggregation-prone region (APR) of gastrointestinal peptides of the glucagon family: xFxxWL. We determine nine polymorphic crystal structures of the APR segments of glucagon-like peptides 1 and 2, and exendin and its derivatives. We follow amyloid formation by CD, FTIR, ThT assays, and AFM. We propose that the pH-dependent changes of the protonation states of glutamate/aspartate residues of APRs initiate switching between the amyloid and the folded, monomeric forms of the hormones. We find that pH sensitivity diminishes in the absence of acidic gatekeepers and amyloid formation progresses over a broad pH range. Our results highlight the dual role of short aggregation core motifs in reversible amyloid formation and receptor binding.

A carbapenem antibiotic inhibiting a mammalian serine protease: structure of the acylaminoacyl peptidase-meropenem complex

The structure of porcine AAP (pAAP) in a covalently bound complex with meropenem was determined by cryo-EM to 2.1 Å resolution, showing the mammalian serine-protease inhibited by a carbapenem antibiotic. AAP is a modulator of the ubiquitin-proteasome degradation system and the site of a drug-drug interaction between the widely used antipsychotic, valproate and carbapenems. The active form of pAAP - a toroidal tetramer - binds four meropenem molecules covalently linked to the catalytic Ser587 of the serine-protease triad, in an acyl-enzyme state. AAP is hindered from fully processing the antibiotic by the displacement and protonation of His707 of the catalytic triad. We show that AAP is made susceptible to the association by its unusually sheltered active pockets and flexible catalytic triads, while the carbapenems possess sufficiently small substituents on their β-lactam rings to fit into the shallow substrate-specificity pocket of the enzyme.

The Importance of Mg2+ -Free State in Nucleotide Exchange of Oncogenic K-Ras Mutants

For efficient targeting of oncogenic K-Ras interaction sites, a mechanistic picture of the Ras-cycle is necessary. Herein, we used NMR relaxation techniques and molecular dynamics simulations to decipher the role of slow dynamics in wild-type and three oncogenic P-loop mutants of K-Ras. Our measurements reveal a dominant two-state conformational exchange on the ms timescale in both GDP- and GTP-bound K-Ras. The identified low-populated higher energy state in GDP-loaded K-Ras has a conformation reminiscent of a nucleotide-bound/Mg²⁺ -free state characterized by shortened β2/β3-strands and a partially released switch-I region preparing K-Ras for the interaction with the incoming nucleotide exchange factor and subsequent reactivation. By providing insight into mutation-specific differences in K-Ras structural dynamics, our systematic analysis improves our understanding of prolonged K-Ras signaling and may aid the development of allosteric inhibitors targeting nucleotide exchange in K-Ras.

Omicron Binding Mode: Contact Analysis and Dynamics of the Omicron Receptor-Binding Domain in Complex with ACE2

On 26 November 2021, the WHO classified the Omicron variant of the SARS-CoV-2 virus (B.1.1.529 lineage) as a variant of concern (VOC) (COVID-19 Variant Data, Department of Health, 2022). The Omicron variant contains as many as 26 unique mutations of effects not yet determined (Venkatakrishnan, A., Open Science Framework, 2021). Out of its total of 34 Spike protein mutations, 15 are located on the receptor-binding domain (S-RBD) (Stanford Coronavirus Antiviral & Resistance Database, 2022) that directly contacts the angiotensin-converting enzyme 2 (ACE2) host receptor and is also a primary target for antibodies. Here, we studied the binding mode of the S-RBD domain of the Spike protein carrying the Omicron mutations and the globular domain of human ACE2 using molecular dynamics (MD) simulations. We identified new and key Omicron-specific interactions such as R493 (of mutation Q493R), which forms salt bridges both with E35 and D38 of ACE2, Y501 (N501Y), which forms an edge-to-face aromatic interaction with Y41, and Y505 (Y505H), which makes an H-bond with E37 and K353. The glycan chains of ACE2 also bind differently in the WT and Omicron variants in response to different charge distributions on the surface of Spike proteins. However, while the Omicron mutations considerably improve the overall electrostatic fit of the two interfaces, the total number of specific and favorable interactions between the two does not increase. The dynamics of the complexes are highly affected too, making the Omicron S-RBD:ACE2 complex more rigid; the two main interaction sites, Patches I and II, isolated in the WT complex, become connected in the Omicron complex through the alternating interaction of R493 and R498 with E35 and D38.

Directed Evolution-Driven Increase of Structural Plasticity Is a Prerequisite for Binding the Complement Lectin Pathway Blocking MASP-Inhibitor Peptides

MASP-1 and MASP-2 are key activator proteases of the complement lectin pathway. The first specific mannose-binding lectin-associated serine protease (MASP) inhibitors had been developed from the 14-amino-acid sunflower trypsin inhibitor (SFTI) peptide by phage display, yielding SFTI-based MASP inhibitors, SFMIs. Here, we present the crystal structure of the MASP-1/SFMI1 complex that we analyzed in comparison to other existing MASP-1/2 structures. Rigidified backbone structure has long been accepted as a structural prerequisite for peptide inhibitors of proteases. We found that a hydrophobic cluster organized around the P2 Thr residue is essential for the structural stability of wild-type SFTI. We also found that the same P2 Thr prevents binding of the rigid SFTI-like peptides to the substrate-binding cleft of both MASPs as the cleft is partially blocked by large gatekeeper enzyme loops. Directed evolution removed this obstacle by replacing the P2 Thr with a Ser, providing the SFMIs with high-degree structural plasticity, which proved to be essential for MASP inhibition. To gain more insight into the structural criteria for SFMI-based MASP-2 inhibition, we systematically modified MASP-2-specific SFMI2 by capping its two termini and by replacing its disulfide bridge with varying length thioether linkers. By doing so, we also aimed to generate a versatile scaffold that is resistant to reducing environment and has increased stability in exopeptidase-containing biological environments. We found that the reduction-resistant disulfide-substituted l-2,3-diaminopropionic acid (Dap) variant possessed near-native potency. As MASP-2 is involved in the life-threatening thrombosis in COVID-19 patients, our synthetic, selective MASP-2 inhibitors could be relevant coronavirus drug candidates.

Cryo-EM structure of acylpeptide hydrolase reveals substrate selection by multimerization and a multi-state serine-protease triad

The first structure of tetrameric mammalian acylaminoacyl peptidase, an enzyme that functions as an upstream regulator of the proteasome through the removal of terminal N-acetylated residues from its protein substrates, was determined by cryo-EM and further elucidated by MD simulations. Self-association results in a toroid-shaped quaternary structure, guided by an amyloidogenic β-edge and unique inserts. With a Pro introduced into its central β-sheet, sufficient conformational freedom is awarded to the segment containing the catalytic Ser587 that the serine protease catalytic triad alternates between active and latent states. Active site flexibility suggests that the dual function of catalysis and substrate selection are fulfilled by a novel mechanism: substrate entrance is regulated by flexible loops creating a double-gated channel system, while binding of the substrate to the active site is required for stabilization of the catalytic apparatus – as a second filter before hydrolysis. The structure not only underlines that within the family of S9 proteases homo-multimerization acts as a crucial tool for substrate selection, but it will also allow drug design targeting of the ubiquitin-proteasome system.

The structure of tetrameric mammalian acylaminoacyl peptidase – a key upstream regulator of the proteasome – was determined by cryo-EM (and elucidated by MD), showing a “shutters-and-channels” substrate selection apparatus created by oligomerization.

Dynamically encoded reactivity of Ras enzymes: opening new frontiers for drug discovery

Decoding molecular flexibility in order to understand and predict biological processes-applying the principles of dynamic-structure-activity relationships (DSAR)-becomes a necessity when attempting to design selective and specific inhibitors of a protein that has overlapping interaction surfaces with its upstream and downstream partners along its signaling cascade. Ras proteins are molecular switches that meet this definition perfectly. The close-lying P-loop and the highly flexible switch I and switch II regions are the site of nucleotide-, assisting-, and effector-protein binding. Oncogenic mutations that also appear in this region do not cause easily characterized overall structural changes, due partly to the inherent conformational heterogeneity and pliability of these segments. In this review, we present an overview of the results obtained using approaches targeting Ras dynamics, such as nuclear magnetic resonance (NMR) measurements and experiment-based modeling calculations (mostly molecular dynamics (MD) simulations). These methodologies were successfully used to decipher the mutant- and isoform-specific nature of certain transient states, far-lying allosteric sites, and the internal interaction networks, as well as the interconnectivity of the catalytic and membrane-binding regions. This opens new therapeutic potential: the discovered interaction hotspots present hitherto not targeted, selective sites for drug design efforts in diverse locations of the protein matrix.

Tumor-Suppressor p53TAD1-60 Forms a Fuzzy Complex with Metastasis-Associated S100A4: Structural Insights and Dynamics by an NMR/MD Approach

Conformationally flexible protein complexes represent a major challenge for structural and dynamical studies. We present herein a method based on a hybrid NMR/MD approach to characterize the complex formed between the disordered p53TAD1-60 and the metastasis-associated S100A4. Disorder-to-order transitions of both TAD1 and TAD2 subdomains upon interaction is detected. Still, p53TAD1-60 remains highly flexible in the bound form, with residues L26, M40, and W53 being anchored to identical hydrophobic pockets of the S100A4 monomer chains. In the resulting "fuzzy" complex, the clamp-like binding of p53TAD1-60 relies on specific hydrophobic anchors and on the existence of extended flexible segments. Our results demonstrate that structural and dynamical NMR parameters (cumulative Δδ, SSP, temperature coefficients, relaxation time, hetNOE) combined with MD simulations can be used to build a structural model even if, due to high flexibility, the classical solution structure calculation is not possible.

Configuration-Controlled Crystal and/or Gel Formation of Protected d-Glucosamines Supported by Promiscuous Interaction Surfaces and a Conformationally Heterogeneous Solution State

The configuration-dependent self-association mode of the two anomers of O-Ac,N-Fmoc-d-glucosamine, a foldamer building block, leading to gel and/or single crystal formation is described. The β-anomer of the sugar amino acid (2) forms a gel from various solvents (confirmed by SEM, rheology measurements, NMR, and ECD spectroscopy), whereas the α-anomer (1) does not form a gel with any solvent tested. Transition from the solution state to a gel is coupled to a concurrent shift of the Fmoc-groups: from a freely rotating (almost symmetrical) to a specific, asymmetric orientation. Whereas the crystal structure of the α-anomer is built as an evenly packed 3D system, the β-anomer forms a looser superstructure of well-packed 2D layers. Modeling indicates that in the lowest energy, but scarcely sampled conformer of the β-anomer, the Fmoc-group bends above the sugar moiety, stabilized by intramolecular CH↔π interactions between the aromatic rings. It is concluded that possessing an extended and promiscuous interaction surface and a conformationally heterogeneous solution state are among the basic requirements of gel formation for a candidate molecule.

Structural impact of GTP binding on downstream KRAS signaling

Oncogenic RAS proteins, involved in ∼30% of human tumors, are molecular switches of various signal transduction pathways. Here we apply a new protocol for the NMR study of KRAS in its (inactive) GDP- and (activated) GTP-bound form, allowing a comprehensive analysis of the backbone dynamics of its WT-, G12C- and G12D variants. We found that Tyr32 shows opposite mobility with respect to the backbone of its surroundings: it is more flexible in the GDP-bound form while more rigid in GTP-complexes (especially in WT- and G12D-GTP). Using the G12C/Y32F double mutant, we showed that the presence of the hydroxyl group of Tyr32 has a marked effect on the G12C-KRAS-GTP system as well. Molecular dynamics simulations indicate that Tyr32 is linked to the γ-phosphate of GTP in the activated states – an arrangement shown, using QM/MM calculations, to support catalysis. Anchoring Tyr32 to the γ-phosphate contributes to the capture of the catalytic waters participating in the intrinsic hydrolysis of GTP and supports a simultaneous triple proton transfer step (catalytic water → assisting water → Tyr32 → O1G of the γ-phosphate) leading to straightforward product formation. The coupled flip of negatively charged residues of switch I toward the inside of the effector binding pocket potentiates ligand recognition, while positioning of Thr35 to enter the coordination sphere of the Mg²⁺ widens the pocket. Position 12 mutations do not disturb the capture of Tyr32 by the γ-phosphate, but (partially) displace Gln61, which opens up the catalytic pocket and destabilizes catalytic water molecules thus impairing intrinsic hydrolysis.

Nucleotide exchange to the physiological, activated, GTP-bound form of KRAS results in the anchoring of Tyr32 within the active site.

Gas Sensing by Bacterial H-NOX Proteins: An MD Study

Gas sensing is crucial for both prokaryotes and eukaryotes and is primarily performed by heme-based sensors, including H-NOX domains. These systems may provide a new, alternative mode for transporting gaseous molecules in higher organisms, but for the development of such systems, a detailed understanding of the ligand-binding properties is required. Here, we focused on ligand migration within the protein matrix: we performed molecular dynamics simulations on three bacterial (Ka, Ns and Cs) H-NOX proteins and studied the kinetics of CO, NO and O₂ diffusion. We compared the response of the protein structure to the presence of ligands, diffusion rate constants, tunnel systems and storage pockets. We found that the rate constant for diffusion decreases in the O₂ > NO > CO order in all proteins, and in the Ns > Ks > Cs order if single-gas is considered. Competition between gases seems to seriously influence the residential time of ligands spent in the distal pocket. The channel system is profoundly determined by the overall fold, but the sidechain pattern has a significant role in blocking certain channels by hydrophobic interactions between bulky groups, cation-π interactions or hydrogen bonding triads. The majority of storage pockets are determined by local sidechain composition, although certain functional cavities, such as the distal and proximal pockets are found in all systems. A major guideline for the design of gas transport systems is the need to chemically bind the gas molecule to the protein, possibly joining several proteins with several heme groups together.

Tissue-Specific Requirement for the GINS Complex During Zebrafish Development

Efficient and accurate DNA replication is particularly critical in stem and progenitor cells for successful proliferation and survival. The replisome, an amalgam of protein complexes, is responsible for binding potential origins of replication, unwinding the double helix, and then synthesizing complimentary strands of DNA. According to current models, the initial steps of DNA unwinding and opening are facilitated by the CMG complex, which is composed of a GINS heterotetramer that connects Cdc45 with the mini-chromosome maintenance (Mcm) helicase. In this work, we provide evidence that in the absence of GINS function DNA replication is cell autonomously impaired, and we also show that gins1 and gins2 mutants exhibit elevated levels of apoptosis restricted to actively proliferating regions of the central nervous system (CNS). Intriguingly, our results also suggest that the rapid cell cycles during early embryonic development in zebrafish may not require the function of the canonical GINS complex as neither zygotic Gins1 nor Gins2 isoforms seem to be present during these stages.

Protein Aggregation in a Nutshell: The Splendid Molecular Architecture of the Dreaded Amyloid Fibrils

The recent high-resolution structures of amyloid fibrils show that the organization of peptide segments into amyloid aggregate architecture is a general process, though the morphology is more complex and intricate than suspected previously. The amyloid fibrils are often cytotoxic, accumulating as intracellular inclusions or extracellular plaques and have the ability to interfere with cellular physiology causing various cellular malfunctions. At the same time, the highly ordered amyloid structures also present an opportunity for nature to store and protect peptide chains under extreme conditions - something that might be used for designing storage, formulation, and delivery of protein medications or for contriving bio-similar materials of great resistance or structure-ordering capacity. Here we summarize amyloid characteristics; discussing the basic morphologies, sequential requirements and 3D-structure that are required for the understanding of this newly (re)discovered protein structure - a prerequisite for developing either inhibitors or promoters of amyloid-forming processes.

The Route from the Folded to the Amyloid State: Exploring the Potential Energy Surface of a Drug-Like Miniprotein

Invited for the cover of this issue is the group of András Perczel at Eötvös Loránd University, Budapest, Hungary and colleagues from Osaka University, Japan. The image depicts the amyloid buildup of an Exenatide derivate miniprotein (E5) monitored on a simplified hyperspace. Read the full text of the article at 10.1002/chem.201903826.

The Route from the Folded to the Amyloid State: Exploring the Potential Energy Surface of a Drug-Like Miniprotein

The amyloid formation of the folded segment of a variant of Exenatide (a marketed drug for type‐2 diabetes mellitus) was studied by electronic circular dichroism (ECD) and NMR spectroscopy. We found that the optimum temperature for E5 protein amyloidosis coincides with body temperature and requires well below physiological salt concentration. Decomposition of the ECD spectra and its barycentric representation on the folded‐unfolded‐amyloid potential energy surface allowed us to monitor the full range of molecular transformation of amyloidogenesis. We identified points of no return (e.g.; T=37 °C, pH 4.1, c_E5=250 μm, c_NaCl=50 mm, t>4–6 h) that will inevitably gravitate into the amyloid state. The strong B‐type far ultraviolet (FUV)‐ECD spectra and an unexpectedly strong near ultraviolet (NUV)‐ECD signal (Θ_{≈275–285  nm}) indicate that the amyloid phase of E5 is built from monomers of quasi‐elongated backbone structure (φ≈−145°, ψ≈+145°) with strong interstrand Tyr↔Trp interaction. Misfolded intermediates and the buildup of “toxic” early‐stage oligomers leading to self‐association were identified and monitored as a function of time. Results indicate that the amyloid transition is triggered by subtle misfolding of the α‐helix, exposing aromatic and hydrophobic side chains that may provide the first centers for an intermolecular reorganization. These initial clusters provide the spatial closeness and sufficient time for a transition to the β‐structured amyloid nucleus, thus the process follows a nucleated growth mechanism.

Assignment of Vibrational Circular Dichroism Cross-Referenced Electronic Circular Dichroism Spectra of Flexible Foldamer Building Blocks: Towards Assigning Pure Chiroptical Properties of Foldamers

Assignment of the most established electronic circular dichroism (ECD) spectra of polypeptides and foldamers is either “evidence based” or relies on the 3D structures of longer oligomers of limited internal dynamics, which are derived from NMR spectroscopy (or X‐ray) data. Critics warn that the use of NMR spectroscopy and ECD side by side has severe limitations for flexible molecules because explicit knowledge of conformational ensembles is a challenge. Herein, an old–new method of comparing ab initio computed and measured vibrational circular dichroism (VCD) data is presented to validate both the structures (conf(i)) and their relative weights (c(i)) that make up the conformational ensemble. Based on the array of {conf(i), c(i)}, the pure ECD spectra, g(i)^conf(i), can be ab initio calculated. The reconstructed spectrum Σc(i)g(i)^conf(i) can thus help to assign any experimental ECD counterparts. Herein, such a protocol is successfully applied to flexible foldamer building blocks of sugar β‐amino acid diamides. The epimeric pair of the model system was selected because these molecules were conformationally tunable by simple chemical modification, and thus, the robustness of the current approach could be probed. The initial hydrogen bond (NH⋅⋅⋅O) eliminated by N‐methylation reorients the amide plain, which influences the chiroptical properties of the foldamer building block; this structural change is successfully monitored by changes to the VCD and ECD transitions, which are now assigned to pure conformers. The current method seems to be general and effective without requiring extensive CPU and spectroscopic resources.

Achieving Functionality Through Modular Build-up: Structure and Size Selection of Serine Oligopeptidases

Enzymes of the prolyl oligopeptidase family (S9 family) recognize their substrates not only by the specificity motif to be cleaved but also by size - they hydrolyze oligopeptides smaller than 30 amino acids. They belong to the serine-protease family, but differ from classical serine-proteases in size (80 kDa), structure (two domains) and regulation system (size selection of substrates). This group of enzymes is an important target for drug design as they are linked to amnesia, schizophrenia, type 2 diabetes, trypanosomiasis, periodontitis and cell growth. By comparing the structure of various members of the family we show that the most important features contributing to selectivity and efficiency are: (i) whether the interactions weaving the two domains together play a role in stabilizing the catalytic triad and thus their absence may provide for its deactivation: these oligopeptidases can screen their substrates by opening up, and (ii) whether the interaction-prone β-edge of the hydrolase domain is accessible and thus can guide a multimerization process that creates shielded entrance or intricate inner channels for the size-based selection of substrates. These cornerstones can be used to estimate the multimeric state and selection strategy of yet undetermined structures.

Hydration shell differentiates folded and disordered states of a Trp-cage miniprotein, allowing characterization of structural heterogeneity by wide-line NMR measurements

Hydration properties of folded and unfolded/disordered miniproteins were monitored in frozen solutions by wide-line 1H-NMR. The amount of mobile water as function of T (-80 °C < T < 0 °C) was found characteristically different for folded (TC5b), semi-folded (pH < 3, TCb5(H+)) and disordered (TC5b_N1R) variants. Comparing results of wide-line 1H-NMR and molecular dynamics simulations we found that both the amount of mobile water surrounding proteins in ice, as well as their thaw profiles differs significantly as function of the compactness and conformational heterogeneity of their structure. We found that (i) at around -50 °C ~50 H2Os/protein melt (ii) if the protein is well-folded then this amount of mobile water remains quasi-constant up to -20 °C, (iii) if disordered then the quantity of the lubricating mobile water increases with T in a constant manner up to ~200 H2Os/protein by reaching -20 °C. Especially in the -55 °C ↔ -15 °C temperature range, wide-line 1H-NMR detects the heterogeneity of protein fold, providing the size of the hydration shell surrounding the accessible conformers at a given temperature. Results indicate that freezing of protein solutions proceeds by the gradual selection of the enthalpically most favored states that also minimize the number of bridging waters.

The mutation-dependent pathogenicity of NPHS2 p.R229Q: A guide for clinical assessment

NPHS2, encoding podocin, is the major gene implicated in steroid-resistant nephrotic syndrome. Its c.686G>A, p.R229Q variant is the first human variant with a mutation-dependent pathogenicity; it is only pathogenic when trans-associated to specific mutations. Secondary to its high allele frequency in the European, South Asian, African, and Latino populations, its benign trans-associations can be accidentally identified in affected patients. Distinguishing pathogenic and benign p.R229Q associations can be challenging. In this paper, we present the currently known pathogenic and benign associations, and show that a rare p.R229Q association can be considered pathogenic if the variant in trans meets the following criteria; it affects the 270-351 residues and alters but does not disrupt the oligomerization, its p.R229Q association is found in a family with slowly progressing focal segmental glomerulosclerosis, but is expected to be rare in the general population (<1:10⁶ ). We show that >15% of the p.R229Q associations identified so far in patients are benign.

C-terminal oligomerization of podocin mediates interallelic interactions

Interallelic interactions of membrane proteins are not taken into account while evaluating the pathogenicity of sequence variants in autosomal recessive disorders. Podocin, a membrane-anchored component of the slit diaphragm, is encoded by NPHS2, the major gene mutated in hereditary podocytopathies. We formerly showed that its R229Q variant is only pathogenic when trans-associated to specific 3' mutations and suggested the causal role of an abnormal C-terminal dimerization. Here we show by FRET analysis and size exclusion chromatography that podocin oligomerization occurs exclusively through the C-terminal tail (residues 283-382): principally through the first C-terminal helical region (H1, 283-313), which forms a coiled coil as shown by circular dichroism spectroscopy, and through the 332-348 region. We show the principal role of the oligomerization sites in mediating interallelic interactions: while the monomer-forming R286Tfs*17 podocin remains membranous irrespective of the coexpressed podocin variant identity, podocin variants with an intact H1 significantly influence each other's localization (r² = 0.68, P = 9.2 × 10^-32). The dominant negative effect resulting in intracellular retention of the pathogenic F344Lfs*4-R229Q heterooligomer occurs in parallel with a reduction in the FRET efficiency, suggesting the causal role of a conformational rearrangement. On the other hand, oligomerization can also promote the membrane localization: it can prevent the endocytosis of F344Lfs*4 or F344* podocin mutants induced by C-terminal truncation. In conclusion, C-terminal oligomerization of podocin can mediate both a dominant negative effect and interallelic complementation. Interallelic interactions of NPHS2 are not restricted to the R229Q variant and have to be considered in compound heterozygous individuals.

First Principles Calculation of the Reaction Rates for Ligand Binding to Myoglobin: The Cases of NO and CO

Ligand binding by proteins is among the most fundamental processes in nature. Among these processes the binding of small gas molecules, such as O₂ , CO and NO to heme proteins has traditionally received vivid interest, which was further boosted by their recently recognized significant role in gas sensing in the body. At the heart of the binding of these ligands to the heme group is the spinforbidden reaction between high-spin iron(II) and the ligand yielding a low-spin adduct. We use computational means to address the complete mechanism of CO and NO binding by myoglobin. Considering that it involves several steps occurring on different time scales, molecular dynamics simulations were performed to address the diffusion of the ligand through the enzyme, and DFT calculations in combination with statistical rate calculation to investigate the spin-forbidden reaction. The calculations yielded rate constants in qualitative agreement with experiments and revealed that the bottleneck of NO and CO binding is different; for NO, diffusion was found to be rate-limiting, whereas for CO, the spin-forbidden step is the slowest.

Predictable Conformational Diversity in Foldamers of Sugar Amino Acids

A systematic conformational search was carried out for monomers and homohexamers of furanoid β-amino acids: cis-(S,R) and trans-(S,S) stereoisomers of aminocyclopentane carboxylic acid (ACPC), two different aminofuranuronic acids (AFU^α and AFU^β), their isopropylidene derivatives (AFU(ip)), and the key intermediate β-aminotetrahydrofurancarboxylic acid (ATFC). The stereochemistry of the building blocks was chosen to match that of the natural sugar amino acid (xylose and ribose) precursors (XylAFU and RibAFU). The results show that hexamers of cis-furanoid β-amino acids show great variability: while hydrophobic cyclopentane (cis-ACPC)₆ and hydrophilic (XylAFU^α/β)₆ foldamers favor two different zigzagged conformation as hexamers, the backbone fold turns into a helix in the case of (cis-ATFC)₆ (10-helix) and (XylAFU(ip))₆ (14-helix). Trans stereochemistry resulted in hexamers exclusively with the right-handed helix conformation, (H₁₂^P)₆, regardless of their polarity. We found that the preferred oligomeric structure of XylAFU^α/β is conformationally compatible with β-pleated sheets, while that of the trans/(S,S) units matches with α-helices of proteins.

Phosphorylation as conformational switch from the native to amyloid state: Trp-cage as a protein aggregation model

The 20 residue long Trp-cage miniprotein is an excellent model for both computational and experimental studies of protein folding and stability. Recently, great attention emerged to study disease-related protein misfolding, aggregation, and amyloid formation, with the aim of revealing their structural and thermodynamic background. Trp-cage is sensitive to both environmental and structure-modifying effects. It aggregates with ease upon structure destabilization, and thus it is suitable for modeling aggregation and amyloid formation. Here, we characterize the amyloid formation of several sequence modified and side-chain phosphorylated Trp-cage variants. We applied NMR, circular dichroism (CD) and Fourier transform infrared (FTIR) spectroscopies, molecular dynamics (MD) simulations, and transmission electron microscopy (TEM) in conjunction with thioflavin-T (ThT) fluorescence measurements to reveal the structural consequences of side-chain phosphorylation. We demonstrate that the native fold is destabilized upon serine phosphorylation, and the resultant highly dynamic structures form amyloid-like ordered aggregates with high intermolecular β-structure content. The only exception is the D9S(P) variant, which follows an alternative aggregation process by forming thin fibrils, presenting a CD spectrum of PPII helix, and showing low ThT binding capability. We propose a complex aggregation model for these Trp-cage miniproteins. This model assumes an additional aggregated state, a collagen triple helical form that can precede amyloid formation. The phosphorylation of a single serine residue serves as a conformational switch, triggering aggregation, otherwise mediated by kinases in cell. We show that Trp-cage miniprotein is indeed a realistic model of larger globular systems of composite folding and aggregation landscapes and helps us to understand the fundamentals of deleterious protein aggregation and amyloid formation.

A self-compartmentalizing hexamer serine protease from Pyrococcus horikoshii: substrate selection achieved through multimerization

Oligopeptidases impose a size limitation on their substrates, the mechanism of which has long been under debate. Here we present the structure of a hexameric serine protease, an oligopeptidase from Pyrococcus horikoshii (PhAAP), revealing a complex, self-compartmentalized inner space, where substrates may access the monomer active sites passing through a double-gated "check-in" system, first passing through a pore on the hexamer surface and then turning to enter through an even smaller opening at the monomers' domain interface. This substrate screening strategy is unique within the family. We found that among oligopeptidases, a residue of the catalytic apparatus is positioned near an amylogenic β-edge, which needs to be protected to prevent aggregation, and we found that different oligopeptidases use different strategies to achieve such an end. We propose that self-assembly within the family results in characteristically different substrate selection mechanisms coupled to different multimerization states.

Direct hydride shift mechanism and stereoselectivity of P450nor confirmed by QM/MM calculations

Nitric oxide reductase (P450(nor)) found in Fusarium oxysporum catalyzes the reduction of nitric oxide to N(2)O in a multistep process. The reducing agent, NADH, is bound in the distal pocket of the enzyme, and direct hydride transfer occurs from NADH to the nitric oxide bound heme enzyme, forming intermediate I. Here we studied the possibility of hydride transfer from NADH to both the nitrogen and oxygen of the heme-bound nitric oxide, using quantum chemical and combined quantum mechanics/molecular mechanics (QM/MM) calculations, on two different protein models, representing both possible stereochemistries, a syn- and an anti-NADH arrangement. All calculations clearly favor hydride transfer to the nitrogen of nitric oxide, and the QM-only barrier and kinetic isotope effects are good agreement with the experimental values of intermediate I formation. We obtained higher barriers in the QM/MM calculations for both pathways, but hydride transfer to the nitrogen of nitric oxide is still clearly favored. The barriers obtained for the syn, Pro-R conformation of NADH are lower and show significantly less variation than the barriers obtained in the case of anti conformation. The effect of basis set and wide range of functionals on the obtained results are also discussed.

Chemotaxis induced by SXWS tetrapeptides in Tetrahymena--overlapping chemotactic effects of SXWS sequences and their identical amino acids

The chemotactic potential of SXWS peptides and the components of the extracellular domain of cytokine receptors were investigated in Tetrahymena as a functional index of substitution with different amino acids in the position 'X' of the tetrapeptide. Data obtained demonstrate that position X plays a special determining role in the ligand, SEWS and STWS possess extremely strong chemoattractant ability, and aromatic amino acids result in chemorepellent ligands. Diverse effects of structurally related molecules, for example, SNWS-SDWS, demonstrate a highly sensitive discrimination potential in the applied model system. Physicochemical characteristics (hydropathy, residue size, and solvent-exposed area) of the amino acids were correlated with the chemotactic activity. Data obtained by computer-assisted conformation analysis of SXWS peptides and the highly overlapping chemotactic effects of the investigated SXWS peptides as well as the presence of the amino acids in the 'X' position indicate that member 'X' of the SXWS sequence performs a special role in interactions with the chemotaxis receptors in the membrane.

Structure and catalysis of acylaminoacyl peptidase: closed and open subunits of a dimer oligopeptidase

Acylaminoacyl peptidase from Aeropyrum pernix is a homodimer that belongs to the prolyl oligopeptidase family. The monomer subunit is composed of one hydrolase and one propeller domain. Previous crystal structure determinations revealed that the propeller domain obstructed the access of substrate to the active site of both subunits. Here we investigated the structure and the kinetics of two mutant enzymes in which the aspartic acid of the catalytic triad was changed to alanine or asparagine. Using different substrates, we have determined the pH dependence of specificity rate constants, the rate-limiting step of catalysis, and the binding of substrates and inhibitors. The catalysis considerably depended both on the kind of mutation and on the nature of the substrate. The results were interpreted in terms of alterations in the position of the catalytic histidine side chain as demonstrated with crystal structure determination of the native and two mutant structures (D524N and D524A). Unexpectedly, in the homodimeric structures, only one subunit displayed the closed form of the enzyme. The other subunit exhibited an open gate to the catalytic site, thus revealing the structural basis that controls the oligopeptidase activity. The open form of the native enzyme displayed the catalytic triad in a distorted, inactive state. The mutations affected the closed, active form of the enzyme, disrupting its catalytic triad. We concluded that the two forms are at equilibrium and the substrates bind by the conformational selection mechanism.

Comparative computational analysis of active and inactive cofactors of nitric oxide synthase

Nitric oxide synthases (NOSs) are heme proteins that catalyze the formation of nitric oxide from L-Arg in the presence of oxygen. Of the two electrons required for the first step of the reaction, the second is primarily donated by the tetrahydrobiopterin (H4B) cofactor bound adjacent to the heme, which is eventually reduced back to resting state by the ultimate electron source of the reaction, the flavins of the NOS reductase domain. Density functional theory calculations were carried out to identify those protonation states of different cofactor molecules that best support radicalization of the cofactor and the coupled increase in the electron density of the heme-bound oxygen molecule. Three cofactor molecules were studied, native H4B, an active analogue, 5-methyl-H4B, and the inactive 4-amino-H4B. Findings support the emerging model where H4B and 5-methyl-H4B are coupled proton/electron sources of NOS catalysis, while 4-amino-H4B is an inhibitor due to its inability to donate the catalytically required proton.

Interaction of antagonists with calmodulin: insights from molecular dynamics simulations

We report results of 12 ns, all-atom molecular dynamics simulation (MDS) and Poisson-Boltzmann free energy calculations (PBFE) on calmodulin (CaM) bound to two molecules of trifluoperazine (TFP) and of N-(3,3, diphenylpropyl)- N'-[1- R-(3,4-bis-butoxyphenyl)-ethyl]-propylenediamine (DPD). X-ray data show very similar structures for the two complexes, yet the antagonists significantly differ with respect to their CaM binding affinities, the neutral DPD is much more potent. The goal of the study was to unravel the reason why TFP is less potent although its positive charge should facilitate binding. The electrostatic energy terms in CHARMM and binding free energy terms of the PBFE approach showed TFP a better antagonist, while inspection of hydrophobic contacts supports DPD binding. Detailed inspection of the amino acid contributions of PBFE calculations unravel that steric reasons oppose the favorable binding of TFP. Structural conditions are given for a successful drug design strategy, which may benefit also from charge-charge interactions.

Binding mode analysis of the NADH cofactor in nitric oxide reductase: a theoretical study

Nitric oxide reductase (Nor), a member of the cytochrome P450 superfamily, takes part in the denitrification process of fungi by reducing NO to N(2)O. Evidence indicates that Nor binds NADH, source of the reducing equivalents of the reaction, within its large hydrophilic ligand binding cavity on the distal side of heme and receives electrons directly from the cofactor. Here we present a binding mode analysis of the structure of the Nor-NO-NADH complex, performed in three steps. The NADH cofactor was first docked into the enzyme interior using the Monte Carlo multiple minimum algorithm, refined by low-mode conformational search and the final arrangement was obtained in a 5ns NPT molecular dynamics simulation. The NADH cofactor, in our results, is positioned--by Arg174, Lys291, Asp393 and several water molecules--within reactive distance of the NO binding spot suggesting a direct hydride shift mechanism between the two. The catalytically required water molecule is captured by NADH and the cofactor not only retains the suggested H-bonded proton transfer pathway between the active site and the solvent, but provides structural restraint for its members. We also found that direct interaction is formed between the cofactor and propionate A of the heme group, which flips from the proximal to the distal side of the heme plane in order to become an H-bonding partner of NADH. The role of Arg64 and Glu71 was suggested to be fixing the residues of the translocated helix B' to their new position.

Protonation state of Asp30 exerts crucial influence over surface loop rearrangements responsible for NO release in nitrophorin 4

pK(a) values of ionizable residues were calculated for the crystal structures describing the pH and NO binding dependant conformations of nitrophorin 4, a pH sensitive NO carrier heme protein. Comparison of resultant H-bonding patterns allowed the identification of the amino acids that take part in signaling pH change. We carried out MD simulations to show that the protonation state of Asp30, buried in the closed conformation, is crucial for maintaining the tight packed conformation of the closed form of the complex - presenting a model for the functional decrease of NO binding affinity of nitrophorins at physiological pH.

Role of proximal His93 in nitric oxide binding to metmyoglobin. Application of continuum solvation in Monte Carlo protein simulations

Monte Carlo protein simulations with continuum solvation were used to explore the conformational mobility of NO within the active site of metmyoglobin. To the best of our knowledge this is the first application of a continuum solvation model for exploring protein binding sites. The usefulness of the Monte Carlo conformational analysis was demonstrated in comparative molecular dynamics simulations. Analysis of conformer populations revealed that Monte Carlo conformational analysis is more effective in mapping the relevant region of the potential surface. Distribution of low-energy conformers obtained for the metmyoglobin-NO complex was found to depend on the orientation of proximal His93. Different charge distributions corresponding to the two experimentally verified possible torsions of this proximal residue result in strong binding of NO or its release to a nearby hydrophobic trap. Conformer populations obtained by Monte Carlo conformational analysis were grouped into three main families: one, with the NO directly bound to the iron, appears when the CA-CB-CG-CD2 torsion of His93 is at its ligand binding value (-113 degrees); and two conformers exist where NO is trapped in a nearby hydrophobic pocket, the same cavity that was determined to be the geminate trap of CO in ferrous Mb as a result of the torsional flip of His93 to its ligand releasing state (-125 degrees). Based on this analysis, we suggest that the electrostatic rearrangement coupled to the conformational fluctuation of the proximal His leads to the geminate trapping of the ligand. Conformational rearrangement of the proximal side would provide the possibility of rebinding of the ligand to Fe.