Characterization of the flexibility of the peripheral stalk of prokaryotic rotary A‐ATPases by atomistic simulations

Rotary ATPases are involved in numerous physiological processes, with the three distinct types (F/A/V‐ATPases) sharing functional properties and structural features. The basic mechanism involves the counter rotation of two motors, a soluble ATP hydrolyzing/synthesizing domain and a membrane‐embedded ion pump connected through a central rotor axle and a stator complex. Within the A/V‐ATPase family conformational flexibility of the EG stators has been shown to accommodate catalytic cycling and is considered to be important to function. For the A‐ATPase three EG structures have been reported, thought to represent conformational states of the stator during different stages of rotary catalysis. Here we use long, detailed atomistic simulations to show that those structures are conformers explored through thermal fluctuations, but do not represent highly populated states of the EG stator in solution. We show that the coiled coil tail domain has a high persistence length (～100 nm), but retains the ability to adapt to different conformational states through the presence of two hinge regions. Moreover, the stator network of the related V‐ATPase has been suggested to adapt to subunit interactions in the collar region in addition to the nucleotide occupancy of the catalytic domain. The MD simulations reported here, reinforce this observation showing that the EG stators have enough flexibility to adapt to significantly different structural re‐arrangements and accommodate structural changes in the catalytic domain whilst resisting the large torque generated by catalytic cycling. These results are important to understand the role the stators play in the rotary‐ATPase mechanism. Proteins 2016; 84:1203–1212. © 2016 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Systematic classification of CDR-L3 in antibodies: implications of the light chain subtypes and the VL-VH interface

Antibody modeling is widely used for the analysis of antibody-antigen interactions and for the design of potent antibody drugs. The antibody combining site is composed of six complementarity determining regions (CDRs). The CDRs, except for CDR-H3, which is the most diverse CDR, form limited numbers of canonical structures, which can be identified from the amino acid sequences. A method to classify the CDR-H3 structure from its amino acid sequence was previously proposed. However, since those CDR structures were classified, many more antibody crystal structures have been determined. We performed systematic analyses of the CDR-L3 structures and found novel canonical structures, and we also classified a previously identified canonical structure into two subtypes. In addition, two differently defined canonical structures in the kappa and lambda subtypes were classified into the same canonical structure. We also identified a key residue in CDR-L3, which determines the conformation of CDR-H3. Several analyses of CDR-L3 loops longer than nine residues were performed. These new findings should be useful for structural modeling and are eventually expected to accelerate the design of antibody drugs.     

Automated docking of monosaccharide substrates and analogues and methyl α-Acarviosinide in the glucoamylase active site

Glucoamylase is an important industrial glucohydrolase with a large specificity range. To investigate its interaction with the monosaccharides D-glucose, D-mannose, and D-galactose and with the substrate analogues 1-deoxynojirimycin, D-glucono-1,5-lactone, and methyl αacarviosinide, MM3(92)-optimized structures were docked into its active site using AutoDock 2.1. The results were compared to structures of glucoamylase complexes obtained by protein crystallography. Charged forms of some substrate analogues were also docked to assess the degree of protonation possessed by glucoamylase inhibitors. Many forms of methyl αa-carviosinide were conformationally mapped by using MM3(92), characterizing the conformational pH dependence found for the acarbose family of glucosidase inhibitors. Their significant conformers, representing the most common states of the inhibitor, were used as initial structures for docking. This constitutes a new approach for the exploration of binding modes of carbohydrate chains. Docking results differ slightly from x-ray crystallographic data, the difference being of the order of the crystallographic error. The estimated energetic interactions, even though agreeing in some cases with experimental binding kinetics, are only qualitative due to the large approximations made by AudoDock force field. © 1997 Wiley-Liss, Inc.     

Computation of 3D structure and distribution of helical parameters for D-period segments of human fibrillar collagen III

The structures of D-period segments of collagen (234 amino residues or ～1/4 of whole length) are established by methods of molecular mechanics and geometry analysis. Each D-period segment proves to have a unique spatial structure. The distributions of local helical parameters along the molecule are calculated. It is found that a second hydrogen bond is formed in every case when the second residue in the tripeptide G-X-Y is an amino acid. With such a combined H-bond network, all the peptide CO groups of glycines and of the third residues in tripeptides have quasi-equivalent positions on the surface of the collagen molecule. The local deformations of the polyproline II helix in the triple complex give rise to the observed modulation of structure at the macromolecular level, which may be important for the mutual orientation of collagen molecules during fibrillogenesis.     

