Genetic Spectrum of Familial Hypercholesterolaemia in the Malaysian Community: Identification of Pathogenic Gene Variants Using Targeted Next-Generation Sequencing.

Familial hypercholesterolaemia (FH) is caused by mutations in lipid metabolism genes, predominantly in low-density lipoprotein receptor (LDLR), apolipoprotein B (APOB), proprotein convertase subtilisin/kexin-type 9 (PCSK9) and LDL receptor adaptor protein 1 (LDLRAP1). The prevalence of genetically confirmed FH and the detection rate of pathogenic variants (PV) amongst clinically diagnosed patients is not well established. Targeted next-generation sequencing of LDLR, APOB, PCSK9 and LDLRAP1 was performed on 372 clinically diagnosed Malaysian FH subjects. Out of 361 variants identified, 40 of them were PV (18 = LDLR, 15 = APOB, 5 = PCSK9 and 2 = LDLRAP1). The majority of the PV were LDLR and APOB, where the frequency of both PV were almost similar. About 39% of clinically diagnosed FH have PV in PCSK9 alone and two novel variants of PCSK9 were identified in this study, which have not been described in Malaysia and globally. The prevalence of genetically confirmed potential FH in the community was 1:427, with a detection rate of PV at 0.2% (12/5130). About one-fourth of clinically diagnosed FH in the Malaysian community can be genetically confirmed. The detection rate of genetic confirmation is similar between potential and possible FH groups, suggesting a need for genetic confirmation in index cases from both groups. Clinical and genetic confirmation of FH index cases in the community may enhance the early detection of affected family members through family cascade screening.

TMDIM: an improved algorithm for the structure prediction of transmembrane domains of bitopic dimers.

[Formula: see text]-Helical transmembrane proteins are the most important drug targets in rational drug development. However, solving the experimental structures of these proteins remains difficult, therefore computational methods to accurately and efficiently predict the structures are in great demand. We present an improved structure prediction method TMDIM based on Park et al. (Proteins 57:577-585, 2004) for predicting bitopic transmembrane protein dimers. Three major algorithmic improvements are introduction of the packing type classification, the multiple-condition decoy filtering, and the cluster-based candidate selection. In a test of predicting nine known bitopic dimers, approximately 78% of our predictions achieved a successful fit (RMSD <2.0 Å) and 78% of the cases are better predicted than the two other methods compared. Our method provides an alternative for modeling TM bitopic dimers of unknown structures for further computational studies. TMDIM is freely available on the web at https://cbbio.cis.umac.mo/TMDIM . Website is implemented in PHP, MySQL and Apache, with all major browsers supported.

Helical integrity and microsolvation of transmembrane domains from Flaviviridae envelope glycoproteins.

Charged and polar amino acids in the transmembrane domains of integral membrane proteins can be crucial for protein function and also promote helix-helix association or protein oligomerization. Yet, our current understanding is still limited on how these hydrophilic amino acids are efficiently translocated from the Sec61/SecY translocon into the cell membrane during the biogenesis of membrane proteins. In hepatitis C virus, the putative transmembrane segments of envelope glycoproteins E1 and E2 were suggested to heterodimerize via a Lys-Asp ion-pair in the host endoplasmic reticulum. Therefore in this work, we carried out molecular dynamic simulations in explicit lipid bilayer and solvent environment to explore the stability of all possible bridging ion-pairs using the model of H-segment helix dimers. We observed that, frequently, several water molecules penetrated from the interface into the membrane core to stabilize the charged and polar pairs. The hydration time and amount of water molecules in the membrane core depended on the position of the charged residues as well as on the type of ion-pairs. Similar microsolvation events were observed in simulations of the putative E1-E2 transmembrane helix dimers. Simulations of helix monomers from other members of the Flaviviridae family suggest that these systems show similar behaviors. Thus this study illustrates the important contribution of water microsolvation to overcome the unfavorable energetic cost of burying charged and polar amino acids in membrane lipid bilayers. Also, it emphasizes the novel role of bridging charged or polar interactions stabilized by water molecules in the hydrophobic lipid bilayer core that has an important biological function for helix dimerization in several envelope glycoproteins from the family of Flaviviridae viruses.

