Biomedical subjects

Louis A Clark

Publications and source records attributed to Louis A Clark.

7 recordsLinked to original sources

PubMed2006

Trends in antibody sequence changes during the somatic hypermutation process.

Probable germline gene sequences from thousands of aligned mature Ab sequences are inferred using simple computational matching to known V(D)J genes. Comparison of the germline to mature sequences in a structural region-dependent fashion allows insights into the methods that nature uses to mature Abs during the somatic hypermutation process. Four factors determine the residue type mutation patterns: biases in the germline, accessibility from single base permutations, location of mutation hotspots, and functional pressures during selection. Germline repertoires at positions that commonly contact the Ag are biased with tyrosine, serine, and tryptophan. These residue types have a high tendency to be present in mutation hotspot motifs, and their abundance is decreased during maturation by a net conversion to other types. The heavy use of tyrosines on mature Ab interfaces is thus a reflection of the germline composition rather than being due to selection during maturation. Potentially stabilizing changes such as increased proline usage and a small number of double cysteine mutations capable of forming disulfide bonds are ascribed to somatic hypermutation. Histidine is the only residue type for which usage increases in each of the interface, core, and surface regions. The net overall effect is a conversion from residue types that could provide nonspecific initial binding into a diversity of types that improve affinity and stability. Average mutation probabilities are approximately 4% for core residues, approximately 5% for surface residues, and approximately 12% for residues in common Ag-contacting positions, excepting the those coded by the D gene.

Amino Acid Sequence↗

PubMed2006

Affinity enhancement of an in vivo matured therapeutic antibody using structure-based computational design.

Improving the affinity of a high-affinity protein-protein interaction is a challenging problem that has practical applications in the development of therapeutic biomolecules. We used a combination of structure-based computational methods to optimize the binding affinity of an antibody fragment to the I-domain of the integrin VLA1. Despite the already high affinity of the antibody (Kd approximately 7 nM) and the moderate resolution (2.8 A) of the starting crystal structure, the affinity was increased by an order of magnitude primarily through a decrease in the dissociation rate. We determined the crystal structure of a high-affinity quadruple mutant complex at 2.2 A. The structure shows that the design makes the predicted contacts. Structural evidence and mutagenesis experiments that probe a hydrogen bond network illustrate the importance of satisfying hydrogen bonding requirements while seeking higher-affinity mutations. The large and diverse set of interface mutations allowed refinement of the mutant binding affinity prediction protocol and improvement of the single-mutant success rate. Our results indicate that structure-based computational design can be successfully applied to further improve the binding of high-affinity antibodies.

Amino Acid Substitution↗

PubMed2005

Simulated adsorption properties and synthesis prospects of homochiral porous solids based on their heterochiral analogs.

Molecular simulations of chiral molecules in porous heterochiral materials were performed to investigate fundamental adsorption properties and possibilities for production of homochiral porous solids. Zeolite BEA polymorph A and zeotype UCSB-7K each provide separated pores of opposite chirality. Single enantiomer and racemic mixture adsorption results are presented and indicate that significant equilibrium enantiomeric excesses of 40-70% in UCSB-7K and 10% in BEA can be achieved. Larger, better-fitting molecules display higher enantiomeric excesses. For dimethylallene, which moves on molecular dynamics time scales in UCSB-7K, self-diffusivities vary by almost an order of magnitude between the two opposite-handed UCSB-7K pores for a given enantiomer. The predicted properties indicate that equilibrium and nonequilibrium strategies using related homochiral materials for separations may be successful. To this end, a discussion of strategies for selectively blocking pores of one chirality on the basis of enantiomer segregation is provided.

Journal Article↗

PubMed2005

Protein aggregation determinants from a simplified model: cooperative folders resist aggregation.