Molecular mechanical properties of short-sequence peptide enzyme mimics

While considerable attempts have been made to recreate the high turnover rates of enzymes using synthetic enzyme mimics, most have failed and only a few have produced minimal reaction rates that can barely be considered catalytic. One particular approach we have focused on is the use of short-sequence peptides that contain key catalytic groups in close proximity. In this study, we designed six different peptides and tested their ability to mimic the catalytic mechanism of the cysteine proteases. Acetylation and deacylation by Ellman’s Reagent trapping experiments showed the importance of having phenylalanine groups surrounding the catalytic sites in order to provide greater proximity between the cysteine, histidine, and aspartate amino acid R-groups. We have also carried out all-atom molecular dynamics simulations to determine the distance between these catalytic groups and the overall mechanical flexibility of the peptides. We found strong correlations between the magnitude of fluctuations in the Cys-His distance, which determines the flexibility and interactions between the cysteine thiol and histidine imidazole groups, and the deacylation rate. We found that, in general, shorter Cys-His distance fluctuations led to a higher deacylation rate constant, implying that greater confinement of the two residues will allow a higher frequency of the acetyl exchange between the cysteine thiol and histidine imidazole R-groups. This may be the key to future design of peptide structures with molecular mechanical properties that lead to viable enzyme mimics.     

Structural homology of spinach acyl carrier protein and Escherichia coli acyl carrier protein based on NMR data

Acyl carrier proteins (ACPs) from spinach and from Escherichia coli have been used to demonstrate the utility of proton NMR for comparison of homologous structures. The structure of E. coli ACP had been previously determined and modeled as a rapid equilibrium among multiple conformational forms (Kim and Prestegard, Biochemistry 28:8792–8797, 1989). Spinach ACP showed two slowly exchanging forms and could be manipulated into one form for structural study. Here we compare this single form to postulated multiple forms of E. coli ACP using the limited amount of NOE data available for the spinach protein. A number of long-range NOE contacts were present between homologous residues in both spinach and E. coli ACP, suggesting tertiary structural homology. To allow a more definitive structural comparison, a method was developed to use spinach ACP NOE constraints to search for regions of structural divergence from two postulated forms of E. coli ACP. The homologous regions of the two protein sequences were aligned, additional distance constraints were extracted from the E. coli structure, and these were mapped onto the spinach sequence. These distance constraints were combined with experimental NOE constraints and a distance geometry simulated annealing protocol was used to test for compatibility of the constraints. All of the experimental spinach NOE constraints could be successfully combined with the E. coli data, confirming the general hypothesis of structural homology. A better fit was obtained with one form, suggesting a preferential stabilization of that form in the spinach case. Proteins 27:131–143 © 1997 Wiley-Liss, Inc.     

Molecular modeling of green fluorescent protein: structural effects of chromophore deprotonation

Molecular dynamics (MD) simulations were carried out to study the conformational rearrangement induced by deprotonation of the fluorescent chromophore in GFP, as well as the associated changes in the hydrogen-bonding network. For both the structures with either a neutral or an anionic chromophore, it was found that the beta-barrel was stable and rigid, and the conformation of the chromophore was consistent with the available x-ray structure. The conformational change in Thr203 due to deprotonation was also found to be consistent with the three-state isomerization model. Although GFP is highly fluorescent, denatured-GFP is nonfluorescent, indicating that the environment of the protein plays an important role in its fluorescence behavior. Our MD simulations, which explore the effect of the protein shell on the conformation of the chromophore, find the flexibility of the central chromophore to be significantly restricted due to the rigid nature of the protein shell. The hydrogen-bonding between the chromophore and neighboring residues was also shown to contribute to the chromophore rigidity. In addition to the MD studies, quantum mechanics/molecular mechanics (QM/MM) ONIOM calculations were carried out to investigate the effect of the beta-barrel on the internal rotation in the chromophore. Along with providing quantitative values for torsional rotation barriers about the bridging bond in the chromophore, the ONIOM calculations also validate our MD force field parameters.     