Chaos-embedded particle swarm optimization approach for protein-ligand docking and virtual screening.

BACKGROUND: Protein-ligand docking programs are routinely used in structure-based drug design to find the optimal binding pose of a ligand in the protein's active site. These programs are also used to identify potential drug candidates by ranking large sets of compounds. As more accurate and efficient docking programs are always desirable, constant efforts focus on developing better docking algorithms or improving the scoring function. Recently, chaotic maps have emerged as a promising approach to improve the search behavior of optimization algorithms in terms of search diversity and convergence speed. However, their effectiveness on docking applications has not been explored. Herein, we integrated five popular chaotic maps-logistic, Singer, sinusoidal, tent, and Zaslavskii maps-into PSOVina[Formula: see text], a recent variant of the popular AutoDock Vina program with enhanced global and local search capabilities, and evaluated their performances in ligand pose prediction and virtual screening using four docking benchmark datasets and two virtual screening datasets. RESULTS: Pose prediction experiments indicate that chaos-embedded algorithms outperform AutoDock Vina and PSOVina in ligand pose RMSD, success rate, and run time. In virtual screening experiments, Singer map-embedded PSOVina[Formula: see text] achieved a very significant five- to sixfold speedup with comparable screening performances to AutoDock Vina in terms of area under the receiver operating characteristic curve and enrichment factor. Therefore, our results suggest that chaos-embedded PSOVina methods might be a better option than AutoDock Vina for docking and virtual screening tasks. The success of chaotic maps in protein-ligand docking reveals their potential for improving optimization algorithms in other search problems, such as protein structure prediction and folding. The Singer map-embedded PSOVina[Formula: see text] which is named PSOVina-2.0 and all testing datasets are publicly available on https://cbbio.cis.umac.mo/software/psovina .

Mutation Study of Malaysian Patients with Ornithine Transcarbamylase Deficiency: Clinical, Molecular, and Bioinformatics Analyses of Two Novel Missense Mutations of the OTC Gene.

Ornithine transcarbamylase deficiency (OTCD), an X-linked disorder that results from mutations in the OTC gene, causes hyperammonemia and leads to various clinical manifestations. Mutations occurring close to the catalytic site of OTCase can cause severe OTCD phenotypes compared with those caused by mutations occurring on the surface of this protein. In this study, we report two novel OTC missense mutations, Q171H and N199H, found in Malaysian patients. Q171H and N199H caused neonatal onset OTCD in a male and late OTCD in a female, respectively. In silico predictions and molecular docking were performed to examine the effect of these novel mutations, and the results were compared with other 30 known OTC mutations. In silico servers predicted that Q171H and N199H, as well as 30 known missense mutations, led to the development of OTCD. Docking analysis indicated that N-(phosphonoacetyl)-L-ornithine (PALO) was bound to the catalytic site of OTCase mutant structure with minimal conformational changes. However, the mutations disrupted interatomic interactions in the catalytic site. Therefore, depending on the severity of disruption occurring at the catalytic site, the mutation may affect the efficiency of mechanism and functions of OTCase.

Molecular Docking and Molecular Dynamics Simulation Studies to Predict Flavonoid Binding on the Surface of DENV2 E Protein.

Dengue infections are currently estimated to be 390 million cases annually. Yet, there is no vaccine or specific therapy available. Envelope glycoprotein E (E protein) of DENV mediates viral attachment and entry into the host cells. Several flavonoids have been shown to inhibit HIV-1 and hepatitis C virus entry during the virus-host membrane fusion. In this work, molecular docking method was employed to predict the binding of nine flavonoids (baicalin, baicalein, EGCG, fisetin, glabranine, hyperoside, ladanein, quercetin and flavone) to the soluble ectodomain of DENV type 2 (DENV2) E protein. Interestingly, eight flavonoids were found to dock into the same binding pocket located between the domain I and domain II of different subunits of E protein. Consistent docking results were observed not only for the E protein structures of the DENV2-Thai and DENV2-Malaysia (a homology model) but also for the E protein structures of tick-borne encephalitis virus and Japanese encephalitis virus. In addition, molecular dynamics simulations were performed to further evaluate the interaction profile of the docked E protein-flavonoid complexes. Ile4, Gly5, Asp98, Gly100 and Val151 residues of the DENV2-My E protein that aligned to the same residues in the DENV2-Thai E protein form consistent hydrogen bond interactions with baicalein, quercetin and EGCG during the simulations. This study demonstrates flavonoids potentially form interactions with the E protein of DENV2.

