Elucidation of the cause for reduced activity of abnormal human plasmin containing an Ala55-Thr mutation: importance of highly conserved Ala55 in serine proteases

Plasminogen missense variants and their involvement in cardiovascular and inflammatory disease

Human plasminogen (PLG), the zymogen of the fibrinolytic protease, plasmin, is a polymorphic protein with two widely distributed codominant alleles, PLG/Asp453 and PLG/Asn453. About 15 other missense or non-synonymous single nucleotide polymorphisms (nsSNPs) of PLG show major, yet different, relative abundances in world populations. Although the existence of these relatively abundant allelic variants is generally acknowledged, they are often overlooked or assumed to be non-pathogenic. In fact, at least half of those major variants are classified as having conflicting pathogenicity, and it is unclear if they contribute to different molecular phenotypes. From those, PLG/K19E and PLG/A601T are examples of two relatively abundant PLG variants that have been associated with PLG deficiencies (PD), but their pathogenic mechanisms are unclear. On the other hand, approximately 50 rare and ultra-rare PLG missense variants have been reported to cause PD as homozygous or compound heterozygous variants, often leading to a debilitating disease known as ligneous conjunctivitis. The true abundance of PD-associated nsSNPs is unknown since they can remain undetected in heterozygous carriers. However, PD variants may also contribute to other diseases. Recently, the ultra-rare autosomal dominant PLG/K311E has been found to be causative of hereditary angioedema (HAE) with normal C1 inhibitor. Two other rare pathogenic PLG missense variants, PLG/R153G and PLG/V709E, appear to affect platelet function and lead to HAE, respectively. Herein, PLG missense variants that are abundant and/or clinically relevant due to association with disease are examined along with their world distribution. Proposed molecular mechanisms are discussed when known or can be reasonably assumed.

Statistical Learning of Protein Elastic Network from Positional Covariance Matrix

Positional fluctuation and covariance during protein dynamics are key observables for understanding the molecular origin of biological functions. A frequently employed potential energy function for describing protein structural variation at the coarse-gained level is elastic network model (ENM). A long-standing issue in biomolecular simulation is thus the parametrization of ENM spring constants from the components of positional covariance matrix (PCM). Based on sensitivity analysis of PCM, the direct-coupling statistics of each spring, which is a specific combination of position fluctuation and covariance, is found to exhibit prominent signal of parameter dependence. This finding provides the basis for devising the objective function and the scheme of running through the effective one-dimensional optimization of every spring by self-consistent iteration. Formal derivation of the positional covariance statistical learning (PCSL) method also motivates the necessary data regularization for stable calculations. Robust convergence of PCSL is achieved in taking an all-atom molecular dynamics trajectory or an ensemble of homologous structures as input data. The PCSL framework can also be generalized with mixed objective functions to capture specific property such as the residue flexibility profile. Such physical chemistry-based statistical learning thus provides a useful platform for integrating the mechanical information encoded in various experimental or computational data.

NBCZone: Universal three-dimensional construction of eleven amino acids near the catalytic nucleophile and base in the superfamily of (chymo)trypsin-like serine fold proteases

(Chymo)trypsin-like serine fold proteases belong to the serine/cysteine proteases found in eukaryotes, prokaryotes, and viruses. Their catalytic activity is carried out using a triad of amino acids, a nucleophile, a base, and an acid. For this superfamily of proteases, we propose the existence of a universal 3D structure comprising 11 amino acids near the catalytic nucleophile and base - Nucleophile-Base Catalytic Zone (NBCZone). The comparison of NBCZones among 169 eukaryotic, prokaryotic, and viral (chymo)trypsin-like proteases suggested the existence of 15 distinct groups determined by the combination of amino acids located at two "key" structure-functional positions 54_T and 55_T near the catalytic base His57_T. Most eukaryotic and prokaryotic proteases fell into two major groups, [ST]A and TN. Usually, proteases of [ST]A group contain a disulfide bond between cysteines Cys42_T and Cys58_T of the NBCZone. In contrast, viral proteases were distributed among seven groups, and lack this disulfide bond. Furthermore, only the [ST]A group of eukaryotic proteases contains glycine at position 43_T, which is instrumental for activation of these enzymes. In contrast, due to the side chains of residues at position 43_T prokaryotic and viral proteases do not have the ability to carry out the structural transition of the eukaryotic zymogen-zyme type.