Using molecular dynamics simulations on crambin to evaluate the suitability of different continuum dielectric and hydrogen atom models for protein simulations

Molecular dynamics simulations of enzymes with enough explicit waters of solvation to realistically account for solute-solvent interactions can burden the computational resources required to perform the simulation by more than two orders of magnitude. Since enzyme simulations even with an implicit solvation model can be imposing for a supercomputer, it is important to assess the suitability of different continuum dielectric models for protein simulations. A series of 100-picosecond molecular dynamics simulations were performed on the X-ray crystal structure of the protein crambin to examine how well computed structures, obtained using seven continuum dielectric and two hydrogen atom models, agreed with the X-ray structure. The best level of agreement between computed and experimental structures was obtained using a constant dielectric of 2 and the all-hydrogen model. Continuum dielectric models of 1, 1r, and 2r also led to computed structures in reasonably good agreement with the X-ray structure. In all cases, the all-hydrogen model gave better agreement than the united atom model, although, in one case, the difference was not significant. Dielectric models of 4, 80, and 4r with either hydrogen model yielded significantly poorer fits. It is especially noteworthy that the observed trends did not semiquantitatively converge until about 50 picoseconds into the simulations, suggesting that validation studies for protein calculations based on energy minimizations or short simulations should be viewed with caution.     

The response of T4 lysozyme to large-to-small substitutions within the core and its relation to the hydrophobic effect.

To further examine the structural and thermodynamic basis of hydrophobic stabilization in proteins, all of the bulky non-polar residues that are buried or largely buried within the core of T4 lysozyme were substituted with alanine. In 25 cases, including eight reported previously, it was possible to determine the crystal structures of the variants. The structures of four variants with double substitutions were also determined. In the majority of cases the "large-to-small" substitutions lead to internal cavities. In other cases declivities or channels open to the surface were formed. In some cases the structural changes were minimal (mainchain shifts < or = 0.3 A); in other cases mainchain atoms moved up to 2 A. In the case of Ile 29 --> Ala the structure collapsed to such a degree that the volume of the putative cavity was zero. Crystallographic analysis suggests that the occupancy of the engineered cavities by solvent is usually low. The mutants Val 149 --> Ala (V149A) and Met 6 --> Ala (M6A), however, are exceptions and have, respectively, one and two well-ordered water molecules within the cavity. The Val 149 --> Ala substitution allows the solvent molecule to hydrogen bond to polar atoms that are occluded in the wild-type molecule. Similarly, the replacement of Met 6 with alanine allows the two solvent molecules to hydrogen bond to each other and to polar atoms on the protein. Except for Val 149 --> Ala the loss of stability of all the cavity mutants can be rationalized as a combination of two terms. The first is a constant for a given class of substitution (e.g., -2.1 kcal/mol for all Leu --> Ala substitutions) and can be considered as the difference between the free energy of transfer of leucine and alanine from solvent to the core of the protein. The second term can be considered as the energy cost of forming the cavity and is consistent with a numerical value of 22 cal mol(-1) A(-3). Physically, this term is due to the loss of van der Waal''s interactions between the bulky sidechain that is removed and the atoms that form the wall of the cavity. The overall results are consistent with the prior rationalization of Leu --> Ala mutants in T4 lysozyme by Eriksson et al. (Eriksson et al., 1992, Science 255:178-183).     

Progress in fold recognition

The prediction experiment reveals that fold recognition has become a powerful tool in structural biology. We applied our fold recognition technique to 13 target sequences. In two cases, replication terminating protein and prosequence of subtilisin, the predicted structures are very similar to the experimentally determined folds. For the first time, in a public blind test, the unknown structures of proteins have been predicted ahead of experiment to an accuracy approaching molecular detail. In two other cases the approximate folds have been predicted correctly. According to the assessors there were 12 recognizable folds among the target proteins. In our postprediction analysis we find that in 7 cases our fold recognition technique is successful. In several of the remaining cases the predicted folds have interesting features in common with the experimental results. We present our procedure, discuss the results, and comment on several fundamental and technical problems encountered in fold recognition. © 1995 Wiley-Liss, Inc.     