Two-chain aggregation simulations using minimalist models of proteins G, L, and mutants were used to investigate the fundamentals of protein aggregation. Mutations were selected to break up repeats of hydrophobic beads in the sequence while maintaining native topology and folding ability. Data are collected under conditions in which all chain types have similar folded populations and after equilibrating the separated chains to minimize competition between folding and aggregation. Folding cooperativity stands out as the best single-chain determinant under these conditions and for these simple models. It can be experimentally measured by the width of the unfolding transition during thermal denaturation and loosely related to population of intermediate-like states during folding. Additional measures of cooperativity and other properties such as radius of gyration fluctuations and patterning of hydrophobic residues are also examined. Initial contact system states with transition-state characteristics can be identified and are more expanded than average initial contact states. Two-chain minimalist model aggregates are considerably less structured than their native states and have minimal domain-swapping features.

Amino Acid Sequence↗

PubMed2005

Influence of denatured and intermediate states of folding on protein aggregation.

We simulate the aggregation thermodynamics and kinetics of proteins L and G, each of which self-assembles to the same alpha/beta [corrected] topology through distinct folding mechanisms. We find that the aggregation kinetics of both proteins at an experimentally relevant concentration exhibit both fast and slow aggregation pathways, although a greater proportion of protein G aggregation events are slow relative to those of found for protein L. These kinetic differences are correlated with the amount and distribution of intrachain contacts formed in the denatured state ensemble (DSE), or an intermediate state ensemble (ISE) if it exists, as well as the folding timescales of the two proteins. Protein G aggregates more slowly than protein L due to its rapidly formed folding intermediate, which exhibits native intrachain contacts spread across the protein, suggesting that certain early folding intermediates may be selected for by evolution due to their protective role against unwanted aggregation. Protein L shows only localized native structure in the DSE with timescales of folding that are commensurate with the aggregation timescale, leaving it vulnerable to domain swapping or nonnative interactions with other chains that increase the aggregation rate. Folding experiments that characterize the structural signatures of the DSE, ISE, or the transition state ensemble (TSE) under nonaggregating conditions should be able to predict regions where interchain contacts will be made in the aggregate, and to predict slower aggregation rates for proteins with contacts that are dispersed across the fold. Since proteins L and G can both form amyloid fibrils, this work also provides mechanistic and structural insight into the formation of prefibrillar species.

Bacterial Proteins↗

PubMed2004

Computational elucidation of the transition state shape selectivity phenomenon.

The most commonly cited example of a transition state shape selective reaction, m-xylene disproportionation in zeolites, is examined to determine if the local spatial environment of a reaction can significantly alter selectivity. In the studied reaction, ZPE-corrected rate limiting energy barriers are 136 kJ/mol for the methoxide-mediated pathway and 109 to 145 kJ/mol for the diphenylmethane-mediated pathway. Both pathways are likely to contribute to selectivity and disfavor one product isomer (1,3,5-trimethylbenzene), but relative selectivity to the other two isomers varies with pore geometry, mechanistic pathway, and inclusion of entropic effects. Most importantly, study of one pathway in three different common zeolite framework types (FAU, MFI, and MOR) allows explicit and practically oriented consideration of pore shape. Variation of the environment shape at the critical transition states is thus shown to affect the course of reaction. Barrier height shifts on the order of 10-20 kJ/mol are achievable. Observed selectivities do not agree with the transition state characteristics calculated here and, hence, are most likely due to product shape selectivity. Further examination of the pathways highlights the importance of mechanistic steps that do not result in isomer-defining bonds and leads to a more robust definition of transition state shape selectivity.

Journal Article↗

PubMed2003

Stable mechanistically-relevant aromatic-based carbenium ions in zeolite catalysts.

The existence of carbenium ion species is assumed in many zeolite catalysis mechanisms. Using computational techniques that include environmental effects, a benzenium-type carbenium ion is identified in zeolite catalysts for the first time. Localization of nearby transition states indicate that this species may play an important role as an intermediate in the bimolecular m-xylene disproportionation reaction. The barrier to back-donation of the proton from the benzenium ion is at least 50 kJ/mol, meaning that this species may be spectroscopically observable. An additional carbenium ion intermediate, formed by abstraction of a hydride from m-xylene, is also predicted. The stability of this second new carbenium ion suggests that aromatic-based carbenium ions are likely to be intermediates in many zeolite-catalyzed reactions. Two types of fundamentally different fully periodic calculations support the stability predictions.

Journal Article↗