Inhibition of bacterial undecaprenyl pyrophosphate synthase by small fungal molecules

Viridicatumtoxin and spirohexaline, small fungal molecules with a tetracyclic scaffold and an additional spirobicyclic ring in common, were found to inhibit bacterial undecaprenyl pyrophosphate (UPP) synthase with IC₅₀ values of 4 and 9 μm, respectively. These molecules showed weak inhibitory activity against catalytically related enzymes such as bacterial octaprenyl pyrophosphate synthase and yeast dehydrodolichyl pyrophosphate synthase, indicating that the compounds preferentially inhibit UPP synthase. They showed antimicrobial activity, particularly against Gram-positive bacteria including methicillin-resistant Staphylococcus aureus (MRSA). Furthermore, molecular modeling strongly suggested that the hydrophobic spirobicyclic ring of viridicatumtoxin interacts with three hydrophobic clefts of the active site in MRSA UPP synthase.

[Mode of action of microbial anti-MRSA agents]

Methicillin-resistant Staphylococcus aureus (MRSA) is known as a major nosocomial pathogen that has also developed resistance to many antibiotics. Moreover, MRSA resistance to a last-resort antibiotic, vancomycin, has been reported. Therefore, new anti-infectious agents to prevent and treat MRSA infection are needed. Based on this background, our group has focused on the discovery of new microbial agents active against MRSA infection. Viridicatumtoxin and spirohexaline, produced by Penicillium sp. FKI-3368, were isolated as inhibitors of undecaprenyl pyrophosphate (UPP) synthase of Staphylococcus aureus, which was involved in cell wall synthesis. Viridicatumtoxin and spirohexaline with a pentacyclic spiro skeleton inhibited UPP synthase activity with an IC(50) value of 4.0 and 9.0 µM, respectively. Actually, the growth of gram-positive bacteria including MRSA was strongly inhibited by the compounds. Our computational modeling experiments indicated that spirohexaline A was inserted into the substrate pocket of UPP synthase and interacted with Glu(88) via a carbamoyl group of the compound, with Ala(76), Met(54) and Asn(35) via three hydroxyl groups, and with certain hydrophobic amino acids via a spiro ring. Cyslabdan, produced by Streptomyces sp. K04-0144, was isolated as a potentiator of β-lactam imipenem activity against MRSA. The compound consisted of a labdan skeleton and an N-acetylcysteine. Cyslabdan potentiated imipenem activity by over 1000 fold, drastically reducing the MIC value of imipenem against MRSA from 16 to 0.03 µg/mL. The binding proteins of cyslabdan were investigated in the lysate of MRSA to identify FemA, which was involved in the formation of the pentaglycine interpeptide bridge in MRSA peptidoglycan.

Protein structure prediction in CASP6 using CHIMERA and FAMS

In CASP6, the CHIMERA-group predicted full-atom models of all targets using SKE-CHIMERA, a Web-user interface system for protein structure prediction that allows human intervention at necessary stages; we used a lot of information from our own data and from publicly available data. Using SKE-CHIMERA, we iterated manual step (template selection and alignment by the in-house program CHIMERA) and automatic step (three-dimensional model building by the in-house program FAMS). The official CASP6 assessment showed that CHIMERA-group was one of the most successful predictors in homology modeling, especially for FR/H (Fold Recognition/Homologous). In this article, we introduce the method of CHIMERA-group and discuss its successes and failures in CASP6.

Evaluation of homology modeling of the severe acute respiratory syndrome (SARS) coronavirus main protease for structure based drug design

To accelerate the development of drugs against severe acute respiratory syndrome (SARS), we constructed a homology model of the SARS coronavirus main protease using our modeling software, FAMS Ligand&Complex, and released it before the X-ray structure was solved. The X-ray structure showed our model as accurately predicted and useful for structure based drug design.

Remarkably different structures and reaction mechanisms of ketoreductases for the opposite stereochemical control in the biosynthesis of BIQ antibiotics