Protein-protein docking by simulating the process of association subject to biochemical constraints

We present a computational procedure for modeling protein-protein association and predicting the structures of protein-protein complexes. The initial sampling stage is based on an efficient Brownian dynamics algorithm that mimics the physical process of diffusional association. Relevant biochemical data can be directly incorporated as distance constraints at this stage. The docked configurations are then grouped with a hierarchical clustering algorithm into ensembles that represent potential protein-protein encounter complexes. Flexible refinement of selected representative structures is done by molecular dynamics simulation. The protein-protein docking procedure was thoroughly tested on 10 structurally and functionally diverse protein-protein complexes. Starting from X-ray crystal structures of the unbound proteins, in 9 out of 10 cases it yields structures of protein-protein complexes close to those determined experimentally with the percentage of correct contacts >30% and interface backbone RMSD <4 A. Detailed examination of all the docking cases gives insights into important determinants of the performance of the computational approach in modeling protein-protein association and predicting of protein-protein complex structures.     

Modeling of inhibitor–metalloenzyme interactions and selectivity using molecular mechanics grounded in quantum chemistry

We investigated the binding properties of the metalloprotease inhibitors hydroxamate, methanethiolate, and methylphosphoramidate to a model coordination site occurring in several Zn2+ metalloproteases, including thermolysin. This was carried out using both the SIBFA (sum of interactions between fragments ab initio-computed) molecular mechanics and the SCF/MP2 procedures for the purpose of evaluating SIBFA as a metalloenzyme modeling tool. The energy-minimized structures were closely similar to the X-ray crystallographic structures of related thermolysin-inhibitor complexes. We found that selectivity between alternative geometries and between inhibitors usually stemmed from multiple interaction components included in SIBFA. The binding strength sequence is hydroxamate > methanethiolate ≥ methylphosphoramidate from multiple interaction components included in SIBFA. The trends in interaction energy components, rankings, and preferences for mono- or bidentate binding were consistent in both computational procedures. We also compared the Zn2+ vs. Mg2+ selectivities in several other polycoordinated sites having various “hard” and “soft” qualities. This included a hexahydrate, a model representing Mg2+/Ca2+ binding sites, a chlorophyll-like structure, and a zinc finger model. The latter three favor Zn2+ over Mg2+ by a greater degree than the hydrated state, but the selectivity varies widely according to the ligand “softness.” SIBFA was able to match the ab initio binding energies by <2%, with the SIBFA terms representing dispersion and charge-transfer contributing the most to Zn2+/Mg2+ selectivity. These results showed this procedure to be a very capable modeling tool for metalloenzyme problems, in this case giving valuable information about details and limitations of “hard” and “soft” selectivity trends. Proteins 31:42–60, 1998. © 1998 Wiley-Liss, Inc.     

Toward better refinement of comparative models: predicting loops in inexact environments

Achieving atomic-level accuracy in comparative protein models is limited by our ability to refine the initial, homolog-derived model closer to the native state. Despite considerable effort, progress in developing a generalized refinement method has been limited. In contrast, methods have been described that can accurately reconstruct loop conformations in native protein structures. We hypothesize that loop refinement in homology models is much more difficult than loop reconstruction in crystal structures, in part, because side-chain, backbone, and other structural inaccuracies surrounding the loop create a challenging sampling problem; the loop cannot be refined without simultaneously refining adjacent portions. In this work, we single out one sampling issue in an artificial but useful test set and examine how loop refinement accuracy is affected by errors in surrounding side-chains. In 80 high-resolution crystal structures, we first perturbed 6-12 residue loops away from the crystal conformation, and placed all protein side chains in non-native but low energy conformations. Even these relatively small perturbations in the surroundings made the loop prediction problem much more challenging. Using a previously published loop prediction method, median backbone (N-Calpha-C-O) RMSD's for groups of 6, 8, 10, and 12 residue loops are 0.3/0.6/0.4/0.6 A, respectively, on native structures and increase to 1.1/2.2/1.5/2.3 A on the perturbed cases. We then augmented our previous loop prediction method to simultaneously optimize the rotamer states of side chains surrounding the loop. Our results show that this augmented loop prediction method can recover the native state in many perturbed structures where the previous method failed; the median RMSD's for the 6, 8, 10, and 12 residue perturbed loops improve to 0.4/0.8/1.1/1.2 A. Finally, we highlight three comparative models from blind tests, in which our new method predicted loops closer to the native conformation than first modeled using the homolog template, a task generally understood to be difficult. Although many challenges remain in refining full comparative models to high accuracy, this work offers a methodical step toward that goal.     