ENRI: A tool for selecting structure-based virtual screening target conformations.

Finding pharmaceutically relevant target conformations from an arbitrary set of protein conformations remains a challenge in structure-based virtual screening (SBVS). The growth in the number of available conformations, either experimentally determined or computationally derived, obscures the situation further. While the inflated conformation space potentially contains viable druggable targets, the increase of conformational complexity, as a consequence, poses a selection problem. To address this challenge, we took advantage of machine learning methods, namely an over-sampling and a binary classification procedure, and present a novel method to select druggable receptor conformations. Specifically, we trained a binary classifier on a set of nuclear receptor conformations, wherein each conformation was labeled with an enrichment measure for a corresponding SBVS. The classifier enabled us to formulate suggestions and identify enriching SBVS targets for six of seven nuclear receptors. Further, the classifier can be extended to other proteins of interest simply by feeding new training data sets to the classifier. Our work, thus, provides a methodology to identify pharmaceutically interesting receptor conformations for nuclear receptors and other drug targets.

Stability, orientation and position preference of the stem region (residues 689-703) in Hepatitis C Virus (HCV) envelope glycoprotein E2: a molecular dynamics study.

Envelope glycoproteins of Hepatitis C Virus (HCV) play an important role in the virus assembly and initial entry into host cells. Conserved charged residues of the E2 transmembrane (TM) domain were shown to be responsible for the heterodimerization with envelope glycoprotein E1. Despite intensive research on both envelope glycoproteins, the structural information is still not fully understood. Recent findings have revealed that the stem (ST) region of E2 also functions in the initial stage of the viral life cycle. We have previously shown the effect of the conserved charged residues on the TM helix monomer of E2. Here, we extended the model of the TM domain by adding the adjacent ST segment. Explicit molecular dynamics simulations were performed for the E2 amphiphilic segment of the ST region connected to the putative TM domain (residues 683-746). Structural conformation and behavior are studied and compared with the nuclear magnetic resonance (NMR)-derived segment of E2 ( 2KQZ.pdb). We observed that the central helix of the ST region (residues 689 - 703) remained stable as a helix in-plane to the lipid bilayer. Furthermore, the TM domain appeared to provide minimal contribution to the structural stability of the amphipathic region. This study also provides insight into the orientation and positional preferences of the ST segment with respect to the membrane lipid-water interface.

Contribution of charged and polar residues for the formation of the E1-E2 heterodimer from Hepatitis C Virus.

The transmembrane domains of the envelope glycoprotein E1 and E2 have crucial multifunctional roles in the biogenesis of hepatitis C virus. We have performed molecular dynamics simulations to investigate a structural model of the transmembrane segments of the E1-E2 heterodimer. The simulations support the key role of the Lys370-Asp728 ion pair for mediating the E1-E2 heterodimerization. In comparison to these two residues, the simulation results also reveal the differential effect of the conserved Arg730 residue that has been observed in experimental studies. Furthermore, we discovered the formation of inter-helical hydrogen bonds via Asn367 that stabilize dimer formation. Simulations of single and double mutants further demonstrate the importance of the ion-pair and polar interactions between the interacting helix monomers. The conformation of the E1 fragment in the simulation of the E1-E2 heterodimer is in close agreement with an NMR structure of the E1 transmembrane segment. The proposed model of the E1-E2 heterodimer supports the postulated cooperative insertion of both helices by the translocon complex into the bilayer.