Two ketoreductases, RED1 and RED2, are involved in the biosynthesis of actinorhodin in Streptomyces coelicolor A3(2) and dihydrogranaticin in S. violaceoruber Tu22, respectively. They are responsible for the stereospecific reductions of the bicyclic intermediate to give (S)- or (R)-DNPA, although there is no similarity between their amino acid sequences. Biotransformation using synthetic analogous substrates revealed that the substrate specificities are quite different. Homology modelling studies and site directed mutagenesis showed remarkable differences in three-dimensional structures and catalytic mechanisms between RED1 and RED2.

Population-based distribution of plasminogen activity and estimated prevalence and relevance to thrombotic diseases of plasminogen deficiency in the Japanese: the Suita Study

Reduced plasminogen activity with a normal level of antigen is commonly observed in Japanese individuals. The first reported patient with plasminogen deficiency was accompanied with deep vein thrombosis. The present study examines whether heterozygous or homozygous deficiency of plasminogen is a risk factor for thrombotic disease. This study measures the plasminogen activity of 4517 individuals in the general population, determines the cut-off to define plasminogen deficiency, and identifies plasminogen deficiencies in the control groups and thrombotic disease groups. In another study, we examined the phenotypes of consecutive patients with homozygous plasminogen deficiency detected in our hospital. We found 173 and two of 4517 individuals to have heterozygous and homozygous deficiency with normal plasminogen antigen level, respectively, and 19 to have heterozygous deficiency with reduced antigen levels. The incidence of plasminogen deficiency in an age- and sex-matched control group (13/324, 4.01% for deep vein thrombosis or 13/330, 3.94% for stroke) selected from the 4517 individuals was not significantly different from those in patients with deep vein thrombosis (3/108, 2.78%) or cardioembolic stroke (6/110, 5.55%). Among 19 patients with homozygous plasminogen deficiency showing about 10% plasminogen activity, none had deep vein thrombosis. These findings indicate that neither heterozygous nor homozygous plasminogen deficiency constitutes a significant risk factor for thrombotic disease.

An automatic homology modeling method consisting of database searches and simulated annealing

We introduce a method of homology modeling consisting of database searches and simulated annealing. All processes involving searches for homologous proteins, alignment, the construction of Calpha atoms, construction of main-chain atoms, and the construction of side-chain atoms are performed automatically. In this method, main-chain conformations are generated from the weighted average of mainchain coordinates in reference proteins. The weight is defined by the local space homology representing the similarity of environmental residues at topologically equivalent positions in reference proteins. Side-chain conformations are generated for constructed main-chain atoms by database searches, and main-chain atoms are optimized for the fixed side-chain conformations. These two processes, i.e., the side-chain generation and main-chain optimization, are repeated several times. This type of construction provides a structure similar to the x-ray structure, in particular, for main-chain and side-chain atoms in the residues belonging to structurally conserved regions (SCRs). The accuracy of our method was evaluated for 14 proteins whose structures are known. The average root mean square deviation between models and x-ray structures was 2.29 A for all atoms, and the percentage of chi1 angles within 30 degrees was 72.6% for SCRs residues. Some models were in good agreement with their respective x-ray structures. Our method, which has the advantage of being automated, gives results similar to, or better than, published results for three widely used test proteins. Our software, FAMS, is available on the World Wide Web.

A homology modeling method of an icosahedral viral capsid: inclusion of surrounding protein structures

A methodological development is presented for homology modeling of an icosahedrally symmetric assembly of proteins. In the method, a main-chain structure of an asymmetric unit of a protein assembly is constructed and structure refinement is performed, taking the surrounding symmetry-related proteins into consideration with rotational symmetry boundary conditions. To test the procedure, three models of a poliovirus capsid were constructed with different modeling conditions based on the X-ray structure of a rhinovirus capsid. Model S and model N were constructed with and without considering surrounding proteins, respectively. Model N2 was obtained by refinement in rotational symmetry boundary conditions of the structure of model N. The three models were compared with the X-ray structure of a poliovirus capsid. Root mean square deviations and C alpha distances indicate that model S is the most accurate. Examination of the intermolecular short contacts indicates that model S and model N2 are superior to model N, because they do not make severe intermolecular short contacts. Symmetric intermolecular interactions are important for both the structural fragment search and energy minimization to predict better loop structures. The programs developed in this study are thus valuable in homology modeling of an icosahedral viral capsid.