Psychotherapy of alcoholics; the concept of psychotherapeutic motivation

The psychiatric patient's compliance with requirements of therapy (attendance, responsiveness, verbalization, etc.) often is adopted by the therapist as the chief criterion of motivation for cure. By these standards alcoholics often are judged to be poorly motivated and unsuited to the therapist's mode of therapy. In such cases, the method may have to be adapted to the case, and a more permissive attitude taken toward the patient's apparent noncooperation.     

Three-dimensional finite element simulations of the mechanical response of the fingertip to static and dynamic compressions

The analysis of the mechanics of the contact interactions of fingers/handle and the stress/strain distributions in the soft tissues in the fingertip is essential to optimize design of tools to reduce many occupation-related hand disorders. In the present study, a three-dimensional (3D) finite element (FE) model for the fingertip is proposed to simulate the nonlinear and time-dependent responses of a fingertip to static and dynamic loadings. The proposed FE model incorporates the essential anatomical structures of a finger: skin layers (outer and inner skins), subcutaneous tissue, bone and nail. The soft tissues (inner skin and subcutaneous tissue) are considered to be nonlinearly viscoelastic, while the hard tissues (outer skin, bone and nail) are considered to be linearly elastic. The proposed model has been used to simulate two loading scenarios: (a) the contact interactions between the fingertip and a flat surface and (b) the indentation of the fingerpad via a sharp wedge. For case (a), the predicted force/displacement relationships and time-dependent force responses are compared with the published experimental data; for case (b), the skin surface deflection profiles were predicted and compared with the published experimental observations. Furthermore, for both cases, the time-dependent stress/strain distributions within the tissues of the fingertip were calculated. The good agreement between the model predictions and the experimental observations indicates that the present model is capable of predicting realistic time-dependent force/displacement responses and stress/strain distributions in the soft tissues for dynamic loading conditions.     

Assessment of the ability to model proteins with leucine-rich repeats in light of the latest structural information

The three-dimensional structures of leucine-rich repeat (LRR)-containing proteins from five different families were previously predicted based on the crystal structure of the ribonuclease inhibitor, using an approach that combined homology-based modeling, structure-based sequence alignment of LRRs, and several rational assumptions. The structural models have been produced based on very limited sequence similarity, which, in general, cannot yield trustworthy predictions. Recently, the protein structures from three of these five families have been determined. In this report we estimate the quality of the modeling approach by comparing the models with the experimentally determined structures. The comparison suggests that the general architecture, curvature, "interior/exterior" orientations of side chains, and backbone conformation of the LRR structures can be predicted correctly. On the other hand, the analysis revealed that, in some cases, it is difficult to predict correctly the twist of the overall super-helical structure. Taking into consideration the conclusions from these comparisons, we identified a new family of bacterial LRR proteins and present its structural model. The reliability of the LRR protein modeling suggests that it would be informative to apply similar modeling approaches to other classes of solenoid proteins.     

Sampling and scoring: A marriage made in heaven

Most structure prediction algorithms consist of initial sampling of the conformational space, followed by rescoring and possibly refinement of a number of selected structures. Here we focus on protein docking, and show that while decoupling sampling and scoring facilitates method development, integration of the two steps can lead to substantial improvements in docking results. Since decoupling is usually achieved by generating a decoy set containing both non‐native and near‐native docked structures, which can be then used for scoring function construction, we first review the roles and potential pitfalls of decoys in protein–protein docking, and show that some type of decoys are better than others for method development. We then describe three case studies showing that complete decoupling of scoring from sampling is not the best choice for solving realistic docking problems. Although some of the examples are based on our own experience, the results of the CAPRI docking and scoring experiments also show that performing both sampling and scoring generally yields better results than scoring the structures generated by all predictors. Next we investigate how the selection of training and decoy sets affects the performance of the scoring functions obtained. Finally, we discuss pathways to better alignment of the two steps, and show some algorithms that achieve a certain level of integration. Although we focus on protein–protein docking, our observations most likely also apply to other conformational search problems, including protein structure prediction and the docking of small molecules to proteins.Proteins 2013; 81:1874–1884. © 2013 Wiley Periodicals, Inc.     

Investigation of the molecular similarity in closely related protein systems: The PrP case study

The amyloid conversion is a massive detrimental modification affecting several proteins upon specific physical or chemical stimuli characterizing a plethora of diseases. In many cases, the amyloidogenic stimuli induce specific structural features to the protein conferring the propensity to misfold and form amyloid deposits. The investigation of mutants, structurally similar to their native isoform but inherently prone to amyloid conversion, may be a viable strategy to elucidate the structural features connected with amyloidogenesis. In this article, we present a computational protocol based on the combination of molecular dynamics (MD) and grid‐based approaches suited for the pairwise comparison of closely related protein structures. This method was applied on the cellular prion protein (PrP^C) as a case study and, in particular, addressed to the quali/quantification of the structural features conferred by either E200K mutations and treatment with CaCl₂, both able to induce the scrapie conversion of PrP. Several schemes of comparison were developed and applied to this case study, and made up suitable of application to other protein systems. At this purpose an in‐house python codes has been implemented that, together with the parallelization of the GRID force fields program, will spread the applicability of the proposed computational procedure. Proteins 2015; 83:1751–1765. © 2015 Wiley Periodicals, Inc.     

不同泪道硅胶引流管留置时间对慢性泪囊炎患者视力、生活质量及复发率的影响

目的:探讨不同泪道硅胶引流管留置时间对慢性泪囊炎患者的生活质量、视力及复发率的影响。方法:回顾性选取2018年1月2019年12月期间于我院就诊的慢性泪囊炎患者91例(128眼),根据术后泪道硅胶引流管留置时间分为A、B两组,其中A组泪道硅胶引流管留置时间6周,44例(61眼),B组泪道硅胶引流管留置时间12周,47例(67眼)。对比两组疗效、视力、并发症发生率、主诉溢泪发生率、生活质量及复发率。结果:A组拔管当天的总有效率高于B组(P<0.05),两组拔管后3个月总有效率对比无明显差异(P>0.05)。拔管后6个月两组社会功能、躯体疼痛、精神健康、生理功能、精力、情感职能、生理职能、总体健康维度评分均较术前升高(P<0.05),但两组组间对比未见统计学差异(P>0.05)。两组患者术前、拔管后3个月视力组间及组内比较均未见统计学差异(P>0.05)。两组患者主诉溢泪发生率比较未见统计学差异(P>0.05)。A组并发症发生率、复发率低于B组(P<0.05)。结论:泪道硅胶引流管留置时间的长短对慢性泪囊炎患者疗效、视力、生活质量、主诉溢泪发生率无明显影响,但留置6周者并发症发生率、复发率低于留置12周者,提示临床应视患者具体情况尽量缩短泪道硅胶引流管留置时间。     

Molecular modelling of the inclusion complexes between β-cyclodextrin and (R)/(S)-methylphenobarbitone and its application to HPLC

Molecular modelling of β-cyclodextrin and optimisation of its potential energy suggests that a favoured conformation is that distorted from a symmetrical torus. The inclusion of water molecules into the torus cavity simulates the increased stability in an aqueous solvent. Complexes of β-cyclodextrin with (R)- and (S)-enantiomers of methylphenobarbitone have been modelled and energetically optimised by the application of molecular mechanics. The simulations suggests that the guest molecules adopt an orientation in which the phenyl ring is projected into the torus cavity, with in each case the plane of the ring parallel to a longer axis of the distorted torus and slightly displaced from the axis through the torus cavity. It is suggested that the asymmetry in the macrocyclic ring contributes to chiral recognition as a result of additional discriminatory binding to the barbiturate ring residue of each enantiomer, which occupy different 3D geometries. The enantiomers form complexes of different minimum potential energies. The resulting difference in complex stability can be related to the behaviour of β-cyclodextrin, as a mobile phase additive in reverse-phase HPLC to effect chiral separation of rac-medthylphenobarbitone during chromatography. © 1994 Wiley-Liss, Inc.